Matlab進行FFT變換

1、代碼

close all;          %先關閉所有圖片
Adc=2;              %直流分量幅度
A1=3;               %頻率F1信號的幅度
A2=1.5;             %頻率F2信號的幅度
F1=50;              %信號1頻率(Hz)
F2=75;              %信號2頻率(Hz)
Fs=256;             %采樣頻率(Hz)，根據奈奎斯特定理：采樣頻率必須大於2倍信號頻率
P1=-30;             %信號1相位(度)
P2=90;              %信號相位(度)
N=256;              %采樣點數，為了方便進行FFT運算，通常N取2的整數次方
t=[0:1/Fs:N/Fs];    %采樣時刻


S=Adc+A1*cos(2*pi*F1*t+pi*P1/180)+A2*cos(2*pi*F2*t+pi*P2/180); %信號
plot(S);            %顯示原始信號
title('原始信號');  

figure;
Y = fft(S,N);       %做FFT變換，結果為N點的復數，每一個點就對應着一個頻率點
Ayy = (abs(Y));     %取模運算，對數值是求絕對值，對復數是求幅值。這里就是復數轉幅度值
plot(Ayy(1:N));     %顯示原始的FFT模值結果
title('FFT 模值');

figure;
Ayy=Ayy/(N/2);      %換算成實際的幅度，針對半頻譜
Ayy(1)=Ayy(1)/2;    %直流和奈奎斯特頻率處還需要/2
F=([1:N]-1)*Fs/N;   %換算成實際的頻率值，采到每個點對應的頻率
plot(F(1:N/2),Ayy(1:N/2));    %顯示換算后的FFT模值結果
title('幅度-頻率曲線圖');

figure;
Pyy=[1:N/2];
for i=1:N/2
Pyy(i)=phase(Y(i));             %計算相位
Pyy(i)=Pyy(i)*180/pi;           %換算為角度
end;
plot(F(1:N/2),Pyy(1:N/2));      %顯示相位圖
title('相位-頻率曲線圖');

2、結果