傅里葉變換（fft）matlab程序三

clear all %清除內(nèi)存所有變量

close all %關(guān)閉所有打開的圖形窗口

%% 執(zhí)行FFT點(diǎn)數(shù)與原信號長度相等（100點(diǎn)）

% 構(gòu)建原信號

N=100; % 信號長度（變量@@@@@@@）

Fs=1; % 采樣頻率

dt=1/Fs; % 采樣間隔

t=［0:N-1］*dt; % 時(shí)間序列

xn=cos（2*pi*0.24*［0:99］）+cos（2*pi*0.26*［0:99］）;

xn=［xn，zeros（1，N-100）］; % 原始信號的值序列

subplot（3，2，1） % 變量@@@@@@@

plot（t，xn） % 繪出原始信號

xlabel（‘時(shí)間/s’），title（‘原始信號（向量長度為100）’） % 變量@@@@@@@

% FFT分析

NN=N; % 執(zhí)行100點(diǎn)FFT

XN=fft（xn，NN）/NN; % 共軛復(fù)數(shù)，具有對稱性

f0=1/（dt*NN）; % 基頻

f=［0:ceil（（NN-1）/2）］*f0; % 頻率序列

A=abs（XN）; % 幅值序列

subplot（3，2，2），stem（f，2*A（1:ceil（（NN-1）/2）+1）），xlabel（‘頻率/Hz’） % 繪制頻譜（變量@@@@@@@）

axis（［0 0.5 0 1.2］） % 調(diào)整坐標(biāo)范圍

title（‘執(zhí)行點(diǎn)數(shù)等于信號長度（單邊譜100執(zhí)行點(diǎn)）’）; % 變量@@@@@@@

%% 執(zhí)行FFT點(diǎn)數(shù)大于原信號長度

% 構(gòu)建原信號

N=100; % 信號長度（變量@@@@@@@）

Fs=1; % 采樣頻率

dt=1/Fs; % 采樣間隔

t=［0:N-1］*dt; % 時(shí)間序列

xn=cos（2*pi*0.24*［0:99］）+cos（2*pi*0.26*［0:99］）;

xn=［xn，zeros（1，N-100）］; % 原始信號的值序列

subplot（3，2，3） % 變量@@@@@@@

plot（t，xn） % 繪出原始信號

xlabel（‘時(shí)間/s’），title（‘原始信號（向量長度為100）’） % 變量@@@@@@@

% FFT分析

NN=120; % 執(zhí)行120點(diǎn)FFT（變量@@@@@@@）

XN=fft（xn，NN）/NN; % 共軛復(fù)數(shù)，具有對稱性

f0=1/（dt*NN）; % 基頻

f=［0:ceil（（NN-1）/2）］*f0; % 頻率序列

A=abs（XN）; % 幅值序列

subplot（3，2，4），stem（f，2*A（1:ceil（（NN-1）/2）+1）），xlabel（‘頻率/Hz’） % 繪制頻譜（變量@@@@@@@）

axis（［0 0.5 0 1.2］） % 調(diào)整坐標(biāo)范圍

title（‘執(zhí)行點(diǎn)數(shù)大于信號長度（單邊譜120執(zhí)行點(diǎn)）’）; % 變量@@@@@@@

%% 執(zhí)行FFT點(diǎn)數(shù)與原信號長度相等（120點(diǎn)）

% 構(gòu)建原信號

N=120; % 信號長度（變量@@@@@@@）

Fs=1; % 采樣頻率

dt=1/Fs; % 采樣間隔

t=［0:N-1］*dt; % 時(shí)間序列

xn=cos（2*pi*0.24*［0:99］）+cos（2*pi*0.26*［0:99］）;

xn=［xn，zeros（1，N-100）］; % 原始信號的值序列

subplot（3，2，5） % 變量@@@@@@@

plot（t，xn） % 繪出原始信號

xlabel（‘時(shí)間/s’），title（‘原始信號（向量長度為120）’） % 變量@@@@@@@

% FFT分析

NN=120; % 執(zhí)行120點(diǎn)FFT（變量@@@@@@@）

XN=fft（xn，NN）/NN; % 共軛復(fù)數(shù)，具有對稱性

f0=1/（dt*NN）; % 基頻

f=［0:ceil（（NN-1）/2）］*f0; % 頻率序列

A=abs（XN）; % 幅值序列

subplot（3，2，6），stem（f，2*A（1:ceil（（NN-1）/2）+1）），xlabel（‘頻率/Hz’） % 繪制頻譜（變量@@@@@@@）

axis（［0 0.5 0 1.2］） % 調(diào)整坐標(biāo)范圍

title（‘執(zhí)行點(diǎn)數(shù)等于信號長度（單邊譜120執(zhí)行點(diǎn)）’）; % 變量@@@@@@@

結(jié)果