如何用matlab编程实现高阶累计量的计算？

上下求索＿quest

本人正在做关于滚动轴承故障诊断的毕业设计 ，现在想计算离散轴承振动信号的高阶累计量，以区分不同的故障特征。
设x(n)为离散时间实值平稳随机过程，其二阶距、三阶距、四阶矩分别定义为

m2x(i)=E［x(n)x(n+i)］　　(4)

m3x(i,j)=E［x(n)x(n+i)x(n+j)］　　(5)

m4x(i,j,k)=E［x(n)x(n+i)x(n+j)x(n+k)］　　(6)

　　若x(n)为零均值平稳随机过程，则其二阶累计量、三阶累计量、四阶累计量分别定义为

c2x(i)=m2x(i)　　(7)

c3x(i,j)=m3x(i,j)　　(8)

c4x(i,j,k)=m4x(i,j,k)-m2x(i)m2x(j-k)-m2x(j)m2x(k-i)-m2x(k)m2x(i-j)　　(9)
请问各位高手，我如何利用matlab实现上述算法，谢谢！