0208106224 发表于 2012-5-16 13:10

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

下面这段是我弄的一个EMD去噪,但不知道为什么在matlab里老是一直处于死循环状态!哪位大侠帮帮忙看看
function =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,:))
页: [1]
查看完整版本: 一个肌电信号如何通过EMD去噪呢?