实验四用FFT作谱分析.doc_第1页
实验四用FFT作谱分析.doc_第2页
实验四用FFT作谱分析.doc_第3页
实验四用FFT作谱分析.doc_第4页
实验四用FFT作谱分析.doc_第5页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

数字信号处理实验报告（四）-用FFT作谱分析 学院： 班级： 姓名： 学号：实验四 用FFT作谱分析一实验目的：1：加深对DFT算法原理及基本性质的理解；2：熟悉FFT算法原理及子程序的应用；3：学习FFT对连续信号的和时域离散信号作谱分析的方法。了解可能出现的的分析误差和原因，以便在实际中正确应用FFT；二实验仪器、材料及软件： 计算机，MATLAB软件。三实验结果及分析： function xk=mydft(xn,N)n=0:1:N-1;k=n;WN=exp(-j*2*pi/N);nk=n*k;WNnk=WN.nk;xk=xn*WNnk;function y=myditfft(x)m=nextpow2(x);N=2m;if length(x)N; x=x,zeros(1,N-length(x);endnxd=bin2dec(fliplr(dec2bin(1:N-1,m)+1;y=x(nxd);for mm=1:m Nmr=2mm; u=1; WN=exp(-i*2*pi/Nmr);for j=1:Nmr/2 for k=j:Nmr:N kp=k+Nmr/2; t=y(kp)*u; y(kp)=y(k)-t; y(k)=y(k)+t; end u=u*WN;endendn=0:1:1023;t=n*0.001;x=6*sin(60*pi*t);y=fft(x);h=abs(y);a=angle(y);plot(h);figure(6);plot(a,k-);clc;clear alln=0:1023;dt=0.001;f2=20*sin(60*pi*n*dt)+5*cos(120*pi*n*dt);tic;a=mydft(f2,1024);toc;subplot(221);plot(a);tic;b=myditfft(f2);toc;subplot(222);plot(b);tic;c=fft(f2);toc;subplot(223);plot(c);clc;clear alln=0:1023;dt=0.001;x=20*sin(60*pi*n*dt)+5*cos(120*pi*n*dt);y=fft(x);h=abs(y);a=
人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论