2022年数字信号处理实验离散系统变换域分析 .pdf

《2022年数字信号处理实验离散系统变换域分析 .pdf》由会员分享，可在线阅读，更多相关《2022年数字信号处理实验离散系统变换域分析 .pdf（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、个人资料整理仅限学习使用实验 4 离散系统的变换域分析一、实验目的1、熟悉对离散系统的频率响应分析方法；2、加深对零、极点分布的概念理解。二、实验原理离散系统的是时域方程为其变换域分析方法如下：频域：系统的频率响应为：Z 域：系统的转移函数为：分解因式：其中和称为零、极点。三、实验要求：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 9 页个人资料整理仅限学习使用编程实现系统参数输入，绘出幅度频率响应和相位曲线和零、极点分布图。三、求系统num=0.0528 0.0797 0.1295 0.1295 0.797 0.0528 。den=

2、1 -1.8107 2.4947 -1.8801 0.9537 -0.2336 。z,p,k=tf2zp(num,den。disp(零点。disp(z。disp(极点。disp(p。disp(增益系数 。disp(k。figure(1 zplane(num,den figure(2 freqz(num,den,128 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 9 页个人资料整理仅限学习使用4.12 求滤波器输出的前十个采样值；num=0.3 -0.25 0.1。den=1 0 0。z,p,k=tf2zp(num,den。m=abs

3、(p。disp(零点。disp(z。disp(极点。disp(p。z,p,k=tf2zp(num,den。m=abs(p。disp(零点。disp(z。disp(极点。disp(p。subplot(2,2,1。zplane(num,den 。title(零 极点图 。k=256。w=0:pi/k:pi 。h=freqz(num,den,w。subplot(2,2,2。subplot(2,2,3。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 9 页个人资料整理仅限学习使用plot(w/pi,abs(h。grid title(幅度谱 xl

4、abel(omega/pi。ylabel(幅值 subplot(2,2,4。plot(w/pi,angle(h。grid title(相位谱 xlabel(omega/pi。ylabel(弧度 clear。B=0.3 -0.25 0.1。A=1。n=0:9。X=2*(n=0-(n=2 。Y=filter(A,B,X 。stem(n,Y,.。title(滤波器输出采样值 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 9 页个人资料整理仅限学习使用6.22 对于传输函数求脉冲响应和阶跃响应num=5。den=1 -1 0.09。subpl

5、ot(2,1,1 impz(num,den,40。h,t=impz(num,den,40。title(脉冲响应 step=ones(1,100 。subplot(2,1,2 y=conv(h,step。n=0:138。stem(n,y,filled title(阶跃响应 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页，共 9 页个人资料整理仅限学习使用7.26 滤波器的传输函数为： H 不计算 ，画出幅度响应曲线。k=256。num=0,0,1,1。den=1,-1,0,0。w=0:pi/k:pi 。h=freqz(num,den,w。su

6、bplot(2,2,1。plot(w/pi,real(h。grid title(实部。xlabel(omega/pi。ylabel(幅度 subplot(2,2,2。plot(w/pi,imag(h 。grid title(虚部。xlabel(omega/pi。ylabel(Amplitude subplot(2,2,3。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 9 页个人资料整理仅限学习使用plot(w/pi,abs(h。grid title(幅度谱 。xlabel(omega/pi。ylabel(幅值 subplot(2,2,

7、4。plot(w/pi,angle(h。grid title(相位谱 。xlabel(omega/pi。ylabel(弧度 7.28 由七项滑动平均滤波器的零极点推断滤波器形状解：由七项滑动滤波器特点，写出其对应传输函数num=1 1 1 1 1 1 1。den=7 0 0 0 0 0 0。z,p,k=tf2zp(num,den。disp(零点。disp(z。disp(极点。disp(z。disp(增益系数 。disp(k。figure(1 zplane(num,den 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页，共 9 页个人资料整理仅限学习使用figure(2 freqz(num,den,128 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页，共 9 页个人资料整理仅限学习使用精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共 9 页