www.gusucode.com > MATLAB仿真技术与应用__的配套例题和习题的matlab源代码 > 第七章/example10.m

    %MATLAB program 7-10
%State space mode
clear all;
clc;
A=[-21 –9.25 –3.28 –1.56;16 0 0 0;0 8 0 0;0 0 2 0];    %输入状态方程中的各矩阵
B=[2;0;0;0];
C=[0 1.56 1.17 1.56]; 
D=0.5;
[num1,den1]=ss2tf(A,B,C,D);           %状态空间模型转换成传递函数模型
G1=tf(num1,den1);
G=feedback(G1,1);                    %反馈联接
[num,den]=tfdata(G,'v')
[aa,bb,cc,dd]=tf2ss(num,den)            %传递函数模型转换成状态空间模型
sys=ss(aa,bb,cc,dd)
t=0:0.01:10;                           %时间
u=[0,ones(1,1000)];                     %阶跃信号
lsim(sys,u,t);                          %仿真
grid;