马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clear;
t=1:500;
t=t*1/2000;
x=sin(2*pi*20*t);
y=0.4*sin(2*pi*40*t+140);
z=x+y;
tspan=0.4;
plot(t,z)
imf=emd(z);
emd_visu(z,1:length(z),imf)
%%%%%%%%%%%%%求边际谱
figure;
fs=1000; %fs为采样频率;
N=1000; %采样点数
[A,fa,tt]=hhspectrum(z);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E);
for k=1:size(E,1)
bjp(k)=sum(E(k,:))*1/fs*1/tspan;
end
H=size(E,1);
f=(0:H-1)/H*(fs/2);
plot(f,bjp);
上面的程序得到的边际谱是否正确呢,我综合论坛的程序做的。谢谢! |