数字信号处理实验离散系统变换域分析研究.docVIP

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 实验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) 版权申明 本文部分内容，包括文字、图片、以及设计等在网上搜集整理.版权为个人所有 This article includes some parts, including text, p