循环统计量，高阶统计量
有人做循环平稳分析的么，涉及到循环双谱等知识的，可以一起来交流下么？ 循环双谱matlab程序function cyclic_bispectrum(x,fs,M1,M2)
%%% 循环双谱程序
%%% x 输入的信号
%%%=====================================================
%%% 第一步：求信号的三次变换,依据为： y(i,n)=x(i)x(i+n)x(i+n) 后两项n取一样，为方便运算。
%%%
M = length(x); % 信号采样点数
if M>20000
M = 20000;
end
% 初始化
N1 = fix(3*M/4); N2 = M - N1; % 必须满足 M=N1+N2，且N1>N2。
y = zeros(N1,N2);
for i=1:N1
for n=1:N2
y(i,n) = x(i)*x(i+n)*x(i+n);
end
end
%%%=====================================================
%%% 第二步：求信号的三次变换,依据为： Rtk(n)=1/N1*sum(y(i,n)*exp(-1i*2*pi*i*k/N1))
%%% M1对应 2pi 弧度，分辨率 2pi/M1
Rtk = zeros(N2,M1);
for n = 1:N2
% Rtk(n,:) = 1/N1*fftshift(fft(y(:,n),M1));
Rtk(n,:) = fft(y(:,n),M1);
end
%%%=====================================================
%%% 第三步：求循环双谱切片谱,依据为： Imk(m)=sum(Rtk(1:N2).*exp(-1i*2*pi*m*(1:N2)/N2))
%%%
Imk = zeros(M1,M2);
for m = 1:M1
% Imk(m,:) = fftshift(fft(Rtk(:,m),M2));
Imk(m,:) = fft(Rtk(:,m),M2);
end
aaa = Imk(1:M1/2,1:M2/2);
x_lab = fs/2*linspace(0,1,M1/2);
y_lab = fs/2*linspace(0,1,M2/2);
% I_abs = abs(Imk);
I_abs = abs(aaa);
%%% 画图
figure
mesh(y_lab,x_lab,I_abs)
% mesh(I_abs)
% figure
% plot(x_lab,I_abs(:,1))
%
figure
plot(x_lab,I_abs)
%
% figure
% plot(x_lab,I_abs(:,round(M1/2)-1)) 循环双谱及在周期平稳类故障中的应用
作者:苏中元； 贾民平； 许飞云； 胡建中；
摘要:论述了未知循环频率的周期平稳信号循环双谱的估计方法;研究了在循环双谱的循环累积量计算中涉及变量的简化存储方法,提出了该变量矩阵是对称阵,通过算法可以三角阵的元素来表述,克服了循环双谱传统估计方法计算量较大的缺陷,提高了运算效率;提出了循环双谱对调相故障信号的分析能力,以及对加性噪声的处理能力,仿真并验证了该方法的有效性,并将其应用于旋转机械状态分析。 amnesia 发表于 2016-9-5 13:35
循环双谱及在周期平稳类故障中的应用
作者:苏中元； 贾民平； 许飞云； 胡建中；
哇 多谢兄台。咱们可以一起交流这个。 楼主研究这个做什么 《基于高阶循环统计量的频谱感知方法》 高阶统计量3,4阶以上呢，是不是阶数越高越好呢，以后有时间好好研究下。
zswseu 发表于 2016-9-12 17:45
高阶统计量3,4阶以上呢，是不是阶数越高越好呢，以后有时间好好研究下。
阶数有什么影响？？ Edinburgh 发表于 2016-9-13 08:48
阶数有什么影响？？
我也不清楚呀，阶数的选取，可能会根据分析的结构合理的选择吧。
zswseu 发表于 2016-9-13 10:39
我也不清楚呀，阶数的选取，可能会根据分析的结构合理的选择吧。
阶数对精度没有影响吗 Edinburgh 发表于 2016-9-13 12:50
阶数对精度没有影响吗
这方面我也只是听说而已，还没系统研究，还等大神来解答吧。
zswseu 发表于 2016-9-13 15:36
这方面我也只是听说而已，还没系统研究，还等大神来解答吧。
搜噶应该是对精度有影响的
页:
[1]