对时间序列数据做频域分析的matlab代码

以下是一个简单的 MATLAB 代码示例，用于对时间序列数据进行频域分析：

% 读取时间序列数据
data = load('timeseries.txt');

% 计算采样频率和采样点数
fs = 100; % 假设采样频率为100Hz
n = length(data);

% 对数据进行傅里叶变换
Y = fft(data);

% 计算频率轴上的频率值
f = (0:n-1)*(fs/n);

% 计算幅值谱和相位谱
amplitude = abs(Y/n);
phase = unwrap(angle(Y))/n;

% 绘制幅值谱和相位谱
subplot(2,1,1);
plot(f,amplitude);
xlabel('Frequency (Hz)');
ylabel('Amplitude');
title('Amplitude Spectrum');
subplot(2,1,2);
plot(f,phase);
xlabel('Frequency (Hz)');
ylabel('Phase');
title('Phase Spectrum');

这个代码示例假设时间序列数据已经存储在一个名为 `timeseries.txt` 的文本文件中。在实际应用中，你需要根据自己的数据格式进行相应的修改。该代码首先计算采样频率和采样点数，然后对数据进行傅里叶变换。接下来，代码计算频率轴上的频率值，并计算幅值谱和相位谱。最后，代码绘制了幅值谱和相位谱两张图像，分别显示频率和幅值之间的关系以及频率和相位之间的关系。