大家帮忙看看边际谱的幅值怎么不对啊?
下面我是对信号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);
= hhspectrum(imf(1:end-1,:));
= 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('幅值');
回复 楼主 的帖子
不好意思,图片传不上去啊!得到10Hz的信号幅值为3,而25Hz的信号的幅值只有1.8.回复 2楼 的帖子
衰减问题,以前的帖子讨论过了,看一下置顶帖关于边际谱的讨论,注意整周期采样问题
页:
[1]