|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
下面我是对信号x=5*sin(2*pi*25*t)+3*sin(2*pi*10*t)+exp(2*t)进行emd分解后求边际谱,但边际谱的幅值怎么不对?请高手指导!
程序:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc
clear;
fs=100; %fs为采样频率;
tspan=2;
t=1/fs:1/fs:tspan;
N=length(t); %采样点数
y1=5*sin(2*pi*25*t);
y2=3*sin(2*pi*10*t);
y3=exp(2*t);
x=y1+y2+y3;
imf=emd(x);
[A,f,tt] = hhspectrum(imf(1:end-1,:));
[im,tt1,ff] = toimage(A,f,tt);
figure
for k=1:size(im,1)
bjp(k)=sum(im(k,:))*1/fs;
end
ff=ff*fs;
plot(ff,bjp);
xlabel('频率 / Hz');
ylabel('幅值'); |
|