北航惯导第一次大作业
《惯性导航原理》第一次大作业
一、 原理分析
惯导系统为指北方位的平台系统，则利用比力方程以及陀螺提供的东、北、天三个比力数据，即可计算得到在每个数据采集点的平台即时速度，再通过经纬度的计算公式，就可以得到每个数据采集点平台的即时经纬度，以每个数据采集点为下一个采集点的起点，即可对速度和经纬度进行累计计算，从而得到平台在运动过程中任意时刻的速度和位置情况。 1.模型公式的推导
载体相对地球运动时，加速度计测得的比力表达式，称为比力方程，方程如下：
g
V V f ep ep ie
ep
-?++=)2(?? （1）
在指北方案中，平台模拟地理坐标系，将上式中平台坐标系用地理坐标系代入得：
t t t et
t
ie
t t g
V f V +?+-=)2(?
?
（2）
系统中测量的是比力分量，将上式写成分量形式
=-
+ （3）
又因为地球的自转角速率为：
(4)
地理坐标系相对于地球坐标系的角速率为：
= (5)
将（4）（5）两个式子带入（3）式，即可得到如下方程组：
（6）
2.速度计算
作业要求只考虑水平通道，因此只需要计算正东、正北两个方向的速度即可。理论上计算得到t x V 、t y V 后，再积分一次可得到速度值，即
??
??
?+=+=
??t t
y t y
t y t
t x t x t x V dt V V V dt V
V 0
00
但在本次计算过程中，三个方向的速度均是从零开始在各时间节点上的累加，并不是t 的函数，因此速度计算可以由以下方程组实现：
（7）
此方程组表示了从第i 个采集点到第（i+1）个采集点的速度递推公式。 方程中Rx 表示卯酉圈的曲率半径，Ry 表示子午圈的曲率半径，计算方法如下：
Rx=Re/(1-esin 2
L);
Ry=Re/(1+2e-3esin 2
L); (8)
由于平台在运动中纬度L 也在不断变化，因此，计算过程中应当追踪两个半径的
变化，正如方程组（7）所示。
另外方程组中g 表征平台所处纬度下的重力加速度：
g=g 0(1+0.0052884sin 2L-0.0000059sin 2L) (9) 为了尽可能减小计算的累积误差，计算过程中可以对g 的变化进行追踪。 3.经纬度计算
载体所在位置的地理纬度L ，经度λ可由下列方程求得：
???
???
?+=+=?
?t
xt
t
x
t
yt t
y
Ldt R V L dt R V L 0
sec λλ
与速度的计算相同，经纬度也不是t 的函数，可以由累加得到：
二、程序流程图
三、导航结果
1.系统位置坐标曲线图
2.系统东向速度随时间变化曲线图
3.系统北向速度随时间变化曲线图
4．系统纬度、经度、东向速度、北向速度的终点值
四、小结
程序运行结果显示本次平台的运动是一个以初始点位置为中心的往返运动，从东向和北向速度的变化情况来看也是如此。
由于时间问题，自己有两种想法没有验证，不知道正确与否。一是把陀螺仪角速率信息ω看做平台相对地球的角速率信息，从而可将方程组（3）直接建立成模型计算。二是在本题目中，天向速度是存在的，我的模型计算结果是-1.5554m/s，因此高度是有变化的，可以加入高度通道进行计算，不知道结果会差多少。
惯性导航是比较有难度的课程，本来觉得很有压力，但是俞老师可以一个台阶一个台阶引导我学习，并激发了我的学习兴趣，非常感谢。我非常喜欢本门课程，对授课方式也非常满意，再次谢谢俞老师！
五、源程序
clear all;
clc;
Vx(1)=0; %初始化变量
Vy(1)=0;
Vz(1)=0;
L(1)=40.162565402/180*pi; %将初始位置经纬度变换成弧度
J(1)=116.343692076/180*pi;
Wie=7.2921E-5; %给公式中的常数赋值
pi=3.141592654;
re=6378245;
e=1/298.3;
h=37.74319;
g0=9.78049;
load('E:\学习\惯性导航\2011秋《惯性导航原理》第一次作业\fw');
%读取文件中的数据
Fx=f(1,:); %提取正东方向比力数据并定义
Fy=f(2,:); %提取正北方向比力数据并定义
Fz=f(3,:); %提取天向比力数据并定义
Wx=w(1,:); %提取陀螺正东方向角速率并定义
Wy=w(2,:); %提取陀螺正北方向角速率并定义
Wz=w(3,:); %提取陀螺天向角速率并定义
i=1;
t=1E-2;
for i=1:120000 %共有120000个数据，利用for程序完成对数据
的循环处理
g=g0*(1+0.0052884*sin(L(i))^2-0.0000059*sin(2*L(i))^2);
%在纬度发生变化的情况下追踪重力加速度的变化，减小计算误差Rx(i)=re/(1-e*sin(L(i))^2); %卯酉圈主曲率半径
Ry(i)=re/(1+2*e-3*e*sin(L(i))^2); %子午圈主曲率半径
Vx(i+1)=(Fx(i)+(2*Wie*sin(L(i))+Vx(i)*tan(L(i))/Rx(i))*Vy(i)-(2*Wie*c os(L(i))+Vx(i)/Rx(i))*Vz(i))*t+Vx(i);%正东方向速度计算
Vy(i+1)=(Fy(i)-(2*Wie*sin(L(i))+Vx(i)*tan(L(i))/Rx(i))*Vx(i)+Vy(i)*Vz (i)/Ry(i))*t+Vy(i); %正北方向速度计算
Vz(i+1)=(Fz(i)+(2*Wie*cos(L(i)+Vx(i))/Rx(i))*Vx(i)+Vy(i)*Vy(i)/Ry(i)-g)*t+Vz(i); %天向速度计算
L(i+1)=t*Vy(i)/Ry(i)+L(i); %纬度计算
J(i+1)=t*Vx(i)/(Rx(i)*cos(L(i)))+J(i);%经度计算
i=i+1; %i与1累加，直至将120000组数据处理完成循环结束
end
figure(1)
plot(L*180/pi,J*180/pi);xlabel('经度'),ylabel('纬度'); %以经度为横轴，纬度为纵轴（单位为：度）作出系统位置坐标曲线图
T=0:0.01:1200;
figure(2)
plot(T,Vy);xlabel('时间/秒'),ylabel('东向速度/m/s'); %以时间为横轴（单位：秒），东向速度为纵轴作出系统速度随时间变化曲线图
figure(3)
plot(T,Vx);xlabel('时间/秒'),ylabel('北向速度/m/s'); %以时间为横轴（单位：秒），北向速度为纵轴作出系统速度随时间变化曲线图
Lzz=L(120001)*180/pi %输出最终位置的纬度
Jzz=J(120001)*180/pi %输出最终位置的经度
Vxzz=Vx(120001) %正东方向的速度最终值
Vyzz=Vy(120001) %正北方向的速度最终值
北航最优化方法大作业参考
北航最优化方法大作业参考
1 流量工程问题 1.1 问题重述 定义一个有向网络G=(N,E)，其中N是节点集，E是弧集。令A是网络G的点弧关联矩阵，即N×E阶矩阵，且第l列与弧里(I,j)对应，仅第i行元素为1，第j行元素为-1，其余元素为0。再令b m=(b m1,…,b mN)T，f m=(f m1,…,f mE)T，则可将等式约束表示成： Af m=b m 本算例为一经典TE算例。算例网络有7个节点和13条弧，每条弧的容量是5个单位。此外有四个需求量均为4个单位的源一目的对，具体的源节点、目的节点信息如图所示。这里为了简单，省区了未用到的弧。此外，弧上的数字表示弧的编号。此时，c=((5,5…,5)1 )T， ×13 根据上述四个约束条件，分别求得四个情况下的最优决策变量x=((x12,x13,…,x75)1× )。 13 图 1 网络拓扑和流量需求
1.2 7节点算例求解 1.2.1 算例1（b1=[4;-4;0;0;0;0;0]T） 转化为线性规划问题： Minimize c T x1 Subject to Ax1=b1 x1>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x1*=[4 0 0 0 0 0 0 0 0 0 0 0 0]T 对应的最优值c T x1=20 1.2.2 算例2（b2=[4;0;-4;0;0;0;0]T） Minimize c T x2 Subject to Ax2=b2 X2>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x2*=[0 4 0 0 0 0 0 0 0 0 0 0 0]T 对应的最优值c T x2=20 1.2.3 算例3（b3=[0;-4;4;0;0;0;0]T） Minimize c T x3 Subject to Ax3=b3 X3>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x3*=[4 0 0 0 4 0 0 0 0 0 0 0 0]T 对应的最优值c T x3=40
惯性导航作业
惯性导航作业
一、数据说明： 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 1.用薄钢板制造一体积5m 3，长度不小于4m ，无上盖的货箱，要求钢板耗量最小。确定货箱的长x 1、宽x 2和高x 3。试列出问题的数学模型。 解：min 32312122x x x x x x z ++= s.t 5321=x x x 41≥x 0,,321≥x x x 2.将下面的线性规划问题表示为标准型并用单纯形法求解 max f=x 1+2x 2+x 3 s ．t ．2x 1+x 2-x 3≤2 -2x 1+x 2-5x 3≥-6 4x 1+x 2+x 3≤6 x i ≥0 i=1，2，3 解：先化标准形： Min 321x x x z -+= 224321=+-+x x x x 6525321=++-x x x x 646321=+++x x x x 列成表格： 1 2 1 610011460105122001112----- 可见此表已具备1°，2°，3°三个特点，可采用单纯形法。首先从底行中选元素-1，由2/2,6/2,6/4最小者决定选第一行第一列的元素2，标以记号，迭代一次得 1 2 1 2102310401162010021212 11-------- 再从底行中选元素-2/3，和第二列正元素1/2，迭代一次得 1 2 12 32 30 210231040116201002121211- ------ 再从底行中选元素-3，和第二列正元素2，迭代一次得 4 2 3 3 410120280114042001112--- 再迭代一次得 10 2 30 2 10 6 221023 1010213000421021013-- 选取最优解： 惯性导航技术综合实验 实验五惯性基组合导航及应用技术实验 惯性/卫星组合导航系统车载实验 一、实验目的 ①掌握捷联惯导/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 的位置及位置误差图： 发动机空燃比控制器 引言：我主要从事自动化相关研究。这里介绍我曾经接触过的发动机空燃比控制器设计中的优化问题。 发动机空燃比控制器设计中的最优化问题 AFR =a f m m && (1) 空燃比由方程（1）定义，在发动机运行过程中如果控制AFR 稳定在14.7可以获 得最好的动力性能和排放性能。如果假设进入气缸的空气流量a m &可以由相关单元检测得到，则可以通过控制进入气缸的燃油流量f m &来实现空燃比的精确控制。由于实际发动机的燃油喷嘴并不是直接对气缸喷燃油，而是通过进气歧管喷燃油，这么做会在进 气歧管壁上液化形成油膜，因此不仅是喷嘴喷出的未液化部分燃油会进入气缸，油膜 蒸发部分燃油也会进入气缸，如方程（2）。这样如何更好的喷射燃油成为了一个问题。 1110101122211ττττ?? ?? -?? ??????????=+????????-????????????-???? ? ??? ?? ????????? ?f f f v X x x u x x X x y =x && (2) 其中12、,==ff fv x m x m &&=f y m &,=fi u m &这里面，表示油膜蒸发量ff m &、fv m &表示为液化部分燃油、fi m &表示喷嘴喷射的燃油，在τf 、τv 、X 都已知的情况下，由现代控制理论知识，根据系统的增广状态空间模型方程（3） 0000001 1 011011114.70ττττ????-?? ??????????=-+-??????????????? ??????????????? ?? ??=?????? f f v v a X X u +q q m y q x x x &&& (3) 其中()0 14.7?t a q = y -m &。由极点配置方法，只要设计控制器方程（4），就可以 使得y 无差的跟踪阶跃输入，那么y 也能较好的跟踪AFR *a m /&。 12-- u =K q K x (4) 这里面的12、K K 确定，可由主导极点概念降维成两个参数12C ,C ，虽然都是最终稳态无差，但是目标是使得瞬态过程中y 和阶跃输入y r 的差异尽可能的小。所以原问 基于运动规划的惯性导航系统动态实验 二零一三年六月十日 实验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 惯性导航基础课程大作业报告（一）光纤陀螺误差建模与分析 班级：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; 航空科学与工程学院 2016年研究生入学考试复试大纲 一、复试方式：笔试＋面试 二、复试组织： 1、笔试：由航空学院统一组织，考试科目及复试大纲另见《航空科学与工程学院2013年考研复试安排》。 2、口试：以学科专业组为单位，由3-5位硕士生导师组成面试小组（组长为教授），每位考生的面试时间为20分钟。 三、复试流程和评分标准： 1）检查并核实考生面试所必备的个人证件和材料；考生可以提供有助于证明自己背景和能力的相关材料，证件和材料完备是面试的必要条件。 2）考生用英语口述个人基本情况、兴趣等，面试小组老师就考生基本情况提问，考生用英文回答问题。 3）考生朗读一段考场指定的专业外语短文，并口头翻译成中文。 4）面试小组老师就基础理论知识提问，学生用中文回答问题。 5）面试小组老师就专业知识提问，学生用中文回答问题。 面试结束后考生退场，在3-5个工作日后见航空学院网站“招生就业”栏目的“研究生招生”，会通知出学院的拟录取名单，在7层的研究生教学橱窗也会公布。 四、考场纪律 考生准时到达指定的复试考场，遵守考场秩序，尊重考试教师。 五、各学科专业组具体复试内容及参考书： 1、飞行力学与飞行安全系2016年硕士研究生入学复试程序 方式： 由3~6位硕士生导师组成面试小组，每位考生的面试时间为20分钟。 范围： 面试范围包括英语口语能力、专业英语阅读理解能力、专业基础理论知识和专业知识。具体环节如下： 1）对考生学习背景、心理、爱好和志愿等基本情况的了解。 2）考察考生的英语阅读和口头表达能力。 3）基础理论和专业知识面试。基础理论包括自动控制原理、理论力学和材料力学。专业知识包括飞行力学、飞行安全、飞行器总体设计、空气动力学等。 参考书： 基础理论可以选用任何一本考生熟悉的《自动控制原理》、《理论力学》、《材料力学》教材。专业课可以参考《飞机飞行动力学》（熊海泉编）或《飞机飞行性能》、《飞机的稳定与控制》等方面的参考书。 面试流程和评分标准： 1）检查并核实考生面试所必备的个人证件和材料；证件和材料完备是面试的必要条件。2）考生用英语口述个人基本情况、兴趣等，面试小组老师就考生基本情况提问，考生回答问题。 3）读一段指定的专业外语，并口头翻译成中文。 4）面试小组老师就基础理论知识提问，学生回答问题。 5）面试小组老师就专业知识提问，学生回答问题。 6）问答结束后，考生退场，面试老师根据考核要求和面试情况，对考生进行评分。 7）所有考生面试结束后，面试老师根据总体情况，对所有考生进行综合评估和比较，给出面试成绩。 2、人机与环境工程/制冷及低温工程2016年硕士研究生入学复试程序 方式： 由3~5位硕士生导师组成面试小组，每位考生的面试时间为20分钟。 范围： 1）英语阅读和口头表达能力。 2）对考生心理、基本情况的了解。 3）基础理论和专业知识面试。基础理论包括：自动控制原理，理论力学，流体力学；专业知识包括工程热力学，传热学，人机工程，低温制冷。考生可以选择其中1门基础理论和1门专业课作为面试内容，或者是综合知识。 参考书： 可以选用任何一本考生熟悉的《自动控制原理》、《理论力学》、《流体力学》教材。专业课可以选用考生熟悉的《工程热力学》，《传热学》，《人机工程》，低温制冷等方面的参考书。 面试流程和评分标准： 1）检查并核实考生面试所必备的个人证件和材料；证件和材料完备是面试的必要条件. 2）考生用英语口述个人基本情况、兴趣等，面试小组老师就考生基本情况提问，考生回答问题。 3）读一段指定的专业外语，并口头翻译成中文。 4）面试小组老师就基础理论知识提问，学生回答问题。 5）面试小组老师就专业知识提问，学生回答问题。 6) 问答结束后，考生退场，面试老师根据考核要求和面试情况，对考生进行评分。 数值分析第三次大作业 一、算法的设计方案： （一）、总体方案设计： 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 目标：使用带双步位移的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 实 验一 陀螺仪关键参数测试与分析实验 加速度计关键参数测试与分析实验 二零一三年五月十二日 实验一陀螺仪关键参数测试与分析实验 一、实验目得 通过在速率转台上得测试实验,增强动手能力与对惯性测试设备得感性认识；通过对陀螺仪测试数据得分析,对陀螺漂移等参数得物理意义有清晰得认识,同时为在实际工程中应用陀螺仪与对陀螺仪进行误差建模与补偿奠定基础。 二、实验内容 利用单轴速率转台,进行陀螺仪标度因数测试、零偏测试、零偏重复性测试、零漂测试实验与陀螺仪标度因数与零偏建模、误差补偿实验。 三、实验系统组成 单轴速率转台、MEMS 陀螺仪(或光纤陀螺仪)、稳压电源、数据采集系统与分析系统。 四、实验原理 1.陀螺仪原理 陀螺仪就是角速率传感器,用来测量载体相对惯性空间得角速度，通常输出与角速率对应得电压信号。也有得陀螺输出频率信号(如激光陀螺）与数字信号（把模拟电压数字化)。以电压表示得陀螺输出信号可表示为: (1-1）式中就是与比力有关得陀螺输出误差项，反映了陀螺输出受比力得影响,本实验不考虑此项误差。因此，式（１－１）简化为 (１-2）由（１－2）式得陀螺输出值所对应得角速度测量值: (1－3) 对于数字输出得陀螺仪,传感器内部已经利用标度因数对陀螺仪模拟输出进行了量化,直接输出角速度值,即： （１－4）就是就是陀螺仪得零偏，物理意义就是输入角速度为零时，陀螺仪输出值所对应得角速度。且 (１-5) 精度受陀螺仪标度因数、随机漂移、陀螺输出信号得检测精度与得影响。通常与表现为有规律性,可通过建模与补偿方法消除,表现为随机特性,可通过信号滤波方法抵制。因此，准确标定与就是实现角速度准确测量得基础。 五、陀螺仪测试实验步骤 1)标度因数与零偏测试实验 ａ、接通电源，预热一定时间; b、陀螺工作稳定后，测量静止情况下陀螺输出并保存数据; 航空科学与工程学院 《飞行力学实验班》课程实验飞机典型模态特性仿真 实验报告 学生姓名:姜南 学号:11051136 专业方向:飞行器设计与工程 指导教师:王维军 (2014年 6 月29日 一、实验目的 飞机运动模态是比较抽象的概念, 是课程教学中的重点和难点。本实验针对这一问题,采用计算机动态仿真和在人-机飞行仿真实验平台上的驾驶员在环仿真实验,让学生身临其境地体会飞机响应与模态特性的关系,加深对飞机运动模态特性的理解。 二、实验内容 1.纵向摸态特性实验 计算某机在某状态下的短周期运动、长周期运动的模态参数;进行时域的非实时或实时仿真实验,操纵升降舵激发长、短周期运动模态,并由结果曲线分析比较模态参数;放宽飞机静稳定性,观察典型操纵响应曲线,并通过驾驶员在环实时仿真体验飞机的模态特性变化。 2.横航向模态特性实验 计算某机在某状态下的滚转、荷兰滚、螺旋模态参数;进行时域仿真计算,操纵副翼或方向舵,激发滚转、荷兰滚等运动模态,并由结果曲线分析比较模态参数。 三、各典型模态理论计算方法及模态参数结果表 1 纵向模态纵向小扰动运动方程 0000 1 00 0e p e p e p u w e u w q p u w q X X u u X X g Z Z w w Z Z Z q q M M M M M δδδδδ δδδθθ????????-???? ????????? ? ???????????=+??????????????????? ?????????????????? A =[ X u X ?w Z u Z w 0?g Z q 0M ?u M ?w0 M q 010] =[?0.01999980.0159027?0.0426897?0.04034850?32.2869.6279 0?0.00005547?0.001893500?0.54005010] A 的特征值方程 |λ+0.0199998?0.01590270.0426897 λ+0.0403485032.2 ?869.627900.000055470.001893500λ+0.540050 ?1λ |=0 特征根λ1,2=?0.290657205979137±1.25842158268078i λ3,4=?0.00954194402086311±0.0377636398212079i 半衰期t 1/2由公式t 1/2=? ln2λ 求得,分别为 t 1/2,1=2.38475828674173s t 1/2,3=72.6421344585972s 振荡频率ω分别为 ω1=1.25842158268078rad/s ω3=0.0377636398212079rad/s 周期T 由公式T = (完整)北航惯性导航作业二. 编辑整理： 尊敬的读者朋友们： 这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（(完整)北航惯性导航作业二.）的内容能够给您的工作和学习带来便利。同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。 本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为(完整)北航惯性导航作业二.的全部内容。 惯性导航作业 一、数据说明： 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。 北京航空航天大学 数值分析大作业八 学院名称自动化 专业方向控制工程 学号 学生姓名许阳 教师孙玉泉 日期2014 年11月26 日 一．题目 关于x , y , t , u , v , w 的方程组(A.3) ???? ?? ?=-+++=-+++=-+++=-+++79 .0sin 5.074.3cos 5.007.1cos sin 5.067.2cos 5.0y w v u t x w v u t y w v u t x w v u t (A.3) 以及关于z , t , u 的二维数表（见表A-1）确定了一个二元函数z =f (x , y )。 表A-1 二维数表 t z u 0 0.4 0.8 1.2 1.6 2 0 -0.5 -0.34 0.14 0.94 2.06 3.5 0.2 -0.42 -0.5 -0.26 0.3 1.18 2.38 0.4 -0.18 -0.5 -0.5 -0.18 0.46 1.42 0.6 0.22 -0.34 -0.58 -0.5 -0.1 0.62 0.8 0.78 -0.02 -0.5 -0.66 -0.5 -0.02 1.0 1.5 0.46 -0.26 -0.66 -0.74 -0.5 1. 试用数值方法求出f (x , y ) 在区域}5.15.0,8.00|), {≤≤≤≤=y x y x D （上的近似表达式 ∑∑===k i k j s r rs y x c y x p 00 ),( 要求p (x , y )以最小的k 值达到以下的精度 ∑∑==-≤-=10020 7210)],(),([i j i i i i y x p y x f σ 其中j y i x i i 05.05.0,08.0+==。 2. 计算),(),,(* ***j i j i y x p y x f (i =1,2,…,8 ; j =1,2,…,5) 的值，以观察p (x , y ) 逼 近f (x , y )的效果，其中j y i x j i 2.05.0,1.0**+==。 《结构优化设计》 大作业报告 实验名称: 拓扑优化计算与分析 1、引言 大型的复杂结构诸如飞机、汽车中的复杂部件及桥梁等大型工程的设计问题，依靠传统的经验和模拟实验的优化设计方法已难以胜任，拓扑优化方法成为解决该问题的关键手段。近年来拓扑优化的研究的热点集中在其工程应用上，如: 用拓扑优化方法进行微型柔性机构的设计，车门设计，飞机加强框设计，机翼前缘肋设计，卫星结构设计等。在其具体的操作实现上有两种方法，一是采用计算机语言编程计算，该方法的优点是能最大限度的控制优化过程，改善优化过程中出现的诸如棋盘格现象等数值不稳定现象，得到较理想的优化结果，其缺点是计算规模过于庞大，计算效率太低;二是借助于商用有限元软件平台。本文基于matlab软件编程研究了不同边界条件平面薄板结构的在各种受力情况下拓扑优化，给出了几种典型结构的算例，并探讨了在实际优化中优化效果随各参数的变化，有助于初学者初涉拓扑优化的读者对拓扑优化有个基础的认识。 2、拓扑优化研究现状 结构拓扑优化是近20年来从结构优化研究中派生出来的新分支，它在计算结构力学中已经被认为是最富挑战性的一类研究工作。目前有关结构拓扑优化的工程应用研究还很不成熟，在国外处在发展的初期，尤其在国内尚属于起步阶段。1904 年Michell在桁架理论中首次提出了拓扑优化的概念。自1964 年Dorn等人提出基结构法，将数值方法引入拓扑优化领域，拓扑优化研究开始活跃。20 世纪80 年代初，程耿东和N. Olhoff在弹性板的最优厚度分布研究中首次将最优拓扑问题转化为尺寸优化问题，他们开创性的工作引起了众多学者的研究兴趣。1988年Bendsoe和Kikuchi发表的基于均匀化理论的结构拓扑优化设计，开创了连续体结构拓扑优化设计研究的新局面。1993年Xie.Y.M和Steven.G.P 提出了渐进结构优化法。1999年Bendsoe和Sigmund证实了变密度法物理意义的存在性。2002 年罗鹰等提出三角网格进化法，该方法在优化过程中实现了退化和进化的统一，提高了优化效率。目前常使用的拓扑优化设计方法可以分为两大类：退化法和进化法。结构拓扑优化设计研究，已被广泛应用于建筑、航天航空、机械、海洋工程、生物医学及船舶制造等领域。 3、拓扑优化建模（SIMP） 结构拓扑优化目前的主要研究对象是连续体结构。优化的基本方法是将设计区域划分为有限单元，依据一定的算法删除部分区域，形成带孔的连续体，实现连续体的拓扑优化。连续体结构拓扑优化方法目前比较成熟的是均匀化方法、变密度方法和渐进结构优化方法。 变密度法以连续变量的密度函数形式显式地表达单元相对密度与材料弹性模量之间的对应关系，这种方法基于各向同性材料，不需要引入微结构和附加的均匀化过程，它以每个单元的相对密度作为设计变量，人为假定相对密度和材料弹性模量之间的某种对应关系，程序实现简单，计算效率高。变密度法中常用的插值模型主要有:固体各向同性惩罚微结构模型(solidisotropic microstructures with penalization，简称SIMP)和材料属性的合理近似模型(rational approximation ofmaterial properties，简称RAMP)。而本文所用即为SIMP插值模型。最优化方法大作业答案
北航惯性导航综合实验五实验报告
最优化方法大作业
北航惯性导航综合实验四实验报告
北航数值分析大作业第一题幂法与反幂法
北航惯性导航大作业
北航飞行器设计与应用力学系.doc
北航数值分析报告第三次大作业
北航惯性导航综合实验四实验报告
北航数值分析大作业第二题精解
北航惯性导航综合实验一实验报告
北航飞行力学实验班飞机典型模态特性仿真实验报告(精)
(完整)北航惯性导航作业二.
北航数值分析大作业第二题
北航数值分析报告大作业第八题
结构优化设计大作业(北航)