test2.m

来自「新的信号处理方法－希尔伯特黄变换的EMD分解MATLAB程序 把程序看懂 对希尔」· M 代码 · 共 27 行

clear
clc
load spo2_0p2.txt;%spo2_3p.txt;%spo2_0p2.txt ;% spo2_3p.txt;spo2_0p4.txt spo2_3p.txt
result = spo2_0p2;
len=50;
t = 0:0.001:0.5;

x = 1:501;

%第二，第三个字节是红光脉搏波数据
rdxn = result([(x-1)*5+4])*256 +result([(x-1)*5+5]);
stop = [0.1,0.1,0.1];
[imf,ort,nbits]=myemd(rdxn,t,stop,2);

imfsize=size(imf);
i=imfsize(1);
figure;
for j=1:i
    subplot(i/2+1,2,j);plot(imf(j,:));
end

figure;
newrdxn=rdxn;
for k=1:i
    newrdxn=newrdxn-imf(k,:);
    subplot(i/2+1,2,k);plot(newrdxn);
end