1.matlab求解系统响应
LTI系统的零状态响应的函数lsim,其调用格式
y=lsim(sys,f,t)
t表示计算系统响应的抽样点向量,f为系统输入信号,sys时LTI系统模型,用来表示微分方程
模型sys要借助tf函数获得,其调用方式为
sys=tf(b,a)
b,a分别为微分方程左端和右端的各项系数
系统的微分方程为
sys=tf([1],[1 22 77]);
t=0:0.01:5;
f=10*sin(2*pi*t);
y=lsim(sys,f,t);
plot(t,y)
xlabel('Time(sec)')
ylabel('y(t)')
2.卷积的matlab求解
matlab处理离散卷积的函数为conv(f1,f2),对序列做卷积运算。
处理连续信号的卷积时,需要对连续信号取相同的卷积步长,结果再乘以实际步长(连续信号的取样间隔)
t11=0;t12=1;t21=0;t22=2;
t1=t11:0.001:t12;
ft1=2*rectpuls(t1-0.5,1);
t2=t21:0.001:t22
ft2=t2;
t3=t11+t21:0.001:t12+t22;
ft3=conv(ft1,ft2);
ft3=ft3*0.001;
plot(t3,ft3)
title('ft1(t)*ft2(t)')