縦の操舵応答

縦の運動方程式　

行列表示に変形

安定軸に対する安定微係数（可変特性研究機: ロッキードP2V-7をベース）

clear; close all;

g=32.174;%ヤード・ポンド法なので9.8ではない
Za=-236;
Zu=-0.227;
Xa=14.7;
Xu=-0.0215;
Madot=-0.28;
Ma=-3.78;
Mu=0;
Mq=-0.992;
Zde=-12.9;
Mde=-2.48;
U0=293.8;
Zq=-5.76;

A=[Xu Xa 0 -g;
    Zu/U0,Za/U0 (U0+Zq)/U0 0;
    Mu+Madot*Zu/U0 Ma+Madot*Za/U0 Mq+Madot*(U0+Zq)/U0 0;
    0 0 1 0];   
B=[0; Zde/U0; Mde+Madot*Zde/U0; 0];
C=[1 0 0 0;
   0 1 0 0;
   0 0 1 0 ;
   0 0 0 1];
D=[0;0;0;0];
sys=ss(A,B,C,D);


opts = stepDataOptions('InputOffset',0,'StepAmplitude',0.1);
step(sys,opts);
xlim([0,50]);
big;