|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
感觉用时频工具箱中的瞬时频率函数instfreq求得的瞬时频率跟直观不太一致。
程序代码如下:
clear all;
clc
load E:\recent_read\matwork\modelsemg\mu3.mat
sigbip = sigbip(1001:2000);
sigbip = sigbip/max(sigbip);
[imf,ort,nbits] = emd(sigbip); % <-------- adapted to current emd.m
%emd_visu(sigbip,1:2000,imf,1);
%画第一个imf
figure
imf1_an = hilbert(imf(1,:));
instfre_imf1 = instfreq(imf1_an');
subplot(3,1,1);
plot(1:1000,imf(1,:));
title('imf1');
subplot(3,1,2);
plot([1:length(instfre_imf1)],instfre_imf1);
title('imf1 inst frequency');
subplot(3,1,3);
plot([1:1000],sigbip)
title('signal');
[ 本帖最后由 dailiangren 于 2008-3-27 04:01 编辑 ] |
|