基于MATLAB的连续时间系统的频域分析【文档知识】

《基于MATLAB的连续时间系统的频域分析【文档知识】》由会员分享，可在线阅读，更多相关《基于MATLAB的连续时间系统的频域分析【文档知识】（21页珍藏版）》请在装配图网上搜索。

1、课程设计任务书学生姓名： 专业班级：电子科学与技术0701班 指导教师： 刘金根 工作单位： 信息工程学院 题 目: 基于MATLAB的连续时间系统的频域分析 初始条件：MATLAB 6.5 微机要求完成的主要任务: 深入研究连续时间信号和系统时域分析的理论知识。利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能，实现连续时间系统频域分析。1.利用MATLAB分析系统的频率特性；2.用MATLAB实现连续时间信号的采样及重构；3. 撰写MATLAB应用实践课程设计说明书。时间安排：学习MATLAB语言的概况 第1天学习MATLAB语言的基本知识 第2、3天学习MATLAB语言的应

2、用环境，调试命令，绘图能力 第4、5天课程设计 第6-9天答辩 第10天指导教师签名： 年 月 日系主任（或责任教师）签名： 年 月 日青草绿2目 录摘 要1.绪 论12.对课题内容的分析22.1连续时间信号概述22.2采样定理 22.3 总体思路23MATLAB的仿真实现33.1利用MATLAB分析系统的频率特性33.1.1低通滤波器的频率特性33.1.2高通滤波器的频率特性43.1.3全通滤波器的频率特性63.1.4帯通滤波器的频率特性63.2用MATLAB实现连续时间信号的采样及重构63.2.1过采样73.2.2等采样73.2.3欠采样84.心得体会105.参考文献11附录12摘要本文介

3、绍了基于MATLAB的连续时间系统的频域分析。首先利用MATLAB分析了系统的频率特性，分别分析了基于连续时间系统的低通滤波器、高通滤波器、带通滤波器、全通滤波器的频域特性，并依次做出了它们的时域冲激响应波形、频域内幅频特性波形、相频特性波形。在编程过程中分别用到了y=abs( ) 、y=angle( ) 、h=freqs(b,a,w ) 等函数。然后用MATLAB实现了连续时间信号的采样及重构，并以f(t)=Sa(t)为例，分别以过采样、等采样、欠采样三种情况,绘出原信号、采样信号、重构信号的时域波形图。关键词：连续时间系统；频特性；采样；重构AbstractThis article int

4、roduced based on the MATLAB run-on time system frequency range analysis.First has analyzed the system frequency characteristic using MATLAB, analyzed separately based on the run-on time system low pass filter, has passed the filter, the bandpass filter high, all passes the filter the frequency range

5、 characteristic, and has in turn made in their time domain impulse response profile, the frequency range the amplitude-frequency characteristic profile, the frequency characteristic profile.Used y=abs separately in the programming process (), y=angle (), h=freqs (b, a, w) and so on the functions.The

6、n has realized the run-on time signal sampling and restructuring with MATLAB, and take f(t)=Sa(t) as the example, respectively by the sampling, and so on the sampling, has owed the sampling three kind of situations, draws the original signal, the sampling signal, the restructuring signal time domain

7、 oscillogram.Key word: Run-on time system; Frequency characteristic; Sampling; Restructuring青草绿21 绪论MATLAB是国际上公认的优秀科技应用软件，它的基本功能是数值计算、符号运算、图形控制，它的出现给“信号与系统”课程的计算机辅助教学带来了福音，使利用计算机辅助学生完成“信号与系统”课程的数值计算、信号与系统分析的可视化建模及仿真调试成为可能。该软件由公司于是1984 年推出，经过十几年的发展与完善，目前已成为科技界最流行的应用软件。它的主要特点是：（）高效的数值计算及符号计算功能，能使用户从繁杂

