北航第二次应用泛函作业

第二次应用泛函作业

从数学的本质来看，最基本的集合有两类：线性空间（有线性结构的集合）、度量空间（有度量结构的集合）。对线性空间而言，主要研究集合的描述，直观地说就是如何清楚地告诉地别人这个集合是什么样子。为了描述清楚，就引入了基（相当于三维空间中的坐标系）的概念，所以对于一个线性空间来说，只要知道其基即可，集合中的元素只要知道其在给定基下的坐标即可。

但线性空间中的元素没有“长度”（相当于三维空间中线段的长度），为了量化线性空间中的元素，所以又在线性空间引入特殊的“长度”，即范数。赋予了范数的线性空间即称为赋范线性空间。但赋范线性空间中两个元素之间没有角度的概念，为了解决该问题，所以在线性空间中又引入了内积的概念。

因为有度量，所以可以在度量空间、赋范线性空间以及内积空间中引入极限，但抽象空间中的极限与实数上的极限有一个很大的不同就是，极限点可能不在原来给定的集合中，所以又引入了完备的概念，完备的内积空间就称为Hilbert空间。

线性赋范空间就是定义了范数的线性空间，所谓范数就是线性空间到数域的一个映射，其满足范数公理（正定性，齐次性，三角不等式），可以理解为线性空间元素的长度。内积空间就是定义了内积的线性空间，而内积可以看成是两个元素作用生成一个数，按一般向量内积理解即可。度量空间是定义了度量的线性空间，也就是两个元素之间的“长度”，满足正定性、对称性、三角不等式。

一般而言，定义了内积可以诱导出范数（也就是与自己做内积再开根号），定义了范数可以诱导出度量（两元素的度量即为元素差的范数），但度量诱导范数需要加一点限制。所谓希尔伯特空间就是定义了内积的线性空间，并且按照内积诱导出的度量是完备的（完备就是柯西序列在内部收敛）特别的，实数域上的有限维希尔伯特空间叫做欧几里得空间；复数域上的有限维希尔伯特空间叫做酉空间。

泛函分析的空间理论一个有重大意义的事实是：n

R中的极限论及基于极限概念的分析理论的许多结果，仅依赖于模长性质，而与模长的定义无关，因而实际上并不依赖于Euclid空间的特殊构造。这就说明了若某个向量空间上定义了一种类似于模长的概念，它具有正定性，齐次性，三角不等式性质，则定义极限之后，就可将Euclid空间中那些仅依赖于性质的概念与结论直接推广于该空间，从而得到经典分析的一个具有广阔发展余地的拓广，从而引向赋范空间概念。

惯性导航作业

惯性导航作业

一、数据说明： 1：惯导系统为指北方位的捷连系统。初始经度为116.344695283度、纬度为39.975172度,高度h为30米。初速度 v0=[-9.993908270;0.000000000;0.348994967]。 2：jlfw中为600秒的数据，陀螺仪和加速度计采样周期分别为为1/100秒和1/100秒。 3：初始姿态角为[2 1 90]（俯仰，横滚，航向,单位为度）,jlfw.mat中保存的为比力信息f_INSc(单位m/s^2)、陀螺仪角速率信息wib_INSc(单位rad/s)，排列顺序为一~三行分别为X、Y、Z向信息. 4: 航向角以逆时针为正。 5:地球椭球长半径re=6378245;地球自转角速度wie=7.292115147e-5;重力加速度g=g0*(1+gk1*c33^2)*(1-2*h/re)/sqrt(1-gk2*c33^2)； g0=9.7803267714;gk1=0.00193185138639;gk2=0.00669437999013;c33=sin(lat纬度); 二、作业要求： 1：可使用MATLAB语言编程，用MATLAB编程时可使用如下形式的语句读取数据：load D:\...文件路径...\jlfw,便可得到比力信息和陀螺仪角速率信息。用角增量法。 2：(1) 以系统经度为横轴，纬度为纵轴（单位均要转换为：度）做出系统位置曲线图； (2) 做出系统东向速度和北向速度随时间变化曲线图（速度单位：m/s，时间单位：s）； (3) 分别做出系统姿态角随时间变化曲线图（俯仰，横滚，航向,单位转换为：度，时间单位：s）； 以上结果均要附在作业报告中。 3：在作业报告中要写出“程序流程图、现阶段学习小结”，写明联系方式。

北航数值分析大作业一

《数值分析B》大作业一 SY1103120 朱舜杰 一．算法设计方案： 1.矩阵A的存储与检索 将带状线性矩阵A[501][501]转存为一个矩阵MatrixC[5][501] . 由于C语言中数组角标都是从0开始的，所以在数组MatrixC[5][501]中检索A的带内元素a ij的方法是： A的带内元素a ij=C中的元素c i-j+2,j 2.求解λ1，λ501，λs ①首先分别使用幂法和反幂法迭代求出矩阵按摸最大和最小的特征值λmax和λmin。λmin即为λs； 如果λmax>0,则λ501=λmax；如果λmax<0,则λ1=λmax。 ②使用带原点平移的幂法（mifa（）函数），令平移量p=λmax，求 出对应的按摸最大的特征值λ，max， 如果λmax>0,则λ1=λ，max+p；如果λmax<0,则λ501=λ，max+p。 3.求解A的与数μk=λ1+k（λ501-λ1）/40的最接近的特征值λik （k=1,2，…，39）。 使用带原点平移的反幂法，令平移量p=μk，即可求出与μk最接近的特征值λik。 4.求解A的（谱范数）条件数cond（A）2和行列式d etA。 ①cond（A）2=|λ1/λn|，其中λ1和λn分别是矩阵A的模最大和 最小特征值。

