北航《算法与数据结构》在线作业一 辅导资料

北航《算法与数据结构》在线作业一 辅导资料
北航《算法与数据结构》在线作业一 辅导资料

北航《算法与数据结构》在线作业一

一、单选题(共 25 道试题,共 100 分。)

1. 在一个单链表HL中,若要在指针q所指结点的后面插入一个由指针p所指向的结点,则执行()。

. q->nxt=p->nxt;p->nxt=q;

. p->nxt=q->nxt;q=p;

. q->nxt=p->nxt;p->nxt=q;

. p->nxt=q->nxt;q->nxt=p;

正确答案:

2. 快速排序的记录移动次数()比较次数,其总执行时间为O(nlog2n)。

. 大于

. 大于等于

. 小于等于

. 小于

正确答案:

3. 下述几种排序方法中,平均查找长度最小的是()

. 插入排序

. 选择排序

. 快速排序

. 归并排序

正确答案:

4. 对顺序表上的插入、删除算法的时间复杂性分析来说,通常以()为标准操作

. 条件判断

. 结点移动

. 算术表达式

. 赋值语句

正确答案:

5.

对于含有n个顶点条边的无向连通图,利用Prim算法生成最小代价生成树其时间复杂度为( )。

. O(log2n)

. O(n2)

. O(n)

. O(log2)

正确答案:

6. 下列关于栈的叙述正确的是()。

. 栈是非线性结构

. 栈是一种树状结构

. 栈具有先进先出的特征

. 栈具有后进先出的特征

正确答案:

7. 堆是一个键值序列{k1,k2,…, kn},对i=1,2,…,|_n/2_|,满足( )

. ki≤k2i≤k2i+1

. ki

. ki≤k2i且ki≤k2i+1(2i+1≤n)

. ki≤k2i 或ki≤k2i+1(2i+1≤n)

正确答案:

8. 如果只想得到1024个元素组成的序列中第5个最小元素之前的部分排序的序列,用()方法最快。

. 起泡排序

. 快速排序

. 简单选择排序

. 堆排序

正确答案:

9. 队列的插入操作是在()进行。

. 队首

. 队尾

. 队前

. 队后

正确答案:

10. 具有65个结点的完全二叉树其深度为()。

. 8

. 7

. 6

. 5

正确答案:

11. 某二叉树结点的前序序列为、、、、、G、,中序遍历为、、、、、、G。该二叉树结点的后序序列为 ( )。

. ,,,,,G,

. ,,,,,G,

. ,G,,,,,

. ,G,,,,,

正确答案:

12. 某程序的时间复杂度为(3n+nlog2n+n 2+8), 其数量级表示为()。

. O(n)

. O(nlog2n)

. O(n 2)

. O(log2n)

正确答案:

13. 单链表的一个存储结点包含()

. 数据域或指针域

. 指针域或链域

. 指针域和链域

. 数据域和链域

正确答案:

14. 若由森林转化得到的二叉树是非空的二叉树,则二叉树形状是()。

. 根结点无右子树的二叉树

. 根结点无左子树的二叉树

. 根结点可能有左二叉树和右二叉树

. 各结点只有一个儿子的二叉树

正确答案:

15. 对于顺序表,以下说法错误的是()

. 顺序表是用一维数组实现的线性表,数组的下标可以看成是元素的绝对地址

. 顺序表的所有存储结点按相应数据元素间的逻辑关系决定的次序依次排列

. 顺序表的特点是:逻辑结构中相邻的结点在存储结构中仍相邻

. 顺序表的特点是:逻辑上相邻的元素,存储在物理位置也相邻的单元中

正确答案:

16. 在一个单链表HL中,若要向表头插入一个由指针p指向的结点,则执行()。

. HL=p;p->nxt=HL;

. p->nxt=HL;HL=p;

. p->nxt=HL;p=HL;

. p->nxt=HL->nxt;HL->nxt=p;

正确答案:

17. 设有向图有n个顶点和条边,采用领接表作为其存储表示,在进行拓扑排序时,总的计算时间为()。

. O(nlog)

. O(n+)

. O(n*)

. O(n的平方)

正确答案:

18. 非空的循环单链表h的尾节点(由p所指向)满足()。

. p->nxt=NULL

. p=NULL

. p->nxt=h

. p=h

正确答案:

19. 对于单链表表示法,以下说法错误的是()

. 数据域用于存储线性表的一个数据元素

