在MATLAB中,可以使用fft函數(shù)來實(shí)現(xiàn)傅里葉變換。下面是一個(gè)示例代碼:
% 生成一個(gè)信號
Fs = 1000; % 采樣頻率
t = 0:1/Fs:1; % 時(shí)間序列
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);
% 進(jìn)行傅里葉變換
N = length(x); % 信號長度
X = fft(x); % 傅里葉變換結(jié)果
f = (0:N-1)*(Fs/N); % 頻率坐標(biāo)
power = abs(X).^2/N; % 頻譜功率
% 繪制頻譜圖
plot(f,power)
xlabel('Frequency (Hz)')
ylabel('Power')
該代碼首先生成了一個(gè)由兩個(gè)正弦信號疊加而成的信號。然后使用fft函數(shù)對信號進(jìn)行傅里葉變換,得到頻譜信息。最后,使用plot函數(shù)繪制出頻譜圖。