离散系统频域分析及matlab实现.docx
- 文档编号:30114880
- 上传时间:2023-08-05
- 格式:DOCX
- 页数:21
- 大小:120.51KB
离散系统频域分析及matlab实现.docx
《离散系统频域分析及matlab实现.docx》由会员分享,可在线阅读,更多相关《离散系统频域分析及matlab实现.docx(21页珍藏版)》请在冰豆网上搜索。
离散系统频域分析及matlab实现
《数字信号处理》
课程设计报告
离散系统的频域分析及matlab实现
专业:
通信工程
班级:
通信11级
组次:
姓名及学号:
姓名及学号:
组员承担任务
指
导
教
师
评
价
意
见
离散系统的频域分析及matlab实现
一、设计目的
1.熟悉并掌握matlab软件的使用;
2.掌握离散系统的频域特性;
3.学会分析离散系统的频域特性的方法;
二、设计任务
1.设计一个系统函数系统的频率响应进行分析;
2.分析系统的频域响应;
3.分析系统的因果稳定性;
4.分析系统的单位脉冲响应;
三、设计原理
1.系统函数
对于离散系统可以利用差分方程,单位脉冲响应,以及系统函数对系统进行描述。
在本文中利用系统函数H(z)进行描述。
若已知一个差分方程为
M
N
Z变换,最终可以得到系统函数的一
y(n)bix(ni)
aiy(n
i),则可以利用双边取
i0
i1
M
i
般式H(z),H(z)
Y(z)
biz
i
0
。
若已知系统的单位脉冲响应,则直接将其进行Z
X(z)
N
aiz
i
i
0
变换就可以得到系统函数H(z)。
系统函数表征系统的复频域特性。
2.系统的频率响应:
利用Z变化分析系统的频率响应:
设系统的初始状态为零,系统对输入为单位脉冲序列(n)的响应输出称为系统的单位脉冲响应h(n)。
对h(n)进行傅里叶变换,得到:
H(ej
n)h(n)ejn|H(ejwn)|ej(w)
-
其中|H(ejwn)|称为系统的幅频特性函数,
()称为系统的相位特性函数。
H(ejw)表示的
是系统对特征序列ejwn的响应特性。
对于一个系统输入信号为
x(n)
ejn,则系统的输出信
号为H(ejw)ejwn。
由上可以知道单频复指数信号
ejwn通过频率响应函数为
H(ejw)后,输
出仍为单频复指数信号,其幅度放大了|H(ejw)
|,相移为
()。
对于系统函数H(z)与H(w)之间,若系统函数
H(z)的收敛域包含单位圆
|z|=1,则有
H(ejw)H(z)|zejw,在MATLAB中可以利用freqz函数计算系统的频率响应。
(1)[h,w]=freqz(b,a,n)
可得到n点频率响应,这n个点均匀地分布在上半单位圆(即
),并将这n点频率记
录在w中,相应的频率响应记录在h中。
n最好能取2的幂次方,如果缺省,则n=512。
(2)[h,w]=freqz(b,a,n,'whole')
在之间均匀选取n个点计算频率响应。
(3)[h,w]=freqz(b,a,n,Fs)
Fs为采样频率(以Hz为单位),在0~Fs/2频率范围内选取n个频率点,计算相应的频率响应。
(4)[h,w]=freqz(b,a,n,'whole',Fs)
在0~Fs之间均匀选取n个点计算频率响应。
(5)freqz(b,a)
可以直接得到系统的幅频和相频特性曲线。
其中幅频特性以分贝的形式给出,频率特性曲线的横轴采用的是归一化频率,即Fs/2=1。
3.系统的因果性和稳定性
3.1因果性
因果系统其单位脉冲响应序列h(n)一定是一个因果序列,其z域的条件是其系统函数H(z)
的收敛域一定包含,即点不是极点,极点分布在某个圆内,收敛域在某个圆外。
3.2稳定性
系统稳定就要求
|h(n)|
,由序列的
H(ejw)
存在条件和
H(ejw)
H(z)|zejw可以知道
-
系统稳定的z域条件就是H(z)的收敛域包含单位圆,即极点全部分布在单位圆内部。
由上3.1和3.2可知,利用系统的零极点分布图可以判断系统的因果性和稳定性。
若在零极点分布图中,若系统的极点都分布在单位圆内,则此系统是因果系统,若有极点
分布在单位圆外,则此系统是非因果系统。
在MATLAB中可以利用zplane函数画出系统的零极点分布图。
系统函数的零极点图的绘制:
zplane(b,a)。
其中b为系统函数的分子,为系统函数的分母。
a
4.系统的单位脉冲响应
设系统的初始状态为零,系统对输入为单位脉冲序列(n)的响应输出称为系统的单位脉冲响应h(n)。
对于离散系统可以利用差分方程,单位脉冲响应,以及系统函数对系统进行描述。
单位脉冲响应是系统的一种描述方法,若已知了系统的系统函数,可以利用系统
得出系统的单位脉冲响应。
在MATLAB中利用impz由函数函数求出单位脉冲响应h(n)。
稳定系统的频率响应就是系统函数在单位圆上的取值,计算系统的频率响应,可将
系统函数中的Z变量用ej代入即可得到。
频率响应是一个复函数,其模叫幅度响应,其相
角叫相位响应,它反映了输入序列的频谱经系统后所发生的变化规律。
从幅频曲线上可直
观看到各频率分量的幅度变化情况,从相频曲线上可直观看到各频率分量的相移情况。
根
据频响曲线分析系统对信号频谱的影响,概念清楚、简单直观,对信号综合也意义重大,
但要将一个较复杂的频率响应复函数转化成幅度响应和相位响应并图示,计算量大且容易
出错,图示结果也不一定精确。
利用MATLAB函数这些问题都迎刃而解。
幅频特性,它指的是当
w从0到
变化时,
|H(jw)|的变化特性,记为
A(w),相频
特性,指的是当
w从
0到
变化时,
H(j
)的变化特性称为相频特性,记为
()。
离散系统的幅频特性曲线和相频特性曲线直观的反应了系统对不同频率的输入序列的处理情况。
四、设计过程
下面已知一个系统的系统函数为
H(z)
0.30.06z1
11.1z1
0.55z2
0.125z3
1.系统的频率响应:
MATLAB源程序为:
clc;closeall;clearall;
a=[1-0.3];
b=[1-1.60.9425];
[H,w]=freqz(a,b,’whole’);
subplot(2,1,1);
plot(w/pi,abs(H));
title(‘幅度谱’);
xlabel(‘\omega/\pi’);
ylabel(‘|H(e^j^\omega)|’);
gridon;
subplot(2,1,2);
plot(w/pi,angle(H));
title(‘相位谱’);
xlabel(‘\omega/\pi’);
ylabel(‘phi(\omega)’);
gridon;
程序运行结果如图1:
图1
由图1分析可以知道此系统为带通;
2.分析系统的稳定性:
MATLAB源程序为:
clc;closeall;clearall;
a=[1-0.3];
b=[1-1.60.9425];
zplane(a,b);
title('零极点分布图');
程序运行结果如图2所示:
图2
由图2分析可以知道系统的所有极点都位于单位圆内,所以此系统为稳定系统;
3.分析系统的因果性:
MATLAB源程序为:
clc;closeall;clearall;
a=[1-0.3];
b=[1-1.60.9425];
impz(a,b,30);
title('系统单位脉冲响应');
axis([-330-22])
程序运行结果如图3所示:
图3
由图3分析可以知道此系统的单位脉冲响应序列为一个因果序列,所以此系统为因果系统;
五、收获与体会
用MATLAB辅助分析离散系统,程序简单易懂,计算准确,绘图方便,只要掌握了系
统分析的概念原理和方法,繁杂的分析计算及图形显示均可用MATLAB实现。
利用MATLAB
的计算功能,不仅使大量的手工计算得以简化,也使得系统分析更为简便和高效;利用
MATLAB的绘图功能,有利于分析结果的直观理解,也有利于深入掌握所学的内容,更能培
养学习的兴趣。
六.参考文献:
(1)数字信号处理丁玉美西安电子科技大学出版社
(2)数字信号处理及MATLAB实现余成波清华大学出版社
(3)数字信号处理教程--matlab释义与实现陈怀琛电子工业出版社
(4)matlab7辅助信号处理技术与应用飞思科技产品研发中心电子工业出版社
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 离散系统 分析 matlab 实现
![提示](https://static.bdocx.com/images/bang_tan.gif)