欢迎来到天天文库
浏览记录
ID:39576366
大小:41.50 KB
页数:5页
时间:2019-07-06
《实验四用FFT作谱分析》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、数字信号处理实验报告(四)----用FFT作谱分析学院:班级:姓名:学号:实验四用FFT作谱分析一.实验目的:1:加深对DFT算法原理及基本性质的理解;2:熟悉FFT算法原理及子程序的应用;3:学习FFT对连续信号的和时域离散信号作谱分析的方法。了解可能出现的的分析误差和原因,以便在实际中正确应用FFT;二.实验仪器、材料及软件:计算机,MATLAB软件。三.实验结果及分析:functionxk=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;functiony=myditfft(x)m=
2、nextpow2(x);N=2^m;iflength(x)3、a=angle(y);plot(h);figure(6);plot(a,'k-');clc;clearalln=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;clearalln=0:1023;dt=0.001;x=20*sin(60*pi*n*dt)+5*cos(4、120*pi*n*dt);y=fft(x);h=abs(y);a=angle(y);subplot(211)plot(h);subplot(212)plot(a);
3、a=angle(y);plot(h);figure(6);plot(a,'k-');clc;clearalln=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;clearalln=0:1023;dt=0.001;x=20*sin(60*pi*n*dt)+5*cos(
4、120*pi*n*dt);y=fft(x);h=abs(y);a=angle(y);subplot(211)plot(h);subplot(212)plot(a);
此文档下载收益归作者所有