在MATLAB中,可以使用fft函數來實現傅里葉變換。下面是一個示例代碼:
% 生成一個信號
Fs = 1000; % 采樣頻率
t = 0:1/Fs:1; % 時間序列
f1 = 50; % 信號頻率1
f2 = 120; % 信號頻率2
A1 = 1; % 信號幅值1
A2 = 0.5; % 信號幅值2
x = A1*sin(2*pi*f1*t) + A2*sin(2*pi*f2*t);
% 進行傅里葉變換
N = length(x); % 信號長度
X = fft(x); % 傅里葉變換結果
f = (0:N-1)*(Fs/N); % 頻率坐標
power = abs(X).^2/N; % 頻譜功率
% 繪制頻譜圖
plot(f,power)
xlabel('Frequency (Hz)')
ylabel('Power')
該代碼首先生成了一個由兩個正弦信號疊加而成的信號。然后使用fft函數對信號進行傅里葉變換,得到頻譜信息。最后,使用plot函數繪制出頻譜圖。