[华工电信数字信号处理实验LTI系统应用.docVIP

一、实验题目 1、用函数y=filter(p,d,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，用y=impz(p,d,N)求系统的冲激响应，再用卷积来计算任意信号作用于系统的响应。 求两个系统 各自的冲激响应，并且比较filter和conv 函数的区别 2、用函数[z，p，K]=tf2zp（num，den）求得有理分式形式的系统转移函数的零、极点，用函数zplane（z，p）绘出零、极点分布图；也可以用函数zplane（num，den）直接绘出有理分式形式的系统转移函数的零、极点分布图。使h=freqz(num,den,w)函数可求系统的频率响应，w是频率的计算点，如w=0:pi/255:pi, h是复数，abs(h)为幅度响应，angle(h)为相位响应。另外，在MATLAB中，可以用函数 [r，p，k]=residuez（num，den）完成部分分式展开计算；可以用函数sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系统的串联。 求如下系统函数的零、极点，并且绘图，求系统的频率响应， 并将它转换成1阶系统的并联，2阶系统的串联 3、设计不同长度的滑动平均滤波器，且分析该系统的频率响应，分析滤波器长度对信号平滑效果的影响，对输入和输出延迟的影响。 二、实验过程 第一题： 对于系统： （1）、实验代码： a1=[1,0.75,0.125]; b1=[1,-1]; n=0:10; y=impz(b1,a1,10); figure; stem(y) title(y1impz) xlabel(x); ylabel(y); x1=[1 zeros(1,10)]; y1filter=filter(b1,a1,x1); figure(2); stem(n,y1filter); title(y1filter); xlabel(x); ylabel(y); x1=[1 zeros(1,10)]; h=impz(b1,a1,10); y1conv=conv(h,x1); n=0:19; figure(3); stem(n,y1conv,filled) title(y1conv) xlabel(x); ylabel(y); （2）、实验效果图： 对于系统： （1）、实验代码： a2=[1]; b2=[0.25,1,1,1]; n=0:10; y=impz(b2,a2,10); figure; stem(y) title(y2impz) xlabel(x); ylabel(y); x1=[1 zeros(1,10)]; y1filter=filter(b1,a1,x1); figure; stem(n,y1filter); title(y2filter); xlabel(x); ylabel(y); x1=[1 zeros(1,10)]; h=impz(b2,a2,10); y1conv=conv(h,x1); n=0:19; figure; stem(n,y1conv,filled) title(y2conv); xlabel(x); ylabel(y); （2）、实验效果图： 结果分析： y=impz(p,d,N)是用来实现冲激响应的，d表示差分方程输出y的系数，p表示输入x的系数，N表示冲激响应输出的序列个数。 y=filter(p,d,x)用来实现差分方程，d表示差分方程输出y的系数，p表示输入x的系数，而x表示输入序列。当输入为单位冲激函数时，输出为系统的冲激响应，输出结果长度数等于x的长度。 而y=conv(x,h)是用来实现卷级的，对x序列和h序列进行卷积，当输入为单位冲激函数时，输出为系统的冲激响应，输出的结果个数等于x的长度与h的长度之和减去1。 第二题： 实验代码 num=[1 -0.1 -0.3 -0.3 -0.2]; den=[1 0.1 0.2 0.2 0.5]; [z,p,k]=tf2zp(num,den); m=abs(p); disp(零点);disp(z); disp(极点);disp(p); disp(增益系数);disp(k); sos=zp2sos(z,p,k); disp(二阶系统串联);disp(real(sos)); disp(一阶系统并联); [r,p,k]=residuez(num,den) figure(1); zplane(num,den) figure(2); k=256; num=[1 -0.1 -0.3 -0.3 -0.2]; den=[1 0.1 0.2 0.2 0.5]; w=0:pi/k:pi; h=freqz(num,den,w); subplot(2,2,1); plot(w/pi,real