飞行力学大作业 SY1105401 白斌
26
飞行力学大作业 SY1105401 白斌
结果分析：
1. 从小扰动和全量方程的阶跃输入图形可以看出，小扰动只是在小范围内的全量方程近似的结果。随着时间的增加，小扰动的近似结果与实际的结果相差越来越大。
2. 全量方程的结果输出图形表明，在升降舵和副翼阶跃输入的条件下，系统是发散的。即使有很小的扰动，系统也不稳定，最终会变得难以控制。而小扰动线化的结果表明，小扰动则对扰动具有较好的稳定性，相对全量方程较容易控制。
3. 由于小扰动做了近似处理，横纵向分开计算，因此导致了与全量方程的不相符。但实际情况是横纵向有着恨复杂的耦合，因此必须寻找更精确地简化近似，来得到更好的结果。
27
飞行力学大作业 SY1105401 白斌
五、总结与心得
经过将近两个月的努力，终于完成了本门课的大作业。由于自己本科学的是机械工程及自动化，没有航空的背景，所以感觉比较吃力。但是自己在这个过程中一直没有消极对待，而是积极主动地去学习，因此收获颇丰。在此，我把在完成本次课程大作业过程中的一些收获和体会总结如下，以期能对以后的学习研究提供很好的经验。
1.勤能补拙。尽管自己没有航空背景，但是通过自己努力付出，最终还是完成了本次的大作业。过程曲折艰难，但收获知识的同时，亦收获了自信。
2.提高了独自分析问题的能力，积累了一些科学研究的方法。本次作业过程当中，我依照老师给的提示，独立认真的圆满完成既定任务。遇到问题，我更倾向于独立思考，不轻易向他人寻求帮助，以锻炼自身的独立解决问题的能力。 3.熟悉了Matlab的操作。Matlab是一款功能强大的软件，熟悉它的相关操作，对我们今后的科研有很大的帮助。
28
飞行力学大作业 SY1105401 白斌
六、附录
1、计算配平迎角和舵偏角的程序。
Cm0=[0.205 0.081 -0.046 -0.174 -0.259 0.168 0.077 -0.02 -0.145 -0.202 0.186 0.107 -0.009 -0.121 -0.184 0.196 0.11 -0.005 -0.127 -0.193 0.213 0.11 -0.006 -0.129 -0.199 0.251 0.141 0.1 -0.102 -0.15 0.245 0.127 0.006 -0.097 -0.16 0.238 0.119 -0.001 -0.113 -0.167 0.252 0.133 0.014 -0.087 -0.104 0.231 0.108 0 -0.084 -0.076 0.198 0.081 -0.013 -0.069 -0.041 0.192 0.093 0.032 -0.006 -0.005 ]';
alpha0=[-10 -5 0 5 10 15 20 25 30 35 40 45 ];
elevator0=[-24 -12 0 12 24 ]';
A0=[0.77 0.241 -0.1 -0.416 -0.731 -1.053 -1.366 -1.646
29
飞行力学大作业 SY1105401 白斌
-1.917 -2.12 -2.248 -2.229 ];
for i=33.29:0.001:33.292 a=i
for j=0.711:0.001:0.732
A=interp1(alpha0,A0,i,'spline'); Cz=A-0.19*j/25.0;
Fz=1/2.0*Cz*0.7356*60^2*27.8709+9298.64*9.8*cos(i*3.14159/180); Cm=griddata(alpha0,elevator0,Cm0,i,j,'v4'); if abs(Fz)<1 && abs(Cm)<0.00001
dlmwrite('C:\\Users\\kotk\\Desktop\\trim2\\alpha.txt',i,'-append') dlmwrite('C:\\Users\\kotk\\Desktop\\trim2\\alpha.txt',' ','-append') dlmwrite('C:\\Users\\kotk\\Desktop\\trim2\\elevator.txt',j,'-append') dlmwrite('C:\\Users\\kotk\\Desktop\\trim2\\elevator.txt',' ','-append') dlmwrite('C:\\Users\\kotk\\Desktop\\trim2\\Cm.txt',Cm,'-append') dlmwrite('C:\\Users\\kotk\\Desktop\\trim2\\Cm.txt',' ','-append') dlmwrite('C:\\Users\\kotk\\Desktop\\trim2\\Cz.txt',Cz,'-append') dlmwrite('C:\\Users\\kotk\\Desktop\\trim2\\Cz.txt',' ','-append') dlmwrite('C:\\Users\\kotk\\Desktop\\trim2\\Fz.txt',Fz,'-append') dlmwrite('C:\\Users\\kotk\\Desktop\\trim2\\Fz.txt',' ','-append') end end end
2、S函数的程序
switch flag, case 0,
[sys,x0,str,ts]=mdlInitializeSizes; case 1,
sys=mdlDerivatives(t,x,u); case 2,
sys=mdlUpdate(t,x,u); case 3,
sys=mdlOutputs(t,x,u); case 4,
sys=mdlGetTimeOfNextVarHit(t,x,u);
30