clf;
fs=100;N=128;???%采樣頻率和數(shù)據(jù)點(diǎn)數(shù)
n=0:N-1;t=n/fs;???%時間序列
x=0.5*sin（2*pi*15*t）+2*sin（2*pi*40*t）;?%信號
y=fft（xN）;????%對信號進(jìn)行快速Fourier變換
mag=abs（y）;?????%求得Fourier變換后的振幅
f=n*fs/N;????%頻率序列
subplot（221）plot（fmag）;???%繪出隨頻率變化的振幅
xlabel（‘頻率/Hz‘）;
ylabel（‘振幅‘）;title（‘N=128‘）;grid?on;
subplot（222）plot（f（1:N/2）mag（1:N/2））;?%繪出Nyquist頻率之前隨頻率變化的振幅
xlabel（‘頻率/Hz‘）;
ylabel（‘振幅‘）;title（‘N=128‘）;grid?on;
%對信號采樣數(shù)據(jù)為1024點(diǎn)的處理
fs=100;N=1024;n=0:N-1;t=n/fs;
x=0.5*sin（2*pi*15*t）+2*sin（2*pi*40*t）;?%信號
y=fft（xN）;???%對信號進(jìn)行快速Fourier變換
mag=abs（y）;???%求取Fourier變換的振幅
f=n*fs/N;
subplot（223）plot（fmag）;?%繪出隨頻率變化的振幅
xlabel（‘頻率/Hz‘）;
ylabel（‘振幅‘）;title（‘N=1024‘）;grid?on;
subplot（224）
plot（f（1:N/2）mag（1:N/2））;?%繪出Nyquist頻率之前隨頻率變化的振幅
xlabel（‘頻率/Hz‘）;
ylabel（‘振幅‘）;title（‘N=1024‘）;grid?on;






clf;fs=100;?%采樣頻率
Ndata=32;?%數(shù)據(jù)長度
N=32;?%FFT的數(shù)據(jù)長度
n=0:Ndata-1;t=n/fs;???%數(shù)據(jù)對應(yīng)的時間序列
x=0.5*sin（2*pi*15*t）+2*sin（2*pi*40*t）;???%時間域信號
y=fft（xN）;???%信號的Fourier變換
mag=abs（y）;????%求取振幅
f=（0:N-1）*fs/N;?%真實(shí)頻率
subplot（221）plot（f（1:N/2）mag（1:N/2）*2/N）;?%繪出Nyquist頻率之前的振幅
xlabel（‘頻率/Hz‘）;ylabel（‘振幅‘）;
title（‘Ndata=32?Nfft=32‘）;grid?on;
Ndata=32;???%數(shù)據(jù)個數(shù)
N=128;?????%FFT采用的數(shù)據(jù)長度
n=0:Ndata-1;t=n/fs;???%時間序列
x=0.5*sin（2*pi*15*t）+2*sin（2*pi*40*t）;
y=fft（xN）;
mag=abs（y）;
f=（0:N-1）*fs/N;?%真實(shí)頻率
subplot（222）plot（f（1:N/2）mag（1:N/2）*2/N）;?%繪出Nyquist頻率之前的振幅
xlabel（‘頻率/Hz‘）;ylabel（‘振幅‘）;
title（‘Ndata=32?Nfft=128‘）;grid?on;
Ndata=136;???%數(shù)據(jù)個數(shù)
N=128;?????%FFT采用的數(shù)據(jù)個數(shù)
n=0:Ndata-1;t=n/fs;?%時間序列
x=0.5*sin（2*pi*15*t）+2*sin（2*pi*40*t）;
y=fft（xN）;
mag=abs（y）;
f=（0:N-1）*fs/N;???%真實(shí)頻率
subplot（223）plot（f（1:N/2）mag（1:N/2）*2/N）;?%繪出Nyquist頻率之前的振幅
xlabel（‘頻率/Hz‘）;ylabel（‘振幅‘）;
title（‘Ndata=136?Nfft=128‘）;grid?on;
Ndata=136;????%數(shù)據(jù)個數(shù)
N=512;????%FFT所用的數(shù)據(jù)個數(shù)
n=0:Ndata-1;t=n/fs;?%時間序列
x=0.5*sin（2*pi*15*t）+2*sin（2*pi*40*t）;
y=fft（xN）;
mag=abs（y）;
f=（0:N-1）*fs/N;???%真實(shí)頻率
subplot（224）plot（f（1:N/2）mag（1:N/2）*2/N）;?%繪出Nyquist頻率之前的振幅
xlabel（‘頻率/Hz‘）;ylabel（‘振幅‘）;
title（‘Ndata=136?Nfft=512‘）;grid?on;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2322??2014-01-09?21:10??FFT.m

-----------?---------??----------?-----??----

?????????????????2322????????????????????1