|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我用G. Rilling 2007年写的版本求出爆破波形的HHT谱图,然后用论坛上破凰写的求边际谱的程序求爆破波形的边际谱图,但是得出的结果是消减的太快,不知是为什么呢,哪位高手可以进来帮忙,帮看看呀!
这是我的程序:
sig=rsac(a);
t=sig(:,1); %信号的时域采样点 =36001
pm=sig(:,2); %幅值
N=length(t);
fs=50; %采样频率
T=N*1/fs;
[imf,ort,nbits] = emd(pm);
emd_visu(pm,t,imf);
[A,f,tt] = hhspectrum(imf(1:end-1,:));
[im,tt,ff] = toimage(A,f,tt);
disp_hhs(im,[],fs);
figure;
for k=1:size(im,1); %size(im)=[400 35999]
bjp1(k)=sum(im(k,:))*1/fs;
end
NN=size(im,1);
f=(0:NN-1)/NN*(fs/3);
plot(f,bjp1,'r');
xlabel('频率HZ');
ylabel('幅值'); |
|