数字信号处理实验精编.doc

  1. 1、本文档共20页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
PAGE  PAGE - 20 - 目 录 实验1 离散时间信号的频域分析2 实验2 FFT算法与应用7 实验3 IIR数字滤波器的设计12 实验4 FIR数字滤波器的设计17 实验1 离散时间信号的频域分析 一.实验目的 信号的频域分析是信号处理中一种有效的工具。在离散信号的时域分析中,通常将信号表示成单位采样序列(n)的线性组合,而在频域中,将信号表示成复变量e或 e的线性组合。通过这样的表示,可以将时域的离散序列映射到频域以便于进一步的处理。 在本实验中,将学习利用MATLAB计算离散时间信号的DTFT和DFT,并加深对其相互关系的理解。 实验原理 (1)DTFT和DFT的定义及其相互关系。序列x(n)DTFT定义为= 它是关于自变量的复函数,且是以2为周期的连续函数。可以表示为,其中,和分别是实部和虚部;还可以表示为 =,其中, 和分别是的幅度函数和相位函数;它们都是的实函数,也是以2为周期的周期函数。 序列的N点DFT定义为,是周期为N的序列。与的关系:是对)在一个周期中的谱的等间隔N点采样,即 ,而可以通过对内插获得,即 (2)使用到的MATLAB命令有基于DTFT离散时间信号分析函数以及求解序列的DFT函数。 1)基于 DTFT离散时间信号分析函数有:freqz,real,imag,abs,angle。函数freqz可以用来计算一个以e的有理分式形式给出的序列的DTFT值。freqz的形式多样,常见的有H=freqz(num,den,w),其中num表示序列有理分式DTFT的分子多项式系数,den表示分母多项式系数(均按z的降幂排列),矢量w表示在0~2中给定的一系列频率点集合。freqz函数的其他形式参见帮助文件。在求出DTFT值后,可以使用函数real,imag,abs和angle分别求出并绘出起实部、虚部和相位谱。例如)=利用函数freqz计算出,然后利用函数abs和angle分别求出幅频特性与相位特性最后利用plot命令绘出曲线。 求解序列DFT的函数有:fft,ifft。函数fft(x)可以计算R点序列的R点DFT值;而fft(x,N)则计算R点序列的N点DFT,若RN,则直接截取R点DFT的前N点,若RN,则x先进行补零扩展为N点序列再求N点DFT。函数ifft(X)可以计算R点的谱序列的R点IDFT值;而ifft(X,N)同fft(x,N)的情况。 三、实验设备 计算机、MATLAB软件 四、实验内容 (1)编程计算并画出下面DTFT的实部,虚部、幅度和相位谱。 程序如下:num=[0.1313 -0.1553 0.1313 0.0518]; den=[1 1.2828 1.0388 0.3418]; w=0:0.001:2*pi; H=freqz(num,den,w); figure subplot(221) plot(w,real(H)); title(实部) grid on; subplot(222) plot(w,imag(H)) title(虚部) grid subplot(223) plot(w,abs(H)) title(幅度) grid subplot(224) plot(w,abs(H)) title(相位谱) grid on; 波形如下: 图1-1 DTFT的实部,虚部、幅度和相位谱 (2)计算32点序列x(n)=cosn,0n31的32点和64点DFT,分别绘出幅度谱图形,并绘出该序列的DTFT图形。 程序如下:n1=0:15; n2=0:31; x1=cos((13*pi*n1)/32) x2=cos((13*pi*n2)/32) X1=fft(x1); X2=fft(x2); subplot(211) plot(n1,abs(X1)) title(32点DFT) grid on; subplot(212) plot(n2,abs(X2)) title(64点DFT) grid on; 波形如下: 图1-2 序列DFT幅度谱图形和DTFT图形 DTFT程序如下: A=1; n=0:31; x=cos((13*pi*n)/32); B=x; w=0:0.01:2*pi*2; [H]=freqz(B,A,w); magH=abs(H); subplot(2,1,1); plot(w,magH);grid; ylable(Magnitude); subplot(2,1,2); plot(w,phaH);grid; xlable(w); ylable(phase); grid on; 波形如下: 图1-3 DTFT图形 实验2 FFT算法与应用 一、实验目的 在理论学习的基础上,通过本次实验。加深对快速傅里叶

文档评论(0)

贪玩蓝月 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档