声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2348|回复: 8

[编程技巧] HHT编程检测谐波,各位大虾进来看看,帮帮忙啊

[复制链接]
发表于 2011-10-25 21:41 | 显示全部楼层 |阅读模式

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

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

x
波形如下所示

fn=50;
x1=sin(2*pi*fn*t);
x2=p.*[0.5*sin(2*pi*3*fn*t)+0.3*sin(2*pi*7*fn*t)];
x=x1+x2;
想对imf1求瞬时频率,出来的结果不对
麻烦各位高手帮忙想想是哪里出错了


U(V46IYR)D7SL590Z0IY549.jpg
回复
分享到:

使用道具 举报

发表于 2011-10-25 23:45 | 显示全部楼层
1.好像不是编程问题吧!
2.好像资讯也没给齐
发表于 2011-10-27 17:17 | 显示全部楼层
请把问题说具体点,最好有程序和数据,才好帮你啊
 楼主| 发表于 2011-10-28 11:10 | 显示全部楼层
回复 3 # 月生才 的帖子

我的问题就是出来的波形怎么不是平稳的……是震荡的曲线,是不是和断点效应问题有关系呢
发表于 2011-10-28 20:05 | 显示全部楼层
回复 4 # o蛋塔o 的帖子

你把完整的程序贴出来,大家分析下嘛。还有你那个图是什么图呢?imf1的瞬时频率?
 楼主| 发表于 2011-10-29 09:07 | 显示全部楼层
回复 5 # 月生才 的帖子

clc;
clear;
t=0:1/3200:0.4;
p=zeros(1,length(t));
aa=find(t>=0.12&t<=0.24);
p(aa)=1;
fn=50;
x1=sin(2*pi*fn*t);
x2=p.*[0.5*sin(2*pi*3*fn*t)+0.3*sin(2*pi*7*fn*t)];
x=x1+x2;
plot(t,x);
imf=emd(x);
figure(2);
cemd_visu(x,1:length(x),imf);
麻烦你了……
发表于 2011-10-29 09:25 | 显示全部楼层
回复 6 # o蛋塔o 的帖子

lz帖子上的那张图是怎么画的啊?代码没有啊!
发表于 2011-10-29 14:27 | 显示全部楼层
你把完整的程序贴出来,大家分析下嘛。还有你那个图是什么图呢?imf1的瞬时频率?
发表于 2011-12-1 18:04 | 显示全部楼层
求瞬时频率可以用F. Auger的instfreq.m函数求,cemd_visu.m的输出是画出原输入信号、分解后得到的imfs和余量,并且相位可调。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-5 16:30 , Processed in 0.075557 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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