. 指针域或链域用于存放一个指向本结点所含数据元素的直接后继所在结点的指针

. 所有数据通过指针的链接而组织成单链表

. NULL称为空指针,它不指向任何结点,只起标志作用

正确答案:

20. 设有1000个无序的元素,希望用最快的速度挑选出其中前10个最大的元素,最好()排序法。

. 起泡排序

. 快速排序

. 堆排序

. 基数排序

正确答案:

21. 下列有关图遍历的说法中不正确的是()。

. 连通图的深度优先搜索是个递增过程

. 图的广度优先搜索中邻接点的寻找具有“先进先出”的特征

. 非连通图不能用深度优先搜索法

. 图的遍历要求每个顶点仅被访问一次

正确答案:

22. 下列图的说法中正确的是()。

. 一个具有 n 个顶点的无向完全图的边数为 n(n-1)

. 连通图的生成树是该图的一个极大连通子图

. 图的广度优先搜索是一个递归过程

. 在非连通图的遍历过程中,每调用一次深度优先搜索算法都得到该图的一个连通分量

正确答案:

23. 在以下栈的基本运算中,不是加工型运算的是().

. lnitStk(S)

. Push(S,X)

. Pop(S)

. mpty(S)

正确答案:

24. 队列操作的原则是()。

. 先进先出

. 后进先出

. 只能进行插入

. 只能进行删除

正确答案:

25. 以下说法错误的是()

. 线性表的元素可以是各种各样的,逻辑上相邻的元素在物理位置上不一定相邻

. 在线性表的顺序存储结构中,逻辑上相邻的两个元素在物理位置上不一定相邻

. 在线性表的链式存储结构中,逻辑上相邻的元素在物理位置上不一定相邻

. 线性表的链式存储结构的特点是用一组任意的存储单元存储线性表的数据元素

正确答案:

北航最优化方法大作业参考

北航最优化方法大作业参考

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

北航数值分析大作业一

《数值分析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

惯性导航作业

惯性导航作业

一、数据说明: 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:在作业报告中要写出“程序流程图、现阶段学习小结”,写明联系方式。

最优化方法大作业答案

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-- 选取最优解:

北航数理统计回归分析大作业

应用数理统计第一次大作业 学号: 姓名: 班级: 2013年12月

国家财政收入的多元线性回归模型 摘 要 本文以多元线性回归为出发点,选取我国自1990至2008年连续19年的财政收入为因变量,初步选取了7个影响因素,并利用统计软件PASW Statistics 17.0对各影响因素进行了筛选,最终确定了能反映财政收入与各因素之间关系的“最优”回归方程: 46?578.4790.1990.733y x x =++ 从而得出了结论,最后我们用2009年的数据进行了验证,得出的结果在误差范围内,表明这个模型可以正确反映影响财政收入的各因素的情况。 关键词:多元线性回归,逐步回归法,财政收入,SPSS 0符号说明 变 量 符号 财政收入 Y 工 业 X 1 农 业 X 2 受灾面积 X 3 建 筑 业 X 4 人 口 X 5 商品销售额 X 6

进出口总额X7

1 引言 中国作为世界第一大发展中国家,要实现中华民族的伟大复兴,必须把发展放在第一位。近年来,随着国家经济水平的飞速进步,人民生活水平日益提高,综合国力日渐强大。经济上的飞速发展并带动了国家财政收入的飞速增加,国家财政的状况对整个社会的发展影响巨大。政府有了强有力的财政保证才能够对全局进行把握和调控,对于整个国家和社会的健康快速发展有着重要的意义。所以对国家财政的收入状况进行研究是十分必要的。 国家财政收入的增长,宏观上必然与整个国家的经济有着必然的关系,但是具体到各个方面的影响因素又有着十分复杂的相关原因。为了研究影响国家财政收入的因素,我们就很有必要对其财政收入和影响财政收入的因素作必要的认识,如果能对他们之间的关系作一下回归,并利用我们所知道的数据建立起回归模型这对我们很有作用。而影响财政收入的因素有很多,如人口状况、引进的外资总额,第一产业的发展情况,第二产业的发展情况,第三产业的发展情况等等。本文从国家统计信息网上选取了1990-2009年这20年间的年度财政收入及主要影响因素的数据,包括工业,农业,建筑业,批发和零售贸易餐饮业,人口总数等。文中主要应用逐步回归的统计方法,对数据进行分析处理,最终得出能够反映各个因素对财政收入影响的最“优”模型。 2解决问题的方法和计算结果 2.1 样本数据的选取与整理 本文在进行统计时,查阅《中国统计年鉴2010》中收录的1990年至2009年连续20年的全国财政收入为因变量,考虑一些与能源消耗关系密切并且直观上

