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;
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

夜神明月

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值