在MATLAB中進行傅里葉變換的一般步驟如下:
下面是一個簡單的示例代碼來說明這些步驟:
% 信號參數
Fs = 1000; % 采樣頻率
T = 1/Fs; % 采樣周期
L = 1000; % 采樣點數
% 創建時間向量和信號向量
t = (0:L-1)*T; % 時間向量
x = sin(2*pi*50*t) + sin(2*pi*120*t); % 信號向量
% 進行傅里葉變換
X = fft(x); % 傅里葉變換結果
P2 = abs(X/L); % 雙邊頻譜
P1 = P2(1:L/2+1); % 單邊頻譜
P1(2:end-1) = 2*P1(2:end-1); % 幅度修正
% 創建頻率向量
f = Fs*(0:(L/2))/L;
% 繪制時域波形和頻域波形
subplot(2,1,1);
plot(t, x);
title('時域波形');
xlabel('時間 (s)');
ylabel('幅度');
subplot(2,1,2);
plot(f, P1);
title('頻域波形');
xlabel('頻率 (Hz)');
ylabel('幅度');
運行以上代碼,將會繪制出信號的時域波形和頻域波形。