实验五 快速Fourier变换(FFT)及应用.docVIP

  1. 1、本文档共6页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验五 快速Fourier变换(FFT)及应用 实验目的 验证频域采样定理。 在理论学习的基础上,通过本实验,加深对FFT的理解,熟悉MATLAB中的有关函数。 应用FFT对典型信号进行频谱分析。 了解应用FFT进行信号频谱分析过程中可能出现的问题,以便在实际中能够正确应用FFT。 应用FFT实现序列的线性卷积。 实验内容 验证频域采样定理。 利用MATLAB产生一个长度为N的三角波序列,并完成以下要求: 计算N=30时的64点DFT,并图示 对在[0,2??上进行??点抽样,得到的32点IDFT,即得到; 绘制出的波形,观察和的关系,并加以说明。 解:MATLAB程序清单如下: M=64; % 指定DFT点数 N=30; % 指定序列长度 n=0:N-1; xn=2*[0:N/2,N/2-1:-1:1]/N; % 产生幅度为1的N点三角波序列 Xk=fft(xn,M); % 计算Xk=DFT[x(n)]; Xk1=Xk(1:2:M); % 对Xk隔点抽取得到Xk1 xn1=ifft(Xk1); % 对Xk1作IDFT得到xn1 n1=0:2*M; xc=xn1(mod(n1,M/2)+1); % 对xn1以M/2为周期进行延拓 subplot(2,2,1);stem(n,xn,.);grid; title([num2str(M/2) 点三角波序列x(n)]); subplot(2,2,2);k=0:M-1;stem(k,abs(Xk),.);grid; axis([0,M,0,max(Xk)]); title([三角波序列x(n)的 num2str(M) 点DFT:X(k)]); subplot(2,2,4);k1=0:M/2-1;stem(k1,abs(Xk1),.);grid; axis([0,M/2,0,max(Xk)]); title([隔点抽取X(k)得到 num2str(M/2) 点DFT:X_1(k)]); subplot(2,2,3);stem(n1,xc,.);grid; axis([0,2*M,0,max(xn1)]); title(序列x_1(n)的周期延拓); 由程序运行结果可以看出,在频域[0,2??上采样点数小于离散序列的长度时,将产生时域混叠,不能由。只有当频域采样点数大于等于序列长度时,才能由频域采样无失真的恢复出原序列,即。 已知x(n) 是N1点序列,序列的有值区间为[0, N1-1],h(n)是N2 点序列,序列的有值区间为[0,N2-1],现将序列右移m 位,即序列的有值区间变为[m, m+N1-1],然后对这序列做N 点圆周卷积得y(n),试问y(n)中哪个n 值的范围对应于x(n)* h(n)的结果。 m=2; N=7; xn=[1 2 2 1 -1 2]; hn=ones(1,4); y1=conv(xn,hn) % xn和hn的线性卷积 xn=[xn zeros(1,N-length(xn))]; hn1=[hn zeros(1,N-length(hn))]; yc1=ifft(fft(xn).*fft(hn1)) % xn和hn的N点圆周卷积 hn2=[zeros(1,m) hn]; xn=[xn zeros(1,N-length(xn))]; hn2=[hn2 zeros(1,N-length(hn2))]; yc2=ifft(fft(xn).*fft(hn2)) % xn和移位后的hn的N点圆周卷积 理解高密度谱和高分辨率频谱的概念。 设 取 将(1)中的补零加长到 增加抽样点的数目,取 ???L??程序清单如下: ?定义序列x(n) function y=xn(n) y=cos(0.48*pi*n)+cos(0.*pi*n)+cos(0.51*pi*n); % 主程序清单 clf xn=inline(cos(0.485*pi*n)+cos(0.51*pi*n)+cos(0.52*pi*n),n); % 定义一个inline局部函数 N1=10; N2=100; n1=0:N1-1;k1=n1; n2=0:N2-1;k2=n2; w1=k1*2*pi/N1; w2=k2*2*pi/N2; w3=w2; x1=xn(n1); Xk1=fft(x1); x2=[x1 zeros(1,N2-N1)]; Xk2=fft(x2); x3=xn(n2); Xk3=fft(x3); subplot(3,2,1) stem(n1,x1,.);grid; title(信号x_1(n),0=n=9); subplot(3,2,2) stem(w1/pi, abs(Xk1),.);grid; title(DFT[x_1

文档评论(0)

bhl0572 + 关注
实名认证
文档贡献者

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

1亿VIP精品文档

相关文档