3.4.3 時間応答 ― MATLAB


 制御対象の状態空間表現の定義 ...... M ファイル plant.m (p.59)
clear
format compact

A = [  0  1
     -10 -2 ];
B = [ 0
      1 ];
C = [ 1  0 ];
D = 0;
sysP = ss(A,B,C,D);

 シミュレーション結果の表示 ...... M ファイル plot_data.m (p.59)
plot(t,y);
grid;
xlabel('t [s]');
ylabel('y(t)');


(c) 任意の時間応答


 任意の時間応答の描画 ...... M ファイル response_lsim.m (p.62)
plant
    
t = 0:0.01:5;
u(1:200)   = t(1:200)/2;
u(201:501) = ones(1,301);
x0 = [ 0.1  0 ]';

y = lsim(sysP,u,t,x0);
plot_data
>> response_lsim

 任意の時間応答の描画 ...... M ファイル response_lsim2.m (本には記載していない)
plant
    
t = 0:0.01:5;
u(1:200)   = t(1:200)/2;
u(201:501) = ones(1,301);
x0 = [ 0.1  0 ]';

lsim(sysP,u,t,x0);
>> response_lsim2

 任意の時間応答の描画(グラフのカスタマイズ) ...... M ファイル response_lsim3.m(本には記載していない)
plant
    
t = 0:0.01:5;
u(1:200)   = t(1:200)/2;
u(201:501) = ones(1,301);
x0 = [ 0.1  0 ]';

y = lsim(sysP,u,t,x0);

% --- グラフのカスタマイズ --------------------------------
figure(1)
plot(t,y,'r','linewidth',2);
grid;

xlim([0 5]);
ylim([-0.1 0.2]);

set(gca,'xtick',[0:1:5]);
set(gca,'ytick',[-0.1:0.1:0.2]);

set(gca,'fontname','times','fontsize',20);

xlabel('time [s]','fontname','times','fontsize',22);
ylabel('{\it{y}}({\it{t}})','fontname','times','fontsize',22);


% --- グラフのカスタマイズ --------------------------------
figure(2)
plot(t,u,'r','linewidth',2);
grid;

xlim([0 5]);
ylim([0 1.5]);

set(gca,'xtick',[0:1:5]);
set(gca,'ytick',[0:0.5:1.5]);

set(gca,'fontname','times','fontsize',20);

xlabel('time [s]','fontname','times','fontsize',22);
ylabel('{\it{u}}({\it{t}})','fontname','times','fontsize',22);
>> response_lsim3