8、的数学运算分析中解脱出来。（）完备的图形处理功能，实现计算结果和编程的可视化。（）友好的用户界面及接近数学表达式的自然化语言，使学者易于学习和掌握。（）功能丰富的应用工具箱（如信号处理工具箱），为用户提供了大量方便实用的处理工具。运用MATLAB 对信号与线性系统进行分析与实现的具体方法和过程，其目的在于：（） 让学生在学习“信号与系统”课程的同时，掌握MATLAB 的应用，对MATLAB语言在低年级学生中的推广应用起到促进作用。（） 学会应用MATLAB 的数值计算功能，将学生从繁琐的数学运算中解脱出来，从而将便多的时间留于对信号与系统的基本分析方法和应用的理解与思考。（） 让学生将课程中的

9、重点、难点及部分课后练习用MATLAB 进行形象、直观的可视化计算机模拟与仿真实现，从而加深对信号与系统基本原理、方法及应用的理解，以培养学生主动获取知识和独立解决问题的能力，为学习后继专业课打下坚实的基础。本文将以MATLAB为工具，对信号与系统在联系时间系统的频域进行分析。2 对课题内容的分析2.1连续时间信号概述在某一时间区间内，除若干个不连续点外，如果任意识可都可给出确定的函数值，则称该信号为连续时间信号，简称为连续信号。从严格意义上讲，MATLAB数值计算的方法并不能处理连续时间信号，然而，可利用连续信号在等时间间隔的取样值来近似表示连续信号，即当取样时间间隔足够小时，这些离散样值能

10、够被MATLAB处理，并且能较好地近似表示连续信号。2.2采样定理对于一个有限频宽信号进行理想采样，当采样频率时，采样值唯一确定；当此采样信号通过截止频率的理想低通滤波器后，原始信号可以完全重建。通常把最低允许的采样频率2称为奈奎斯特频率，把最大允许的采样间隔称为奈奎斯特间隔。2.3总体思路利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能，实现连续时间系统频域分析。其中分析系统的频率特性时用到了y=abs( ) 、y=angle( ) 、h=freqs(b,a,w ) 等函数。实现连续时间信号的采样及重构时把采样分为了过采样、等采样、欠采样三种情况，分别做它们的原信号、采样信号

11、、重构信号的时域波形图，其中用到Sa(t)=sinc(t/pi)、fa=f*Ts*wc/pi*sinc(wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t)等函数。3 设计内容3.1利用MATLAB分析系统的频率特性3.1.1低通滤波器的频率特性低通滤波器模型： 取k=70,b=10,c=90可得图3.1.1图3.1.1低通滤波器时域冲激响应波形、频域内幅频特性、相频特性。对于低通滤波器，可以从图3.1.1看出，当wH(s)0.8,而其相频特性为随着w的增加而减小并逐渐趋近于-。3.1.2高通滤波器的频率特性高通滤波器模型：以下解法只供参考，方法

12、不唯一：若令1）当，可以得出：2）考虑一种情况:当，由以上两条件可得=-90-157j, =-90+157j,k=1,并由此可通过matlab得图3.1.2图3.1.2高通滤波器时域冲激响应波形、频域内幅频特性、相频特性。对于高通滤波器，可以从图3.1.2看出，当w50时 1.2H(s)0.8,而其相频特性为随着w的增加而减小并逐渐趋近于0。3.1.3全通通滤波器的频率特性全通滤波器模型：可以取=4，通过matlab可得图3.1.3图3.1.3全通通滤波器时域冲激响应波形、频域内幅频特性、相频特性。当w 从0增大时，H( jw)的幅频特性是一条数值为1的水平线，即对输入信号的各频率分量都进行等

13、值传输；而j(w)从0开始下降，最终趋于- 。这种网络称为全通网络，在传输系统中常用来进行相位校正，如作相位均衡器或移相器。3.1.4带通滤波器的频率特性带通滤波器模型：可取=25,b=40, =18, =306由matlab可得图3.1.4图3.1.4带通滤波器时域冲激响应波形、频域内幅频特性、相频特性。对于带通滤波器，可以从图3.1.4看出，当40w10时 1.2H(s)0.8,而其相频特性为随着w的增加先增加而后减小并逐渐趋近于-/2。3.2用MATLAB实现连续时间信号的采样及重构我们选取信号f (t ) = Sa(t )作为被采样的信号，是因为：第一， f (t)是一个带限信号，其=

