MATLAB数字信号处理课程设计.docx
《MATLAB数字信号处理课程设计.docx》由会员分享,可在线阅读,更多相关《MATLAB数字信号处理课程设计.docx(24页珍藏版)》请在沃文网上搜索。
1、目录目录2摘要3一设计目的和要求4二设计原理及方法42.1设计原理42.2设计方法7三实验内容73.1信号的产生73.2滤波器的设计103.3信号滤波分析13参考文献19学习心得20摘要几乎所有的工程技术领域都要涉及到信号处理问题,信号处理一般是包括数据采集以及对信号进行分析、变换、综合、估计与识别。对于数字信号来说,数字信号的幅度和时间都是离散值,数字信号处理是采用数值计算的方法完成对信号的处理。而待处理的的信号往往夹带着噪声。这就需要数字滤波器对信号进行滤波处理,滤除其中的噪声,得到想要的信号。所谓数字滤波器,就是输入、输出都是数字信号的,通过数值计算处理改变输入信号所含频率成分的相对比例
2、,或者滤除某些频率成分的数字器件或程序。常用的经典滤波器有低通、高通、带通、带阻。关键字:信号产生 数字滤波器噪声 频谱分析一设计目的和要求1、产生一个连续信号,包含低频,中频,高频分量。2、对产生的信号进行采样,进行频谱分析。3、分别设计高通、低通、带通滤波器对信号进行滤波处理,观察滤波后信号的频谱。二设计原理及方法2.1设计原理理论上信号的采样要符合奈奎斯特采样定律,就是采样频率要高一点,一般为被采信号最高频率的2倍,只有这样,才能保证频域不混叠,也就是采样出来数字信号中包含了被采信号的所有信息,而且没有引入干扰。这就是信号的时域采样。频谱分析是指对信号进行频域谱的分析,观察其频域的各个分
3、量的功率大小,其理论基础是傅立叶变换,现在一般采用数字的方法,也就是将时域信号数字化后做FFT,可以得到频域的波形。数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。可以设计系统的频率响应,让它满足一定的要求,从而对通过该系统的信号的某些特定的频率成分进行过滤,这就是滤波器的基本原理。如果系统是一个连续系统,则滤波器称为模拟滤波器。如果系统是一个离散系统,则滤波器称为数字滤波器。信号通过线性系统后,其输出就是输入信号和系统冲激响应的卷积。除了外,的波形将不同于输入波形。从频域分析来看,信号通过线性系统后,输出信号的频谱将是输入信号的频谱与系统传递函
4、数的乘积。除非为常数,否则输出信号的频谱将不同于输入信号的频谱,某些频率成分较大的模,因此,中这些频率成分将得到加强,而另外一些频率成分的模很小甚至为零,中这部分频率分量将被削弱或消失。因此,系统的作用相当于对输入信号的频谱进行加权。IIR滤波器的设计原理:IIR数字滤波器的设计一般是利用目前已经很成熟的模拟滤波器的设计方法来进行设计,通常采用模拟滤波器原型有butterworth函数、chebyshev函数、bessel函数、椭圆滤波器函数等。IIR数字滤波器的设计步骤:(1)按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标;(2)根据模拟滤波器技术指标设计为响应的模拟低通滤
5、波器;(3)很据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器;(4)如果要设计的滤波器是高通、带通或带阻滤波器,则首先把它们的技术指标转化为模拟低通滤波器的技术指标,设计为数字低通滤波器,最后通过频率转换的方法来得到所要的滤波器。采样是连续信号数字处理的第一个关键环节。对一个连续信号进行理想采样的过程可用(1.1)式表示。 其中为的理想采样,为周期冲激脉冲,即的傅里叶变换为将(1.2)式代入(1.1)式并进行傅里叶变换,式中的就是采样后得到的序列, 即的傅里叶变换为 比较(1.5)和(1.4)可知为了在数字计算机上观察分析各种序列的频域特性,通常对在上进行M点采样来观察分析。对长度
6、为N的有限长序列,有其中 滤波器的设计,以巴特沃斯滤波器为例以s替换,将幅度平方函数写成s函数 复变量:s=,上式表示有2N个极点,极点用下面公式表示:=K=0,1,2,32N-1。2.2设计方法首先利用MATLAB分别产生低频中频高频信号,然后进行叠加得到连续时间信号;对所产生的连续时间信号进行采样,得到数字信号;对信号进行FFT频谱分析,绘制其频谱图;根据信号频谱分析的结果,分别设计高通,低通,带通滤波器,得到滤波器的幅频及相频特性。三实验内容3.1信号的产生在MATLAB里面利用正玄函数产生中、低高频的合成信号y=sin(2*pi*50*x)+sin(2*pi*200*x)+sin(2*
7、pi*600*x);产生的信号由50HZ、200HZ、600HZ合成。图1信号的产生对合成信号进行频谱分析利用MATLAB内部函数:z=fft(y);,对信号进行频谱分析图2产生信号的频谱通过对信号的频谱分析我们可以清楚的看到产生的信号的频谱,分为高、低、中三种频率。实验的源程序为:Fs=5000; x=0:1/Fs:20*pi;y=sin(2*pi*50*x)+sin(2*pi*200*x)+sin(2*pi*600*x);%产生低频、中频、高频信号subplot(2,1,1);plot(x(1:50),y(1:50);title(产生的信号);ylabel(幅度);xlabel(时间)z=
8、fft(y,2048);%傅里叶变换频谱分析subplot(2,1,2);plot(0:1023)*2500/1024,abs(z(1:1024);title(信号频谱分析);ylabel(幅度);xlabel(频率单位(HZ))axis(0,800,0,1500)3.2滤波器的设计低通滤波器的设计wp=2*pi*100;ws=2*pi*150;Rp=0.1;As=50;N,wc=buttord(wp,ws,Rp,As,s);B,A=butter(N,wc,s);k=0:511;fk=0:150/512:150;wk=2*pi*fk;HK=freqs(B,A,wk);plot(fk,20*lo
9、g(abs(HK);grid ontitle(低通滤波器);ylabel(幅度/db);xlabel(频率单位(HZ));axis(0,200,-1,1);通过上述程序的运行得到下述波形:图3低通滤波器结果分析:图3为设计的低通滤波器,通带截止频率为100HZ,阻带截止频率为150HZ。图示得到的低通滤波器满足实验要求。带通滤波器:图4带通滤波器带通滤波器代码:fpl=150;fpu=250;fsl=100;fsu=500;fs=5000;wp=2*fpl/fs,2*fpu/fs;ws=2*fsl/fs,2*fsu/fs;rp=1;rs=40;N,wpo=ellipord(wp,ws,rp,r
10、s);B,A=ellip(N,rp,rs,wpo);HK,w=freqz(B,A,512);plot(w/(2*pi)*fs,HK);grid ontitle(带通滤波器);ylabel(幅度/db);xlabel(频率单位(HZ));axis(0,500,-1,1);实验结果:通过对带通滤波器的设计得到图4所示的带通波形,通带截止频率为100HZ、300HZ,阻带截止频率为150HZ、250HZ。上述参数经过试验满足实验要求。高通滤波器的设计高通滤波器代码:wp=2*pi*500;ws=2*pi*400;Rp=0.01;As=50;N,wc=buttord(wp,ws,Rp,As,s);BH
11、,AH=butter(N,wc,high,s);k=0:511;fk=100:900/4096:900;wk=2*pi*fk;HK=freqs(BH,AH,wk);plot(fk,20*log(abs(HK);grid ontitle(高通滤波器);ylabel(幅度/db);xlabel(频率单位(HZ));图5高通滤波器实验结果分析:通过对高通滤波器的设计得到图5所示的高通滤波器,通带截止频率为400HZ阻带截止频率为500HZ。经试验验证满足设计指标。5.信号处理把产生的信号分别通过设计的高、低、带通滤波器观察实验波形。3.3信号滤波分析低通滤波后的信号分析:图6低通滤波后的信号低通滤波
12、器实验代码:%信号产生Fs=5000; x=0:1/Fs:2*pi;y=sin(2*pi*50*x)+sin(2*pi*200*x)+sin(2*pi*600*x);%产生低频、中频、高频信号subplot(5,1,1);plot(x(1:50),y(1:50);grid ontitle(产生的信号);ylabel(幅度);xlabel(时间(t))%傅里叶变换频谱分析z=fft(y,2048);subplot(5,1,2);plot(0:1023)*2500/1024,abs(z(1:1024);grid ontitle(产生信号的频谱分析);ylabel(幅度);xlabel(频率单位(H
- 1.请仔细阅读文档,确保文档完整性,对于不预览、不比对内容而直接下载带来的问题本站不予受理。
- 2.下载的文档,不会出现我们的网址水印。
- 3、该文档所得收入(下载+内容+预览)归上传者、原创作者;如果您是本文档原作者,请点此认领!既往收益都归您。
下载文档到电脑,查找使用更方便
20 积分
下载 | 加入VIP,下载更划算! |
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 数字信号 处理 课程设计