clear;
clc;
close all;
A = [-2,-2.5,-0.5;
1,0,0;
0,1,0];
B = [1;0;0];
C = [0 1.5 1];
D = [0];
sys = ss(A,B,C,D);
x0 = [1,0,2];
t= 0:0.01:20;
b1=0; %信号发生突变的时刻
b2=2; %信号发生突变的时刻
u=(2*stepfun(t,b1)-1.5*stepfun(t,b2))';%产生单位的阶跃信号
t = t';
figure(1);
step(sys)
figure(2);
impulse(sys)
[y,t,x] = lsim(sys,u,t,x0);
figure(3);
grid on
plot(t,x(:,1),'r')
hold on
plot(t,x(:,2),'b')
hold on
plot(t,x(:,3),'g')
xlabel("幅值");
ylabel("时间/秒");
title('状态全响应');
figure 使用为figure(数字),lsim函数使用