%---------------------------------------------------------------%
%---Plot time history of the satellite attitude and its affix---%
%---------------------------------------------------------------%
if(exist('w_sat','var')==1)
%-----------------卫星姿态变化曲线---------------------------%
% figure('Name','卫星姿态变化曲线','position',[100 100 600 400],'NumberTitle','off');
% 
% subplot(2,1,1);plot(w_sat(:,1),w_sat(:,2),'-',...
%                     w_sat(:,1),w_sat(:,3),'--',...
%                     w_sat(:,1),w_sat(:,4),':','LineWidth',1.5);
% legend('滚动轴','俯仰轴','偏航轴');
% grid on;
% ylabel('星体角速度\omega_s(deg/s)');
% 
% subplot(2,1,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