②矩阵A的行列式可先对矩阵A进行LU分解后，detA等于U所有对角线上元素的乘积。 二．源程序 #include #include #include #include #include #include #include #define E 1.0e-12 /*定义全局变量相对误差限*/ int max2(int a,int b) /*求两个整型数最大值的子程序*/ { if(a>b) return a; else return b; } int min2(int a,int b) /*求两个整型数最小值的子程序*/ { if(a>b) return b; else return a; } int max3(int a,int b,int c) /*求三整型数最大值的子程序*/ { int t; if(a>b) t=a; else t=b; if(t

机械原理大作业

机械原理大作业 This model paper was revised by the Standardization Office on December 10, 2020

机械原理大作业三 课程名称：机械原理 设计题目：齿轮传动设计 院系： 班级： 设计者： 学号： 指导教师： 设计时间： 1、设计题目 机构运动简图 机械传动系统原始参数

2、传动比的分配计算 电动机转速min /745r n =，输出转速m in /1201r n =，min /1702r n =， min /2303r n ，带传动的最大传动比5.2max =p i ,滑移齿轮传动的最大传动比4m ax =v i ，定轴齿轮传动的最大传动比4m ax =d i 。 根据传动系统的原始参数可知，传动系统的总传动比为： 传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三部分实现。设带传动的传动比为5.2max =p i ，滑移齿轮的传动比为321v v v i i i 、、，定轴齿轮传动的传动比为f i ，则总传动比 令 4max 1==v v i i 则可得定轴齿轮传动部分的传动比为 滑移齿轮传动的传动比为 设定轴齿轮传动由3对齿轮传动组成，则每对齿轮的传动比为 3、齿轮齿数的确定 根据滑移齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮5、6、7、8、9和10为角度变位齿轮，其齿数： 35,18,39,14,43,111098765======z z z z z z ；它们的齿顶高系数1=* a h ，径向间 隙系数25.0=*c ，分度圆压力角020=α，实际中心距mm a 51'=。

北航自动控制元件复习提纲

自控元件复习提纲 一、关于考试（从学长处获得信息，仅供大家参考）： 1、简答：eg：1）、为什么空载实验可以测r m ，x m？ 2）、变压器中是否必须有无功功率？ 2、论述：eg：1）、变压器工作原理。 2）、电机基本特点。 论述时必要时要作图说明。如要论述两相伺服电机为何有无自转特性时应该作出单相工作时两相伺服电机的机械特性曲线来说明。 3、计算：两个计算，一个直流，一个交流。 二、各章重点内容概述（参考往届笔记）： 第1章直流磁路及其计算 重点章节：1-1，1-2，1-3 磁路总是闭合的；磁路计算的正反两类任务；等效磁路的画法第2章直流电磁铁及其典型应用 重点章节：2-1，2-2，另外第三节中的继电器的主要技术指标大家也要注意第3章直流电机的一般问题 重点章节：3-1，3-2，3-4，3-5，3-6 发电机、电动机的识别，电角的概念，电枢反应的概念，直流电机的电枢电动 势和电磁转矩的计算方法，直流电机的电势平衡关系、转矩平衡关系、功率平 衡关系，电枢绕组的具体原理如果不明白大家不必深究，只要知道概念就可以， 这里不是重点。 PS：本章是比较重要的一章，计算题可能出在这一章。 第4章直流测速发电机和直流伺服电动机 重点章节：4-1，4-2 第三节要掌握直流力矩电动机的特点，知道其应用场合以 及为什么要用在这些场合。 准确理解直流伺服电动机的工作原理、四种工作状态，准确掌握直流测速发电 机和直流伺服电动机的工作原理及特性（输入、输出等） 第5章变压器 重点章节：5-2，5-3，5-4，5-5，5-6 空载电压平衡式和等效电路及相量图，负载运行的电压平衡式和磁动势平衡式 及相量图和等效电路P101的图5-12大家仔细看。参数、额定数据和特性。磁 场问题转化到电路问题。 第6章异步电动机 重点章节：6-1，6-3，6-4，6-5 转差率S，功率传递，相电动势，第三节可能考计算，注意一下习题中的最后 两道计算题6.12，6.15 另外，第六章和第五章联系很紧密，大家可以结合起来复习效果更好。 第7章两相电机 重点章节：7-1，7-2 分解磁场，无自转的分析，稳定运行范围大，第一节的第（四）部分其它大家 可以不用管。异步测速发电机的输出电压的特点，原理，频率与转速无关，只 与电源有关。 第8章同步电动机

北航惯性导航综合实验五实验报告

惯性导航技术综合实验 实验五惯性基组合导航及应用技术实验

惯性/卫星组合导航系统车载实验 一、实验目的 ①掌握捷联惯导/GPS组合导航系统的构成和基本工作原理； ②掌握采用卡尔曼滤波方法进行捷联惯导/GPS组合的基本原理； ③掌握捷联惯导 /GPS组合导航系统静态性能； ④掌握动态情况下捷联惯导 /GPS组合导航系统的性能。 二、实验内容 ①复习卡尔曼滤波的基本原理（参考《卡尔曼滤波与组合导航原理》第二、五章）； ②复习捷联惯导/GPS组合导航系统的基本工作原理（参考以光衢编著的《惯性导航原理》第七章）； 三、实验系统组成 ①捷联惯导/GPS组合导航实验系统一套； ②监控计算机一台。 ③差分 GPS接收机一套； ④实验车一辆； ⑤车载大理石平台； ⑥车载电源系统。 四、实验内容 1）实验准备 ①将IMU紧固在车载大理石减振平台上，确认IMU的安装基准面紧靠实验平台； ②将IMU与导航计算机、导航计算机与车载电源、导航计算机与监控计算