北航惯性导航大作业

惯性导航基础课程大作业报告(一)光纤陀螺误差建模与分析 班级: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;

北航博士研究生培养方案

交通科学与工程学院 道路与铁道工程(082301) 博士研究生培养方案 一、适用学科 道路与铁道工程(081401) 二、培养目标 1.坚持党的基本路线,热爱祖国,遵纪守法,品行端正,诚实守信,身心健康,具有良好的科研道德和敬业精神。 2.适应科技进步和社会发展的需要,在本学科上掌握坚实宽广的基础理论和系统深入的专门知识;熟练掌握一门外语;具有独立从事科学研究的能力;具有良好的综合素质。 3.在科学或专门技术上做出创造性的成果。 三、培养方向 1.道路与铁道工程的检测与加固; 2.土木工程结构分析与设计理论; 3.岩土本构理论及工程应用; 4.土木工程施工技术与材料; 5.工程结构仿真。 四、学制 学历博士研究生学制为3年。 博士研究生一般在入学后1年内完成课程学习,应在文献综述与开题报告前完成课程学分,应在博士论文答辩前完成全部学分和培养要求的有关环节。 鼓励博士研究生从入学开始就进行学位论文研究工作;文献综述与开题报告至申请学位论文答辩的时间间隔不得少于1年。 五、知识结构、课程设置与学分要求 1.知识结构要求 (1)基础理论与专业基础知识 高等工程数学与数学基础(数值分析、数理统计、矩阵理论、最优化理论与算法、数理方程、常微分方程、数学试验),专业基础知识(变分与有限元素法原理、高等混凝土结构、高等土力学、高等土木工程材料学、高等结构动力学、工程结构可靠度、工程塑性力学)。 (2)专业综合知识 混凝土结构非线性分析,高等钢结构,混凝土徐变力学,基础工程学,建设项目管理,高等岩石力学,建筑结构健康诊治,混凝土结构试验,岩土工程测试技术,建筑结构无损检测技术,土动力学,建筑结构有限元分析与应用,组合结构,城市地下工程,理论土力学与现代岩石测试技术,道路与铁道工程学科综合课。 (3)学科前沿与交叉学科知识 现代工程结构进展,材料科学进展,空间数据处理,科技信息检索与利用,科学

最优化方法大作业

发动机空燃比控制器 引言:我主要从事自动化相关研究。这里介绍我曾经接触过的发动机空燃比控制器设计中的优化问题。 发动机空燃比控制器设计中的最优化问题 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 的差异尽可能的小。所以原问

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

《数值分析》计算实习题目 第一题: 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;

北航卡尔曼滤波课程-捷联惯导静基座初始对准实验

卡尔曼滤波实验报告 捷联惯导静基座初始对准实验 一、实验目的 ①掌握捷联惯导的构成和基本工作原理; ②掌握捷联惯导静基座对准的基本工作原理; ③了解捷联惯导静基座对准时的每个系统状态的可观测性; ④了解双位置对准时系统状态的可观测性的变化。 二、实验原理 选取状态变量为:[]T E N E N U x y x y z X V V δδεεε=ψψψ??,其

中导航坐标系选为东北天坐标系,E V δ为东向速度误差,N V δ为北向速度误差,E ψ为东向姿态误差角,N ψ为北向姿态误差角,U ψ为天向姿态误差角,x ?为东向加速度偏置,y ?为北向加速度偏置,x ε为东向陀螺漂移,y ε为北向陀螺漂移,z ε为天向陀螺漂移。则系统的状态模型为: X AX W =+ (1) 其中, 1112212211 12 1321222331323302sin 000002sin 000000000sin cos 0000sin 000000cos 0000000000000000000000000000000000000000000000000000 0L g C C L g C C L L C C C L C C C L C C C A Ω-? ? ??-Ω????Ω-Ω? ?-Ω????Ω=? ?????? ?????????? ? [00000]E N E N U T V V W W W W W W δδψψψ=,E D V W W δψ 为零均值高斯 白噪声,分别为加速度计误差和陀螺漂移的噪声成分,Ω为地球自转角速度,ij C 为姿态矩 阵n b C 中的元素,L 为当地纬度。 量测量选取两个水平速度误差:[ ]T E N Z V V δδ=,则量测方程为: 10000000000100000000E E N N V X V δηδη???? ??=+???????????? (2) 即Z HX η=+ 其中,H 为量测矩阵,[]T E N ηηη=为量测方程的随机噪声状态矢量,为零均值高 斯白噪声。 要利用基本卡尔曼滤波方程进行状态估计,需要将状态方程和量测方程进行离散化。 系统转移矩阵为: 2323/1111102!3!! n n k k k k k k n T T T I TA A A A n ∞ -----=Φ=++++=∑ (3)

