北航惯性导航大作业

惯性导航基础课程大作业报告（一）光纤陀螺误差建模与分析

班级：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，εz对北向速度误差δVy均有地球自转周期，舒勒周期分量的影响。其中，εx对δVy的影响较大;εy，εz产生的影响几乎相近。

3.εx，εy，εz对东向速度误差δL的影响

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;

a=Ws^2-Wie^2;

b=sin(L);

c=cos(L);

d=sin(Wie*t);

e=cos(Wie*t);

f=sin(Ws*t);

h=cos(Ws*t);

mcL1=Ws^2*0.1/a*(d/Wie-f/Ws);

mcL2=(Ws^2*Wie*b/a*(h/Ws^2-e/Wie^2)+b/Wie)*0.1;

mcL3=(Ws^2*c*e/Wie/a-Wie*c*h/a-c/Wie)*0.1;

plot(t,[mcL1',mcL2',mcL3']);

title('Ex,Ey,Ez对mcL的影响');

xlabel('时间t');

ylabel('mcL(t)');

legend('Ex-mcL1','Ey-mcL2','Ez-mcL3');

grid;

分析：εx，εy，εz对纬度误差δL均有地球自转周期的影响，εx还会有舒勒周期分量的影响。

4.εx，εy，εz对东向速度误差δλ的影响

clc;clear all;

t=1:0.01:50;

pi=3.14;

g=9.8;

L=pi/180*39;

Ws=2*pi/84.4*60;

Wie=2*pi/24;

R=g/(Ws)^2;

a=Ws^2-Wie^2;

b=sin(L);

c=cos(L);

d=sin(Wie*t);

e=cos(Wie*t);

f=sin(Ws*t);

h=cos(Ws*t);

mcLONG1=(tan(L)/Wie*(1-e)-Wie*c*h/a)*0.1;

mcLONG2=(sec(L)*(Ws^2-Wie^2*c^2)*f/Ws/a-Ws^2*tan(L)*b*d/Wie/a-t*c)*0.1;

mcLONG3=(Ws^2*b*d/Wie/a-Wie^2*b*f/a/Ws-t*b)*0.1;

plot(t,[mcLONG1',mcLONG2',mcLONG3']);

title('Ex,Ey,Ez对mcLONG.的影响');

xlabel('时间t');

ylabel('mcLONG.(t)');

legend('Ex-mcLONG.1','Ey-mcLONG.2','Ez-mcLONG.3');

grid;

axis square;

分析：εx，εy，εz对经度误差δλ均有地球自转周期的影响，εy还会有舒勒周期分量的影响，其中，εy，εz还产生了随时间累积的分量。

5.εx，εy，εz对东向速度误差δφX的影响

clc;clear all;

t=1:0.01:25;

pi=3.14;

g=9.8;

L=pi/180*39;

Ws=2*pi/84.4*60;

Wie=2*pi/24;

R=g/(Ws)^2;

a=Ws^2-Wie^2;

b=sin(L);

c=cos(L);

d=sin(Wie*t);

e=cos(Wie*t);

f=sin(Ws*t);

h=cos(Ws*t);

mcAngle1=(Ws*f-Wie*d)/a*0.1;

mcAngle2=Wie*b*(e-h)/a*0.1;

mcAngle3=Wie*c*(h-e)/a*0.1;

plot(t,[mcAngle1',mcAngle2',mcAngle3']);

title('Ex,Ey,Ez对mcAngle的影响');

xlabel('时间t');

ylabel('mcAngle(t)');

legend('Ex-mcAngle1','Ey-mcAngle2','Ez-mcAngle3');

grid;

axis square;

分析：εx，εy，εz对水平方位误差φx均有地球自转周期，舒勒周期分量的影响，其中，εx对φx产生的影响最大。

5.εx，εy，εz对东向速度误差δφY的影响

clc;clear all;

t=1:0.01:25;

pi=3.14;

g=9.8;

L=pi/180*39;

Ws=2*pi/84.4*60;

Wie=2*pi/24;

R=g/(Ws)^2;

a=Ws^2-Wie^2;

b=sin(L);

c=cos(L);

d=sin(Wie*t);

e=cos(Wie*t);

f=sin(Ws*t);

h=cos(Ws*t);

mcAngle1=Wie*b*(h-e)/a*0.1;

mcAngle2=((Ws^2-Wie^2*c^2)/Ws/a*f-Wie*b^2/a*d)*0.1;

mcAngle3=Wie*b*c/a*(d-Wie/Ws*f)*0.1;

plot(t,[mcAngle1',mcAngle2',mcAngle3']);

title('Ex,Ey,Ez对mcAngle y的影响');

xlabel('时间t');

ylabel('mcAngle y(t)');

legend('Ex-mcAngle1','Ey-mcAngle2','Ez-mcAngle3');

grid;

axis square;

分析：εx，εy，εz对水平方位误差φy均有地球自转周期的影响，而εx，εy还产生了舒勒周期分量的影响，其中，εy对φy 产生的影响最大。

5.εx，εy，εz对东向速度误差δφZ的影响

clc;clear all;

t=1:0.01:25;

pi=3.14;

g=9.8;

L=pi/180*39;

Ws=2*pi/84.4*60;

Wie=2*pi/24;

R=g/(Ws)^2;

a=Ws^2-Wie^2;

b=sin(L);

c=cos(L);

d=sin(Wie*t);

e=cos(Wie*t);

f=sin(Ws*t);

h=cos(Ws*t);

mcAngle1=(sec(L)/Wie*(1-e)+Wie*b*b/c*(h-e)/a)*0.1;

mcAngle2=(Wie^2*b*c-Ws^2*b/c)/a*(d/Wie-f/Ws)*0.1;

mcAngle3=((Ws^2-Wie^2*c^2)*d/Wie/a-Wie^2*b^2*f/Ws/a)*0.1;

plot(t,[mcAngle1',mcAngle2',mcAngle3']);

title('Ex,Ey,Ez对mcAngle z的影响');

xlabel('时间t');

ylabel('mcAngle z(t)');

legend('Ex-mcAngle1','Ey-mcAngle2','Ez-mcAngle3');

grid;

分析：εx，εy，εz对方位姿态误差φz均有地球自转周期的影响，而εy还产生了舒勒周期分量的影响，其中，εx对φz产生的影响最大。

(二) 加速度计零偏引起的系统误差

1. Δx对φy,φz以及Δy对φx的影响

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;

a=Ws^2-Wie^2;

b=sin(L);

c=cos(L);

d=sin(Wie*t);

e=cos(Wie*t);

f=sin(Ws*t);

h=cos(Ws*t);

amc=0.0001*g;

mcAnglez=b/c/g*(1-h)*amc*180/pi*3600;

mcAngley=(1-h)/g*amc*180/pi *3600;

mcAnglex=-(1-h)/g*amc*180/pi *3600;

subplot(311);plot(t,mcAnglex,'r-');xlabel('时间t');ylabel('mcAnglex (t)');

legend('amcy-mcAnglex');grid;

subplot(312);plot(t,mcAngley,'g-');xlabel('时间t');ylabel('mcAngley (t)');

legend('amcx-mcAngley');grid;

subplot(313);plot(t,mcAnglez,'b-');xlabel('时间t');ylabel('mcAnglez (t)');

legend('amcx-mcAnglez');grid;

分析：Δx对φy,φz以及Δy对φx的影响包含了常值分量和舒乐振荡分量。

2.Δx对δVx, Δy对δVy, Δy对δL, Δx对δλ的影响

clc;clear all;

t=1:0.01:25;

pi=3.14;

g=9.8;

L=pi/180*39;

Ws=2*pi/84.4*60;

Wie=2*pi/24;

R=g/(Ws)^2;

a=Ws^2-Wie^2;

b=sin(L);

c=cos(L);

d=sin(Wie*t);

e=cos(Wie*t);

f=sin(Ws*t);

h=cos(Ws*t);

amc=0.0001*g;

mcVx=f/Ws*amc;

mcVy=f/Ws*amc;

mcL=(1-h)/g*amc*180/pi*3600;

mcLONG=sec(L)*amc/g*(1-h)*180/pi*3600;

subplot(221);plot(t,mcVx,'r:');xlabel('时间t');ylabel('mcVx (t)');grid;

subplot(222);plot(t,mcVy,'g:');xlabel('时间t');ylabel('mcVy (t)');grid;

subplot(223);plot(t,mcL,'b:');xlabel('时间t');ylabel('mcL (t)');grid;

subplot(224);plot(t,mcLONG,'b:');xlabel('时间t');ylabel('mcLONG (t)');grid;

分析：Δx对δVx, Δy对δVy, Δy对δL, Δx对δλ,均包含了常值分量和舒乐振荡分量的影响。

（三）起始误差对系统误差的影响

1. δVx0,δVy0,δL0,φx0,φy0,φz0对δVx的影响

syms s;

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;

a=s^2+Ws^2;

b=s^2+Wie^2;

c=sin(L);

d=cos(L);

c11=s/a;

c12=0;

c13=s*g*Wie*c/a/b;

c14=c13;

c15=-g*(s^2+Wie^2*d^2)/a/b;

c16=-g*Wie^2*c*d/a/b;

mcVx1=ilaplace(c11*0.1);

mcVx3=ilaplace(c13*0.0005*pi/180);

mcVx4=ilaplace(c14*20/3600*pi/180);

mcVx5=ilaplace(c15*20/3600*pi/180);

mcVx6=ilaplace(c16*5/60*pi/180);

mcVx1t=subs(mcVx1);

mcVx3t=subs(mcVx3);

mcVx4t=subs(mcVx4);

mcVx5t=subs(mcVx5);

mcVx6t=subs(mcVx6);

plot(t,mcVx1t);title('mcVx0对mcVx的影响');xlabel('时间t');ylabel('mcVx(t)');grid; figure(2);

plot(t,mcVx3t);title('mcL0对mcVx的影响');xlabel('时间t');ylabel('mcVx(t)');grid; figure(3);

plot(t,mcVx4t,'g*',t,mcVx5t,'g+',t,mcVx6t,'g-');

legend('anglex0-mcVx','angley0-mcVx','anglez0-mcVx');

title('初始水平姿态角，方位姿态角对mcVx的影响');xlabel('时间t');ylabel('mcVx(t)');grid;

分析：δVx0对δVx只有舒乐周期震荡的影响；δL0对δVx有地球自转周期和舒乐周期震荡两个方面的影响；φx0,φy0,φz0对δVx有地球自转周期和舒乐振荡周期的影响，且φy0产生的影响最大。

2. δVx0,δVy0,δL0,φx0,φy0,φz0对δVy的影响

syms s;

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;

a=s^2+Ws^2;

b=s^2+Wie^2;

c=sin(L);

d=cos(L);

c21=0;

c22=s/a;

c23=-g*Wie^2/a/b;

c24=s^2*g/a/b;

c25=s*g*Wie*c/a/b;

c26=-s*g*Wie*d/a/b;

mcVy2=ilaplace(c22*0.1);

mcVy3=ilaplace(c23*0.0005*pi/180);

mcVy4=ilaplace(c24*20/3600*pi/180);

mcVy5=ilaplace(c25*20/3600*pi/180);

mcVy6=ilaplace(c26*5/60*pi/180);

mcVy1t=subs(mcVy1);

mcVy3t=subs(mcVy3);

mcVy4t=subs(mcVy4);

mcVy5t=subs(mcVy5);

mcVy6t=subs(mcVy6);

plot(t,mcVy1t);title('mcVy0对mcVy的影响');xlabel('时间t');ylabel('mcVy(t)');grid; figure(2);

plot(t,mcVy3t);title('mcL0对mcVy的影响');xlabel('时间t');ylabel('mcVy(t)');grid; figure(3);

plot(t,mcVy4t,'.',t,mcVy5t,'+',t,mcVy6t,'-');

legend('anglex0-mcVy','angley0-mcVy','anglez0-mcVy');

title('初始水平姿态角，方位姿态角对mcVy的影响');xlabel('时间t');ylabel('mcVy(t)');grid;

分析：δVy0对δVy只有舒乐周期震荡的影响；δL0对δVy有地球自转周期和舒乐周期震荡两个方面的影响；φx0,φy0,φz0对δVy有地球自转周期和舒乐振荡周期的影响，且φx0,φz0产生的影响较大。

3. δVx0,δVy0,δL0,φx0,φy0,φz0对δL0的影响

syms s;

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;

a=s^2+Ws^2;

b=s^2+Wie^2;

c=sin(L);

d=cos(L);

c31=0;

c32=1/a/R;

c33=s/a;

c34=s*Ws^2/a/b;

c35=Ws^2*Wie*c/a/b;

c36=-Ws^2*Wie*d/a/b;

mcL2=ilaplace(c32*0.1);

mcL3=ilaplace(c33*0.0005);

mcL4=ilaplace(c34*20/3600);

mcL5=ilaplace(c35*20/3600);

mcL6=ilaplace(c36*5/60);

mcL2t=subs(mcL2);

mcL3t=subs(mcL3);

mcL4t=subs(mcL4);

mcL5t=subs(mcL5);

mcL6t=subs(mcL6);

plot(t,mcL2t);title('mcVy0对mcL的影响');xlabel('时间t');ylabel('mcL(t)');grid; figure(2);

plot(t,mcVy3t);title('mcL0对mcL的影响');xlabel('时间t');ylabel('mcL(t)');grid; figure(3);

plot(t,mcL4t,'.',t,mcL5t,'+',t,mcL6t,'-');

legend('angleX0-mcL','angley0-mcL','anglez0-mcL');

title('初始水平姿态角，方位姿态角对mcL的影响');xlabel('时间t');ylabel('mcL(t)');grid;

分析：δVy0对δL只有舒乐周期震荡的影响；δL0对δL有地球自转周期和舒乐周期震荡两个方面的影响；φx0,φy0,φz0对δL 有地球自转周期和舒乐振荡周期的影响，且φz0产生的影响较大。

4. δVx0,δVy0,δL0,φx0,φy0,φz0对φx的影响

syms s;

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;

a=s^2+Ws^2;

b=s^2+Wie^2;

c=sin(L);

d=cos(L);

c41=0;

c42=-1/a/R;

c43=-s*Wie^2/a/b;

c44=s^3/a/b;

c45=s^2*Wie*c/a/b;

c46=-s^2*Wie*d/a/b;

anglex2=ilaplace(c42*0.1);

anglex3=ilaplace(c43*0.0005);

anglex4=ilaplace(c44*20/3600);

anglex5=ilaplace(c45*20/3600);

anglex6=ilaplace(c46*5/60);

anglex2t=subs(anglex2);

anglex3t=subs(anglex3);

anglex4t=subs(anglex4);

anglex5t=subs(anglex5);

anglex6t=subs(anglex6);

plot(t,anglex2t);title('mcVy0对abglex的影响');xlabel('时间t');ylabel('abglex');grid;

figure(2);

plot(t,anglex3t);title('mcL0对anglex的影响');xlabel('时间t');ylabel('anglex');grid;

figure(3);

plot(t,anglex4t,'.',t,anglex5t,'+',t,anglex6t,'-');

legend('angleX0-anglex','angley0-anglex','anglez0-anglex');

title('初始水平姿态角，方位姿态角对anglex的影响');xlabel('时间t');ylabel('anglex');grid;

惯性导航作业

惯性导航作业

一、数据说明： 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

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

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

惯性/卫星组合导航系统车载实验 一、实验目的 ①掌握捷联惯导/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 的位置及位置误差图：

最新北航数理统计大作业-多元线性回归

北航数理统计大作业-多元线性回归

应用数理统计多元线性回归分析 （第一次作业） 学院： 姓名： 学号： 2013年12月

交通运输业产值的多元线性回归分析 摘要：本文基于《中国统计年鉴》（2012年版）统计数据，寻找影响交通运输业发展的因素，包括工农业发展水平、能源生产水平、进出口贸易交流以及居民消费水平等，利用统计软件SPSS对各因素进行了筛选分析，采用逐步回归法得到最优多元线性回归模型，并对模型的回归显著性、拟合度以及随机误差的正态性进行了检验，最后可以利用有效的最优回归模型对将来进行预测。 关键字：多元线性回归，逐步回归，交通运输产值，工业产值，进出口总额1，引言 交通运输业指国民经济中专门从事运送货物和旅客的社会生产部门，包括铁路、公路、水运、航空等运输部门。它是国民经济的重要组成部分，是保证人们在政治、经济、文化、军事等方面联系交往的手段，也是衔接生产和消费的一个重要环节。交通运输业在现代社会的各个方面起着十分重要的作用，因此研究交通运输业发展水平与各个影响因素间的关系显得十分重要，建立有效的数学相关模型对于预测交通运输业的发展，制定相关政策方案提供依据。根据经验交通运输业的发展受到工农业发展、能源生产、进出口贸易以及居民消费水平等众因素的影响，故建立一个完整精确的数学模型在理论上基本无法实现，并且在实际运用中也没有必要，一种简单有效的方式就是寻找主要影响因素，分析其与指标变量的相关性，建立多元线性回归模型就是一种有效的方式。 变量与变量之间的关系分为确定性关系和非确定性关系，函数表达确定性关系。研究变量间的非确定性关系，构造变量间经验公式的数理统计方法称为

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

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

实验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;

北航航空工程大型通用软件应用大作业样本

航空科学与工程学院 《航空工程大型通用软件应用》大作业 机翼结构设计与分析 组号第3组 小组成员11051090 赵雅甜 11051093 廉佳 11051100 王守财 11051108 刘哲 11051135 张雄健 11051136 姜南 6月

目录 一 CATIA部分....................................... 错误!未定义书签。( 一) 作业要求..................................... 错误!未定义书签。( 二) 作业报告..................................... 错误!未定义书签。 1、三维模型图................................... 错误!未定义书签。 2、工程图....................................... 错误!未定义书签。 二 FLUENT部分...................................... 错误!未定义书签。( 一) 作业要求..................................... 错误!未定义书签。( 二) 作业报告..................................... 错误!未定义书签。 1、计算方法和流程............................... 错误!未定义书签。 2、网格分布图................................... 错误!未定义书签。 3、气动力系数................................... 错误!未定义书签。 4、翼型表面压力曲线............................. 错误!未定义书签。 5、翼型周围压力云图............................. 错误!未定义书签。 6、翼型周围x方向速度云图....................... 错误!未定义书签。 7、翼型周围y方向速度云图....................... 错误!未定义书签。 8、翼型周围x方向速度矢量图..................... 错误!未定义书签。 9、翼型周围y方向速度矢量图..................... 错误!未定义书签。 10、流线图...................................... 错误!未定义书签。 三 ANSYS部分....................................... 错误!未定义书签。( 一) 作业要求..................................... 错误!未定义书签。( 二) 作业报告..................................... 错误!未定义书签。 1、机翼按第一强度理论计算的应力云图............. 错误!未定义书签。 2、机翼按第二强度理论计算的应力云图............. 错误!未定义书签。 3、机翼按第三强度理论计算的应力云图............. 错误!未定义书签。 4、机翼按第四强度理论计算的应力云图............. 错误!未定义书签。

北航惯性导航大作业

惯性导航基础课程大作业报告（一）光纤陀螺误差建模与分析 班级：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 =-，

北航涡轮泵大作业

学号姓名成绩 《冲击式涡轮和反力式涡轮的设计计算》 总结：对冲击式涡轮和反力式涡轮进行设计计算，得到计算结果，具体见表1 和表2。 表1 反力式涡轮的计算结果 表2 冲击式涡轮的计算结果

根据计算结果，我们对比可以得到冲击式涡轮和反力式涡轮的相同点 是： 冲击式涡轮和反力式涡轮在计算功率时，均由泵的功率决定，由 T P N N =∑ 计算。 不同点具体见表3. 表3 反力式涡轮和冲击式涡轮的比较 1. 冲击式涡轮出口压力值取决于涡轮排气是直接排入周围环境还是导入辅助喷管，但两种情况下出口压强和反力式相比均很小。而反力式涡轮通常用于补燃式的液体火箭发动机中的涡轮泵中，所

以在不记喷注器压降的条件下，涡轮的出口压力等于燃烧室的压力。 2.在计算反力式涡轮的参数时，由于反力度容易确定，在分析过程 中广泛采用热力反力度。

反力式涡轮的设计计算 一.反力式涡轮参数的选择 在具有冲击式涡轮的供应系统(无补燃发动机系统)中，由燃气发生器产生的富燃燃气驱动涡轮，涡轮不冷却，富燃燃气的温度在1000～1200K 的范围内，比富氧燃气的允许温度(600～800K)高得多。另外，富燃燃气的气体常数比富氧燃气的气体常数大一些，这些都有利于减小需通过涡轮的燃气流量。 涡轮流量m t q 是具有冲击式涡轮的供应系统的主要参数之一。m t q 值越小，发动机的比冲就越高。涡轮流量m t q 可由泵和涡轮的功率平衡： T Pf Po N N N =+ 泵的需用功率降低，可减小通过涡轮的燃气流量，因此应尽量提高泵的效率。选定泵的结构并确定其效率后，可根据功率平衡求出所需的涡轮燃气流量，由此确定涡轮的效率。 涡轮入口压力(燃气发生器压力)取决于氧化剂泵的出口压力。当用燃料冷却推力室时，燃料泵出口压力比氧化剂泵的出口压力高。 涡轮出口压力之值取决于涡轮排气是直接排入周围环境还是导入辅助喷管。 冲击式涡轮计算的原始数据为： (1)涡轮的设计功率：涡轮功率T N 由泵所需的功率决定，由涡轮泵装置设计任务给定： 其中，T N —涡轮的设计功率，又称涡轮的轴功率； Pf N —燃料泵的轴功率； Pf N —氧化剂泵的轴功率。 (2)涡轮的设计角速度：涡轮的设计转速ω由泵不发生汽蚀时允许的最大角速度确定； (3)涡轮工质的物理常数和温度：涡轮进口总压*0P 、进口总温*0T 、和出口静压2P ；涡轮工质的绝热指数k 和气体常数R 。 二.反力式涡轮参数的选择

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

数理统计第二次大作业材料行业股票的聚类分析与判别分析 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

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

基于运动规划的惯性导航系统动态实验 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

北航_现代控制理论结课大作业

1. 控制系统任务的物理描述 为了满足飞机品质的要求，飞机的纵向运动和横侧向运动都需要有能够连续工作的阻尼器，以用来调整飞机的飞行姿态，避免其出现不必要的俯仰和倾斜。维持飞机纵向运动的阻尼器称为俯仰阻尼器，维持飞机横侧向运动的阻尼器称为偏航阻尼器。本次课程大作业旨在通过运用Matlab 的经典控制系统设计工具对某型飞机偏航阻尼器进行控制系统的设计。 2. 控制系统对象的数学模型 巡航状态下，某型飞机侧向运动的状态空间模型为： 111 12131411122212223242122131 3233343132234142434441424()1()()()()2()()()3()()4t x t a a a a b b t x t a a a a b b u t a a a a b b u t x t t a a a a b b x t t x x x x ??????????????????????????????????=+???????????????????????? ?????????? 111121314122122 2324234()()()()()()x t c c c c y t x t c c c c y t x t x t ??????????=?????????????? 式中： 1()x t ：侧滑角（单位为rad ） 2()x t ：偏航角速度（单位为/rad s ） 3()x t ：滚转角速度（单位为/rad s ） 4()x t ：倾斜角（单位为rad ） 输入向量及输出向量分别为： 1()u t ：方向舵偏角（单位为rad ） 2()u t ：副翼偏角（单位为rad ）

北航数理统计期末考试题

材料学院研究生会 学术部 2011年12月 2007-2008学年第一学期期末试卷 一、（6分，A 班不做）设x 1，x 2，…，x n 是来自正态总体2(,)N μσ的样本，令 )x x T -= ， 试证明T 服从t -分布t （2） 二、（6分，B 班不做）统计量F-F(n,m)分布，证明 111(,)F F n m αααα-的（0<<1）的分位点x 是。 三、（8分）设总体X 的密度函数为 其中1α>-，是位置参数。x 1，x 2，…，x n 是来自总体X 的简单样本，试求参数α的矩估计和极大似然估计。 四、（12分）设总体X 的密度函数为 1x exp x (;) 0 , p x μμσσσ??-? -≥??? =????? ，其它， 其中,0,μμσσ-∞<<+∞>已知，是未知参数。x 1，x 2，…，x n 是来自总体X 的简单样本。 （1）试求参数σ的一致最小方差无偏估计σ∧ ； （2）σ∧ 是否为σ的有效估计？证明你的结论。

五、（6分，A 班不做）设x 1，x 2，…，x n 是来自正态总体211(,)N μσ的简单样本，y 1，y 2，…，y n 是来自正态总体222(,)N μσ的简单样本，且两样本相互独立，其中221122,,,μσμσ是未知参数，2212σσ≠。为检验假设012112:, :,H H μμμμ=≠可令12, 1,2,..., , ,i i i z x y i n μμμ=-==-则上述假设检验问题等价于0111:0, :0,H H μμ=≠这样双样本检验问题就变为单检验问题。基于变换后样本z 1，z 2，…，z n ，在显著性水平α下，试构造检验上述问题的t-检验统计量及相应的拒绝域。 六、（6分，B 班不做）设x 1，x 2，…，x n 是来自正态总体20(,)N μσ的简单样本，0μ已知，2σ未知，试求假设检验问题 22220010:, :H H σσσσ≥<的水平为α 的UMPT 。 七、（6分）根据大作业情况，试简述你在应用线性回归分析解决实际问题时应该注意哪些方面？ 八、（6分）设方差分析模型为 总离差平方和 试求A E(S )，并根据直观分析给出检验假设012:...0P H ααα====的拒绝域形式。 九、（8分）某个四因素二水平试验，除考察因子A 、B 、C 、D 外，还需考察A B ?，B C ?。今选用表78(2)L ，表头设计及试验数据如表所示。试用极差分析指出因子的主次顺序和较优工艺条件。

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

实 验一 陀螺仪关键参数测试与分析实验 加速度计关键参数测试与分析实验 二零一三年五月十二日 实验一陀螺仪关键参数测试与分析实验 一、实验目得 通过在速率转台上得测试实验,增强动手能力与对惯性测试设备得感性认识；通过对陀螺仪测试数据得分析,对陀螺漂移等参数得物理意义有清晰得认识,同时为在实际工程中应用陀螺仪与对陀螺仪进行误差建模与补偿奠定基础。 二、实验内容 利用单轴速率转台,进行陀螺仪标度因数测试、零偏测试、零偏重复性测试、零漂测试实验与陀螺仪标度因数与零偏建模、误差补偿实验。 三、实验系统组成 单轴速率转台、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、引言 大型的复杂结构诸如飞机、汽车中的复杂部件及桥梁等大型工程的设计问题，依靠传统的经验和模拟实验的优化设计方法已难以胜任，拓扑优化方法成为解决该问题的关键手段。近年来拓扑优化的研究的热点集中在其工程应用上，如: 用拓扑优化方法进行微型柔性机构的设计，车门设计，飞机加强框设计，机翼前缘肋设计，卫星结构设计等。在其具体的操作实现上有两种方法，一是采用计算机语言编程计算，该方法的优点是能最大限度的控制优化过程，改善优化过程中出现的诸如棋盘格现象等数值不稳定现象，得到较理想的优化结果，其缺点是计算规模过于庞大，计算效率太低;二是借助于商用有限元软件平台。本文基于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插值模型。