机、GPS 接收机与导航计算机、GPS 天线与GPS 接收机、GPS 接收机与GPS 电池之间的连接线正确连接； ③ 打开GPS 接收机电源，确认可以接收到4颗以上卫星； ④ 打开电源，启动实验系统。 2） 捷联惯导/GPS 组合导航实验 ① 进入捷联惯导初始对准状态，记录IMU 的原始输出，注意5分钟内严禁移动实验车和IMU ； ② 实验系统经过5分钟初始对准之后，进入导航状态； ③ 移动实验车，按设计实验路线行驶； ④ 利用监控计算机中的导航软件进行导航解算，并显示导航结果。 五、 实验结果及分析 (一) 理论推导捷联惯导短时段（1分钟）位置误差，并用1分钟惯导实验数据验证。 1、一分钟惯导位置误差理论推导： 短时段内（t<5min ），忽略地球自转0ie ω=，运动轨迹近似为平面1/0R =，此时的位置误差分析可简化为： （1） 加速度计零偏?引起的位置误差：2 10.88022t x δ?==m （2） 失准角0φ引起的误差：2 02 0.92182g t x φδ==m （3） 陀螺漂移ε引起的误差：3 30.01376 g t x εδ==m 可得1min 后的位置误差值123 1.8157m x x x x δδδδ=++= 2、一分钟惯导实验数据验证结果： （1）纯惯导解算1min 的位置及位置误差图：

专业点题北航机械原理

一、齿轮传动的基本概念 渐开线齿轮的啮合特点：(1)渐开线齿廓能够保证定传动比;(2)渐开线齿廓之间的正压力方向不变;(3)渐开线齿廓传动具有可分性。 齿轮机构的特点是：传动平稳、适用范围广、效率高、结构紧凑、工作可靠、寿命长。但制造和安装精度高、制造费用大，且不适合于距离较远的两轴之间的传动。齿轮传动可以用来传递任意轴间的运动和动力。 齿轮传动按照一对齿轮传递的相对运动分为平面齿轮传动和空间齿轮传动，平面齿轮传动又分为直齿圆柱齿轮传动、斜齿圆柱齿轮传动和人字齿轮传动;按照工作条件可以分为开式传动、半开式传动和闭式传动。 齿轮传动的基本要求是：传动准确、平稳;承载能力强。 二、齿轮传动的设计与计算 齿廓曲线与齿廓啮合基本定律：在啮合传动的任一瞬时，两轮齿廓曲线在相应接触点的功法线必须通过按给定传动比确定的该瞬时的节点。 渐开线齿轮啮合的正确条件：啮合轮齿的工作侧齿廓的啮合点必须总是在啮合线上，即两齿轮的模数和压力角应该分别相等。 齿轮传动的无侧隙啮合及标准齿轮的安装：一个齿轮节圆上的齿厚等于另一个齿轮节圆上的齿槽宽是无侧隙啮合的条件;外啮合齿轮的标准中心距为，内啮合是标准中心距为。

齿轮及其变位的相关计算：相关参数为齿数、模数、分度圆压力角、齿顶高系数和顶隙系数及标准直齿轮的几何尺寸计算，包括分度圆直径、齿顶高、齿根高、齿全高、齿顶圆直径、齿根圆直径、基圆直径、齿距、齿厚、齿槽宽、中心距、顶隙以及变位齿轮的变位系数等。 渐开线齿轮的根切现象：用展成法加工齿轮式，若刀具的齿顶线或齿顶圆与啮合线的焦点超过被切齿轮的极限点，则刀具的齿顶会将被切齿轮的齿根的渐开线齿廓切去了一部 分。避免根切的最小齿数，用标准齿条刀具切制标准齿轮时，因为 ，最少齿数为17。 三、机构的组成 构件指独立的运动单元，两个构件直接接触组成仍能产生某些相对运动的连接叫运动副。运动副按照相对运动的范围可以分为平面运动副和空间运动副;按运动副元素分为：低副-面接触、应力低;高副-点接触或线接触，应力高。其中运动副元素是只形成运动副的组建之间直接接触的部分。 四、机构自由度的计算 机构相对于机架所具有的独立运动的数目，叫机构的自由度。设一个平面机构由N个构件组成，其中必定有一个构件为机架，其活动构件数为n=N-1.设机构共有个低副、 个高副，因为在平面机构中每个低副和高副分别限制两个自由度和一个自由度，故平面机构的自由度为。在计算平面机构的自由度时，应该注意三种特殊情况：(1)复合铰链：三个或更多的构件在同一处联接成同轴线的两个或更多个转动副，就构成了复合铰链，计算自由度时应该按照两个或更多个运动副计算。(2)局部自由度：在有些机构中，为了其他一些非运动的原因，设置了附加机构，这种附加机构的运动是完全独立的，对整个

电气工程专业考研专业课初试科目及复试内容汇总