14、 1 ；第二，它是一个典型的信号，是分析其他信号的基础，因此完全有必要对它信号特征详加了解。此外，应该指出的是，实际信号中，绝大多数都不是严格意义上的带限信号，这时根据实际精度要求来确定信号的带宽 。对于一个有限频宽信号进行理想采样，当采样频率时，采样值唯一确定；当此采样信号通过截止频率的理想低通滤波器后，原始信号可以完全重建。通常把最低允许的采样频率2称为奈奎斯特频率，把最大允许的采样间隔称为奈奎斯特间隔。3.2.1过采样过采样是使用远大于奈奎斯特采样频率的频率对输入信号进行采样。对于f(t)=Sa(t)，它的带宽=1，为了由f (t) 的采样信号不失真的重构f(t)，由时域采样定理知采样间

15、隔，取，即为过采样。利用MATLAB 中的抽样函数Sinc(t)=sin()/来表示Sa(t),有Sa(t)=Sinc(t/)。图3.2.1过采样原信号、采样信号、重构信号的时域波形图3.2.2等采样 由采样定理知f(t)=Sa(t)的采样间隔，取=，即为等采样图3.2.2等采样原信号、采样信号、重构信号的时域波形图3.2.3欠采样欠采样的信号频率(大于fs/2),由采样定理知f(t)=Sa(t)的采样间隔，取=1.5，即为欠采样。图3.2.3欠采样原信号、采样信号、重构信号的时域波形图4.心得体会MATLAB对自己以后的工作和学习很有帮助。在刚开始设计时，有点茫然，连设计的题意都不太明白，在

16、逐步的深入学习、了解之后慢慢有所了解。开始写程序时也是无从下手，是通过查阅一些图书资料和网上资料后开始进入正题。在这次课程设计中，我通过多方面地搜集资料，成功地用MATLAB编写出低通、高通、全通、帯通滤波器时域冲激响应波形、频域内幅频特性、相频特性的仿真实现，及过采样、等采样、欠采样的原信号、采样信号、重构信号的时域波形图。通过这次的实践，我明白了要将理论与实际相结合的道理，尽管这个过程会有一些辛苦，但通过努力实现后，就能大大深化我对知识的理解程度，增长实践经验。这表现在我对连续时间系统的频域的理解的加深、MATLAB的功能特性都有了进一步的认识。总之，我在本次课程设计中学到了很多关于MAT

