% 已知系统y''(t)+4y'(t)+3y(t)=2f'(t)+f(t),初始值为y(0_)=1,y'(0_)=2,
% 计算在输入信号为f(t)=u(t)时的系统零输入,零状态以及全响应。
clear all;
t=0:0.01:10;
a=[1 4 3];
b=[2 1];
sys=tf(b,a);
f=heaviside(t+0.01);
eq='D2y+4*Dy+3*y=0'; %定义符号微分方程表达式
cond='y(0)=1,Dy(0)=2'; %初始条件
z=dsolve(eq,cond);
z=simplify(z);
subplot(3,1,1);
ezplot(z,t);
title('零输入响应');
grid on;
y=lsim(sys,f,t);
subplot(3,1,2);
plot(t,y);
xlabel('t');
title('零状态响应');
grid on;
[A,B,C,D]=tf2ss(b,a);
sys1=ss(A,B,C,D);
zi=[1,2];
y1=lsim(sys1,f,t,zi);
subplot(3,1,3);
plot(t,y1);
xlabel('t');
title('全响应');
零输入和零状态的图形加起来应该是全响应。。为什么图形不是。。
% 计算在输入信号为f(t)=u(t)时的系统零输入,零状态以及全响应。
clear all;
t=0:0.01:10;
a=[1 4 3];
b=[2 1];
sys=tf(b,a);
f=heaviside(t+0.01);
eq='D2y+4*Dy+3*y=0'; %定义符号微分方程表达式
cond='y(0)=1,Dy(0)=2'; %初始条件
z=dsolve(eq,cond);
z=simplify(z);
subplot(3,1,1);
ezplot(z,t);
title('零输入响应');
grid on;
y=lsim(sys,f,t);
subplot(3,1,2);
plot(t,y);
xlabel('t');
title('零状态响应');
grid on;
[A,B,C,D]=tf2ss(b,a);
sys1=ss(A,B,C,D);
zi=[1,2];
y1=lsim(sys1,f,t,zi);
subplot(3,1,3);
plot(t,y1);
xlabel('t');
title('全响应');
零输入和零状态的图形加起来应该是全响应。。为什么图形不是。。