《电气工程专业考研专业课初试科目及复试内容汇总》 自动化专业的考研方向 自动化专业方向很广，考的时候还分双控，模式，电力电机等等方面，你可以参看学校是否在这个方面有无国家重点实验室，是不是国家重点学科来比较。 1. 清华， 2.中科院， 3.上海交大， 4.浙大，5华工，北航，东南，东北大学，西安交大，哈 尔滨工业大学，中国科技大学，华北电力，天津大学，东南大学，华中科技，武汉大学天津大学自动化 一般说来，初试的分数是最重要的，特别是考外校。当然，你的动手能力也是很重要的，还有你的英语口语，考研复试都是要考虑的。例如上海交大的复试，双控353的复试线，有380的被刷下来，就是英语口语已经专业课不是很扎实的。考外校的话依据学校而定是否要找导师 动手能力强，参加电子设计大赛都是作为你考研复试的参考，还是好好的准备初试的考试吧，毕竟它是个门槛。 【电气工程及其自动化】 北京工业大学 421自动控制原理 复试：1、电子技术2、计算机原理 北京航空航天大学 [双控] 432控制理论综合或433控制工程综合 [检测] 433控制工程综合或436检测技术综合 [系统] 431自动控制原理或451材料力学或841概率与数理统计 [模式] （自动化学院）433控制工程综合或436检测技术综合、（宇航学院）423信息类专业综合或431自动控制原理或461计算机专业综合 [导航] （自动化学院）432控制理论综合或433控制工程综合、（宇航学院）431自动控制原理 复试：无笔试。1) 外语口语与听力考核；2) 专业基础理论与知识考核；3) 大学阶段学习成绩、科研活动以及工作业绩考核；4) 综合素质与能力考核 北京化工大学 440电路原理 复试：综合1(含自动控制原理和过程控制系统及工程)、综合2(含自动检测技术装置和传感器原理及应用)、综合3(含信号与系统和数字信号处理) 注：数学可选择301数学一或666数学（单） 北京交通大学 [双控/检测]404控制理论 [模式]405通信系统原理或409数字信号处理 复试： [电子信息工程学院双控]常微分方程 [机械与电子控制工程学院检测]综合复试（单片机、自动控制原理） [计算机与信息技术学院模式] 信号与系统或操作系统

北航惯性导航综合实验四实验报告

基于运动规划的惯性导航系统动态实验 二零一三年六月十日

实验4.1 惯性导航系统运动轨迹规划与设计实验 一、实验目的 为进行动态下简化惯性导航算法的实验研究，进行路径和运动状态规划，以验证不同运动状态下惯导系统的性能。通过实验掌握步进电机控制方法，并产生不同运动路径和运动状态。 二、实验内容 学习利用6045B 控制板对步进电机进行控制的方法，并控制电机使运动滑轨产生定长运动和不同加速度下的定长运动。 三、实验系统组成 USB_PCL6045B 控制板（评估板）、运动滑轨和控制计算机组成。 四、实验原理 IMU安装误差系数的计算方法 USB_PCL6045B 控制板采用了USB 串行总线接口通信方式，不必拆卸计算机箱就可以在台式机或笔记本电脑上进行运动控制芯片PCL6045B 的学习和评估。 USB_PCL6045B 评估板采用USB 串行总线方式实现评估板同计算机的数据交换，由评估板的FIFO 控制回路完成步进电机以及伺服电机的高速脉冲控制，任意2 轴的圆弧插补，2-4 轴的直线插补等运动控制功能。USB_PCL6045B 评估板上配置了全部PCL6045B 芯片的外部信号接口和增量编码器信号输入接口。由 USB_PCL6045B 评估测试软件可以进行PCL6045B 芯片的主要功能的评估测试。

图4-1-1USB_PCL6045B 评估板原理框图 如图4-1-1 所示，CN11 接口主要用于外部电源连接，可以选择DC5V 单一电源或DC5V/24V 电源。CN12 接口是USB 信号接口，用于USB_PCL6045B 评估板同计算机的数据交换。 USB_PCL6045B 评估板已经完成对PCL6045B 芯片的底层程序开发和硬件资源与端口的驱动，并封装成156 个API 接口函数。用户可直接在VC 环境下利用API 接口函数进行编程。 五、实验内容 1、操作步骤 1）检查电机驱动电源（24V） 2）检查USB_PCL6045B 控制板与上位机及电机驱动器间的连接电缆 3）启动USB_PCL6045B 控制板评估测试系统检查系统是否正常工作。 4）运行编写的定长运动程序，并比较实际位移与设定位移。

北航数值分析大作业第一题幂法与反幂法

《数值分析》计算实习题目 第一题： 1. 算法设计方案 （1）1λ，501λ和s λ的值。 1)首先通过幂法求出按模最大的特征值λt1，然后根据λt1进行原点平移求出另一特征值λt2，比较两值大小，数值小的为所求最小特征值λ1，数值大的为是所求最大特征值λ501。 2)使用反幂法求λs ，其中需要解线性方程组。因为A 为带状线性方程组，此处采用LU 分解法解带状方程组。 （2）与140k λλμλ-5011=+k 最接近的特征值λik 。 通过带有原点平移的反幂法求出与数k μ最接近的特征值 λik 。 （3）2cond(A)和det A 。 1）1=n λλ2cond(A)，其中1λ和n λ分别是按模最大和最小特征值。 2）利用步骤（1）中分解矩阵A 得出的LU 矩阵，L 为单位下三角阵,U 为上三角阵，其中U 矩阵的主对角线元素之积即为det A 。 由于A 的元素零元素较多，为节省储存量，将A 的元素存为6×501的数组中，程序中采用get_an_element()函数来从小数组中取出A 中的元素。 2.全部源程序 #include #include void init_a();//初始化A double get_an_element(int,int);//取A 中的元素函数 double powermethod(double);//原点平移的幂法 double inversepowermethod(double);//原点平移的反幂法 int presolve(double);//三角LU 分解 int solve(double [],double []);//解方程组 int max(int,int); int min(int,int); double (*u)[502]=new double[502][502];//上三角U 数组 double (*l)[502]=new double[502][502];//单位下三角L 数组 double a[6][502];//矩阵A int main() { int i,k; double lambdat1,lambdat2,lambda1,lambda501,lambdas,mu[40],det;

