matlab simulink 飞机飞行状态控制

1、内容简介

略
474-可以交流、咨询、答疑

2、内容说明

略

3、仿真分析

close all;
figure(1);

plot(t,gangti(:,3)*180/3.14,'k-',t,ditong(:,1)*180/3.14,'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
xlabel('{\itt}/s')
ylabel('{\it\gamma}/(\circ)')
legend('\it{\gamma}','{\it\gamma}_{d}')

figure(2);
plot(t,gangti(:,4)*180/3.14+gangti(:,3)*180/3.14,'k-',t,ditong(:,2)*180/3.14,'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
ylabel('{\it\theta} /(\circ)');
legend('\it\theta','{\it\theta}_{d}');
xlabel('{\itt}/s')

figure(3);
plot(t,gangti(:,5)*180/3.14,'k-',t,ditong(:,3)*180/3.14,'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
ylabel('{\itQ}/({\circ}/s)');
legend('\itQ','{\itQ}_{d}');
xlabel('{\itt}/s')

figure(4);
plot(t,gangti(:,4)*180/3.14,'k-','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
ylabel('{\itα}/(\circ)');
xlabel('{\itt}/s')

figure(5);
plot(t,cankao(:,2),'k-',t,gangti(:,1),'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
legend('{\itV}_{ref}','{\itV}');
ylabel('{\itV}/(ft/s)');xlabel('{\itt}/s')

figure(6);
plot(t,gangti(:,1)-cankao(:,2),'k-','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
xlabel('{\itt}/s');ylabel('Velocity tracking error /(ft/s)');

figure(7);
plot(t,cankao(:,1),'k-',t,gangti(:,2),'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
legend('{\ith}_{ref}','\ith');
ylabel('{\ith}/m');xlabel('{\itt}/s')

figure(8);
plot(t,gangti(:,2)-cankao(:,1),'k-','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
xlabel('{\itt}/s');ylabel('Altitude tracking error /ft');

figure(9);
plot(t,u(:,1),'k-','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
ylabel('\itΦ');xlabel('{\itt}/s')

figure(10);
plot(t,u(:,2)*180/3.14,'k-','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
xlabel('{\itt}/s');ylabel('{\itδ}_e/(\circ)');

figure(11);
plot(t,eat(:,1)*0.3048,'k-','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
ylabel('{\it\eta}_1');xlabel('{\itt}/s')

figure(12);
plot(t,eat(:,2)*0.3048,'k-','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
ylabel('{\it\eta}_2');xlabel('{\itt}/s')

figure(13);
subplot(211);
plot(t,jianmo(:,1),'k-',t,gangti(:,1),'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
legend('{\itV}_{估计}','{\itV}');
ylabel('{\itV}/(ft/s)');xlabel('{\itt}/s')
subplot(212);
plot(t,jianmo(:,1)-gangti(:,1),'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
legend('{\itV}_{估计误差}');
ylabel('{\itV}/(ft/s)');xlabel('{\itt}/s')

figure(14);
subplot(211);
plot(t,jianmo(:,2),'k-',t,gangti(:,3),'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
legend('{\it\gamma}_{估计}','{\it\gamma}');
ylabel('{\it\gamma}');xlabel('{\itt}/s')
subplot(212);
plot(t,jianmo(:,2)-gangti(:,3),'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
legend('{\it\gamma}_{估计误差}');
ylabel('{\it\gamma}');xlabel('{\itt}/s')

figure(15);
subplot(211);
plot(t,jianmo(:,3),'k-',t,gangti(:,5),'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
legend('{\itQ}_{估计}','{\itQ}');
ylabel('{\itQ}');xlabel('{\itt}/s')
subplot(212);
plot(t,jianmo(:,3)-gangti(:,5),'r-.','linewidth',2);
set(gca,'FontSize',12,'FontName','Times New Roman')
legend('{\itQ}_{估计误差}');
ylabel('{\itQ}');xlabel('{\itt}/s')

 

zhuang 

4、参考论文

略