全日制工程硕士研究生培养方案-北航研究生院-北京航空航天大学

大型飞机高级人才培养班 航空工程全日制工程硕士研究生培养方案 一、适用类别或领域 航空工程(085232) 二、培养目标 材料工程、电子与通信工程、控制工程、航空工程领域全日制工程硕士 (以下简称航空工程等领域全日制工程硕士)是与以上各工程领域任职资格相联系的专业学位,主要为国民经济和国防建设等领域培养应用型、复合型高层次工程技术和工程管理人才。大飞机班旨在探索一条“以国家大型项目人才需求为索引,培养具有献身精神、团结协作精神、开拓创新精神的设计型和复合型人才”的研究生培养新模式,是北航研究生培养体系的一部分。 航空工程等领域全日制工程硕士培养的基本要求是: 1、坚持党的基本路线,热爱祖国、遵纪守法、品行端正、诚实守信、身心健康,具有良好的科研道德和敬业精神。 2、在本领域掌握坚实的基础理论和系统的专门知识,有较宽的知识面和较强的自立能力,具有大飞机设计、制造、运营、管理等领域需求的创造能力和工程实践能力。 3、掌握一门外国语。 三、培养模式及学习年限 1.航空工程等领域全日制工程硕士研究生培养实行导师负责制,或以导师为主的指导小组制,负责制订硕士研究生个人培养计划,选课、组织开题报告、论文中期检查、指导科学研究和学位论文,并与中国商飞、第一飞机设计研究院、西飞公司等航空企业联合培养,实行导师组指导。 2.硕士研究生一般用1学年完成课程学习,课程学习实行学分制,具体学习、考核及管理工作执行《北京航空航天大学研究生院关于研究生课程学习管理规定》。 3.专业实习是全日制工程硕士研究生培养中的重要环节,全日制工程硕士研究生在学期间,应保证不少于0.5年的工程实践。 4.学位论文选题应来源于航空工程等领域工程技术背景。鼓励实行双导师制,其中第一导师为校内导师,校外导师应是与本工程领域相关的专家,也可以根据学生的论文

北航数理统计第二次大作业-数据分析模板

数理统计第二次大作业材料行业股票的聚类分析与判别分析 2015年12月26日

材料行业股票的聚类分析与判别分析摘要

1 引言 2 数据采集及标准化处理 2.1 数据采集 本文选取的数据来自大智慧软件的股票基本资料分析数据,从材料行业的股票中选取了30支股票2015年1月至9月的7项财务指标作为分类的自变量,分别是每股收益(单位:元)、净资产收益率(单位:%)、每股经营现金流(单位:元)、主营业务收入同比增长率(单位:%)、净利润同比增长率(单位:%)、流通股本(单位:万股)、每股净资产(单位:元)。各变量的符号说明见表2.1,整理后的数据如表2.2。 表2.1 各变量的符号说明 自变量符号 每股收益(单位:元)X1 净资产收益率(单位:%)X2 每股经营现金流(单位:元)X3 主营业务收入同比增长率(单位:%)X4 净利润同比增长率(单位:%)X5 流通股本(单位:万股)X6 每股净资产(单位:元)X7 表2.2 30支股票的财务指标 股票代码X1 X2 X3 X4 X5 X6 X7 武钢股份600005-0.0990-2.81-0.0237-35.21-200.231009377.98 3.4444宝钢股份6000190.1400 1.980.9351-14.90-55.011642427.88 6.9197山东钢铁600022-0.11650.060.0938-20.5421.76643629.58 1.8734北方稀土6001110.0830 3.640.652218.33-24.02221920.48 2.2856