哈工大机械原理大作业凸轮机构第四题

Harbin Institute of Technology 机械原理大作业二 课程名称：机械原理 设计题目：凸轮机构设计 姓名：李清蔚 学号：1140810304 班级：1408103 指导教师：林琳

一.设计题目 设计直动从动件盘形凸轮机构，其原始参数见表 1 表一：凸轮机构原始参数 升程（mm ) 升程 运动 角（o） 升程 运动 规律 升程 许用 压力 角（o） 回程 运动 角（o） 回程 运动 规律 回程 许用 压力 角（o） 远休 止角 （o） 近休 止角 （o） 40 90 等加 等减 速30 50 4-5-6- 7多 项式 60 100 120

二.凸轮推杆运动规律 （1）推程运动规律（等加速等减速运动） 推程F0=90° ①位移方程如下： ②速度方程如下： ③加速度方程如下： （2）回程运动规律（4-5-6-7多项式） 回程，F0=90°，F s=100°，F0’=50°其中回程过程的位移方程，速度方程，加速度方程如下：

三.运动线图及凸轮线图 本题目采用Matlab编程，写出凸轮每一段的运动方程，运用Matlab模拟将凸轮的运动曲线以及凸轮形状表现出来。代码见报告的结尾。 1、程序流程框图 开始 输入凸轮推程回 程的运动方程 输入凸轮基圆偏 距等基本参数 输出ds,dv,da图像 输出压力角、曲率半径图像 输出凸轮的构件形状 结束

2、运动规律ds图像如下： 速度规律dv图像如下： 加速度da规律如下图：

3.凸轮的基圆半径和偏距 以ds/dfψ-s图为基础，可分别作出三条限制线（推程许用压力角的切界限D t d t,回程许用压力角的限制线D t'd t'，起始点压力角许用线B0d'')，以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件。 得图如下：得最小基圆对应的坐标位置O点坐标大约为（13，-50）经计算取偏距e=13mm，r0=51.67mm.

北航惯性导航大作业

惯性导航基础课程大作业报告（一）光纤陀螺误差建模与分析 班级：111514 姓名： 学号 2014年5月26日

一.系统误差原理图 二.系统误差的分析 （一）漂移引起的系统误差 1. εx ，εy ，εz 对东向速度误差δVx 的影响 clc;clear all; t=1:0.01:25; g=9.8; L=pi/180*39; Ws=2*pi/84.4*60; Wie=2*pi/24; R=g/(Ws)^2; e=0.1*180/pi; mcVx1=e*g*sin(L)/(Ws^2-Wie^2)*(sin(Wie*t)-Wie*sin(Ws*t)/Ws); mcVx2=e*((Ws^2-(Wie^2)*((cos(L))^2))/(Ws^2-Wie^2)*cos(Ws*t)-(Ws^2)*((sin(L))^2)*cos(Wi e*t)/(Ws^2-Wie^2)-(cos(L))^2); mcVx3=(sin(L))*(cos(L))*R*e*((Ws^2)*cos(Wie*t)/(Ws^2-Wie^2)-(Wie^2)*cos(Ws*t)/(Ws^2-Wi e^2)-1); plot(t,[mcVx1',mcVx2',mcVx3']); title('Ex,Ey,Ez 对Vx 的影响'); xlabel('时间t'); ylabel('Vx(t)'); 0,δλδL ,v v δδ

legend('Ex-mcVx1','Ey-mcVx2','Ez-mcVx3'); grid; axis square; 分析：εx，εy，εz对东向速度误差δVx均有地球自转周期的影响，εx，εy还会有舒勒周期分量的影响，其中，εy对δVx的影响较大。 2.εx，εy，εz对东向速度误差δVy的影响 clc;clear all; t=1:0.01:25; g=9.8; L=pi/180*39; Ws=2*pi/84.4*60; Wie=2*pi/24; R=g/(Ws)^2; e=0.1*180/pi; mcVy1=e*g*(cos(Wie*t)-cos(Ws*t))/(Ws^2-Wie^2); mcVy2=g*sin(L)*e/(Ws^2-Wie^2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t)); mcVy3=g*cos(L)*e/(Ws^2-Wie^2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t)); plot(t,[mcVy1',mcVy2',mcVy3']); title('Ex,Ey,Ez对Vy的影响'); xlabel('时间t'); ylabel('Vy(t)'); legend('Ex-mcVy1','Ey-mcVy2','Ez-mcVy3'); grid; axis square;

北航数值分析报告第三次大作业

数值分析第三次大作业 一、算法的设计方案： （一）、总体方案设计： x y当作已知量代入题目给定的非线性方程组，求（1）解非线性方程组。将给定的(,) i i

得与(,)i i x y 相对应的数组t[i][j],u[i][j]。 （2）分片二次代数插值。通过分片二次代数插值运算，得到与数组t[11][21],u[11][21]]对应的数组z[11][21]，得到二元函数z=(,)i i f x y 。 （3）曲面拟合。利用x[i],y[j],z[11][21]建立二维函数表，再根据精度的要求选择适当k 值，并得到曲面拟合的系数矩阵C[r][s]。 （4）观察和(,)i i p x y 的逼近效果。观察逼近效果只需要重复上面（1）和（2）的过程，得到与新的插值节点(,)i i x y 对应的(,)i i f x y ，再与对应的(,)i i p x y 比较即可，这里求解 (,)i i p x y 可以直接使用（3）中的C[r][s]和k 。 （二）具体算法设计： （1）解非线性方程组 牛顿法解方程组()0F x =的解* x ，可采用如下算法： 1）在* x 附近选取(0) x D ∈，给定精度水平0ε>和最大迭代次数M 。 2）对于0,1, k M =执行 ① 计算() ()k F x 和()()k F x '。 ② 求解关于() k x ?的线性方程组 () ()()()()k k k F x x F x '?=- ③ 若() () k k x x ε∞∞ ?≤，则取*()k x x ≈，并停止计算；否则转④。 ④ 计算(1) ()()k k k x x x +=+?。 ⑤ 若k M <，则继续，否则，输出M 次迭代不成功的信息，并停止计算。 （2）分片双二次插值 给定已知数表以及需要插值的节点，进行分片二次插值的算法： 设已知数表中的点为： 00(0,1,,) (0,1,,)i j x x ih i n y y j j m τ=+=???=+=?? ，需要插值的节点为(,)x y 。 1) 根据(,)x y 选择插值节点(,)i j x y ： 若12h x x ≤+ 或12 n h x x ->-，插值节点对应取1i =或1i n =-，

