实验4 MATLAB在信号与系统中的应用
一、连续信号和系统
【例6.1】连续信号的MA TLAB描述
·程序
>> clear
>> t0=0;tf=5;dt=0.05;t1=1;
>> t=[t0:dt:tf];
>> st=length(t);
%(1)单位冲激函数
>> n1=floor((t1-t0)/dt);
>> x1=zeros(1,st);
>> x1(n1)=1/dt;
>> subplot(2,2,1),stairs(t,x1),grid on
>> axis([0,5,0,1.1])
%(2)单位阶跃函数
>> x2=[zeros(1,n1-1),ones(1,st-n1+1)];
>> subplot(2,2,3),stairs(t,x2),grid on
>> axis([0,5,0,1.1])
%(3)复指数函数
>> alpha=-0.5;w=10;x3=exp((alpha+j*w)*t);
>> subplot(2,2,2),plot(t,real(x3)),grid on
>> subplot(2,2,4),plot(t,imag(x3)),grid on
·总结
1、单位冲激函数无法直接描述,可以把其看成宽度是dt的,幅度是1/dt的矩形脉冲。
2、用stairs命令,可以显示连续信号波形中的不连续点;用plot命令,使波形更光滑。
3、复指数信号可以分解为余弦和正弦信号,它们分别是复数信号的实部和虚部。
【例6.2】LTI系统的零输入响应
a=input('输入分母系数向量a=[a1,a2,...]=');
n=length(a)-1;
Y0=input('输入初始条件向量Y0=[y0,Dy0,D2y0,...]=');
p=roots(a);
V=rot90(vander(p));
c=V\Y0';
dt=input('dt=');
tf=input('tf=');
t=0:dt:tf;
y=zeros(1,length(t));
for k=1:n y=y+c(k)*exp(p(k)*t);end
plot(t,y),grid on