MATLAB-FFT
e小白網(wǎng)址:www.e-xiaobai.com
%FFT實(shí)驗(yàn)數(shù)據(jù)為對連續(xù)信號4*cos(2*pi*35*t)+5*sin(2*pi*60*t)采樣所得?
%采樣時(shí)間間隔為0.005(根據(jù)奈奎斯特采樣定理而定)，共采集1001個(gè)數(shù)據(jù)?
%采樣得到的數(shù)據(jù)保存在FFT實(shí)驗(yàn)數(shù)據(jù).txt文件內(nèi)?
%?t=0.005:0.005:5;?
%?y=4*cos(2*pi*35*t)+5*sin(2*pi*60*t);
filename='FFT實(shí)驗(yàn)數(shù)據(jù).txt';
fileID=fopen(filename);
C=textscan(fileID,'%f %f','HeaderLines',1);
fclose(fileID);
t=C{1,1};
y=C{1,2};
N=length(t);%數(shù)據(jù)個(gè)數(shù)
T=t(2,1)-t(1,1);%時(shí)間間隔
fs=1/T;%采樣頻率
figure(1)
subplot(2,1,1)
plot(t,y)
title('時(shí)域')
xlabel('秒')
%FFT
%NFFT=N;
NFFT=2^nextpow2(N);
y1=fft(y,NFFT);
f=(0:NFFT-1)/NFFT*fs;
y1(1,1)=abs(y1(1,1))/2;%直流量除以2
subplot(2,1,2)
plot(f(1:NFFT/2+1),abs(y1(1:NFFT/2+1))*2/NFFT)
xlabel('赫茲')
ylabel('幅度')
title('幅頻譜')
注:FFT實(shí)驗(yàn)數(shù)據(jù).txt文件可在e小白官網(wǎng)《MATLAB-FFT》文章中下載。
標(biāo)簽：
