声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 14342|回复: 47

[综合] 计算阶比跟踪与谐波成分提取整理贴

  [复制链接]
发表于 2013-1-11 00:36 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 impulse 于 2013-1-11 11:26 编辑

模拟信号由matlab模拟产生,为一个频率随时间线性增加的信号,基频频率从5Hz变化到30Hz,同时含有0.5倍频和2倍频。
代码如下:


  1. %模拟产生变转速振动信号和键相信号
  2. clc;
  3. clear;
  4. Ls=1024*64*4;   %样本数
  5. fs=4096*8;      %采样频率
  6. t0 =0:1/fs:(Ls-1)/fs;
  7. t1=t0(length(t0));
  8. f0=5;     %起始频率
  9. f1=30;    %终止频率
  10. x_1x = 1*chirp(t0,f0,t1,f1,[],280);             %一倍频信号
  11. x_5x = 0.65*chirp(t0,0.5*f0,t1,0.5*f1,[],115);  %0.5倍频信号
  12. x_2x = 0.75*chirp(t0,2*f0,t1,2*f1,[],20);       %2倍频信号
  13. x_noise=0.15*rand(1,length(x_1x));              %噪声
  14. x=x_1x+x_2x+x_5x;%+x_noise;                     %合成信号

  15. [maxtab, mintab] = peakdet(x_1x, 0.5);  %通过基频信号(相当于转轴信号)检波,模拟后面的键相信号,函数附在附件中
  16. t=maxtab(:,1)';
  17. tacho=zeros(1,Ls);
  18. tacho(t)=1;                             %模拟键相信号,假定1为键相脉冲位置
  19. tacho(t+1)=1;
  20. subplot(211)
  21. plot(t0,x);
  22. title('振动信号')
  23. subplot(212)
  24. plot(t0,tacho)
  25. title('转速(键相)信号')
复制代码
1、原始信号如下
rawwave.PNG
2、阶比跟踪重采样后的信号
evenwave.PNG
3、原始信号频谱与重采样后信号阶比谱比较
spect.PNG
4、Gabor变换提取谐波波形
gabor.PNG
5、恒带宽零相位滤波器提取谐波波形
cbf.PNG
6、Vold-Kalman滤波提取谐波波形
vkf.PNG
7、上述方法比较,VKF方法堪称完美。
peakdet.rar (844 Bytes, 下载次数: 269)

点评

赞成: 5.0
赞成: 5
perfect!  发表于 2013-1-11 00:44

评分

2

查看全部评分

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2013-1-16 11:39 | 显示全部楼层
本帖最后由 ufo158 于 2013-1-16 11:43 编辑

顶贴...VKF在计算速度、在临近 交叠阶次的提取上, 效果也很好.
端点上提取结果, 从仿真上看, 优于其他两种, 并且相位想没有偏移。lz VKF带宽怎么取的??常带宽还是百分比?
计算速度上,若用PCG解方程,要慢一些的哈。

发表于 2013-1-19 00:51 | 显示全部楼层
小生问题总是很多:
在文献里面看到有一个叫“Constant Percentage Bandwidth”,麻烦解释一下,谢谢。
 楼主| 发表于 2013-1-19 19:37 | 显示全部楼层

恒百分比带宽,滤波器的带宽与中心频率百分比一定。
发表于 2013-1-22 12:02 | 显示全部楼层
impulse 发表于 2013-1-19 19:37
恒百分比带宽,滤波器的带宽与中心频率百分比一定。

非常谢谢你的解释,但是我还是不懂;
一个老师告诉我们滤波器的截止频率就是分析分析频率的最大值。"滤波器的带宽"是不是就是0~fs/2呢,现在的滤波器主要都是数字滤波器了吧?这个带宽就是由那奎斯特频率决定的吧,那中心频率又是啥子呢?

最好提供相关书籍,我自己看去。
发表于 2013-1-24 11:44 | 显示全部楼层
mxlzhenzhu 发表于 2013-1-22 12:02
非常谢谢你的解释,但是我还是不懂;
一个老师告诉我们滤波器的截止频率就是分析分析频率的最大值。"滤波 ...

这是个带通滤波器, 中心频率也就是需要滤出的那部分分量的频率, 滤波带宽即在此中心频率附近也将被滤出的分量的频率范围.
发表于 2013-2-21 15:07 | 显示全部楼层
thank you
发表于 2013-5-4 06:28 | 显示全部楼层
LZ你好,小弟最近在也在学习cot,但是对怎样实现不是很了解。能不能请你分享一下MATLAB程序,让我学习一下阶比跟踪重采样的方法
发表于 2013-5-6 15:18 | 显示全部楼层
我也想要matlab程序,学习一下阶比跟踪重采样的方法邮箱wutaolww@163.com
 楼主| 发表于 2013-5-6 15:33 | 显示全部楼层
我没有COT的Matlab程序,我的是自己开发的C++商业版本,matlab版本论坛上有一个,虽然效果不是很好,但是基本过程都有
发表于 2013-5-23 11:15 | 显示全部楼层
还需要不断学习啊~~~
发表于 2013-7-19 20:23 | 显示全部楼层
楼主V5,跟楼主学习太多知识
发表于 2013-10-3 11:15 | 显示全部楼层
精华帖。。。学习了
发表于 2013-10-22 19:23 | 显示全部楼层
我又来了,我想问下Vold-Kalman filtering的vkf(y,fs,f,p,bw,multiorder)f bw参数应该怎么取?你能给我推荐些这方面的文章吗?谢谢了。
你的这个数据我用VKF是没有提取到你文中的那样效果,而且差距很远.
提取一阶的参数
p = 1;
bw = 8;
f0=5;     %起始频率
f1=30;    %终止频率
inc=(f1-f0)/Ls;
fdn=[f0+inc:inc:f1];
vkf(y,fs,fdn,p,bw);
能给指点下不?
发表于 2014-4-9 15:56 | 显示全部楼层
主任可否谈一下:阶比分析及阶比跟踪的具体意义?本人很迷惑,3Q
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-3-29 20:32 , Processed in 0.118236 second(s), 25 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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