杭钢股份600126-0.4900-13.190.4184-36.59-8191.0283893.88 3.4497抚顺特钢6003990.219310.080.1703-14.26714.18112962.28 1.4667盛和资源6003920.0247 1.84-0.2141-5.96-19.3739150.00 1.2796宁夏建材6004490.04000.510.3795-22.15-92.3447818.108.7321宝钛股份600456-0.2090-2.53-0.3313-14.81-6070.2043026.578.1497山东药玻6005290.4404 5.26 1.2013 6.5016.7825738.018.5230国睿科技6005620.410011.53-0.2949 3.3018.9416817.86 3.6765海螺水泥600585 1.15169.05 1.1960-13.06-25.33399970.2612.9100华建集团6006290.224012.75-0.57877.90-6.4034799.98 1.8421福耀玻璃6006600.790014.250.9015 3.6017.27200298.63 6.2419宁波富邦600768-0.2200-35.02-0.5129 3.1217.8813374.720.5188马钢股份600808-0.3344-11.710.3939-21.85-689.22596775.12 2.6854亚泰集团6008810.02000.600.1400-23.63-68.16189473.21 4.5127博闻科技6008830.503516.71-0.1010-10.992612.8023608.80 3.0126新疆众和6008880.0523 1.04-0.910662.64162.0464122.59 5.0385西部黄金6010690.0969 3.940.115115.5125.5712600.00 2.4965中国铝业601600-0.0700-2.920.2066-9.0882.79958052.19 2.3811明泰铝业6016770.2688 4.66-1.09040.8227.8640770.247.4850金隅股份6019920.1989 3.390.3310-10.05-39.01311140.26 6.7772松发股份6032680.35007.00-0.3195-4.43-9.622200.00 6.0244方大集团0000550.0950 5.66-0.480939.2920.6742017.94 1.6961铜陵有色0006300.0200 1.220.6132 3.23-30.74956045.21 1.5443鞍钢股份000898-0.1230-1.870.7067-27.32-196.21614893.17 6.4932中钢国际0009280.572714.45-0.4048-14.33410.2441286.57 4.2449中材科技0020800.684610.27 1.219547.69282.1740000.00 6.8936中南重工0024450.1100 4.300.340518.8445.0950155.00 2.7030 2.2 数据的标准化处理 由于不同的变量之间存在着较大的数量级的差别,因此要对数据变量进行标准化处理。本文采用Z得分值法标准化的方法进行标准化,用x的值减去x的均值再除以样本的方差。也就是把个案转换为样本均值为0、标准差为1的样本。如果不同变量的变量值数值相差太大,会导致计算个案间距离时,由于绝对值较小的数值权数较小,个案距离的大小几乎由大数值决定,标准化过程可以解决此类问题,使不同变量的数值具有同等的重要性。经Z标准化输出结果见表 2.2。 表2.2 经Z标准化后的数据 ZX1ZX2ZX3ZX4ZX5ZX6ZX7

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

数值分析第三次大作业 一、算法的设计方案: (一)、总体方案设计: 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 =-,

波音737-800建模大作业

波音737—800飞机飞行模型建立实验 学院:航空自动化 专业:导航制导与控制

1 实验目的 根据飞机所提供的QAR数据,把飞机的飞行过程分为几个阶段,通过受力分析计算得出飞机在各阶段的各个时刻的地速以及飞机当时所处的地球经纬度。这之后,再把计算出来的这些数据与QAR里面的相对应的数据进行比较,得出数据误差。使我们对飞机各阶段的机体受力分析得到验证,最后确定飞机的整个飞行过程的模型。 2 实验内容 分析所得的QAR数据,根据QAR数据对飞机的飞行过程进行分阶段处理。然后查找相关资料,对飞机在飞行各阶段过程中进行受力分析。进而用MATLAB软件编写程序,计算出飞机各个阶段的地速和地球经纬度。最后把计算出来的数据和QAR里相应的数据作比较,用MA TLAB画出比较曲线图,得出计算误差,建立起飞机的飞行过程模型。在整个实验过程中要修学的课程有:《大气数据应用分析》、《导航原理与系统》、《飞机的飞行性能》、《惯性导航原理》、《MATLAB应用与编程》等等。

