马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clear;clc;close all
fs=1000; %fs为采样频率;
N=2000; %采样点数
tspan=2;
t=1/fs:1/fs:tspan;
y1=2*sin(60*pi*t);
y2=5*sin(90*pi*t);
y=[y1;y2;zeros(size(y1))]; %IMF集
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);
[E,tt1]=toimage(A,fa,tt,length(tt));
for k=1:size(E,1)
bjp(k)=sum(E(k,:).*E(k,:))*1/fs*1/tspan;
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
y=y1+y2;
Nsec=512;
window1=hanning(256);
wm1=xcorr(window1,window1); wn1=sum(wm1)/(max(wm1)*511);
[Pyy,f]=pwelch(y,window1,128,Nsec,fs);
figure(2)
plot(f,Pyy*fs/256/wn1,'r');
这样求出的hilbert能量谱能和功率谱做比较么?这样能证明hilbert能量谱正确么?请各位指点!我的目的就是想精确求出能量或功率谱。 |