一个肌电信号如何通过EMD去噪呢？

0208106224

下面这段是我弄的一个EMD去噪，但不知道为什么在matlab里老是一直处于死循环状态！哪位大侠帮帮忙看看
function [mmg1 mmg2]=EmgSmgopen
fid1=fopen('yj01Emg1.emg','r');
mmg1=fread(fid1,'float');
fid2=fopen('yj01Emg2.emg','r');
mmg2=fread(fid2,'float');

fclose(fid1);
fclose(fid2);

%title('肌电信号');

%title('肌音信号');
for i=1:(length(mmg1))/2
    mmgs(i)=mmg1(2*i-1);
    mmgd(i)=mmg1(2*i);
end
%emgTrial=emgTrial(1:51200);
x=mmgs;
figure;
plot(x);
imf=emd_n(x,5);
figure;
subplot(5,1,1);
plot(imf(1,:))