3 实验步骤 3.1 QAR数据分析 QAR数据分析 数据英文数据意义和用途所用仪表备注 1 东经Present Position Longitude 由0°本初子午线向东、西递增到180°导航仪 2 北纬Present Position Latitude 赤道向北递增到90°导航仪 3 磁航向Heading Magnetic 飞机纵轴在地平面上的投影,与磁子午线的 夹角(磁北顺时针转的夹角)。磁偏角:地 球表面任一点的磁子午圈同地理子午圈的夹 角。 磁罗盘上 有罗差修 正器,已经 抵消罗差, 所以磁罗 盘测的基 本就是磁 航向。 4 标准气压高度ALTITUDE 飞机到标准气压平面的垂直距离气压式高度表 5 左无线电高度RADIO HEIGHT Left 飞机到地面的垂直距离 无线电高 度表 6 机场标高AIR/GROUND 机场与海平面的垂直高度 7 左主起落架Left main gear air/end 起落架用于在地面停放及滑行时支撑飞机并 使飞机在地面上灵活运动,并吸收飞机运动 时产生的撞击载荷。主要用来判断飞机是否 起飞。 8 右主起落架Right main gear air/end 9 真空速Computed airspeed 飞机相对于空气的运动速度,根据空速可计 算地速,从而确定已飞距离和待飞时间。 空速表0.5~1.0 10 马赫数MACH 真空速与飞机所在高度的音速之比,当飞机 的M数超过临界M数时,飞机的空气动力特 马赫数表0.5~1.0

北航经济管理复习纲要(From xx_buaa)

固定资产:使用期限较长,单位价值在规定标准以上,在生产过程中为多个生产周期服务,在使用过程中保持原来物质形态的资产。 流动资产:可以在一年或虽然超过一年但仍然是一个生产经营周期内变现或耗用的资产。 无形资产:指没有物质实体而以某种特殊权利和技术知识等资源形态存在并发挥作用的资产。 递延资产:只不能全部计入当期损益,需要分期摊销计入成本的各项费用。 折旧:固定资产由于其价值在多个时期内损耗降低的部分 固定资产折旧:固定资产由于其价值在多个时期内损耗降低的部分。 资金的时间价值:资金在使用中随时间推移所发生的增值。 边际收益:当影响收益的产量或投入要素增加一个单位所增的收益。 边际成本:边际成本指的是每一单位新增生产的产品带来到总成本的增量。 边际利润:单位产量所增加的销售单价扣除边际成本的值。 机会成本:在有限资源及该资源多用途条件下,将该资源用于某种用途而放弃的可能用于其它用途形成的最大代价(付出)。 价值工程:以最低寿命周期成本,可靠地实现必要功能,以功能分析为核心,以提高产品或作业价值为目的的有组织的技术经济活动。 并行工程:是对产品及其相关过程,包括制造过程和支持过程,进行并行、一体化设计的一种系统化方法,目标是降低成本、提高生产率、加快上市速度。 4P(营销组合):市场营销中指产品、价格、渠道与促销。 系统:(钱学森)系统是由相互作用和相互依赖的若干组成部分(要素)结合而成的具有特定功能的有机整体。 市场经济:商品在市场上的价格完全由供需双方决定,没有任何一方(例如政府)加以干涉。 简述全面质量管理的内涵 质量管理仅靠数理统计方法是不够的,还需要一系列的组织管理工作;质量管理活动必须对质量、价格、交货期和服务进行综合考虑,而不仅仅只考虑质量;产品质量的产生、形成和实现过程包括了从市场研究到销售和服务的螺旋上升的循环过程,所以质量管理必须是全过程的管理;产品质量必须同成本联系起来考虑 试说明价格下降使需求量增加的原因 (1)价格降低后,消费者可以用同样的钱买到比此前更多的东西。这相当于消费者实际收入的提高,因而使需求量有所增加。这是由于价格变化所产生的“收入效应”而引起的需求量的增加。 (2)价格降低后,人们会把对替代品的需求转移到这种商品上来,因而使这种商品的需求量增加,这是由于价格变化所产生的“替代效应”引起的。 试述市场均衡价格是怎样形成的 如果市场价格高于均衡价格,,则供给量>均衡产量,此时,卖者找不到足够的买主,就会降低价格;如果市场价格低于均衡价格,,则供给量小于均衡产量,,此时,买者不能如数买到想要的东西,就会抬高价格。如果市场价格等于均衡价格,供给量等于需求量,买者想买的量等于卖者想卖得量,市场达到均衡。 试述系统工程的基本观点 系统整体性观点不着重强调系统单个元素的最优,而是强调整个系统就其功能而言效果最优。 相关与制约观点元素之间存在关系,并且这种关系可以表达。强调尽量地定量或用图表描述出各元素之间或各子系统之间的关系。 系统模拟观点系统可以建立模型,模型是原系统的简化系统,一般要求它具有原系统的主要性能。建模是分析、研究的基础。 系统优化观点 简述开展价值工程工作的六个主要步骤 运用[价值工程]方法开发产品需要按六个步骤(阶段)进行,其分别是:信息收集、创意构想、评估判断、细部发展、汇报审批和追踪实践。 第一步骤的信息收集,包括了设计理念(含功能、条件、标准…等)、成本估价资料、现场状况…等,尽量列出可能的范围,再透过机能(Function)定义和评估,找出标的物中的主要机能(必须是具备的机能),和次要机能(非绝对必要,是用来辅助主要机能)。也就是借着了解问题和机能分析,去筛选和找出问题所在(高成本或成本不合理的项目)。第二步骤是创意构想阶段,这个阶段是在小组成员都对问题充份了解之后针对主要机能开始做脑力激荡,这时候大家仅提构想(方案),不对构想做任何批评,也不考量方案的可行性,大家完全拋开传统模式的思考,让思想任意遨游,经由这个阶段,经常能产生一些具创新性的构想。 第三步骤是评估判断阶段,是对上阶段所提出的各项构想(方案)加以评估分析,首先可删除那些不可行的方案,再对剩余的可行方案做优缺点分析,并依节省成本的潜力及机能的改善做评估,及排列先后次序,然后取其优者,进入下一步的细部发展。 第四步骤,细部发展阶段,对选取之替代方案,就成本、可行性、节省之成本(或提升之机能)做详细完整的叙述。第五步骤,汇报审批阶段,将上阶段所做的报告书对业主做口头报告,这时候业主的接受与否决定了建议方案的是否执行。 第六步骤,追踪与实践,业主接受建议之后,下一个阶段就是落实该建议的执行。因此,这阶段的工作是要追踪确认接受的替代方案已纳入设计中,并协助业主消除替代方案执行的可能障碍。

