循环统计量，高阶统计量

fearless

有人做循环平稳分析的么，涉及到循环双谱等知识的，可以一起来交流下么？

amnesia

循环双谱matlab程序

1. function cyclic_bispectrum(x,fs,M1,M2)
   
2. %%% 循环双谱程序
   
3. %%% x 输入的信号
   
4. 
   
5. %%%=====================================================
   
6. %%% 第一步：求信号的三次变换,依据为： y(i,n)=x(i)x(i+n)x(i+n) 后两项n取一样，为方便运算。
   
7. %%%
   
8. M = length(x); % 信号采样点数
   
9. if M>20000
   
10.     M = 20000;
    
11. end
    
12. 
    
13. % 初始化
    
14. N1 = fix(3*M/4); N2 = M - N1; % 必须满足 M=N1+N2，且N1>N2。
    
15. y = zeros(N1,N2);
    
16. for i=1:N1
    
17.     for n=1:N2
    
18.         y(i,n) = x(i)*x(i+n)*x(i+n);
    
19.     end
    
20. end
    
21. 
    
22. %%%=====================================================
    
23. %%% 第二步：求信号的三次变换,依据为： Rtk(n)=1/N1*sum(y(i,n)*exp(-1i*2*pi*i*k/N1))
    
24. %%% M1对应 2pi 弧度，分辨率 2pi/M1
    
25. Rtk = zeros(N2,M1);
    
26. for n = 1:N2
    
27. %     Rtk(n,:) = 1/N1*fftshift(fft(y(:,n),M1));
    
28.     Rtk(n,:) = fft(y(:,n),M1);
    
29. end
    
30. 
    
31. %%%=====================================================
    
32. %%% 第三步：求循环双谱切片谱,依据为： Imk(m)=sum(Rtk(1:N2).*exp(-1i*2*pi*m*(1:N2)/N2))
    
33. %%%
    
34. Imk = zeros(M1,M2);
    
35. for m = 1:M1
    
36. %     Imk(m,:) = fftshift(fft(Rtk(:,m),M2));
    
37.     Imk(m,:) = fft(Rtk(:,m),M2);
    
38. end
    
39. aaa = Imk(1:M1/2,1:M2/2);
    
40. 
    
41. x_lab = fs/2*linspace(0,1,M1/2);
    
42. y_lab = fs/2*linspace(0,1,M2/2);
    
43. 
    
44. % I_abs = abs(Imk);
    
45. I_abs = abs(aaa);
    
46. %%% 画图
    
47. figure
    
48. mesh(y_lab,x_lab,I_abs)
    
49. % mesh(I_abs)
    
50. 
    
51. % figure
    
52. % plot(x_lab,I_abs(:,1))
    
53. %
    
54. figure
    
55. plot(x_lab,I_abs)
    
56. %
    
57. % figure
    
58. % plot(x_lab,I_abs(:,round(M1/2)-1))

复制代码

amnesia

循环双谱及在周期平稳类故障中的应用

作者:苏中元； 贾民平； 许飞云； 胡建中；

摘要:论述了未知循环频率的周期平稳信号循环双谱的估计方法;研究了在循环双谱的循环累积量计算中涉及变量的简化存储方法,提出了该变量矩阵是对称阵,通过算法可以三角阵的元素来表述,克服了循环双谱传统估计方法计算量较大的缺陷,提高了运算效率;提出了循环双谱对调相故障信号的分析能力,以及对加性噪声的处理能力,仿真并验证了该方法的有效性,并将其应用于旋转机械状态分析。

fearless

amnesia 发表于 2016-9-5 13:35
循环双谱及在周期平稳类故障中的应用

作者:苏中元； 贾民平； 许飞云； 胡建中；

哇 多谢兄台。咱们可以一起交流这个。

Pparis

楼主研究这个做什么

sd

《基于高阶循环统计量的频谱感知方法》

zswseu

高阶统计量3,4阶以上呢，是不是阶数越高越好呢，以后有时间好好研究下。

Edinburgh

zswseu 发表于 2016-9-12 17:45
高阶统计量3,4阶以上呢，是不是阶数越高越好呢，以后有时间好好研究下。

阶数有什么影响？？

zswseu

Edinburgh 发表于 2016-9-13 08:48
阶数有什么影响？？

我也不清楚呀，阶数的选取，可能会根据分析的结构合理的选择吧。

Edinburgh

zswseu 发表于 2016-9-13 10:39
我也不清楚呀，阶数的选取，可能会根据分析的结构合理的选择吧。

阶数对精度没有影响吗

zswseu

Edinburgh 发表于 2016-9-13 12:50
阶数对精度没有影响吗

这方面我也只是听说而已，还没系统研究，还等大神来解答吧。

Edinburgh

zswseu 发表于 2016-9-13 15:36
这方面我也只是听说而已，还没系统研究，还等大神来解答吧。

搜噶  应该是对精度有影响的