数字信号处理20104200王大元实验4案例.docVIP

实验4 离散系统的变换域分析 实验目的 熟悉对离散系统的频率响应分析方法； 加深对零、极点分布的概念理解。 实验原理 离散系统的是时域方程为 其变换域分析方法如下： 频域： 系统的频率响应为： Z域： 系统的转移函数为： 分解因式： 其中 和 称为零、极点。 三、实验要求： 编程实现系统参数输入，绘出幅度频率响应和相位曲线和零、极点分布图。 求系统 num [0.0528 0.0797 0.1295 0.1295 0.797 0.0528]; den [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 4.12 求滤波器输出的前十个采样值； num [0.3 -0.25 0.1]; den [1 0 0]; [z,p,k] tf2zp num,den ; m abs 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 ; plot w/pi,abs h ;grid title 幅度谱 xlabel \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 滤波器输出采样值 6.22 对于传输函数 求脉冲响应和阶跃响应 num [5]; den [1 -1 0.09]; subplot 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 阶跃响应 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 ; subplot 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 ; plot w/pi,abs h ;grid title 幅度谱 ; xlabel \omega/\pi ;ylabel 幅值 subplot 2,2,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 figure 2 freqz num,den,128