数字信号系统的时域、频率特性
分析
定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析
与验证课程
设计
领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计
.
燕山大学 课 程 设 计 说 明 书
设计题目:数字信号系统的时域、频率特性分析与验证 系 别:电气工程系 年级专业:09级工业自动化仪表 学 号: 学生姓名: 指导教师: 教师职称:
.
.
燕山大学《课程设计》任务书
课程名称: 信号处理原理与应用 基层教学单位:电气工程系 指导教师:
(专业)班学号 学生姓名 09仪表1班 级 设计题数字信号系统的时域、频率特性分析与验证 目
设
计
1、 ynynxnxn,,,,,0.51()0.5(1),,,,技 1111
术 2、 ynynxn()2(1)(),,,222
参
数
设 利用Matlab软件稳定和非稳定的数字信号系统进行时域分析(包括单位
计 脉冲响应,单位阶跃响应)和频率特性分析(包括幅频特性,相频特性),
要 其中涉及的参数根据需要自主选择
求
参 [1]谢平,王娜.信号处理原理及应用.北京:机械工业出版社,2008,10 考 [2]楼顺天.MATLAB 5.x程序设计语言.西安:西安电子科技大学出版社 资 [3]薛年喜.MATLAB在数字信号处理中的应用(第二版).清华大学出版社 料
周次 第一周
应 完成全部
方案
气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载
设计: 完 周一、二:分析题目,查阅相关资料,熟悉MATLAB程序设计方法。 成 周二至周五:方案设计,编程,调试,运行,最后完善方案。 内 周六、日:完成课程设计
报告
软件系统测试报告下载sgs报告如何下载关于路面塌陷情况报告535n,sgs报告怎么下载竣工报告下载
书,进行答辩。 容
基层教学指导教 单位主任师签字 签字
说明
关于失联党员情况说明岗位说明总经理岗位说明书会计岗位说明书行政主管岗位说明书
:1、此表一式三份,系、学生各一份,报送院教务科一份。
2、学生那份任务书要求装订到课程设计报告前面。
电气工程学院 教务科
.
.
燕山大学课程设计评审意见表
指导教师评语:
正确完善完善
认真较为合理合理
较认真理论分析方法设计工作态度一般一般
不认真较差较差 成绩:
指导教师:
年 月 日
答辩小组评语:
正确清晰
基本掌握设计结论基本正确
原理了解
不正确不清楚 成绩:
评阅人:
年 月 日 课程设计总成绩:
答辩小组成员签字:
年 月 日
.
.
目 录
第一章 摘要……………………………………1 第二章 引言……………………………………2 第三章 设计数字信号处理系统………………3 第四章 数字信号系统的时域分析……………3
4.1 稳定的数字信号系统分析……………3
4.2 不稳定的数字信号系统分析…………6 第五章 数字信号系统的频率特性分析 ……8
5.1 稳定的数字信号系统分析……………8
5.2 不稳定的数字信号系统分析…………9 第六章 心得体会………………………………11 第七章 参考文献………………………………12
.
.
数字信号系统的时域、频率特性分析与验证
一、摘要
对于系统的时域、频率特性分析与验证问题,本文先设计了两个离散的数字信号处理系统,分别为稳定系统和不稳定系统,用差分方程来描述这两个系统。
1、时域分析:
对于系统的时域分析,本文从差分方程的单位阶跃响应和单位脉冲响应进行分析。运用MATLAB,分别绘制出相应的响应曲线,观察稳定系统与不稳定系统的区别。
2、频率特性:
对于系统的频率特性,本文首先分别对稳定系统和不稳定系统的差分方程进行Z变换,求出它们的传递函数,推到出有理分式,运用MATLAB,绘制出相应的幅频特性曲线和相频特,,Hjw
性曲线,再对稳定系统和不稳定系统进行比较和分析。
关键词:差分方程,阶跃,脉冲,Z变换,幅频特性,相频特性
.
.
二、引言
MATLAB是矩阵实验室的简称,是美国MathWorks公司出品的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技算法、创建用户界面、连接其他编程语言的程序等,主要应用于工程计算、控术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。MATLAB可以进行矩阵运算、绘制函数和数据、实现制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。
数字信号处理系统处理的直接对象是数字信号,处理的方式是数值运算的方式,相对于模拟信号的处理,具有灵活性、精度高、稳定性高的、便于大规模集成等优点。
本文分别从数字信号系统的时域和频率两方面进行分析,并结合实例,运用MATLAB,绘制出相应的特性曲线,简单明了,为今后的学习打下坚实的基础。
.
.
三、设计数字信号处理系统
对于数字信号处理系统,我们用差分方程来描述,设和,,xn
N是阶系统的输入序列和输出序列,则系统的差分方程可表,,yn
示为: MN
(1) ,,,,,,yn,bxn,i,ayn,i,,ii,0,0ii或者 NM
(2) ,,,,ayn,i,bxn,i,a,1,,0ii,,00ii当系统为稳定系统时,满足条件:
,1、有限输入产生有限输出
,,, (3) hn2、(),,,,i=-,,
,Hz3、的收敛域包含单位圆,,,
所以我们可以根据条件设计出稳定系统1和不稳定系统2的差分方程分别如下列式(4),式(5)所示:
ynynxnxn,,,,,0.51()0.5(1) (4) ,,,,1111
(5) ynynxn()2(1)(),,,222
四、数字信号系统的时域分析
时域分析是通过对离散时间系统的数学模型——差分方程进行单位阶跃响应和单位脉冲响应两方面进行分析。 4.1、稳定的数字信号系统分析
4.1.1、单位阶跃响应
当给稳定的系统一个阶跃信号,结合公式(4)构造的un()
系统1,可用递推的方法求出系统的单位节约响应。或者,我们可以直接运用MATLAB来进行求解,绘制出其单位节约响应曲线。
首先我们构造一个阶跃信号stepseq:
function [x,n]=stepseq(n0,n1,n2)
n=[n1:n2];
x=[(n-n0)>=0];
.
.
x=double(x);
将该信号带入式(4)构造的系统中,便可得到系统的阶跃响应
曲线,程序如下:
clear
clc
a=[1,-0.5];
b=[1,0.5];
n=[-10:90];
[x,n]=stepseq(0,-10,90);
h=filter(b,a,x);
stem(n,h);
title('阶跃响应');
ylabel('h(n)');
axis([-10,90,-0.2,3.5]);
xlabel('n');
运用MATLAB,结果如下:
阶跃响应?3.5
3
2.5
2
h(n)1.5
1
0.5
0
-100102030405060708090n
图1
4.1.2、单位脉冲响应
同理,构造一个脉冲信号impseq,如下所示:
function [x,n]=impseq(n0,n1,n2)
n=[-20:120];
x=[(n-0)==0];
.
.
x=double(x);
并将该信号带入系统1中,便可得到其脉冲响应曲线,程序如下
所示:
clear
clc
a=[1,-0.5];
b=[1,0.5];
n=[-20:120];
[x,n]=impseq(0,-20,120);
h=filter(b,a,x);
stem(n,h);
title('脉冲响应');
ylabel('h(n)');
axis([-20,120,-1.1,1.1]);
xlabel('n');
经MATLAB运行后,得到去脉冲响应曲线为:
脉冲响应
1
0.8
0.6
0.4
0.2
0h(n)
-0.2
-0.4
-0.6
-0.8
-1
-20020406080100120n
图2
n由图1,我们可以发现对于稳定的系统,当大于一定值之后
n时,其阶跃响应曲线基本不变,趋于3。而由图2知,当大于一
.
. 定值时,其脉冲响应曲线响应趋于0. 4.2、不稳定的数字信号系统分析
4.2.1单位阶跃响应
对于该系统的单位节约响应,我们可以运用4.1.1中的程序进
行运行,只要把响应的系统2参数改为,改变一下ab,,,[1,2],[1]
的取值范围及坐标轴的范围定义,如下所示: n
clear
clc
a=[1,-2];
b=[1];
n=[-10:90];
[x,n]=stepseq(0,-10,90);
h=filter(b,a,x);
stem(n,h);
title('阶跃响应');
ylabel('h(n)');
xlabel('n');
得到相应的响应曲线为:
27阶跃响应x 102.5
2
1.5
h(n)
1
0.5
0-100102030405060708090n 图3
n当的值改为由0到290时,可得其阶跃响应曲线为:
.
.
87阶跃响应x 104
3.5
3
2.5
2h(n)
1.5
1
0.5
0-50050100150200250300n 图4
由图3和图4可知,不稳定系统的阶跃响应曲线在趋于时,n,
的。 其响应曲线也是趋于,
4.2.2单位脉冲响应
同理,运用4.1.2中的程序,将系统参数相应的改为
,改变相应的坐标得: ab,,,[1,2],[1]
clear
clc
a=[1,-2];
b=[1];
n=[-20:290];
[x,n]=impseq(0,-20,290);
h=filter(b,a,x);
stem(n,h);
title('脉冲响应');
ylabel('h(n)');
xlabel('n')
便可得相应的单位脉冲响应曲线为:
.
.
35脉冲响应x 1014
12
10
8
h(n)6
4
2
0-20020406080100120n
图5
由图3、4和图5,我们可以发现,对于不稳定系统,它的脉冲响应和阶跃响应曲线一样,是趋于的。 ,
五、系统的频率特性分析
离散系统的频域分析就是研究离散信号通过离散系统以后在频谱结构上所产生的变化,因此我们通过绘制出系统的幅频特性和相频特性,实现对系统的频率特性分析。
5.1 稳定的数字信号系统的分析
5.1.1、幅频特性与相频特性
运用MATLAB分析系统时,我们首先将系统的差分方程进行Z变换,得到系统的传递函数为:
Yzz()21, (6) Hz,,,,1Xzz()21,
HjwHz()(),令可以得到: jw11ze,
jw21e,Hjw, (7) ,,1jw21e,
.
.
利用分式,结合以下程序: Hjw()1
clear all
clc
w=linspace(-5,5,200);%-5到5之间取200个点
a=[2,-1];
b=[2,1];%系统参数
H=freqz(b,a,w);%Z变换的频率特性
subplot(2,1,1);%分割窗口
plot(w,abs(H));
grid
title('幅频特性');
xlabel('\omega');
ylabel('|H(j\omega|');
subplot(2,1,2);
plot(w,angle(H));
grid
title('相频特性');
xlabel('\omega');
ylabel('\phi(\omega)'); ab为分式中分母多项式的系数,为分子多项式的系数其中。
得到相应的特性曲线如下图所示:
幅频特性3
2|,
|H(j1
0-5-4-3-2-1012345,相频特性1
0.5
),0(,
-0.5
-1-5-4-3-2-1012345, 图6
.
.
5.2、不稳定的数字信号系统频析
5.2.1、幅频特性与相频特性
对于不稳定系统的频率特性,我们运用MATLAB进行绘制时,只需将5.1.1程序中的ab与换成相应的不稳定系统的参数即可,如,则可得其相应曲线如图所示: ab,,,[1,0],[1,2]
幅频特性
3
2.5
|,2|H(j
1.5
1-5-4-3-2-1012345
,
相频特性
4
2
),0(,
-2
-4-5-4-3-2-1012345
,
图7
但对于稳定系统和不稳定系统的幅频和相频特性曲线,我们是无法直接从其曲线上判断出系统的稳定稳定性的。
.
.
六、心得体会
通过这次,数字信号课程设计,使我加深了对数字信号处理的理解,并且掌握了一些基本的MATLAB编程知识。在设计过程中,由于课上所学知识并不能满足这次设计的需要,我们必须查阅各种资料,学习各种指令的用法,因此加深了我对课本知识的理解。通过这一个星期的课题研究及设计,我对稳定系统与不稳定系统的时域分析和率特性分析原理有了进一步的理解,对MATLAB有了一定的掌握,亲身体会到了信号处理的一些应用。在历时一个星期的课设过程里,我体会了失败的苦涩,也尝到了成功的喜悦。通过这次课程设计,我懂得了理论必须与实践相结合,只有理论是远远不够的,只有把所学知识与实践相结合起来,从理论中得出结论并从实践中检验结论,才能真正提高自己的独立思考和实际动手的能力,往后,一定要注意理论与实践的相结合。总之,这次的实训给予了我不同的学习方法和体验,让我深切的认识到实践的重要性,使我受益匪浅。在以后的学习过程中,我会更加注重自己的操作能力和应变能力,提高自己对MATLAN的 运用能力,全面学习,全面发展~
.
.
七、参考文献
[1]、谢平,王娜.信号处理原理及应用.北京:机械工业出版社,2008,10
[2]、楼顺天.MATLAB 5.x程序设计语言.西安:西安电子科技大学出版社
[3]、薛年喜.MATLAB在数字信号处理中的应用(第二版).清华大学出版社
[4]、鄢景华.自动控制原理(修订版).哈尔滨:哈尔滨工业大学出版社,2006,3
[5]、牛力勇,高庆.数字信号处理实验指导书..
.
第七章附录
7.1材料清单
表3-5材料清单综合布线系统设备配置清单及报价
.
.
7.2平面布置图
(见图纸分册)
.