搜索
信号与系统实验报告实验六：离散信号与系统的时域分析一、实验目的1．学会用MATLAB表示常用离散信号的方法；2．学会用MATLAB实现离散信号卷积的方法；3．学会用MATLAB求解离散系统的单位响应；4．学会用MATLAB求解离散系统的零状态响应；二、实验原理1．离散信号的MATLAB表示表示离散时间信号f(k)需要两个行向量，一个是表示序号k=[]，一个是表示相应函数值f=[]，画图命令是stem。2.离散信号的卷积和两个有限长序列f1，f2卷积可调用MATLAB函数conv，调用格式是f=conv(f1,f2),f是卷积结果，但不显示时间序号，可自编一个函数dconv给出f和k，并画图。function[f,k]=dconv(f1,f2,k1,k2)3．离散系统的单位响应MATLAB提供画系统单位响应函数impz，调用格式是impz(b,a)式中b和a是表示离散系统的行向量；impz(b,a,n)式中b和a是表示离散系统的行向量，时间范围是0~n；impz(b,a,n1,n2)时间范围是n1~n2；y=impz(b,a,n1,n2)由y给出数值序列；4．离散系统的零状态响应MATLAB提供求离散系统零状态响应数值解函数filter，调用格式为filter(b,a,x),式中b和a是表示离散系统的向量，x是输入序列非零样值点行向量，输出向量序号同x一样。三实验验证（截图）2.3.4.5.6.7.8.解答：代码a=[2,-2,1];b=[1,3,2];impz(b,a)impz(b,a,60)impz(b,a,-10:40)图形如下：2.已知y(k)+y(k-1)+0.25*y(k-2)=f(k),输入f(t)=e(k),画输出波形，范围0~15。解答：代码：a=[110.25];b=[1];t=0:20;x=heaviside(t);y=filter(b,a,x)subplot(2,1,1)stem(t,x)title('输入序列')subplot(2,1,2)stem(t,y)title('响应序列')图形如下：