matlab simulink 飞机飞行状态控制

Posted 2022-08-08 studyer_domi

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、参考论文

略