You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
100 lines
4.4 KiB
Plaintext
100 lines
4.4 KiB
Plaintext
%---------------------------------------------------------------%
|
|
%---Plot time history of the satellite attitude and its affix---%
|
|
%---------------------------------------------------------------%
|
|
if(exist('a_sat','var')==1)
|
|
%-----------------卫星姿态变化曲线---------------------------%
|
|
figure('Name','卫星姿态变化曲线','position',[100 100 600 400],'NumberTitle','off');
|
|
|
|
subplot(2,1,1);plot(a_sat(:,1),a_sat(:,2),'-',...
|
|
a_sat(:,1),a_sat(:,3),'--',...
|
|
a_sat(:,1),a_sat(:,4),':','LineWidth',1.5);
|
|
legend('滚动轴','俯仰轴','偏航轴');
|
|
grid on;
|
|
ylabel('星体角速度\omega_s(deg/s)');
|
|
|
|
subplot(2,2,2);plot(a_sat(:,1),a_sat(:,2),'-',...
|
|
a_sat(:,1),a_sat(:,3),'--',...
|
|
a_sat(:,1),a_sat(:,4),':','LineWidth',1.5);
|
|
legend('滚动轴','俯仰轴','偏航轴');
|
|
grid on;
|
|
ylabel('星体姿态角(deg)');
|
|
xlabel('t(s)','fontsize',10,'fontname','华文中宋');
|
|
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %----------------挠性模态及力矩变化曲线---------------------%
|
|
% figure('Name','帆板1挠性模态及力矩曲线','position',[600 100 600 400],'NumberTitle','off');
|
|
%
|
|
% subplot(2,1,1);plot(eta1(:,1),eta1(:,2),'-',...
|
|
% eta1(:,1),eta1(:,3),'--',...
|
|
% eta1(:,1),eta1(:,4),':','LineWidth',1.5);
|
|
% hold on;
|
|
% plot(eta1(:,1),eta1(:,5),'-',...
|
|
% eta1(:,1),eta1(:,6),'--','LineWidth',1);
|
|
% ylabel('+Y帆板模态\eta_1');
|
|
% legend('1阶','2阶','3阶','4阶','5阶');
|
|
%
|
|
% subplot(2,1,2);plot(T_f1(:,1),T_f1(:,2),'-',...
|
|
% T_f1(:,1),T_f1(:,3),'--',...
|
|
% T_f1(:,1),T_f1(:,4),':','LineWidth',1.5);
|
|
% ylabel('+Y帆板振动耦合力矩T_f_1(Nm)');
|
|
% legend('滚动轴','俯仰轴','偏航轴');
|
|
% xlabel('t(s)','fontsize',13,'fontname','华文中宋');
|
|
% %----------------挠性模态及力矩变化曲线---------------------%
|
|
% figure('Name','帆板2挠性模态及力矩曲线','position',[600 100 600 400],'NumberTitle','off');
|
|
% subplot(2,1,1);plot(eta2(:,1),eta2(:,2),'-',...
|
|
% eta2(:,1),eta2(:,3),'--',...
|
|
% eta2(:,1),eta2(:,4),':','LineWidth',1.5);
|
|
% hold on;
|
|
% plot(eta2(:,1),eta2(:,5),'-',...
|
|
% eta2(:,1),eta2(:,6),'--','LineWidth',1);
|
|
% ylabel('-Y帆板模态\eta_2');
|
|
% legend('1阶','2阶','3阶','4阶','5阶');
|
|
%
|
|
% subplot(2,1,2);plot(T_f2(:,1),T_f2(:,2),'-',...
|
|
% T_f2(:,1),T_f2(:,3),'--',...
|
|
% T_f2(:,1),T_f2(:,4),':','LineWidth',1.5);
|
|
% ylabel('-Y帆板振动耦合力矩T_f_2(Nm)');
|
|
% legend('滚动轴','俯仰轴','偏航轴');
|
|
% xlabel('t(s)','fontsize',13,'fontname','华文中宋');
|
|
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
% %----------------帆板工况及转动耦合力矩---------------------%
|
|
% figure('Name','帆板1帆板工况及转动耦合力矩','position',[600 100 600 400],'NumberTitle','off');
|
|
% subplot(2,1,1);plot(w_sp(:,1),w_sp(:,2),'-','LineWidth',1.5);
|
|
% ylabel('+Y帆板工况\omega_a_1(deg/s)');
|
|
% legend('+Y Solar panel');
|
|
%
|
|
% subplot(2,1,2);plot(T_cou1(:,1),T_cou1(:,2),'-',...
|
|
% T_cou1(:,1),T_cou1(:,3),'--',...
|
|
% T_cou1(:,1),T_cou1(:,4),':','LineWidth',1.5);
|
|
% ylabel('+Y帆板转动耦合力矩T_c_o_u_1(Nm)');
|
|
% legend('滚动轴','俯仰轴','偏航轴');
|
|
% xlabel('t(s)','fontsize',13,'fontname','华文中宋');
|
|
% %----------------帆板工况及转动耦合力矩---------------------%
|
|
% figure('Name','帆板2帆板工况及转动耦合力矩','position',[600 100 600 400],'NumberTitle','off');
|
|
% subplot(2,1,1);plot(w_sp(:,1),w_sp(:,3),'-','LineWidth',1.5);
|
|
% ylabel('-Y帆板工况\omega_a_2(deg/s)');
|
|
% legend('-Y Solar panel');
|
|
%
|
|
% subplot(2,1,2);plot(T_cou2(:,1),T_cou2(:,2),'-',...
|
|
% T_cou2(:,1),T_cou2(:,3),'--',...
|
|
% T_cou2(:,1),T_cou2(:,4),':','LineWidth',1.5);
|
|
% ylabel('-Y帆板转动耦合力矩T_c_o_u_2(Nm)');
|
|
% legend('滚动轴','俯仰轴','偏航轴');
|
|
% xlabel('t(s)','fontsize',13,'fontname','华文中宋');
|
|
%
|
|
% %-----------------力矩变化曲线---------------------------%
|
|
% figure('Name','力矩变化曲线','position',[600 100 600 400],'NumberTitle','off');
|
|
%
|
|
% subplot(2,1,1);plot(T_c(:,1),T_c(:,2),'-',...
|
|
% T_c(:,1),T_c(:,3),'--',...
|
|
% T_c(:,1),T_c(:,4),':','LineWidth',1.5);
|
|
% legend('滚动轴','俯仰轴','偏航轴');
|
|
% ylabel('控制力矩T_c(Nm)');
|
|
%
|
|
%
|
|
% subplot(2,1,2);plot(T_p1(:,1),T_p1(:,3),'LineWidth',1.5);
|
|
% ylabel('+Y帆板驱动力矩T_p_1(Nm)');
|
|
% xlabel('t(s)','fontsize',10,'fontname','华文中宋');
|
|
|
|
|
|
else
|
|
warndlg('^_^画图前别忘了仿真^_^','错误');
|
|
end |