原创力文档- 1、本文档共9页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
实验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)
文档评论（0）