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;