matlab积分求idtft函数,MATLAB在信号与系统中的应用

实验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