华南理工大学数字信号处理实验报告3(曹老师)..docVIP

一、实验目的 加深对LTI系统的理解以及分析。 二、实验原理 系统输入、输出关系可用以下差分方程描述： 系统响应为如下的卷积计算式： 当时，h[n]是有限长度的（n：[0，M]），称系统为FIR系统；反之，称系统为IIR系统。 系统的转移函数为 三、实验内容 1、用函数y=filter(p,d,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，用y=impz(p,d,N)求系统的冲激响应，再用卷积来计算任意信号作用于系统的响应。 求两个系统 各自的冲激响应，并且比较filter和conv函数的区别 实验代码如下： clear %离散时间序列x[n] n = 0:9; x = 5*exp(-n); subplot(4,2,1); stem(n,x) title(离散时间序列x[n]); %用filter函数滤波 a1 = [1 , 0.75 , 0.125]; b1 = [1 , -1]; y1 = filter(b1,a1,x); subplot(4,2,3); stem(n,y1) title(filter滤波1); a2 = [1]; b2 = [0 , 0.25 , 0.25 , 0.25 , 0.25]; y2 = filter(b2,a2,x); subplot(4,2,4); stem(n,y2) title(filter滤波2); %求系统的冲激响应 h1 = impz(b1,a1,10); subplot(4,2,5); stem(n,h1) title(冲激响应1); h2 = impz(b2,a2,10); subplot(4,2,6); stem(n,h2) title(冲激响应2); %用conv函数计算卷积 y3 = conv(x,h1); subplot(4,2,7); stem(y3) title(卷积1); y4 = conv(x,h2); subplot(4,2,8); stem(y4) title(卷积2); 实验结果如下： 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阶系统的串联。 求如下系统函数的零、极点，并将它转换成 实验代码如下： clear %系数向量 num = [1 , -0.1 , -0.3 , -0.3 , -0.2]; den = [1 , 0.1 , 0.2 , 0.2 , 0.5]; %求系统转移函数的零极点图 zplane(num,den) title(系统转移函数的零极点图); pause %求系统的频率响应 w = 0:pi/255:pi; h = freqz(num,den,w); subplot(2,1,1); plot(abs(h)) title(幅频响应); subplot(2,1,2); plot(angle(h)) title(相频响应); pause %求系统转换成1阶系统的并联 [r,p,k] = residuez(num,den); hn1 = freqz([r(1)],[1,-p(1)],w); hn2 = freqz([r(2)],[1,-p(2)],w); hn3 = freqz([r(3)],[1,-p(3)],w); hn4 = freqz([r(4)],[1,-p(4)],w); hn = hn1+hn2+hn3+hn4+k; subplot(2,1,1); plot(abs(hn)) title(1阶系统的并联的幅频响应); subplot(2,1,2); plot(angle(hn)) title(1阶系统的并联的相频响应); pause %求系统转换成2阶系统的串联 [z,p,K] = tf2zp(num,den); [sos,g] = zp2sos(z,p,K); hm1 = freqz([sos(1,1),sos(1,2),sos(1,3)],[1,sos(1,5),sos(1,6)],w); hm2 = freqz([sos(2,1),sos(2,2),sos(2,3)],[1,sos(2