北航惯性导航综合实验四实验报告

基于运动规划的惯性导航系统动态实验 GAGGAGAGGAFFFFAFAF

二零一三年六月十日 实验4.1 惯性导航系统运动轨迹规划与设计实验一、实验目的 为进行动态下简化惯性导航算法的实验研究，进行路径和运动状态规划，以验证不同运动状态下惯导系统的性能。通过实验掌握步进电机控制方法，并产生不同运动路径和运动状态。 二、实验内容 学习利用6045B 控制板对步进电机进行控制的方法，并控制电机使运动滑轨产生定长运动和不同加速度下的定长运动。 三、实验系统组成 USB_PCL6045B 控制板（评估板）、运动滑轨和控制计算机组成。 四、实验原理 IMU安装误差系数的计算方法 GAGGAGAGGAFFFFAFAF

USB_PCL6045B 控制板采用了USB 串行总线接口通信方式，不必拆卸计算机箱就可以在台式机或笔记本电脑上进行运动控制芯片PCL6045B 的学习和评估。 USB_PCL6045B 评估板采用USB 串行总线方式实现评估板同计算机的数据交换，由评估板的FIFO 控制回路完成步进电机以及伺服电机的高速脉冲控制，任意 2 轴的圆弧插补，2-4 轴的直线插补等运动控制功能。USB_PCL6045B 评估板上配置了全部PCL6045B 芯片的外部信号接口和增量编码器信号输入接口。由 USB_PCL6045B 评估测试软件可以进行PCL6045B 芯片的主要功能的评估测试。 GAGGAGAGGAFFFFAFAF

图4-1-1USB_PCL6045B 评估板原理框图如图4-1-1 所示，CN11 接口主要用于外部电源连接，可以选择DC5V 单一电源或DC5V/24V 电源。CN12 接口是USB 信号接口，用于USB_PCL6045B 评估板同计算机的数据交换。 USB_PCL6045B 评估板已经完成对PCL6045B 芯片的底层程序开发和硬件资源与端口的驱动，并封装成156 个API 接口函数。用户可直接在VC 环境下利用API 接口函数进行编程。 五、实验内容 GAGGAGAGGAFFFFAFAF

北航惯性导航综合实验一实验报告

实 验一 陀螺仪关键参数测试与分析实验 加速度计关键参数测试与分析实验 二零一三年五月十二日 实验一陀螺仪关键参数测试与分析实验 一、实验目得 通过在速率转台上得测试实验,增强动手能力与对惯性测试设备得感性认识；通过对陀螺仪测试数据得分析,对陀螺漂移等参数得物理意义有清晰得认识,同时为在实际工程中应用陀螺仪与对陀螺仪进行误差建模与补偿奠定基础。 二、实验内容 利用单轴速率转台,进行陀螺仪标度因数测试、零偏测试、零偏重复性测试、零漂测试实验与陀螺仪标度因数与零偏建模、误差补偿实验。 三、实验系统组成 单轴速率转台、MEMS 陀螺仪(或光纤陀螺仪)、稳压电源、数据采集系统与分析系统。

