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

impulse

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

1. 
   
2. 
   
3. %模拟产生变转速振动信号和键相信号
   
4. clc;
   
5. clear;
   
6. Ls=1024*64*4;   %样本数
   
7. fs=4096*8;      %采样频率
   
8. t0 =0:1/fs:(Ls-1)/fs;
   
9. t1=t0(length(t0));
   
10. f0=5;     %起始频率
    
11. f1=30;    %终止频率
    
12. x_1x = 1*chirp(t0,f0,t1,f1,[],280);             %一倍频信号
    
13. x_5x = 0.65*chirp(t0,0.5*f0,t1,0.5*f1,[],115);  %0.5倍频信号
    
14. x_2x = 0.75*chirp(t0,2*f0,t1,2*f1,[],20);       %2倍频信号
    
15. x_noise=0.15*rand(1,length(x_1x));              %噪声
    
16. x=x_1x+x_2x+x_5x;%+x_noise;                     %合成信号
    
17. 
    
18. [maxtab, mintab] = peakdet(x_1x, 0.5);  %通过基频信号（相当于转轴信号）检波，模拟后面的键相信号，函数附在附件中
    
19. t=maxtab(:,1)';
    
20. tacho=zeros(1,Ls);
    
21. tacho(t)=1;                             %模拟键相信号，假定1为键相脉冲位置
    
22. tacho(t+1)=1;
    
23. subplot(211)
    
24. plot(t0,x);
    
25. title('振动信号')
    
26. subplot(212)
    
27. plot(t0,tacho)
    
28. title('转速（键相）信号')
    

复制代码

1、原始信号如下

2、阶比跟踪重采样后的信号

3、原始信号频谱与重采样后信号阶比谱比较

4、Gabor变换提取谐波波形

5、恒带宽零相位滤波器提取谐波波形

6、Vold-Kalman滤波提取谐波波形

7、上述方法比较，VKF方法堪称完美。
peakdet.rar (844 Bytes, 下载次数: 272)

2013-1-11 11:25 上传

点击文件名下载附件

ufo158

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

mxlzhenzhu

小生问题总是很多：
在文献里面看到有一个叫“Constant Percentage Bandwidth”，麻烦解释一下，谢谢。

impulse

mxlzhenzhu 发表于 2013-1-19 00:51

                               
登录/注册后可看大图

小生问题总是很多：
在文献里面看到有一个叫“Constant Percentage Bandwidth”，麻烦解释一下，谢谢。

恒百分比带宽，滤波器的带宽与中心频率百分比一定。

mxlzhenzhu

impulse 发表于 2013-1-19 19:37

                               
登录/注册后可看大图

恒百分比带宽，滤波器的带宽与中心频率百分比一定。

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

最好提供相关书籍，我自己看去。

ufo158

mxlzhenzhu 发表于 2013-1-22 12:02

                               
登录/注册后可看大图

非常谢谢你的解释，但是我还是不懂；
一个老师告诉我们滤波器的截止频率就是分析分析频率的最大值。"滤波 ...

这是个带通滤波器, 中心频率也就是需要滤出的那部分分量的频率, 滤波带宽即在此中心频率附近也将被滤出的分量的频率范围.

zfch

thank you

solamy

LZ你好，小弟最近在也在学习cot，但是对怎样实现不是很了解。能不能请你分享一下MATLAB程序，让我学习一下阶比跟踪重采样的方法

吴涛

我也想要matlab程序，学习一下阶比跟踪重采样的方法邮箱wutaolww@163.com

impulse

我没有COT的Matlab程序，我的是自己开发的C++商业版本，matlab版本论坛上有一个，虽然效果不是很好，但是基本过程都有

ylhlly

还需要不断学习啊~~~

和丹

楼主V5，跟楼主学习太多知识

gxgphjl

精华帖。。。学习了

和丹

我又来了，我想问下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);
能给指点下不？

猫头鹰先生

主任可否谈一下：阶比分析及阶比跟踪的具体意义？本人很迷惑，3Q