在分析信号时,可以先进行傅里叶变换(FFT)分析,在频率方面进行观察
在FFT_Function.m里进行函数的定义
function [FFTCoeff ,f] = FFT_Function(data,fs)
L = length(data); % 信号长度
Y = fft(data);
P2 = abs(Y/L);
FFTCoeff = P2(1:floor(L/2)+1);%floor是取比它小的整数
f = fs*(0:floor(L/2))/L;
在FFT_Analysis.m里进行信号的FFT分析
clc;clear all;
load('N118.mat');
data = a(3,1:1440);
fs = 360;
figure(1);
plot(data);
xlabel('点数');ylabel('幅值');
[FFTCoeff ,f] = FFT_Function(data,fs);%调用函数分析
figure(2);
plot(f,FFTCoeff);
title('心电信号的单边幅度谱');
xlabel('频率 (Hz)');
ylabel('幅值')
结果如下:
嘻嘻!