离散系统的频域分析与零极点分布 1.docx

离散系统的频域分析与零极点分布 1.docx

《离散系统的频域分析与零极点分布 1.docx》由会员分享，可在线阅读，更多相关《离散系统的频域分析与零极点分布 1.docx（18页珍藏版）》请在冰点文库上搜索。

离散系统的频域分析与零极点分布1

课程设计报告

课程名称数字信号课程设计

系别：

工程技术系

专业班级：

电子信息工程0901

学号：

**

姓名：

**

课程题目：

离散系统的频域分析与零极点分布

完成日期：

2012年6月29日

指导老师：

**

2012年6月29日

课程设计目的

（1）掌握通过系统函数求解系统单位脉冲响应电脑方法；

（2）掌握系统函数零极点分布对频率特性的影响；

（3）通过系统函数零极点发布判断系统的稳定性。

课程设计要求

（1）简述通过系统函数零极点分布判断系统因果稳定性的方法；

（2）简述系统函数零极点分布对系统频率响应特性的影响；

（3）完成以上设计实验，并对结果进行分析和解释；

（4）打印程序清单和要求画出的信号波形；

（5）写出本次课程设计的收获和体会。

课程设计注意事项

遵守机房相关规定，按老师要求进行实习，不在实习过程中做与实习不相关的事情。

积极配合老师将实习课程完成，不在实习期间迟到，早退，旷课，严格按照要求和规定来做好这次实习。

课程设计内容

编制Matlab程序，完成以下功能，根据系统函数求出系统的零极点分布图并求解系统的单位脉冲响应；根据零极点分布图判断系统的稳定性；比较不同零极点发布对系统频率响应特性的影响；绘制相关信号的波形。

具体要求如下：

下面四种二阶网络的系统函数具有相同的极点发布：

（1）分别画出各系统的零极点分布图，并判断系统的稳定性；

（2）分别画出系统的幅频特性和相频特性曲线；

（3）分别求出系统的单位脉冲响应，并画出其波形。

课程设计简要操作步骤

（1）打开matlab2006b软件，创建一个M文件；

（2）用zplane函数求系统的零极点，输入MATLAB源程序，点击运行，绘制出零极点分布图，根据分布图判断系统的稳定性；

（3）利用函数freqz计算出

，然后利用函数abs和angle分别求出幅频特性与相频特性，最后利用plot命令绘出曲线；

（4）利用impz函数绘出该系统的单位取样响应，输入MATLAB源程序，绘制出系统单位取样响应的波形图；

（5）通过系统函数零极点分布判断系统因果稳定性，分析系统函数零极点分布对系统频率响应特性的影响。

课程设计心得体会

通过这次课程设计，我学到了很多的东西，不仅巩固了所学的知识，而且学到了很读在书本上学不到的东西。

通过这次课程设计使我懂得了理论与实际相结合的重要性，只有理论知识是远远不够的，只有把所学的理论知识与实际相结合起来，从理论中得出结论，才是真正的知识，才能提高自己的实际动手动手能力和独立思考的能力。

在设计的过程中遇到了各种各样的问题，同时在设计的过程中发现了自己的不足之处，对以前所学的知识理解不够深刻，掌握的不够牢固，通过这次课程设计，把以前所学过的知识重新温故，巩固了所学的知识。

课程设计评语及成绩

评语

成绩

指导教师

（签名）

年月日

离散系统的频域分析与零极点分布

摘要

本课题主要是根据系统函数求出系统的零极点分布图并且求解系统的单位脉冲响应,利用MATLAB软件绘制出系统零极点的分布图,根据零极点在单位圆的分布,判断因果系统的稳定性.再比较不同零极点对系统频率响应特性的影响。

从课题研究和设计过程当中对系统稳定性的判断有了清楚的认识,既极点在单位圆内,则该系统稳定,极点在单位圆外，则该系统为非稳定系统。

同时也对系统函数零极点分布对系统频率响应特性的影响有了深入的了解。

既极点位置主要影响频率响应的峰值及尖锐程度,零点位置主要影响频率响应的谷点位置及形状。

本次课题也对系统的幅频特性曲线和相频特性曲线进行了绘制,并求出了系统的单位脉冲响应以及绘制出了波形图。

关键字：

离散系统，频域分析，零极点分布

一、绪论

编制Matlab程序，完成以下功能，根据系统函数求出系统的零极点分布图，并求解系统的单位脉冲响应；根据零极点分布图判断系统的稳定性；比较不同零极点发布对系统频率响应特性的影响；绘制相关信号的波形。

具体要求如下：

下面四种二阶网络的系统函数具有相同的极点发布：

（1）分别画出各系统的零极点分布图，并判断系统的稳定性；

（2）分别画出系统的幅频特性和相频特性曲线；

（3）分别求出系统的单位脉冲响应，并画出其波形。

二、方案

实验原理

离散系统的时域方程为

其变换域分析方法如下：

频域

系统的频率响应为

Z域

系统的转移函数为

分解因式

，

其中

和

称为零、极点。

