数字信号处理实验五

合集下载

数字信号处理实验实验五汇总

数字信号处理实验实验五汇总

数字信号处理实验报告实验名称:应用FFT实现信号频谱分析学生姓名:z学生学号:学生班级:上课时间:周二上午指导老师:一、 实验目的(1) 能够熟练掌握快速离散傅里叶变换的原理及应用FFT 进行频谱分析的基本方法。

(2) 对离散傅里叶变换的主要性质及FFT 在数字信号处理中的重要作用有进一步的了解。

二、 实验原理1、离散傅里叶变换(DFT )及其主要性质DFT 表示离散信号的离散频谱,DFT 的主要性质中有奇偶对称性、虚实特性等。

通过实验可以加深理解。

实序列的DFT 具有偶对称的实部和奇对称的虚部,这可以证明如下: 由定义,可得X(k)=∑-=1)(N n kn N W n x=)2sin()()2cos()(110kn N n x j kn N n x N n N n ∑∑-=-=-ππX(N-k)=∑-=-10)()(N n n k N NWn x =∑-=-1)(N n kn nNnWWn x =∑-=-1)(N n kn N W n x=)2sin()()2cos()(110kn N n x j kn N n x N n N n ∑∑-=-=+ππ所以,X(k)=X *(N-k)实序列DFT 的这个特性,在本实验中可以通过实指数序列及三角序列看出来。

对于单一频率的三角序列来说,它的DFT 谱线也是单一的,这个物理意义可以从实验中得到验证,在理论上可以推导如下: 设)()2s i n ()(n R n Nn x N π= 其DFT 为 X(k)=∑-=-102)(N n kn Nen x π=kn Nj N n e n N ππ21)2sin(--=∑=kn N j n N j N n N j e e e j πππ22102)(21---=-∑=)(21)1(210)1(2+--=---∑k n N j N n k n N j e e j ππ从而,X(0)=0)(212102=---=∑n N j N n Nj e e j ππX(1)=22)1(21104Nj j N e j N n n N j -==-∑-=-πX(2)=0 …… X(N-2)=0X(N-1)=22)(21210)2(2Nj j N e e j n j N n n N N j =-=-∑-=--ππ以上这串式中X(0)反映了x(n)的直流分量,X(1)是x(n)的一次谐波,又根据虚实特性X *(N-1)=X(1),而其他分量均为零。

数字信号处理实验报告 (5)

数字信号处理实验报告 (5)

实验一信号、系统及系统响应一、实验目的1、熟悉连续信号经理想采样前后的频谱变化关系,加深对时域采样定的理解。

2、熟悉时域离散系统的时域特性。

3、利用卷积方法观察分析系统的时域特性。

4、掌握序列傅里叶变换的计算机实现方法,利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。

二、实验原理采样的的过程既是连续信号离散化的过程。

采用单位冲击串进行采样,为使采样信号能不失真的还原为采样前的信号,根据奈奎斯特采样率,采样频率应该大于信号最高频率的2倍。

因为时域的采样既是对时域的离散化处理,时域离散频域会进行周期延拓,为了防止频域频谱混叠,必须满足奈奎斯特采样定律。

线性卷积的过程为:反褶,移位,相乘,相加。

设一个N1点的序列与一个N2的序列进行卷积则得到N1+N2-1点的序列。

时域卷积,对应频域的相乘。

序列的傅里叶变换即DTFT 。

具有的性质有: 线性,移位性,对偶性,等等。

三、实验内容及步骤1)分析采样序列的特性。

产生采样序列()a x n ,A 444.128=,a =,0Ω=。

a 、 取采样频率s f 1kHz =,即T 1ms =。

观察所采样()a x n 的幅频特性()j X e ω和)(t x a 的幅频特性()X j Ω在折叠频率处有无明显差别。

应当注意,实验中所得频谱是用序列的傅立叶变换公式求得的,所以在频率量度上存在关系:T ω=Ω。

b 、改变采样频率,s f 300Hz =,观察()j X eω的变化并做记录。

c 、 进一步降低采样频率,s f 200Hz =,观察频谱混叠是否明显存在,说明原因,并记录()j X e ω的幅频曲线。

上图是采用不同采样频率时所得到的序列及其对应的傅里叶变换,从图中可以看到,当采样频率比较低时,频谱会发生混叠,且频率越低,混叠现象越明显。

增大采样频率可以有效地防止混叠。

2) 离散信号、系统和系统响应分析。

a 、观察信号()b x n 和系统h ()b n 的时域和频域持性;利用线形卷积求信号()b x n 通过系统h ()b n 的响应y(n),比较所求响应y(n)和h ()b n 的时域及频域特性,注意它们之间有无差异,绘图说明,并用所学结论解释所得结果。

数字信号处理实验五.FIR数字滤波器设计与软件实现

数字信号处理实验五.FIR数字滤波器设计与软件实现

实验五:FIR数字滤波器设计与软件实现一、实验指导1.实验目的(1)掌握用窗函数法设计FIR数字滤波器的原理和方法。

(2)掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。

(3)掌握FIR滤波器的快速卷积实现原理。

(4)学会调用MATLAB函数设计与实现FIR滤波器。

2.实验内容及步骤(1)认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理;(2)调用信号产生函数xtg产生具有加性噪声的信号xt,并自动显示xt及其频谱,如图1所示;图1 具有加性噪声的信号x(t)及其频谱如图(3)请设计低通滤波器,从高频噪声中提取xt中的单频调幅信号,要求信号幅频失真小于0.1dB,将噪声频谱衰减60dB。

先观察xt的频谱,确定滤波器指标参数。

(4)根据滤波器指标选择合适的窗函数,计算窗函数的长度N,调用MATLAB函数fir1设计一个FIR低通滤波器。

并编写程序,调用MATLAB快速卷积函数fftfilt实现对xt的滤波。

绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。

(4)重复(3),滤波器指标不变,但改用等波纹最佳逼近法,调用MATLAB函数remezord和remez设计FIR数字滤波器。

并比较两种设计方法设计的滤波器阶数。

提示:○1MATLAB函数fir1的功能及其调用格式请查阅教材;○2采样频率Fs=1000Hz,采样周期T=1/Fs;○3根据图1(b)和实验要求,可选择滤波器指标参数:通带截止频率fp=120Hz,阻带截至频率fs=150Hz,换算成数字频率,通带截止频率p 20.24pfωπ=T=π,通带最大衰为0.1dB,阻带截至频率s 20.3sfωπ=T=π,阻带最小衰为60dB。

○4实验程序框图如图2所示,供读者参考。

图2 实验程序框图4.思考题(1)如果给定通带截止频率和阻带截止频率以及阻带最小衰减,如何用窗函数法设计线性相位低通滤波器?请写出设计步骤.(2)如果要求用窗函数法设计带通滤波器,且给定通带上、下截止频率为pl ω和pu ω,阻带上、下截止频率为sl ω和su ω,试求理想带通滤波器的截止频率cl cu ωω和。

数字信号处理实验五

数字信号处理实验五

实验五:抽样定理一、实验目的1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。

2、进一步加深对时域、频域抽样定理的基本原理的理解。

3、观察信号抽样与恢复的图形,掌握采样频率的确定方法和内插公式的编程方法。

二、实验内容及步骤1、阅读并输入实验原理中介绍的例题程序,观察输出的数据和图形,结合基本原理理解每一条语句的含义。

2、已知一个连续时间信号f(t)=sinc(t),取最高有限带宽频率f m =1Hz 。

