ex7_3.m

来自「华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 」· M 代码 · 共 19 行

% ex7_3
% 计算例7-2的传递函数矩阵及状态和输出的响应
figure('pos',[50,50,200,150],'color','w');
a=[0 1;-2 -3];
b=[0;1];
c=[1.5 0.5];
d=0;
syms s t;
% G1=inv(s*diag(diag(ones(size(a))))-a)*b;
G1=inv(s*eye(size(a))-a)*b
x=ilaplace(G1/s),
y=ilaplace(c*G1/s),
for I=1:61;
    tt=0.1*(I-1);
    xv(:,I)=subs(x(:),'t',tt);
    yv(I)=subs(y,'t',tt);
end;
plot(0:60,[xv;yv],'k');