17、LAB的知识，获益良多。5.参考文献1 孙祥，徐流美，吴清.MATLAB 7.0基础教程.北京：清华大学出版社，20062 刘泉，姜雪梅.信号与系统.北京：高等教育出版社，20063 唐向宏，岳恒立，邓雪峰.MATLAB及在电子信息类课程中的应用.北京：电子工业出版社，20064 赵静，张瑾，高新科.基于MATLAB的通信系统仿真.北京：北京航空航天大学出版社.20075 梁虹.信号与线性系统分析 基于MATLAB的方法与实现.北京：电子工业出版社，2006.5附录低通b=90;a=1 10 90; sys=tf(b,a);t=0:0.1:5;h=impulse(sys,t);subplot(

18、221);plot(h);gridxlabel(t);ylabel(h(t);title(h(t);h,w=freqs(b,a,100);h1=abs(h);h2=angle(h);subplot(222);plot(w,h1);gridxlabel(角频率(w);ylabel(幅度);title(H(jw)的幅频特性);subplot(223);plot(w,h2*180/pi);gridxlabel(角频率(w);ylabel(相位（度）);title(H(jw)的相频特性);高通k=1;b=1 0 0*k;a=conv(1 90+157j,1 90-157j);sys=tf(b,a);t

19、=0:0.1:10;h=impulse(sys,t);subplot(221);plot(h);gridxlabel(t);ylabel(h(t);title(h(t);h,w=freqs(b,a,100);h1=abs(h);h2=angle(h);subplot(222);plot(w,h1);gridxlabel(角频率(w);ylabel(幅度);title(H(jw)的幅频特性);subplot(223);plot(w,h2*180/pi);gridxlabel(角频率(w);ylabel(相位（度）);title(H(jw)的相频特性);全通b=-1 4;a=1 4;sys=tf(

20、b,a);t=0:0.1:10;h=impulse(sys,t);subplot(221);plot(h);gridxlabel(t);ylabel(h(t);title(h(t);h,w=freqs(b,a,100);h1=abs(h);h2=angle(h);subplot(222);plot(w,h1);gridaxis(0,100,0,1.5)xlabel(角频率(w);ylabel(幅度);title(H(jw)的幅频特性);subplot(223);plot(w,h2*180/pi);gridxlabel(角频率(w);ylabel(相位（度）);title(H(jw)的相频特性)

21、;帯通k=1;b=0 25 40*k;a=conv(1 9+15j,1 9-15j);sys=tf(b,a);t=0:0.1:10;h=impulse(sys,t);subplot(221);plot(h);gridxlabel(t);ylabel(h(t);title(h(t);h,w=freqs(b,a,100);h1=abs(h);h2=angle(h);subplot(222);plot(w,h1);gridxlabel(角频率(w);ylabel(幅度);title(H(jw)的幅频特性);subplot(223);plot(w,h2*180/pi);gridxlabel(角频率(w

22、);ylabel(相位（度）);title(H(jw)的相频特性);过采样t=-15:0.01:15;f=sin(t)./t;subplot(221);plot(t,f);xlabel(t);ylabel(f(t);title(f(t)=sin(t)/t的过采样原信号);gridwm=1;wc=wm;Ts=0.02*pi./wm; ws=2*pi./Ts; n=-100:100; nTs=n*Ts f=sinc(nTs/pi);Dt=0.005;t=-15:Dt:15;fa=f*Ts*wc/pi*sinc(wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,len

23、gth(t);t1=-15:0.5:15;f1=sinc(t1/pi);subplot(222);stem(t1,f1);xlabel(kTs);ylabel(f(kTs);title(sa(t)=sinc(t/pi)的过采样采样信号);subplot(223);plot(t,fa)xlabel(t);ylabel(fa(t);title(由sa(t)=sinc(t/pi)的过采样信号重构sa(t);grid;等采样t=-15:0.01:15;f=sin(t)./t;subplot(221);plot(t,f);xlabel(t);ylabel(f(t);title(f(t)=sin(t)/t

24、的等采样信号);gridwm=1;wc=wm;Ts=0.2*pi./wm; ws=2*pi./Ts; n=-100:100; nTs=n*Ts f=sinc(nTs/pi);Dt=0.005;t=-15:Dt:15;fa=f*Ts*wc/pi*sinc(wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t);t1=-15:0.5:15;f1=sinc(t1/pi);subplot(222);stem(t1,f1);xlabel(kTs);ylabel(f(kTs);title(sa(t)=sinc(t/pi)的等采样信号);subplot(223)

25、;plot(t,fa)xlabel(t);ylabel(fa(t);title(由sa(t)=sinc(t/pi)的等采样信号重构sa(t);grid;欠采样t=-15:0.01:15;f=sin(t)./t;subplot(221);plot(t,f);xlabel(t);ylabel(f(t);title(f(t)=sin(t)/t的欠采样原信号);gridwm=1;wc=wm;Ts=1.5*pi./wm; ws=2*pi./Ts; n=-100:100; nTs=n*Ts f=sinc(nTs/pi);Dt=0.005;t=-15:Dt:15;fa=f*Ts*wc/pi*sinc(wc/

26、pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t);t1=-15:0.5:15;f1=sinc(t1/pi);subplot(222);stem(t1,f1);xlabel(kTs);ylabel(f(kTs);title(sa(t)=sinc(t/pi)的欠采样采样信号);subplot(223);plot(t,fa)xlabel(t);ylabel(fa(t);title(由sa(t)=sinc(t/pi)的欠采样信号重构sa(t);grid;本科生课程设计成绩评定表姓 名陈讲重性 别 男专业、班级电子科学与技术0701课程设计题目： 基于MATLAB的连续时间系统的频域分析课程设计答辩或质疑记录：成绩评定依据：最终评定成绩（以优、良、中、及格、不及格评定）指导教师签字： 年 月 日