导航系统大作业

导航系统

1.简述捷联惯性系统中地理系到机体系的姿态阵b g C 其含义及其功能。 答:含义:导航坐标系g g g O x y z -到机体坐标系b b b O x y z -的一组欧拉角为,,θγψ,导航坐 标系经过3次转动到机体坐标系。g g g x y z 依次沿g O z -、' b O x -、'' b O y -旋转角度-ψ、θ、γ后到b b b x y z 。姿态矩阵中包含了机体的姿态角方位角ψ、俯仰角θ和横滚角γ。 功能:机体陀螺仪输出的角速度信息经过补偿后,积分得到机体坐标系与导航坐标系的姿态信 息和姿态转移矩阵。捷联惯导系统中,加速度计与载体固连,利用姿态阵完成加速度计输出信息从机体坐标到导航坐标的转换。转换后的加速度计信息经过积分可得到机体在导航坐标系下的速度和位置。 2.画出并用式表达速度三角形(地速、控速、风速)及航迹角、航向角与偏流角之间的关系。 答:风速:空气相对于地面的运动速度;空速:飞机相对于空气运动的速度;地速:飞机相对 于地面的运动速度。=+v v v 风地空 航向角:机头在水平面投影与真北方向的夹角?;偏流角:空速矢量和地速矢量之间的夹角,用 δ表示;航迹角:飞机速度矢量在水平面投影与真北方向的夹角。航向角?加上偏流角δ等于地 速v 地的方位角α。 3.简述惯性导航系统、卫星导航系统、多普勒导航、塔康、VOR/DME 、天文导航其各自的基本工作原理、特点及误差特性。 答:一、惯性导航系统 (1)工作原理 以牛顿力学定律为基础,以陀螺仪和加速度计为敏感器件进行导航参数解算。系统根据陀螺

