武汉工程大学matlab实验四lst离散系统的频域分析文档格式.docx
- 文档编号:6283130
- 上传时间:2023-05-06
- 格式:DOCX
- 页数:13
- 大小:226.27KB
武汉工程大学matlab实验四lst离散系统的频域分析文档格式.docx
《武汉工程大学matlab实验四lst离散系统的频域分析文档格式.docx》由会员分享,可在线阅读,更多相关《武汉工程大学matlab实验四lst离散系统的频域分析文档格式.docx(13页珍藏版)》请在冰点文库上搜索。
1、实验目的
1.通过在频域中仿真LSI离散时间系统,理解离散时间系统对输入信号或延迟信号进行频域处理的特性。
2.理解LSI离散时间系统的传输函数和频率响应的概念。
3.理解LSI离散时间系统的滤波特性及滤波器的相关特性。
4.理解并仿真LSI离散时间系统的零、极点分布表征及特性关系。
二、实验设备
计算机,MATLAB语言环境
三、实验基础理论
LSI离散时间系统可用差分方程描述如下:
对应的传输函数和频率响应分别为:
分别有零点和极点。
4、实验内容与步骤
1.传输函数和频率响应分析
按以下的传输函数分别编程计算
和
计算当
时因果LSI离散时间系统的频率响应,并求出它们的群时延及冲激响应的开始部分(前100个值)。
b=[0.15,0,-0.15];
a=[1,-0.5,0.7];
b1=[0.15,0,-0.15];
a1=[0.7,-0.5,1];
figure
(1);
subplot(1,2,1);
plot(grpdelay(b,a));
grid;
title('
系统1的群延时'
)
subplot(1,2,2);
plot(grpdelay(b1,a1));
grid,title('
系统2的群延时'
);
w=0:
0.1:
pi;
z=exp(j*w);
[H1,w]=freqz(b,a);
[H2,w]=freqz(b1,a1);
figure
(2);
subplot(2,2,1);
plot(w,abs(H1)),title('
系统1的频率幅度响应'
subplot(2,2,2);
plot(w,abs(H2));
系统2的频率幅度响应'
subplot(2,2,3);
plot(w,angle(H1));
系统1的频率相位响应'
subplot(2,2,4);
plot(w,angle(H2));
系统2的频率相位响应'
xn=[1,zeros(1,99)];
hn1=filter(b,a,xn);
hn2=filter(b1,a1,xn);
n1=0:
length(hn1)-1;
n2=0:
length(hn2)-1;
figure(3);
stem(n1,hn1),xlabel('
n'
),ylabel('
h1(n)'
),title('
冲激响应hn1'
);
figure(4);
stem(n2,hn2),xlabel('
h2(n)'
冲激响应hn2'
2.画出上面两个LSI离散时间系统对应的零、极点图。
b=[0.15,0,-0.15];
subplot(2,1,1);
zplane(b,a);
h1的零极点分布'
subplot(2,1,2);
zplane(b1,a1);
h2的零极点分布'
;
3.滤波器仿真和特性实验
设计实现一个在0.45
处具有3dB截止角频率
的一阶无限冲激响应低通滤波器和一个无限冲激响应高通滤波器,计算并画出他们各自的增益响应,并证明它们是全通互补和功率互补的。
[b,a]=butter(1,0.45);
[h,w]=freqz(b,a,512);
plot(w/pi,20*log10(abs(h)));
xlabel('
频率/Hz'
ylabel('
幅值/dB'
IIR低通滤波器'
b1=[0.45];
a1=[11];
[BA]=lp2hp(b1,a1,0.45);
[h1,w]=freqz(B,A);
plot(w/pi,abs(h1));
IIR高通滤波器'
4.作图分析
分别生成如下两个因果系统传输函数的零、极点图,并对该图进行研究,推断它们的稳定性。
b=[100];
a=[1-1.8480.85];
b1=[100];
a1=[1-1.8510.85];
figure
(1);
zplane(b,a);
系统1零极点分布图'
plot(w/pi,abs(H1));
axis([01-100600]);
频率响应'
w'
subplot(2,1,1);
zplane(b1,a1);
系统2零极点分布图'
subplot(2,1,2);
plot(w/pi,abs(H2));
频率'
若一个线性时不变系统的冲激响应是绝对可和,则此系统就是BIBO的稳定系统。
由此,无限冲激响应线性时不变系统稳定的必要条件是,随着输入序列点的增加,冲激响应衰减到零。
由上面两个系统频率响应的仿真图可知这两个系统稳定。
五、实验扩展与思考
1.在传输函数和频率响应分析实验中,选择了哪些类型的滤波器?
这些滤波器有何区别?
你会选择哪一种滤波器来滤波,为什么?
答:
在传输函数和频率响应分析实验中,选择了IIR滤波器和FIR滤波器。
两者的区别
阶数
时延
稳定性
运算量
IIR数字滤波器
低
小
存在稳定问题(需先判断系统是否稳定)
较小
FIR数字滤波器
高
大
不存在稳定问题
较大
2.常见的四种理想零相位数字滤波器有哪些?
分别画出它们的频率响应。
四种理想零相位数字滤波器为数字低通滤波器,数字高通滤波器,数字带通滤波器,数字带阻滤波器。
以实现一个在0.45
的滤波器为例,
[h1,omega]=freqz(b,a,512);
plot(omega/pi,20*log10(abs(h1)));
数字低通滤波器'
[h2,w]=freqz(B,A);
plot(w/pi,abs(h2));
数字高通滤波器'
[c,d]=bilinear(b,a,512);
[h3,omega]=freqz(b,a,512);
plot(omega/pi*256,20*log10(abs(h3)));
数字带通滤波器'
[C,D]=bilinear(b1,a1,0.45);
[h4,omega]=freqz(b1,a1,512);
subplot(2,2,4)
plot(omega/pi*256,20*log10(abs(h4)));
数字带阻滤波器'
3.从以上实验中可见,因果无限冲激响应滤波器稳定的条件是什么?
无限冲激响应系统滤波器稳定的条件是,随着输入序列点的增加,冲激响应衰减到零,频率响应衰减到零。
6、总结实验心得体会
通过对lst离散系统的频域分析的实验使我对传输函数和频率函数有了更深刻的记忆。
收获最大的是对两种数字滤波器,IIR数字滤波器和FIR数滤波器的设计理解。
明白了两者的区别,设计要求和设计方法。
了解了数字低通滤波器与数字高通滤波器,数字带通滤波器与数字带阻滤波器两两之间的关系。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 武汉 工程 大学 matlab 实验 lst 离散系统 分析