声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1371|回复: 2

[HHT] 请教高手啊!主程序改怎么写???

[复制链接]
发表于 2009-5-18 23:06 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
谐波信号的 HHT 分析
   设间谐波信号为
                     2
            s(t) =  ∑ Ai sin(2πmi ft)           (13)
                    i=1
式中 A1=1V,A2=0.3V,m1=1,m2=5.5,f=50Hz。
采样频率为3200Hz。
要对这个信号进行HHT分析,我用的是下面三个自定义程序,请教高手给怎么写主程序调用呢??

plot_hht.m

965 Bytes, 下载次数: 18

emd.m

893 Bytes, 下载次数: 17

findpeaks.m

146 Bytes, 下载次数: 16

回复
分享到:

使用道具 举报

 楼主| 发表于 2009-5-18 23:09 | 显示全部楼层

回复 楼主 sxd1168 的帖子

我用的上面几个附件里有人说还是有小错误,希望高手指点一下,我在写毕业论文,好恼火啊。。。
发表于 2009-5-28 19:01 | 显示全部楼层
这是我在论坛里找的,我现在也用的是这个,你用试试吧,里面的几个程序要下工具箱才能用
t=1:500;
t=t*1/400;
x1=sin(2*pi*20*t);
x2=0.4*sin(2*pi*40*t+140);
f=400;
y=linspace(1,2,500);
z=x1+x2+y
%在后面添加一个趋势向可以把三个分量完全分解出来
imf=emd(z);
emd_visu(z,1:length(z),imf,1)
figure (4)
dt=1/f;
h = nspab(imf',500,0,500,dt);
h(1:2,:)=0;


ms=mspc(h,t)

subplot(224)
plot((1:length(ms)),ms)
xlabel('频率(Hz)','FontSize',12);
ylabel('幅值','FontSize',12);
subplot(211)
surf(h(1:200,1:end-20))
shading interp
xlabel('时间(点数)','FontSize',12);
Ylabel('频率(Hz)','FontSize',12);
zlabel('幅值','FontSize',12);
title('a','FontSize',12)


view([-75,25])
yt=subplot(223)
imagesc(h(1:200,:))
Ylabel('频率(Hz)','FontSize',12);
xlabel('时间(点数)','FontSize',12);
set(yt,'ydir','nor')
title('b','FontSize',12);title('c','FontSize',12);

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 20:55 , Processed in 0.083134 second(s), 30 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表