仪的输出建立导航坐标系,根据加速度计输出解算出运载体的速度和位置,从而实现姿态和航向解算。 (2)特点 惯性导航系统不需要任何外来信息,也不会向外辐射任何信息,仅依靠惯性器件就能全天候,全球性的自主三维定位和三维定向,同时具备自主性、隐蔽性和信息的完备性。 (3)误差特性 误差随时间积累,短时间导航精度较高。 二、卫星导航系统 (1)工作原理 以卫星和用户接收机天线之间的距离观测量为基准,根据已知的卫星的瞬时坐标(轨道根数),来确定用户观测点的经纬度和高程信息。 (2)特点 卫星导航系统具有全天候、高精度、自动化、高效益、性能好,应用广的特点,是一种被动式的导航系统。但需要地面站支持,电波易受干扰。 (3)误差特性 在卫星导航系统中,影响测量结果的误差因素有与卫星有关的误差,与观测有关的误差,和与观测站有关的误差。包括卫星时钟、星历误差,也受电离层、对流层和周围环境事物遮挡等影响。长时间导航精度较高。 三、多普勒导航系统 (1)工作原理 多普勒导航系统是一种自助式推算导航系统。机载多普勒雷达向地面发射电波和接收地面的回波,通过测量地面回波的多普勒频移,通过定位解算,即可得到飞行器的位置信息。 (2)特点 多普勒导航系统不需要有地面或卫星发射台,发射的波束窄,角度陡,难以被监测,自主性强,测速精度高,不需要初始对准。 (3)误差特性 影响多普勒导航系统的误差有测速误差和飞机的角度敏感误差。系统的定位误差发散,随时间推移而增大。 四、塔康导航系统 (1)工作原理 塔康导航系统是由塔康地面设备(塔康信标)和机载设备组成。其采用极坐标体制定位,飞机定时向地面台发送和接收信号,机载设备与塔康信标配合连续解算出飞机所在点相对于信标的方位角和距离。 (2)特点

北航飞行器多学科设计优化复习题

飞行器多学科设计优化复习题 1.优化设计问题的三要素是什么?给出一个优化设计问题的例子,分别说明三个要素的具体内容。 三要素分别是设计变量,约束条件和目标函数。 以结构优化设计为例,设计变量可能是蒙皮厚度,前后翼梁缘条厚度,前后翼梁腹板厚度等结构参数;约束条件是机翼强度要求、刚度要求等目标函数是最小化结构重量。 2.飞行器设计一般分哪几个阶段?飞行器多学科优化设计有什么意义? 飞行器设计分三个阶段:概念设计、初步设计、详细设计。 飞行器MDO的意义为: (1)MDO符合系统工程的思想。能有效提高飞行器的设计质量 (2)MDO为飞行器设计提供了一种并行设计模式。 (3)MDO的设计模式与飞行器设计组织体制一致,能够实现更高程度的自动化。 (4)MDO的模块化结构使飞行器设计过程具有很强的灵活性。 3.在飞行器设计过程中,多学科设计优化方法与传统设计方法之间有哪些相同和不同点。 传统的飞行器设计优化中,采取的是一种串行的设计模式,往往首先进行性能设计优化,然后进行结构、操纵和控制系统设计优化,最后进行工艺装备设计。在传统的方法中,各个学科任务成了实现系统设计的最基本单元,影响飞机性能的气动、推进、结构和控制等学科被人为地割裂开来,各学科之间相互耦合所产生的协同效应并未被充分考虑进去,这可能导致失去系统的整体最优解,串行的模式也使得设计时间周期和成本大大增加。 而多学科优化设计技术是一种并行设计模式,它以各子系统、学科的优化设计为基础,在飞行器各个阶段力求各学科的平衡,充分考虑哥们学科之间的相互影响和耦合作用,应用有效的设计/优化策略和分布式计算机网络系统,来组织和管理整个系统的设计过程,通过充分利用各个学科之间的相互作用所产生的协同效应,以获得系统的整体最优解。 相同点在于都有对于子学科的分解,但是MDO更注重子学科间的协同。 4.给出MDO的三种定义,根据你的理解,MDO该如何定义? Definition1:MDO是一种通过充分探索和利用系统中相互作用的协同机制来设计复杂系统和子系统的方法论。 Definition2:MDO是指在复杂工程系统的设计过程中,必须对学科(子系统)之间的相互作用进行分析,并且充分利用这些相互作用进行系统优化合成的方法。 Definition3:多学科设计优化就是进行复杂系统的设计过程中,结合系统的多学科本质,充分利用各种多学科设计与多学科分析工具,最终达到基于多学科优化的方法论。 My Definition:当设计中每个因素都影响另外的所有因素时,确定该改变哪个因素以及改变到什么程度的一种设计方法。 5.多学科设计优化中,什么是学科分析?什么是系统分析? 学科分析:也成为子系统分析或子空间分析,以某一学科设计变量,其他学科对该学科的耦合状态变量和系统的参数为输入,根据某一学科满足的物理规律确定其物理特性的过程 系统分析:对整个系统,给定一组设计变量X,通过求解系统的状态方程得到系统状态变量的过程。 6.什么是多学科设计优化的状态变量?学科状态变量和耦合状态变量之间有什么区别?

相关文档
最新文档