傅里叶变换频谱图原始matlab代码
时间: 2023-05-11 20:01:27 浏览: 168
傅里叶变换是一种非常重要的信号分析工具,可以将一个信号从时域(即时间轴)转换到频域(即频率轴),使我们可以了解信号中各种频率的成分和它们的相对大小。
在Matlab中,可以使用“fft”函数进行快速傅里叶变换,该函数返回的是信号在频域中的幅度与相位信息,可以通过取模操作得到信号的正半频域幅度谱,进而展现信号频域的特征之一 —— 频谱(即幅度谱)。
以下是基于Matlab语言的原始代码,展示如何通过fft函数获得原始信号的频谱。
```matlab
% 定义信号并将其傅里叶变换
t = linspace(0, 2*pi, 5000); % 时间轴
fs = 1/(t(2)-t(1)); % 采样率
x = 0.5*sin(50*t) + 0.7*cos(80*t); % 原始信号
fft_x = fft(x); % 傅里叶变换
% 计算频谱并作图
L = length(x); % 信号的长度
P2 = abs(fft_x/L); % 正半频域幅度谱
P1 = P2(1:L/2+1); % 单侧频谱
P1(2:end-1) = 2*P1(2:end-1); % 计算双侧频谱
f = fs*(0:(L/2))/L; % 频率轴
plot(f,P1) % 绘制频谱
title('频谱图')
xlabel('频率 (赫兹)')
ylabel('幅度')
```
这段代码首先生成一个由sine和cosine函数组合而成的信号,使用f = fs \* (0:(L/2))/L计算频率轴,然后使用fft函数计算这个信号的傅里叶变换,最后使用P2 = abs(fft_x/L)计算信号的正半频域幅度谱,并绘制出在频域各频率点对应的频谱幅值,即可得到原始信号的频谱图。
阅读全文