信号与系统实验（2）

实验内容

a=[1 0 4];% 微分方程左侧的系数向量
n=length(a)-1;  %微分方程阶数
b=[1 1];
p=roots(a);  %特征方程求根
V=rot90(vander(p));  %生成范德蒙特矩阵
C=V\b'; %对应于各特征根系数
dt=0.1;tf=10;
t=0:dt:tf;y=zeros(1,length(t));
for k=1:n y=y+C(k)*exp(p(k)*t); end %对于有多个根用for循环方便，这个题也可以用上个方法
plot(t,y),title('零输入响应');xlabel('t'),ylabel('real(y)'),grid on

a=[1 0 4];  % 微分方程左侧的系数向量  
 b=[1];  %微分方程右侧的系数向量右
 [r,p]=residue(b,a);%residue函数求特征根p1、p2
 h=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t);%叠加各根分量
 dt=0.1; t=0:dt:10;
 x=cos(t);%任意输入信号
 y=conv(x,h)*dt;%求x和h的卷积，长度为2*length(t)-1
 plot(t,y(1:length(t))),title('零状态响应')
 %y(1:length(t))) 保证t和y长度相同

b=[1 6 4];  %微分方程右边系数向量
a=[1,7,10]; %微分方程左边系数向量
subplot(1,2,1)
impulse(b,a,0:0.1:11);% 绘出由向量a和b所定义的连续系统在t1 ~ t2时间范围内，并且以时间间隔p均匀取样的冲激响应的时域波形。
title('冲激响应');
subplot(1,2,2)
step(b,a,0:0.1:11);% 绘出由向量a和b所定义的连续系统在t1 ~ t2时间范围内，并且以时间间隔p均匀取样的阶跃响应的时域波形。
title('阶跃响应');