(1)分别显示原连续信号波形和F s =f m 、F s =2f m 、F s =3f m 三种情况下抽样信号的波形;dt=0.1;f0=1;T0=1/f0; fm=1;Tm=1/fm; t=-2:dt:2; f=sinc(t);subplot(4,1,1);plot(t,f);axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3;fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sinc(n);subplot(4,1,i+1);stem(n,f,'filled');axis([min(n),max(n),1.1*min(f),1.1*max(f)]);课程名称 数字信号处理 实验成绩 指导教师实 验 报 告院系 信息工程学院 班级 学号 姓名 日期end-2-1.5-1-0.50.511.5200.51原连续信号和抽样信号(2)求解原连续信号和抽样信号的幅度谱; dt=0.1;f0=1;T0=1/f0; fm=1;Tm=1/fm; t=-2:dt:2; N=length(t); f=sinc(t); wm=2*pi*fm; k=0:N-1; w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1));axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3;if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2; N=length(n); f=sinc(n); wm=2*pi*fs; k=0:N-1; w=k*wm/N; F=f*exp(-j*n'*w)*Ts;subplot(4,1,i+1);plot(w/(2*pi),abs(F));axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end00.511.522.533.540.20.40.60.811.200.511.522.533.54012(3)用时域卷积的方法(内插公式)重建信号。

数字信号处理实验(1-7)原始实验内容文档(含代码)

数字信号处理实验(1-7)原始实验内容文档(含代码)

实验要求1.每个实验进行之前须充分预习准备,实验完成后一周内提交实验报告;2.填写实验报告时,分为实验题目、实验目的、实验内容、实验结果、实验小结五项;3.实验报告要求:实验题目、实验目的、实验内容、实验结果四项都可打印;但每次实验的实验内容中的重要代码(或关键函数)后面要用手工解释其作用。

实验小结必须手写!(针对以前同学书写实验报告时候抄写代码太费时间的现象,本期实验报告进行以上改革)。

实验一信号、系统及系统响应实验目的:1. 掌握使用MATLAB进行函数、子程序、文件编辑等基本操作;2. 编写一些数字信号处理中常用序列的3. 掌握函数调用的方法。

实验内容:1.在数字信号处理的基本理论和MATLAB信号处理工具箱函数的基础上,可以自己编写一些子程序以便调用。

(1)单位抽样序列δ(n-n0)的生成函数impseq.m(2)单位阶跃序列u(n-n0)的生成函数stepseq.m(3)两个信号相加的生成函数sigadd.m(4)两个信号相乘的生成函数sigmult.m(5)序列移位y(n)=x(n-n0)的生成函数sigshift.m(6)序列翻褶y(n)=x(-n)生成函数sigfold.m(7)奇偶综合函数evenodd.m(8)求卷积和2.产生系列序列,并绘出离散图。

(1) x1(n)=3δ(n-2)-δ(n+4) -5≤n≤5(2) x3(n)=cos(0.04πn)+0.2w(n) 0≤n≤50其中:w(n)是均值为0,方差为1 的白噪声序列。

3.设线性移不变系统的抽样响应h(n)=(0.9)^n*u(n),输入序列x(n)=u(n)-u(n-10),求系统的输出y(n).实验二 系统响应及系统稳定性1.实验目的(1)掌握 求系统响应的方法。

(2)掌握时域离散系统的时域特性。

(3)分析、观察及检验系统的稳定性。

2.实验原理与方法在时域中,描写系统特性的方法是差分方程和单位脉冲响应,在频域可以用系统函数描述系统特性。

数字信号管理方案计划实验报告实验五

数字信号管理方案计划实验报告实验五

物理与电子信息工程学院实验报告实验课程名称:数字信号处理实验名称:FIR数字滤波器设计与软件实现班级:1012341姓名:严娅学号:101234153成绩:_______实验时间:2012年12月20 日一、实验目的(1)掌握用窗函数法设计FIR 数字滤波器的原理和方法。

(2)掌握用等波纹最佳逼近法设计FIR 数字滤波器的原理和方法。

(3)掌握FIR 滤波器的快速卷积实现原理。

(4)学会调用MATLAB 函数设计与实现FIR 滤波器。

二、实验原理1、用窗函数法设计FIR 数字滤波器的原理和方法。

如果所希望的滤波器的理想频率响应函数为 )(ωj d e H ,则其对应的单位脉冲响应为)(n h d =π21ωωωππd e e H j j d )(⎰- (2-1)窗函数设计法的基本原理是用有限长单位脉冲响应序列)(n h 逼近)(n h d 。

由于)(n h d 往往是无限长序列,且是非因果的,所以用窗函数)(n ω将)(n h d 截断,并进行加权处理,得到:)(n h =)(n h d )(n ω (2-2))(n h 就作为实际设计的FIR 数字滤波器的单位脉冲响应序列,其频率响应函数)(ωj d e H 为:)(ωj d e H =∑-=-1)(N n j e n h ω (2-3) 式中,N 为所选窗函数)(n ω的长度。

由第七章可知,用窗函数法设计的滤波器性能取决于窗函数)(n ω的类型及窗口长度N 的取值。

设计过程中,要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N 。

各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见第七章。

这样选定窗函数类型和长度N 后,求出单位脉冲响应)(n h =)(n h d ·)(n ω,并按式(2-3)求出)(ωj e H 。

)(ωj e H 是否满足要求,要进行验算。

一般在)(n h 尾部加零使长度满足于2的整数次幂,以便用FFT 计算)(ωj e H 。

数字信号处理上机实验 作业结果与说明 实验三、四、五

数字信号处理上机实验 作业结果与说明 实验三、四、五

上机频谱分析过程及结果图 上机实验三:IIR 低通数字滤波器的设计姓名:赵晓磊 学号:赵晓磊 班级:02311301 科目:数字信号处理B一、实验目的1、熟悉冲激响应不变法、双线性变换法设计IIR 数字滤波器的方法。

2、观察对实际正弦组合信号的滤波作用。

二、实验内容及要求1、分别编制采用冲激响应不变法、双线性变换法设计巴特沃思、切贝雪夫I 型,切贝雪夫II 型低通IIR 数字滤波器的程序。

要求的指标如下:通带内幅度特性在低于πω3.0=的频率衰减在1dB 内,阻带在πω6.0=到π之间的频率上衰减至少为20dB 。

抽样频率为2KHz ,求出滤波器的单位取样响应,幅频和相频响应,绘出它们的图,并比较滤波性能。

(1)巴特沃斯,双线性变换法Ideal And Designed Lowpass Filter Magnitude Responsefrequency in Hz|H [e x p (j w )]|frequency in pi units|H [ex p (j w )]|Designed Lowpass Filter Phase Response in radians frequency in pi unitsa r g (H [e x p (j w )](2)巴特沃斯,冲激响应不变法(3)切贝雪夫I 型,双线性变换法(4)切贝雪夫Ⅱ型,双线性变换法综合以上实验结果,可以看出,使用不同的模拟滤波器数字化方法时,滤波器的性能可能产生如下差异:使用冲击响应不变法时,使得数字滤波器的冲激响应完全模仿模拟滤波器的冲激响应,也就是时域逼急良好,而且模拟频率和数字频率之间呈线性关系;但频率响应有混叠效应。

frequency in Hz|H [e x p (j w )]|Designed Lowpass Filter Magnitude Response in dBfrequency in pi units|H [e x p (j w )]|frequency in pi unitsa r g (H [e x p (j w )]Ideal And Designed Lowpass Filter Magnitude Responsefrequency in Hz|H [e x p (j w )]|frequency in pi units|H [e xp (j w )]|frequency in pi unitsa r g (H [e x p (j w )]Ideal And Designed Lowpass Filter Magnitude Responsefrequency in Hz|H [e x p (j w )]|frequency in pi units|H [ex p (j w )]|Designed Lowpass Filter Phase Response in radiansfrequency in pi unitsa r g (H [e x p (j w )]使用双线性变换法时,克服了多值映射的关系,避免了频率响应的混叠现象;在零频率附近,频率关系接近于线性关系,高频处有较大的非线性失真。

数字信号处理实验五谱分析

数字信号处理实验五谱分析

用FFT 对信号作频谱分析1.实验目的学习用FFT 对连续信号和时域离散信号进行谱分析的方法,了解可能出现的分析 误差及其原因,以便正确应用FFT 。

2. 实验原理用FFT 对信号作频谱分析是学习数字信号处理的重要内容。

经常需要进行谱分析的信号是模拟信号和时域离散信号。

对信号进行谱分析的重要问题是频谱分辨率D 和分析误差。

频谱分辨率直接和FFT 的变换区间N 有关,因为FFT 能够实现的频率分辨率是N /2π,因此要求D N ≤/2π。

可以根据此式选择FFT 的变换区间N 。

误差主要来自于用FFT 作频谱分析时,得到的是离散谱,而信号(周期信号除外)是连续谱,只有当N 较大时离散谱的包络才能逼近于连续谱,因此N 要适当选择大一些。

周期信号的频谱是离散谱,只有用整数倍周期的长度作FFT ,得到的离散谱才能代表周期信号的频谱。

如果不知道信号周期,可以尽量选择信号的观察时间长一些。

对模拟信号进行谱分析时,首先要按照采样定理将其变成时域离散信号。

如果是模拟周期信号,也应该选取整数倍周期的长度,经过采样后形成周期序列,按照周期序列的谱分析进行。

3.实验步骤及内容(1)对以下序列进行谱分析。

⎪⎩⎪⎨⎧≤≤-≤≤-=⎪⎩⎪⎨⎧≤≤-≤≤+==其它nn n n n n x 其它n n n n n n x n R n x ,074,330,4)(,074,830,1)()()(3241 选择FFT 的变换区间N 为8和16 两种情况进行频谱分析。

分别打印其幅频特性曲线。

并进行对比、分析和讨论。

(2)对以下周期序列进行谱分析。

4()cos 4x n n π=5()cos(/4)cos(/8)x n n n ππ=+选择FFT 的变换区间N 为8和16 两种情况分别对以上序列进行频谱分析。

分别打印其幅频特性曲线。

并进行对比、分析和讨论。

(3)对模拟周期信号进行谱分析6()cos8cos16cos20x t t t t πππ=++选择 采样频率Hz F s 64=,变换区间N=16,32,64 三种情况进行谱分析。

数字信号处理实验报告1-5

数字信号处理实验报告1-5

实验一时域离散信号的产生及时域处理实验目的:了解Matlab软件数字信号处理工具箱的初步使用方法。

掌握其简单的Matlab语言进行简单的时域信号分析。

实验内容:[1.1]已知两序列x1=[0,1,2,3,4,3,2,1,0];n1=[-2:6];x2=[2,2,0,0,0,-2,-2],n2=[2:8].求他们的和ya及乘积yp. 程序如下:x1=[0,1,2,3,4,3,2,1,0];ns1=-2;x2=[2,2,0,0,0,-2,-2];ns2=2;nf1=ns1+length(x1)-1;nf2=ns2+length(x2)-1;ny=min(ns1,ns2):max(nf1,nf2);xa1=zeros(1,length(ny));xa2=xa1;xa1(find((ny>=ns1)&(ny<=nf1)==1))=x1;xa2(find((ny>=ns2)&(ny<=nf2)==1))=x2;ya=xa1+xa2yp=xa1.*xa2subplot(4,4,1),stem(ny,xa1,'.')subplot(4,1,2),stem(ny,xa2,'.')line([ny(1),ny(end)],[0,0])subplot(4,1,3),stem(ny,ya,'.')line([ny(1),ny(end)],[0,0])subplot(4,1,4),stem(ny,yp,'.')line([ny(1),ny(end)],[0,0])[1.2]编写产生矩形序列的程序。

并用它截取一个复正弦序列,最后画出波形。

程序如下:clear;close alln0=input('输入序列起点:n0=');N=input('输入序列长度:N=');n1=input('输入位移:n1=');n=n0:n1+N+5;u=[(n-n1)>=0];x1=[(n-n1)>=0]-[(n-n1-N)>=0];x2=[(n>=n1)&(n<(N+n1))];x3=exp(j*n*pi/8).*x2;subplot(2,2,1);stem(n,x1,'.');xlabel('n');ylabel('x1(n)');axis([n0,max(n),0,1]);subplot(2,2,3);stem(n,x2,'.');xlabel('n');ylabel('x2(n)');axis([n0,max(n),0,1]);subplot(2,2,2);stem(n,real(x3),'.'); xlabel('n');ylabel('x3(n)的实部');line([n0,max(n)],[0,0]);axis([n0,max(n),-1,1]);subplot(2,2,4);stem(n,imag(x3),'.'); xlabel('n');ylabel('x3(n)的虚部');line([n0,max(n)],[0,0]);axis([n0,max(n),-1,1]);[1.3]利用已知条件,利用MATLAB生成图形。

数字信号处理第三版高西全实验

数字信号处理第三版高西全实验

数字信号处理第三版高西全实验数字信号处理第三版高西全实验》是一本旨在介绍数字信号处理的实际应用的教材。

本文档旨在概述该教材的目的和内容。

该教材的目的是通过实验教学的方式,帮助学生更好地理解数字信号处理的原理和技术,并将其应用到实际问题中。

它旨在培养学生的实践能力和解决问题的能力,使他们能够熟练地进行数字信号处理的实际操作。

该教材内容包括许多实验,涵盖了数字信号处理的各个方面。

每个实验都介绍了一个特定的概念或技术,并通过实际的操作和实验数据展示了其应用方式和效果。

学生通过完成实验,可以深入了解数字信号处理的各种算法和方法,研究如何使用相关工具和软件进行信号处理,以及如何分析和评估处理结果。

通过研究《数字信号处理第三版高西全实验》,学生将能够掌握数字信号处理的基本概念和技术,并能够独立地应用这些知识解决实际问题。

这将有助于他们在工程、通信、音视频处理等领域中的职业发展,也为进一步深入研究数字信号处理奠定了坚实的基础。

希望该教材能够对学生们的研究和实践有所帮助,使他们能够更好地理解和运用数字信号处理的方法和技术。

实验目标:本实验旨在介绍数字信号的采样和重构过程,并加深对这两个概念的理解。

实验目标:本实验旨在介绍数字信号的采样和重构过程,并加深对这两个概念的理解。

实验步骤:实验步骤:准备实验所需的信号发生器和示波器设备。

设置信号发生器,产生模拟信号,例如正弦波。

调整示波器参数,将模拟信号接入示波器进行显示。

使用采样器采样模拟信号,并记录采样得到的数字信号。

对采样得到的数字信号进行重构,恢复为原始模拟信号。

使用示波器将重构后的信号进行显示,并比较与原始信号的差异。

实验要点:了解采样和重构的基本概念和原理。

熟悉信号发生器和示波器的操作。

掌握采样器的使用方法。

理解数字信号与模拟信号的差异及其影响。

请参考实验指导书中的详细步骤和注意事项进行实验操作,并记录实验数据和结果。

本文档旨在解释《数字信号处理第三版高西全实验》中的实验二内容。

数字信号处理实验报告完整版[5篇模版]

数字信号处理实验报告完整版[5篇模版]

数字信号处理实验报告完整版[5篇模版]第一篇:数字信号处理实验报告完整版实验 1利用 T DFT 分析信号频谱一、实验目的1.加深对 DFT 原理的理解。

2.应用 DFT 分析信号的频谱。

3.深刻理解利用DFT 分析信号频谱的原理,分析实现过程中出现的现象及解决方法。

二、实验设备与环境计算机、MATLAB 软件环境三、实验基础理论T 1.DFT 与与 T DTFT 的关系有限长序列的离散时间傅里叶变换在频率区间的N 个等间隔分布的点上的 N 个取样值可以由下式表示:212 /0()|()()0 1Nj knjNk NkX e x n e X k k Nπωωπ--====≤≤-∑由上式可知,序列的 N 点 DFT ,实际上就是序列的 DTFT 在 N 个等间隔频率点上样本。

2.利用 T DFT 求求 DTFT方法 1 1:由恢复出的方法如下:由图 2.1 所示流程可知:101()()()Nj j n kn j nNn n kX e x n e X k W eNωωω∞∞----=-∞=-∞=⎡⎤==⎢⎥⎣⎦∑∑∑由上式可以得到:IDFT DTFT第二篇:数字信号处理实验报告JIANGSUUNIVERSITY OF TECHNOLOGY数字信号处理实验报告学院名称:电气信息工程学院专业:班级:姓名:学号:指导老师:张维玺(教授)2013年12月20日实验一离散时间信号的产生一、实验目的数字信号处理系统中的信号都是以离散时间形态存在的,所以对离散时间信号的研究是数字信号的基本所在。

而要研究离散时间信号,首先需要产生出各种离散时间信号。

使用MATLAB软件可以很方便地产生各种常见的离散时间信号,而且它还具有强大绘图功能,便于用户直观地处理输出结果。

通过本实验,学生将学习如何用MATLAB产生一些常见的离散时间信号,实现信号的卷积运算,并通过MATLAB中的绘图工具对产生的信号进行观察,加深对常用离散信号和信号卷积和运算的理解。

数字信号处理实验1--5含代码

数字信号处理实验1--5含代码

数字信号处理实验1--5含代码实验一离散时间信号的时域分析 1. 在MATLAB中利用逻辑关系式n,,0来实现序列,显示范围。

(产生如下,,,n,nn,n,n012图所示的单位脉冲信号的函数为impseq(n0,n1,n2),程序如示例所示),3,n,10并利用impseq函数实现序列:; ,,,,,,yn,2,n,3,,n,6,,xn1nnnn120源代码:impseq.mfunction y=impseq(n0,n1,n2)n=[n1:n2]y=[(n-n0)==0]exp01-1.mfunction impseq(n0,n1,n2)n=-3:1:10y=2*impseq(3,-3,10)+impseq(6,-3,10);stem(n,y)n,,0,,2. 在MATLAB中利用逻辑关系式来实现序列,显示范围。

(自己编写un,nn,n,n012产生单位阶跃信号的函数,函数命名为stepseq(n0,n1,n2)) 并利用编写的stepseq函数实现序列: ,,,,,,yn,un,2,un,2,5,n,10源代码:stepseq.mfunction y=stepseq(n0,n1,n2)n=n1:1:n2y=[(n-n0)>=0]exp01-2.mfunction stepseq(n0,n1,n2)n=-5:1:20y=stepseq(-2,-5,20)+stepseq(2,-5,20)stem(n,y)3. 在MATLAB中利用数组运算符“.^”来实现一个实指数序列。

如: n ,,,,xn,0.30,n,15源代码:n=0:1:15;x=0.3.^nstem(n,x)4. 在MATLAB中调用函数sin或cos产生正余弦序列,如:π,, ,,,,xn,3sin0.4πn,,5cos0.3πn0,n,20,,5,,源代码:n=0:1:20x=11*sin(0.3*pi*n+pi/5)+5*cos(0.3*pi*n)stem(n,x)思考题:1.在MATLAB环境下产生单位脉冲序列和单位阶跃序列各有几种方法,如何使用,2.在MATLAB环境下进行序列的相乘运算时应注意什么问题,实验二离散时间系统的时域分析1. 在MATLAB中利用内部函数conv来计算两个有限长序列的卷积。

数字信号处理实验报告11-12-10

数字信号处理实验报告11-12-10

《数字信号处理》实验报告专业学号姓名实验一 利用FFT 实现快速卷积一、实验目的1.加深理解FFT 在实现数字滤波(或快速卷积)中的重要作用,更好的利用FFT 进行数字信号处理。

2.掌握循环卷积和线性卷积两者之间的关系。

二、实验原理用FFT 来快速计算有限长度序列的线性卷积。

这种方法就是先将输入信号x(n)通过FFT 变换为它的频谱采样值()x k ,然后再和FIR 滤波器的频响采样值H(k)相乘,H(k)可事先存放在存储器中,最后再将乘积H(k)X(k)通过快速傅里叶变换(简称IFFT )还原为时域序列,即得到输出y(n)。

现以FFT 求有限长序列的卷积及求有限长度序列与较长序列间的卷积为例来讨论FFT 的快速卷积方法。

序列x(n)和h(n)的长差不多。

设x(n)的长为N 1,h(n)的长为N 2,要求∑-=-=⊗=1)()()()()(N m m n x m h n y n x n y用FFT 完成这一卷积的具体步骤如下:①为使两有限长序列的线性卷积可用其循环卷积代替而不发生混叠,必须选择循环卷积长度121-+≥N N N ,若采用基2-FFT 完成卷积运算,要求m N 2=(m 为整数)。

②用补零方法使x(n)和h(n)变成列长为N 的序列。

1122()01()01()01()01x n n N x n N n N h n n N h n N n N ≤≤-⎧=⎨≤≤-⎩≤≤-⎧=⎨≤≤-⎩③用FFT 计算x(n)和h(n)的N 点离散傅里叶变换。

④完成X(k)和H(k)乘积,)()()(k H k x k Y = ⑤用FFT 计算 ()Y k 的离散傅里叶反变换得*10*10)(1)(1)(⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=∑∑-=--=N k nk N nk N N k W k Y N W k Y N n y三、主要实验仪器及材料微型计算机、Matlab6.5教学版。

四、实验内容1.数字滤波器的脉冲响应为()22()1/2(),8nN h n R n N ==。

数字信号处理实验五报告

数字信号处理实验五报告

实验5F I R滤波器的设计一、实验目的1.掌握用窗函数法,频率采样法及优化设计法设计FIR 滤波器的原理及方法。

2.熟悉线性相位FIR 滤波器的幅频特性和相频特性。

3.了解各种不同窗函数对滤波器性能的影响。

二、实验内容a)N=45,计算并画出矩形窗、汉明窗、布莱克曼窗的归一化的幅度谱,并比较各自的主要特点。

各自特点:矩形窗函数具有最窄的主瓣宽度,但有最大的旁瓣峰值;汉明窗函数的主瓣稍宽,而旁瓣较小;布莱克曼窗函数主瓣最宽,旁瓣最小。

矩形窗设计的滤波器过渡带最窄,但是阻带最小衰减也最差;布莱克曼窗设计的滤波器阻带衰减最好,过渡带最宽,约为矩形窗设计的三倍。

汉明窗设计的滤波器处于以上二者之间。

b)N=15,带通滤波器的两个通带边界分别是ω1=π,ω2=π。

用汉宁窗设计此线性相位带通滤波器,观察它的实际3dB 和20dB 带宽。

N=45,重复这一设计,观察幅频和相位特性的变化,注意长度N 变化的影响。

N增加,3db带宽和20db带宽分别减小,滤波器特性变好,过渡带变陡,幅频曲线显示其通带较平缓,波动小,阻带衰减大,相频特性曲线显示其相位随频率变化也变大。

c)分别改用矩形窗和布莱克曼窗,设计(2)中的带通滤波器,观察并记录窗函数对滤波器幅频特性的影响,比较三种窗的特点。

矩形窗设计的滤波器过渡带最窄,但阻带最小衰减也最差;汉宁窗设计的滤波器过渡带稍宽,但有较好的阻带衰减;布莱克曼窗设计的滤波器阻带衰减最好,但过渡带最宽。

当使用同种窗设计滤波器时,N越大,主瓣宽度越窄,通带越平坦,过渡带宽越小。

对于同一个N值,当用不同窗设计时,矩形窗的过渡带最窄,但阻带衰减最差;布莱克曼窗的阻带衰减最好,但过渡带最宽;汉明窗的两种特性介于前两者之间。

d)用Kaiser 窗设计一专用线性相位滤波器,N=40,当β=4、6、10 时,分别设计、比较它们的幅频和相频特性,注意β取不同值时的影响。

由图中可以看出,β越大,则窗越窄,过渡带宽越大,主瓣的宽度也相应增加,而频谱的旁瓣越小,阻带最小衰减也越大。

数字信号处理实验报告

数字信号处理实验报告

数字信号处理实验报告数字信号处理实验报告一、实验目的本实验旨在通过数字信号处理的方法,对给定的信号进行滤波、频域分析和采样率转换等操作,深入理解数字信号处理的基本原理和技术。

二、实验原理数字信号处理(DSP)是一种利用计算机、数字电路或其他数字设备对信号进行各种处理的技术。

其主要内容包括采样、量化、滤波、变换分析、重建等。

其中,滤波器是数字信号处理中最重要的元件之一,它可以用来提取信号的特征,抑制噪声,增强信号的清晰度。

频域分析是指将时域信号转化为频域信号,从而更好地理解信号的频率特性。

采样率转换则是在不同采样率之间对信号进行转换,以满足不同应用的需求。

三、实验步骤1.信号采集:首先,我们使用实验室的信号采集设备对给定的信号进行采集。

采集的信号包括噪声信号、含有正弦波和方波的混合信号等。

2.数据量化:采集到的信号需要进行量化处理,即将连续的模拟信号转化为离散的数字信号。

这一步通常通过ADC(模数转换器)实现。

3.滤波处理:将量化后的数字信号输入到数字滤波器中。

我们使用不同的滤波器,如低通、高通、带通等,对信号进行滤波处理,以观察不同滤波器对信号的影响。

4.频域分析:将经过滤波处理的信号进行FFT(快速傅里叶变换)处理,将时域信号转化为频域信号,从而可以对其频率特性进行分析。

5.采样率转换:在进行上述处理后,我们还需要对信号进行采样率转换。

我们使用了不同的采样率对信号进行转换,并观察采样率对信号处理结果的影响。

四、实验结果及分析1.滤波处理:经过不同类型滤波器处理后,我们发现低通滤波器可以有效抑制噪声,高通滤波器可以突出高频信号的特征,带通滤波器则可以提取特定频率范围的信号。

这表明不同类型的滤波器在处理不同类型的信号时具有不同的效果。

2.频域分析:通过FFT处理,我们将时域信号转化为频域信号。

在频域分析中,我们可以更清楚地看到信号的频率特性。

例如,对于噪声信号,我们可以看到其频率分布较为均匀;对于含有正弦波和方波的混合信号,我们可以看到其包含了不同频率的分量。

数字信号处理实验报告MATLAB

数字信号处理实验报告MATLAB

数字信号处理实验报告姓名:班级:09电信一班学号:2)]得下图二,图二图一3.将如下文件另存为:sigadd.m文件function [y,n] = sigadd(x1,n1,x2,n2)% 实现y(n) = x1(n)+x2(n)% -----------------------------% [y,n] = sigadd(x1,n1,x2,n2)% y = 在包含n1 和n2 的n点上求序列和,% x1 = 在n1上的第一序列% x2 = 在n2上的第二序列(n2可与n1不等)n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的长度y1 = zeros(1,length(n)); y2 = y1; % 初始化y1(find((n>=min(n1))&(n<=max(n1))==1))=x1; % 具有y的长度的x1y2(find((n>=min(n2))&(n<=max(n2))==1))=x2; % 具有y的长度的x2y = y1+y2;在命令窗口输入:x1=[1,0.5,0.3,0.4];n1=-1:2;x2=[0.2,0.3,0.4,0.5,0.8,1];n2=-2:3; [y,n] = sigadd(x1,n1,x2,n2)得:y =n=-1:10;x=sin(0.4*pi*n);y=fliplr(x);n1=-fliplr(n);subplot(2,1,1),stem(n,x) subplot(2,1,2),stem(n1,y在命令窗口键入:n=-1:10; x=sin(0.4*pi*n);n (samples)实验结果:1.(1)在命令窗口输入:tic; [am,pha]=dft1(x)N=length(x);w=exp(-j*2*pi/N);for k=1:Nsum=0;for n=1:Nsum=sum+x(n)*w^((k-1)*(n-1));endam(k)=abs(sum);pha(k)=angle(sum);end;toc得到如下结果:am =Columns 1 through 11120.0000 41.0066 20.9050 14.3996 11.3137 9.6215 8.6591 8.1567 8.0000 8.1567 8.6591Columns 12 through 169.6215 11.3137 14.3996 20.9050 41.0066pha =Columns 1 through 110 1.7671 1.9635 2.1598 2.3562 2.5525 2.7489 2.9452 3.1416 -2.9452 -2.7489Columns 12 through 16-2.5525 -2.3562 -2.1598 -1.9635 -1.7671Elapsed time is 0.047000 seconds.(2)在命令窗口输入:tic;[am,pha]=dft2(x)N=length(x);n=[0:N-1];k=[0:N-1];w=exp(-j*2*pi/N);nk=n’*k;wnk=w.^(nk); Xk=x*wnk; am= abs(Xk); pha=angle(Xk); toc得到下图:figure(1)00.10.20.30.40.50.60.70.80.91signal x(n), 0 <= n <= 99(2)在命令窗口键入:n3=[0:1:99];y3=[x(1:1:10) zeros(1,90)]; %添90个零。

数字信号处理实验报告_五个实验

数字信号处理实验报告_五个实验

实验一 信号、系统及系统响应一、 实验目的1、熟悉连续信号经理想采样前后的频谱变化关系,加深对时域采样定理的理解;2、熟悉时域离散系统的时域特性;3、利用卷积方法观察分析系统的时域特性;4、掌握序列傅立叶变换的计算机实现方法,利用序列的傅立叶变换对连续信号、离散信号及系统响应进行频域分析。

二、 实验原理及方法采样是连续信号数字处理的第一个关键环节。

对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件,而且可以加深对傅立叶变换、Z 变换和序列傅立叶变换之间关系式的理解。

对一个连续信号)(t x a 进行理想采样的过程可用下式表示:)()()(^t p t t x x aa其中)(^t x a 为)(t x a 的理想采样,p(t)为周期脉冲,即∑∞-∞=-=m nT t t p )()(δ)(^t x a的傅立叶变换为)]([1)(^s m a m j X T j a X Ω-Ω=Ω∑∞-∞= 上式表明^)(Ωj Xa为)(Ωj Xa的周期延拓。

其延拓周期为采样角频率(T /2π=Ω)。

只有满足采样定理时,才不会发生频率混叠失真。

在实验时可以用序列的傅立叶变换来计算^)(Ωj X a 。

公式如下:Tw jwae X j X Ω==Ω|)()(^离散信号和系统在时域均可用序列来表示。

为了在实验中观察分析各种序列的频域特性,通常对)(jw e X 在[0,2π]上进行M 点采样来观察分析。

对长度为N 的有限长序列x(n),有:n jw N n jw k ke m x eX--=∑=)()(1其中,k Mk πω2=,k=0,1,……M-1 时域离散线性非移变系统的输入/输出关系为∑∞-∞=-==m m n h m x n h n x n y )()()(*)()(上述卷积运算也可在频域实现)()()(ωωωj j j e H e X eY =三、 实验程序s=yesinput(Please Select The Step Of Experiment:\n 一.(1时域采样序列分析 s=str2num(s); close all;Xb=impseq(0,0,1); Ha=stepseq(1,1,10);Hb=impseq(0,0,3)+2.5*impseq(1,0,3)+2.2*impseq(2,0,3)+impseq(3,0,3); i=0;while(s);%时域采样序列分析 if(s==1)k=0;while(1)if(k==0)A=yesinput('please input the Amplitude:\n',...444.128,[100,1000]); a=yesinput('please input the Attenuation Coefficient:\n',...222.144,[100,600]);w=yesinput('please input the Angle Frequence(rad/s):\n',...222.144,[100,600]);endk=k+1;fs=yesinput('please input the sample frequence:\n',...1000,[100,1200]);Xa=FF(A,a,w,fs);i=i+1;string+['fs=',num2str(fs)];figure(i)DFT(Xa,50,string);1=yesinput1=str2num(1);end%系统和响应分析else if(s==2)kk=str2num(kk);while(kk)if(kk==1)m=conv(Xb,Hb);N=5;i=i+1;figure(i)string=('hb(n)');Hs=DFT(Hb,4,string);i=i+1;figure(i)string('xb(n)');DFT(Xb,2,string);string=('y(n)=xb(n)*hb(n)');else if (kk==2)m=conv(Ha,Ha);N=19;string=('y(n)=ha(n)*(ha(n)');else if (kk==3)Xc=stepseq(1,1,5);m=conv(Xc,Ha);string=('y(n)=xc(n)*ha(n)');endendendi=i+1;figure(i)DFT(m,N,string);kk=yesinputkk=str2num(kk);end卷积定理的验证else if(s==3)A=1;a=0.5;w=2,0734;fs=1;Xal=FF(A,a,w,fs);i=i+1;figure(i)string=('The xal(n)(A=1,a=0.4,T=1)'); [Xa,w]DFT(Xal,50,string);i=i+1;figure(i)string =('hb(n)');Hs=DFT(Hb,4,string);Ys=Xs.*Hs;y=conv(Xal,Hb);N=53;i=i+1;figure(i)string=('y(n)=xa(n)*hb(n)');[yy,w]=DFT(y,N,string);i=i+1;figure(i)subplot(2,2,1)plot(w/pi,abs(yy));axis([-2 2 0 2]);xlabel('w/pi');ylabel('|Ys(jw)|');title(FT[x(n)*h(n)]');subplot(2,2,3)plot(w/pi,abs(Ys));axis([-2 2 0 2]);xlabel('w/pi');ylabel('|Ys(jw)|');title('FT[xs(n)].FT[h(n)]'); end end end子函数:离散傅立叶变换及X(n),FT[x(n)]的绘图函数 function[c,l]=DFT(x,N,str) n=0:N-1; k=-200:200; w=(pi/100)*k; l=w; c=x*Xc=stepseq(1,1,5); 子函数:产生信号function c=FF(A,a,w,fs) n=o:50-1;c=A*exp((-a)*n/fs).*sin(w*n/fs).*stepseq(0,0,49); 子函数:产生脉冲信号function [x,n]=impseq(n0,n1,n2) n=[n1:n2];x=[(n-n0)==0];子函数:产生矩形框信号function [x,n]=stepseq(n0,n1,n2) n=[n1:n2];x=[(n-n0>=0)];四、 实验内容及步骤1、认真复习采样理论,离散信号与系统,线性卷积,序列的傅立叶变换及性质等有关内容,阅读本实验原理与方法。

数字信号处理实验五用DFT(FFT)对信号进行频谱分析

数字信号处理实验五用DFT(FFT)对信号进行频谱分析

开课学院及实验室:电子楼3172018年 4月 29 日3()x n :用14()()x n R n =以8为周期进行周期性延拓形成地周期序列.(1> 分别以变换区间N =8,16,32,对14()()x n R n =进行DFT(FFT>,画出相应地幅频特性曲线;(2> 分别以变换区间N =4,8,16,对x 2(n >分别进行DFT(FFT>,画出相应地幅频特性曲线; (3> 对x 3(n >进行频谱分析,并选择变换区间,画出幅频特性曲线.<二)连续信号 1. 实验信号:1()()x t R t τ=选择 1.5ms τ=,式中()R t τ地波形以及幅度特性如图7.1所示.2()sin(2/8)x t ft ππ=+式中频率f 自己选择.3()cos8cos16cos 20x t t t t πππ=++2. 分别对三种模拟信号选择采样频率和采样点数.对1()x t ()R t τ=,选择采样频率4s f kHz =,8kHz ,16kHz ,采样点数用τ.s f 计算.对2()sin(2/8)x t ft ππ=+,周期1/T f =,频率f 自己选择,采样频率4s f f =,观测时间0.5p T T =,T ,2T ,采样点数用p s T f 计算.图5.1 R(t>地波形及其幅度特性对3()cos8cos16cos 20x t t t t πππ=++,选择采用频率64s f Hz =,采样点数为16,32,64. 3. 分别对它们转换成序列,按顺序用123(),(),()x n x n x n 表示.4. 分别对它们进行FFT.如果采样点数不满足2地整数幂,可以通过序列尾部加0满足.5. 计算幅度特性并进行打印.五、实验过程原始记录<数据、图表、计算等)(一> 离散信号%14()()x n R n = n=0:1:10。

数字信号处理--实验五-用DFT(FFT)对信号进行频谱分析

数字信号处理--实验五-用DFT(FFT)对信号进行频谱分析

学生实验报告开课学院及实验室:电子楼3172013年4月29日、实验目的学习DFT 的基本性质及对时域离散信号进行频谱分析的方法,进一步加深对频域概念和数字频率的理解,掌握 MATLAB 函数中FFT 函数的应用。

二、实验原理离散傅里叶变换(DFT)对有限长时域离散信号的频谱进行等间隔采样,频域函数被离散化了, 便于信号的计算机处理。

设x(n)是一个长度为 M 的有限长序列,x(n)的N 点傅立叶变换:X(k)N 1j 三 knDFT[x(n)]N x(n)e N0 k N 1n 0其中WNe.2 jN,它的反变换定义为:1X(n)NkN 1nkX(k)W N0 令z W N k,X(zz WN k则有:N 1x( n)Wj kn 0可以得到,X(k)X(Z)Z WN kZ W N*是Z 平面单位圆上幅角为2kN 的点,就是将单位圆进行N 等分以后第 K 个点。

所以, X(K)是Z 变换在单位圆上的等距采样,或者说是序列傅立叶变换的等距采样。

时域采样在满足Nyquist 定理时,就不会发生频谱混叠。

DFT 是对序列傅立叶变换的等距采样,因此可以用于序列的频谱分析。

如果用FFT 对模拟信号进行谱分析,首先要把模拟信号转换成数字信号,转换时要求知道模拟 信号的最高截至频率,以便选择满足采样定理的采样频率。

般选择采样频率是模拟信号中最高频率的3~4倍。

另外要选择对模拟信号的观测时间,如果采样频率和观测时间确定,则采样点数也确定 了。

这里观测时间和对模拟信号进行谱分析的分辨率有关,最小的观测时间和分辨率成倒数关系。

最小的采样点数用教材相关公式确定。

要求选择的采样点数和观测时间大于它的最小值。

如果要进行谱分析的模拟信号是周期信号,最好选择观测时间是信号周期的整数倍。

如果不知道■ 厂1*1IE向i1A I1f Ii i 0r 1 疋0Jfb-4W0 70000图5.1 R(t)的波形及其幅度特性xn=[on es(1,4),zeros(1,7)];%输入时域序列向量 xn=R4( n)%计算xn 的8点DFTXk16=fft(x n,16);%计算xn 的16点DFTXk32=fft(x n,32); %计算xn 的32点DFTk=0:7;wk=2*k/8;对 x 3(t) cos8 t cos16 t cos20 t ,选择采用频率 f s 64Hz ,采样点数为 16 , 32 , 64。

数字信号处理实验五 龙志鹏

数字信号处理实验五   龙志鹏

5.1 clf;t=0:0.0005:1; f=13;xa=cos(2*pi*f*t); subplot(2,1,1) plot(t,xa);gridxlabel('时间,mesc');ylabel('振幅'); title('连续时间信号x_{a}(t)'); axis([0 1 -1.2 1.2]) subplot(2,1,2); T=0.1; n=0:T:1;xs=cos(2*pi*f*n); k=0:length(n)-1; stem(k,xs);gridxlabel('时间序号n');ylabel('振幅'); title('离散时间序号x[n]'); axis([0 (length(n)-1) -1.2 1.2])0.10.20.30.40.50.60.70.80.91-1-0.500.51时间,mesc 振幅连续时间信号x a (t)12345678910时间序号n振幅离散时间序号x[n]5.4低于抽样周期时: (1)clf; t=0:0.0005:1; f=13;xa=cos(2*pi*f*t); subplot(2,1,1) plot(t,xa);gridxlabel('时间,mesc');ylabel('振幅'); title('连续时间信号x_{a}(t)'); axis([0 1 -1.2 1.2]) subplot(2,1,2); T=0.01; n=0:T:1;xs=cos(2*pi*f*n); k=0:length(n)-1; stem(k,xs);gridxlabel('时间序号n');ylabel('振幅'); title('离散时间序号x[n]'); axis([0 (length(n)-1) -1.2 1.2])0.10.20.30.40.50.60.70.80.91-1-0.500.51时间,mesc 振幅连续时间信号x a (t)102030405060708090100时间序号n振幅离散时间序号x[n](2)clf; t=0:0.0005:1; f=13;xa=cos(2*pi*f*t); subplot(2,1,1) plot(t,xa);gridxlabel('时间,mesc');ylabel('振幅'); title('连续时间信号x_{a}(t)'); axis([0 1 -1.2 1.2])subplot(2,1,2); T=0.05; n=0:T:1;xs=cos(2*pi*f*n); k=0:length(n)-1; stem(k,xs);gridxlabel('时间序号n');ylabel('振幅'); title('离散时间序号x[n]'); axis([0 (length(n)-1) -1.2 1.2])0.10.20.30.40.50.60.70.80.91-1-0.500.51时间,mesc 振幅连续时间信号x a (t)2468101214161820时间序号n振幅离散时间序号x[n]高于抽样周期时: (3)clf; t=0:0.0005:1; f=13;xa=cos(2*pi*f*t); subplot(2,1,1) plot(t,xa);gridxlabel('时间,mesc');ylabel('振幅'); title('连续时间信号x_{a}(t)'); axis([0 1 -1.2 1.2]) subplot(2,1,2); T=0.5; n=0:T:1;xs=cos(2*pi*f*n); k=0:length(n)-1;stem(k,xs);gridxlabel('时间序号n');ylabel('振幅'); title('离散时间序号x[n]'); axis([0 (length(n)-1) -1.2 1.2])0.10.20.30.40.50.60.70.80.91-1-0.500.51时间,mesc 振幅连续时间信号x a (t)0.20.40.60.81 1.2 1.41.61.82时间序号n振幅离散时间序号x[n](4)clf; t=0:0.0005:1; f=13;xa=cos(2*pi*f*t); subplot(2,1,1) plot(t,xa);gridxlabel('时间,mesc');ylabel('振幅'); title('连续时间信号x_{a}(t)'); axis([0 1 -1.2 1.2]) subplot(2,1,2); T=0.3; n=0:T:1;xs=cos(2*pi*f*n); k=0:length(n)-1; stem(k,xs);gridxlabel('时间序号n');ylabel('振幅'); title('离散时间序号x[n]'); axis([0 (length(n)-1) -1.2 1.2])0.10.20.30.40.50.60.70.80.91-1-0.500.51时间,mesc 振幅连续时间信号x a (t)0.511.522.53时间序号n振幅离散时间序号x[n]5.6 clf; T=0.1; f=13;n=(0:T:1)';xs=cos(2*pi*f*n);t= linspace(-0.5,1.5,500)';ya=sinc((1/T)*t(:,ones(size(n)))-(1/T)*n(:,ones(size(t)))')*xs; plot(n,xs,'o',t,ya);grid;xlabel('时间,mesc');ylabel('振幅'); title('重构的时间序号'y_{a}(t)'); axis([0 1 -1.2 1.2])-0.50.51 1.5-1.5-1-0.50.511.5时间,mesc振幅5.10 lf;t=0:0.005:10; xa=2*t.*exp(-t); subplot(2,2,1) plot(t,xa);gridxlabel('时间,mesc');ylabel('振幅'); title('连续时间信号x_{a}(t)'); subplot(2,2,2) wa=0:10/511:10;ha=freqs(2,[1 2 1],wa);plot(wa/(2*pi),abs(ha));grid; xlabel('频率,kHz');ylabel('振幅'); title('|X_{a}(j/Omega)|'); axis([0 5/pi 0 2]); subplot(2,2,3) T=1; n=0:T:10;xs=2*n.*exp(-n); k=0:length(n)-1; stem(k,xs);grid;xlabel('时间序号n');ylabel('振幅'); title('离散时间序号x[n]'); subplot(2,2,4)wd=0:pi/225:pi; hd=freqs(xs,1,wd);plot(wd/(T*pi),T*abs(hd));grid; xlabel('频率,kHz');ylabel('振幅'); title('|X(e^{j\omega})|'); axis([0 1/T 0 2]); 5.145100.20.40.60.8时间,mesc 振幅连续时间信号x a (t)0.51 1.50.511.52频率,kHz振幅|X a (j/Omega)|时间序号n振幅离散时间序号x[n]0.5100.511.52频率,kHz振幅|X (e j )|5.14 clf;Fp=3500;Fs=4500;Wp=2*pi*Fp;Ws=2*pi*Fs;[N,Wn]=buttord(Wp,Ws,0.5,30,'s'); [b,a]=butter(N,Wn,'s'); wa=0:(3*Ws)/511:3*Ws; h=freqs(b,a,wa);plot(wa/(2*pi),20*log10(abs(h)));grid xlabel('频率,Hz');ylabel('增益,dB'); title('增益响应'); axis([0 3*Fs -60 5]);020004000600080001000012000-60-50-40-30-20-10频率,Hz增益,d B增益响应。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
b)采样频率 必须大于等于模拟信号最高频率的两倍以上,才能使采样信号的
频谱不产生频谱混叠。
利用计算机计算上式并不方便,下面我们导出另外一个公式,以便用计算机上进行实验。
理想采样信号 和模拟信号 之间的关系为:
对上式进行傅立叶变换,得到:
在上式的积分号内只有当 时,才有非零值,因此:
上式中,在数值上 = ,再将 代入,得到:
上式的右边就是序列的傅立叶变换 ,即
上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到,只要将自变量ω用 代替即可。
频域采样定理的要点是:
a)对信号x(n)的频谱函数X(ejω)在[0,2π]上等间隔采样N点,得到
则N点IDFT[ ]得到的序列就是原序列x(n)以N为周期进行周期延拓后的主值区序列,公式为:
(4)学会调用MATLAB函数设计与实现FIR滤波器
实验要求
(1)对两种设计FIR滤波器的方法(窗函数法和等波纹最佳逼近法)进行分析比较,简述其优缺点。
(2)附程序清单、打印实验内容要求绘图显示的曲线图。
(3)分析总结实验结果。
(4)简要回答思考题。
实验原理
a)对模拟信号 以间隔T进行时域等间隔理想采样,形成的采样信号的频谱 是原模拟信号频谱 以采样角频率 ( )为周期进行周期延拓。公式为:
③用等波纹最佳逼近法设计的滤波器,其通带和阻带均为等指标均匀分布,没有资源浪费,所以期阶数低得多。
指导教师意见
签名:
年月日
[Ne,fo,mo,W]=remezord(fb,m,dev,Fs);% 确定remez函数所需参数
hn=remez(Ne,fo,mo,W);% 调用remez函数进行设计
Hw=abs(fft(hn,1024));% 求设计的滤波器频率特性
yet=fftfilt(hn,xt,N); % 调用函数fftfilt对xt滤波
①用窗函数法设计的滤波器,如果在阻带截止频率附近刚好满足,则离开阻带截止频率越远,阻带衰减富裕量越大,即存在资源浪费;
②几种常用的典型窗函数的通带最大衰减和阻带最小衰减固定,且差别较大,又不能分别控制。所以设计的滤波器的通带最大衰减和阻带最小衰减通常都存在较大富裕。如本实验所选的blackman窗函数,其阻带最小衰减为74dB,而指标仅为60dB。
(4)根据滤波器指标选择合适的窗函数,计算窗函数的长度N,调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序,调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。
(4)重复(3),滤波器指标不变,但改用等波纹最佳逼近法,调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。
在数字信号处理的应用中,只要涉及时域或者频域采样,都必须服从这两个采样理论的要点。
对比上面叙述的时域采样原理和频域采样原理,得到一个有用的结论,这两个采样理论具有对偶性:“时域采样频谱周期延拓,频域采样时域信号周期延拓”。因此放在一起进行实验。
实验仪器
PC机,MATLAB软件
实验步骤
(1)认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理;
%以下为用等波纹设计法的绘图部分(滤波器损耗函数,滤波器输出信号yw(nT)波形)
%省略
实验数据
实验内容1:
实验内容2:
实验总结
用窗函数法设计滤波器,滤波器长度 Nb=184。滤波器损耗函数和滤波器输出yw(nT)分别如图(a)和(b)所示。
用等波纹最佳逼近法设计滤波器,滤波器长度 Ne=83。滤波器损耗函数和滤波器输出ye(nT)分别如图(c)和(d)所示。
图10.5.2实验程序框图
实验内容
%《数字信号处理(第三版)学习指导》第10章实验5程序
% FIR数字滤波器设计及软件实现
clear all;close all;
%==调用xtg产生信号xt, xt长度N=1000,并显示xt及其频谱,=========
N=1000;xt=xtg(N);
fp=120; fs=150;Rp=;As=60;Fs=1000;% 输入给定指标
(2)调用信号产生函数xtg产生具有加性噪声的信号xt,并自动显示xt及其频谱,如图10.5.1所示;
图10.5.1具有加性噪声的信号x(t)及其频谱如图
(3)请设计低通滤波器,从高频噪声中提取xt中的单频调幅信号,要求信号幅频失真小于,将噪声频谱衰减60dB。先观察xt的频谱,确定滤波器指标参数。
b)由上式可知,频域采样点数N必须大于等于时域离散信号的长度M(即N≥M),才能使时域不产生混叠,则N点IDFT[ ]得到的序列 就是原序列x(n),即 =x(n)。如果N>M, 比原序列尾部多N-M个零点;如果N<M,z则 =IDFT[ ]发生了时域混叠失真,而且 的长度N也比x(n)的长度M短,因此。 与x(n)不相同。
附录1:
贵州大学实验报告
学院: 电气工程学院 专业: 测控技术与仪器 班级:测仪131
姓名
杨凯
学号
实验组
实验时间
指导教师
成绩
实验项目名称
FIR数字滤波器设计与软件实现
实验目的
(1)掌握用窗函数法设计FIR数字滤波器的原理和方法。
(2)掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。
(3)掌握FIR滤波器的快速卷积实现原理。
采样频率Fs=1000Hz,采样周期T=1/Fs;
根据图10.6.1(b)和实验要求,可选择滤波器指标参数:通带截止频率fp=120Hz,阻带截至频率fs=150Hz,换算成数字频率,通带截止频率 ,通带最大衰为,阻带截至频率 ,阻带最小衰为60dB。]
实验程序框图如图10.5.2所示,供读者参考。
% (1) 用窗函数法设计滤波器
wc=(fp+fs)/Fs; %理想低通滤波器截止频率(关于pi归一化)
B=2*pi*(fs-fp)/Fs; %过渡带宽度指标
Nb=ceil(11*pi/B); %blackman窗的长度N
hn=fir1(Nb-1,wc,blackman(Nb));
Hw=abs(fft(hn,1024));% 求设计的滤波器频率特性
ywt=fftfilt(hn,xt,N); %调用函数fftfilt对xt滤波
%以下为用窗函数法设计法的绘图部分(滤波器损耗函数,滤波器输出信号波形)
%省略
% (2) 用等波纹最佳逼近法设计滤波器
fb=[fp,fs];m=[1,0];% 确定remezord函数所需参数f,m,dev
dev=[(10^(Rp/20)-1)/(10^(Rp/20)+1),10^(-As/20)];
两种方法设计的滤波器都能有效地从噪声中提取信号,但等波纹最佳逼近法设计的滤波器阶数低得多,当然滤波实现的运算量以及时延也小得多,从图(b)和(d)可以直观地看出时延差别。
(1)用窗函数法设计线性相位低通滤波器的设计步骤教材中有详细的介绍.
(2)希望逼近的理想带通滤波器的截止频率 分别为:
(3)解释为什么对同样的技术指标,用等波纹最佳逼近法设计的滤波器阶数低?
相关文档
最新文档