四、实验原理 1.陀螺仪原理 陀螺仪就是角速率传感器,用来测量载体相对惯性空间得角速度，通常输出与角速率对应得电压信号。也有得陀螺输出频率信号(如激光陀螺）与数字信号（把模拟电压数字化)。以电压表示得陀螺输出信号可表示为: (1-1）式中就是与比力有关得陀螺输出误差项，反映了陀螺输出受比力得影响,本实验不考虑此项误差。因此，式（１－１）简化为 (１-2）由（１－2）式得陀螺输出值所对应得角速度测量值: (1－3) 对于数字输出得陀螺仪,传感器内部已经利用标度因数对陀螺仪模拟输出进行了量化,直接输出角速度值,即： （１－4）就是就是陀螺仪得零偏，物理意义就是输入角速度为零时，陀螺仪输出值所对应得角速度。且 (１-5) 精度受陀螺仪标度因数、随机漂移、陀螺输出信号得检测精度与得影响。通常与表现为有规律性,可通过建模与补偿方法消除,表现为随机特性,可通过信号滤波方法抵制。因此，准确标定与就是实现角速度准确测量得基础。 五、陀螺仪测试实验步骤 1)标度因数与零偏测试实验 ａ、接通电源，预热一定时间; b、陀螺工作稳定后，测量静止情况下陀螺输出并保存数据;

北航数值分析大作业第二题精解

目标：使用带双步位移的QR 分解法求矩阵10*10[]ij A a =的全部特征值，并对其中的每一个实特征值求相应的特征向量。已知：sin(0.50.2)() 1.5cos( 1.2)(){i j i j ij i j i j a +≠+== (i,j=1,2, (10) 算法： 以上是程序运作的逻辑，其中具体的函数的算法，大部分都是数值分析课本上的逻辑，在这里特别写出矩阵A 的实特征值对应的一个特征向量的求法： ()[]()() []()[]()111111I 00000 i n n n B A I gause i n Q A I u Bu u λλ-?-?-=-?-?? ?-=????→=??????→= ?? ? 选主元的消元 检查知无重特征值 由于=0i A I λ- ，因此在经过选主元的高斯消元以后，i A I λ- 即B 的最后一行必然为零，左上方变 为n-1阶单位矩阵[]()()11I n n -?-，右上方变为n-1阶向量[]()11n Q ?-，然后令n u 1=-，则 ()1,2,,1j j u Q j n ==???-。

这样即求出所有A所有实特征值对应的一个特征向量。 #include #include #include #define N 10 #define E 1.0e-12 #define MAX 10000 //以下是符号函数 double sgn(double a) { double z; if(a>E) z=1; else z=-1; return z; } //以下是矩阵的拟三角分解 void nishangsanjiaodiv(double A[N][N]) { int i,j,k; int m=0; double d,c,h,t; double u[N],p[N],q[N],w[N]; for(i=0;i

北航七系机械学院机械原理大作业

机械原理课程机构设计 实验报告 题目：建筑垃圾破碎机的设计与分析小组成员与学号： 班级： 第1页

建筑垃圾破碎机的设计与分析 摘要 本文简单介绍了建筑垃圾回收再利用的重要性，与工艺性，并自主设计了将颚式破碎机与反击式破碎机相结合的建筑垃圾破碎机。通过solidworks软件对设计机构进行建模，用adams进行仿真分析，验证所设计的机构均达到设计需要与可行性。 关键词：建筑垃圾破碎机、连杆机构、凸轮廓线设计 第2页

目录 1．机构的引出 (4) 1.1 建筑垃圾及其回收利用价值 (4) 1.2颚式破碎机和反击式破碎机各自的利弊分析 (4) 1.3设计新的建筑垃圾破碎机 (6) 2．机构的结构、功能介绍及建模 (7) 2.1 机构设计简图及各部分功能 (7) 2.2尺寸设计及建模 (8) 2.2.1主动轮和各从动轮的传动比 (8) 2.2.2凸轮廓线设计与挡板行程 ................................... 错误！未定义书签。 3．机构的仿真分析 (12) 3.1颚式破碎机的急回特性 (12) 3.2颚式破碎机的传动角验证 (14) 3.3停歇运动导杆机构所带动的下挡板往复运动的间歇性 (14) 4．总结 (17) 第3页

第4页 1． 机构的引出 1.1 建筑垃圾及其回收利用价值 二十一世纪是一个飞速发展的时代，随着城市人口的增加、新农村建设以及城市地铁的大规模扩建，建筑行业的新陈代谢全面加速，建筑垃圾的排放量也随之增加。然而，传统的方法处理建筑垃圾是将建筑垃圾运往乡村或郊外，露天堆放或掩埋。这样不仅破坏植被，降低土壤的生产能力，而且会让建筑垃圾中的有害物质渗入地下水层，污染环境，给人们的生活带来困扰。因此，如何实现建筑垃圾的高效、环保循环利用成为当今人们所面临的一个难题。 建筑垃圾的主要组成部分是废弃混凝土和砖块，而它们都是由水泥和天然砂石拌合而成的，这些都是砖块等建筑材料的重要组成部分。为了最大程度的利用建筑垃圾，首先应该解决的问题就是对其中的大块物料进行破碎，只有这样，破碎后的小快物料才能很好的还原天然砂石的性能，实现建筑垃圾的循环利用。 1.2颚式破碎机和反击式破碎机各自的利弊分析 目前应用较广的破碎机有颚式破碎机与反击式破碎机两种。 颚式破碎机的主体构造如图 1 图 1 颚式破碎机的主体构造 其工作原理为：轮①通过皮带和电机上的主动轮相连，①的转动带动杆②进而带动构件③的摆动（构件③的上端和机架铰接）。构件③通过摆动将体积较大

(完整)北航惯性导航作业二.

(完整)北航惯性导航作业二. 编辑整理： 尊敬的读者朋友们： 这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（(完整)北航惯性导航作业二.）的内容能够给您的工作和学习带来便利。同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。 本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为(完整)北航惯性导航作业二.的全部内容。

惯性导航作业

一、数据说明： 1：惯导系统为指北方位的捷连系统.初始经度为116。344695283度、纬度为 39.975172度,高度h为30米。初速度v0=［—9。993908270;0.000000000； 0.348994967］。 2：jlfw中为600秒的数据，陀螺仪和加速度计采样周期分别为为1/100秒和1/100秒。 3：初始姿态角为［2 1 90］（俯仰，横滚，航向,单位为度)，jlfw。mat中保存的 为比力信息f_INSc(单位m/s^2)、陀螺仪角速率信息wib_INSc(单位rad/s)，排列 顺序为一～三行分别为X、Y、Z向信息. 4：航向角以逆时针为正. 5：地球椭球长半径re=6378245;地球自转角速度wie=7。292115147e-5;重力加速度 g=g0*（1+gk1＊c33^2)＊（1-2*h/re)/sqrt(1—gk2*c33^2）;g0=9.7803267714； gk1=0。00193185138639；gk2=0。00669437999013；c33=sin（lat纬度）; 二、作业要求： 1:可使用 MATLAB语言编程，用MATLAB编程时可使用如下形式的语句读取数据： load D:\..。文件路径。。.\jlfw,便可得到比力信息和陀螺仪角速率信息。用角增 量法。 2:（1) 以系统经度为横轴,纬度为纵轴（单位均要转换为：度）做出系统位置曲线图； （2）做出系统东向速度和北向速度随时间变化曲线图（速度单位:m/s，时间单位:s)； (3) 分别做出系统姿态角随时间变化曲线图（俯仰，横滚,航向,单位转换为：度，时间单位：s)； 以上结果均要附在作业报告中. 3：在作业报告中要写出“程序流程图、现阶段学习小结”,写明联系方式。 （注意程序流程图不是课本上的惯导解算流程，而是你程序分为哪几个模块、是 怎样一步步执行的,什么位置循环等，让别人根据该流程图能够编出相应程序）

