边际谱是否正确呢?
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; %采样点数
=hhspectrum(z);
=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);
上面的程序得到的边际谱是否正确呢,我综合论坛的程序做的。谢谢! 做出结果来了吗 把效果图贴出来
页:
[1]