Matlab实现对信号的傅里叶变换，并画出单边频谱和双边频谱

clear;
clc;
load('CH1.mat'); %加载信号的幅值数据
load('CH2.mat');
% 采样率Fs
Fs = 125000;          % 采样率，单位Hz
% 先对CH1做数据清洗
CH1 = removeOutliers(CH1);
% 获取信号长度
num = length(CH1);
disp(num);
% 生成时间向量
t = linspace(0,1,num);       % 时间向量，从0到1秒，步长为1/Fs
% FFT变换
Y = fftshift(fft(CH1));   
% 计算频率向量
f = Fs*(0:(num/2))/num;

% 绘制时域图像
subplot(1,3,1);
plot(t,CH1);
title("时域图像");
xlabel('时间t/s');
ylabel('幅度');
grid on;
% 绘制双边频谱
subplot(1,3,2);
plot((-num/2:num/2-1)*(Fs/num), abs(Y));
title('双边频谱');
xlabel('频率 (Hz)');
ylabel('幅度');
xlim([-Fs/2 Fs/2]); % 设置x轴范围
grid on;
% 绘制单边频谱
subplot(1,3,3);
plot(f, 2*abs(Y(1:num/2+1))/num); % 单边频谱幅度乘以2（除了直流分量）
title('单边频谱');
xlabel('频率 (Hz)');
ylabel('幅度');
xlim([0 Fs/2]); % 设置x轴范围
grid on;