惯性导航作业

惯性导航系统的基本原理、特点及在现代生活中的应用 惯性导航系统的基本原理 惯性导航系统也称作惯性参考系统，是一种不依赖于外部信息、也不向外部辐射能量(如无线电导航那样)的自主式导航系统。其工作环境不仅包括空中、地面，还可以在水下。惯性导航的基本工作原理是以牛顿力学定律为基础，通过测量载体在惯性参考系的加速度，将它对时间进行积分，且把它变换到导航坐标系中，就能够得到在导航坐标系中的速度、偏航角和位置等信息。 惯性导航系统属于推算导航方式，即从一已知点的位置根据连续测得的运动体航向角和速度推算出其下一点的位置，因而可连续测出运动体的当前位置。惯性导航系统中的陀螺仪用来形成一个导航坐标系，使加速度计的测量轴稳定在该坐标系中，并给出航向和姿态角;加速度计用来测量运动体的加速度，经过对时间的一次积分得到速度，速度再经过对时间的一次积分即可得到距离。 惯性导航技术的理论基础是牛顿力学基本定律。惯性导航系统是以陀螺和加速度计为敏感器件的导航参数解算系统，该系统根据陀螺的输出建立导航坐标系，根据加速度计输出解算出运载体在导航坐标系中的速度和位置。惯性导航系统分成平台式惯性导航系统和捷联式惯性导航系统两大类。平台式惯性导航系统将惯性测量元件安装在惯性平台上，惯性平台稳定在预定的坐标系内，为加速度计提供一个测量基准，并使惯性测量元件体角运动的影响。导航计算机根据加速度计的输出和初始条件进行导航解算，得出载的位置、速度等导航参数。捷联式惯性导航系统将惯性测量元件直接固联在载体上，测量沿载体坐标系的角速度和角加速度，计算机则利用陀螺的输出，进行坐标变换，求解载体的即时速度、位置等导航参数。惯性导航仅依靠惯性装置本身就能在载体内部独立地完成导航任务，不需要与外界发生任何信号联系，具有高度的自主性。这在战略和战术应用上具有重要的意义。但惯性导航的定位误差会随时间逐步增加，必须不断地进行误差修正，才能保证达到要求的精度。 陀螺仪 陀螺仪通常是指安装在万向支架中高度旋转的转子，转子同时可绕垂直于自转轴的一根轴或两根轴进动，前者称单自由度陀螺仪，后者称二自由度陀螺仪。

北航数值分析大作业第二题

数值分析第二次大作业 史立峰 SY1505327

一、 方案 （1）利用循环结构将sin(0.50.2)() 1.5cos( 1.2)() {i j i j ij i j i j a +≠+==(i,j=1,2,……,10)进行赋值，得到需要变换的 矩阵A ； （2）然后，对矩阵A 利用Householder 矩阵进行相似变换，把A 化为上三角矩阵A (n-1)。 对A 拟上三角化，得到拟上三角矩阵A (n-1)，具体算法如下： 记A(1)=A ，并记A(r)的第r 列至第n 列的元素为()n r r j n i a r ij ,,1,;,,2,1) ( +==。 对于2,,2,1-=n r 执行 1. 若 ()n r r i a r ir ,,3,2) ( ++=全为零，则令A(r+1) =A(r),转5；否则转2。 2. 计算 () ∑+== n r i r ir r a d 1 2 )( ()( )r r r r r r r r r r d c a d a c ==-=++则取,0sgn ) (,1)(,1若 )(,12r r r r r r a c c h +-= 3. 令 () n T r nr r r r r r r r r R a a c a u ∈-=++) ()(,2)(,1,,,,0,,0 。 4. 计算 r r T r r h u A p /)(= r r r r h u A q /)(= r r T r r h u p t /= r r r r u t q -=ω T r r T r r r r p u u A A --=+ω)()1( 5. 继续。 （3）使用带双步位移的QR 方法计算矩阵A (n-1)的全部特征值，也是A 的全部特征值，具体算法如下： 1. 给定精度水平0>ε和迭代最大次数L 。 2. 记n n ij n a A A ?-==][) 1()1()1(，令n m k ==,1。