在MATLAB中，可以用函数[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阶系统的串联。

三、过程论述及结果分析

1.分别画出各系统的零极点分布图，并判断系统的稳定性

用zplane函数求系统的零极点，MATLAB源程序为：

B=[1]；

A=[1，-1.6，0.9425]；

Zplane（B,A）,gridon

legend（’零点’，‘极点’）

title（’零极点分布图’）

零极点分布图为：

由上图可知，极点都在圆内，所以该系统稳定。

用zplane函数求系统的零极点，MATLAB源程序为：

B=[1,-0.3]；

A=[1，-1.6，0.9425]；

Zplane（B,A）,gridon

legend（’零点’，‘极点’）

title（’零极点分布图’）

零极点分布图为：

由上图可知，极点都在圆内，所以该系统稳定。

用zplane函数求系统的零极点，MATLAB源程序为：

B=[1,-0.8]；

A=[1，-1.6，0.9425]；

Zplane（B,A）,gridon

legend（’零点’，‘极点’）

title（’零极点分布图’）

零极点分布图为：

由上图可知，极点都在圆内，所以该系统稳定。

用zplane函数求系统的零极点，MATLAB源程序为：

B=[1,-1.6，0.8]；

A=[1，-1.6，0.9425]；

Zplane（B,A）,gridon

legend（’零点’，‘极点’）

title（’零极点分布图’）

零极点分布图为：

由上图可知，极点都在圆内，所以该系统稳定。

2.分别画出系统的幅频特性和相频特性曲线

MATLAB源程序为：

b=[1];

a=[1-1.60.9425];

[HW]=freqz（b,a,400,’whole’）;

Hm=abs（H）;

Hp=angle（H）;

subplot=（211）

plot（w,Hm）,gridon

xlabel（’\omega（rad/s）’）,ylabel（’Magnitude’）

title（’离散系统幅频特性曲线’）

subplot（212）

plot（w,Hp）,gridon

xlabel（’\omega（rad/s）’）,ylabel（’phase’）

title（’离散系统相频特性曲线’）

波形图：

MATLAB源程序为：

b=[1-0.3];

a=[1-1.60.9425];

[HW]=freqz（b,a,400,’whole’）;

Hm=abs（H）;

Hp=angle（H）;

subplot=（211）

plot（w,Hm）,gridon

xlabel（’\omega（rad/s）’）,ylabel（’Magnitude’）

title（’离散系统幅频特性曲线’）

subplot（212）

plot（w,Hp）,gridon

xlabel（’\omega（rad/s）’）,ylabel（’phase’）

title（’离散系统相频特性曲线’）

波形图为：

MATLAB源程序为：

b=[1-0.8];

a=[1-1.60.9425];

[HW]=freqz（b,a,400,’whole’）;

Hm=abs（H）;

Hp=angle（H）;

subplot=（211）

plot（w,Hm）,gridon

xlabel（’\omega（rad/s）’）,ylabel（’Magnitude’）

title（’离散系统幅频特性曲线’）

subplot（212）

plot（w,Hp）,gridon

xlabel（’\omega（rad/s）’）,ylabel（’phase’）

title（’离散系统相频特性曲线’）

波形图为：

MATLAB源程序为：

b=[1-1.60.8];

a=[1-1.60.9425];

[HW]=freqz（b,a,400,’whole’）;

Hm=abs（H）;

Hp=angle（H）;

subplot=（211）

plot（w,Hm）,gridon

xlabel（’\omega（rad/s）’）,ylabel（’Magnitude’）

title（’离散系统幅频特性曲线’）

subplot（212）

plot（w,Hp）,gridon

xlabel（’\omega（rad/s）’）,ylabel（’phase’）

title（’离散系统相频特性曲线’）

波形图为：

3.分别求出系统的单位脉冲响应，并画出其波形

MATLAB源程序为：

a=[1-1.60.9425];

b=[1];

n=0:

30;

impz（b,a,30）,gridon

title（’系统单位取样响应’）

波形图为：

MATLAB源程序为：

a=[1-1.60.9425];

b=[1-0.3];

n=0:

30;

impz（b,a,30）,gridon

title（’系统单位取样响应’）

波形图为：

MATLAB源程序为：

a=[1-1.60.9425];

b=[1-0.8];

n=0:

30;

impz（b,a,30）,gridon

title（’系统单位取样响应’）

波形图为：

MATLAB源程序为：

a=[1-1.60.9425];

b=[1-1.60.8];

n=0:

30;

impz（b,a,30）,gridon

title（’系统单位取样响应’）

波形图为：

四、结论

根据系统函数的极点在图上的分布来判断系统的稳定性,如果极点在单位圆内，那该系统为稳定系统,如果极点在单位圆外，那该系统为非稳定系统。

要获得系统函数H（z）的零极点分布图，可直接应用zplane函数，其语句格式为zplane（A,B）。

其中，B与A分别表示H（z）的分子和分母多项式的系数向量。

它的作用是在Z平面上画出单位圆、零点与极点。

分析信号的方法有两种，一种是时域分析法，一种是频域分析法。

频域分析法是研究控制系统的一种经典方法，是在频域内应用图解分析法评价系统性能的一种工程方法。

该方法是以输入信号的频率为变量，对系统的性能在频率域内进行研究的一种方法。

频域分析法不必直接求解系统的微分方程，而是间接地揭示系统的时域性能，它能方便的显示出系统参数对系统性能的影响，并可以进一步指明如何设计校正。

这种分析方法有利于系统设计，能够估计到影响系统性能的频率范围。

幅频特性,它指的是当ω由0到∞变化时，|G（jω）|的变化特性，记为A（ω）。

相频特性,它指的是当ω由0到∞变化时，∠G（jω）的变化特性称为相频特性，记为（ω）。

离散系统的幅频特性曲线和相频特性曲线直观的反映了系统对不同频率的输入序列的处理情况。

致谢

感谢指导老师杨亚东老师在这次课程设计过程当中对我的悉心指导和讲解，使我在这次课程设计过程中学到了很多课本上学不到的知识，使我对数字信号处理这门课有了更加深入的了解和掌握，相信通过这次课程设计，我的专业课知识会更上一层楼，真心感谢杨老师！

参考文献

[1]张威.MATLAB基础与编程入门（第二版）.西安电子科技大学出版社

[2]高西全.丁玉美.数字信号处理（第三版本）西安电子科技大学出版社

[3]杨亚东.数字信号处理课程设计指导书