matlab求峭度，用公式和自己编写的代码分别求，结果不一样

2658744907

（1）用公式：
x=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32];
y=kurtosis(x)
结果：
y =

    1.7977
（2）用代码：（按matlab的公式写的）
x=[ 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32];
a=mean(x);
b=(x-a).*(x-a).*(x-a).*(x-a);
f=sum(b);
c=f/32;
e=std(x).*std(x).*std(x).*std(x);
d=c/e
结果：
d =

    1.6871这是为什么呢？

ChaChing

直觉是std中n与n-1之间的差异,所以好奇查了下kurtosis的算法(可以edit kurtosis查),不难发生其中是以mean的方式求解
也就是说kurtosis的std是採用n的方式

类似matlab内存函数方式

1. xm=x-mean(x); k=mean(xm.^4)/(mean(xm.^2))^2

复制代码

直接公式方式

1. xm=x-mean(x); k=mean(xm.^4)/(std(x,1))^4

复制代码

ps:可以的话，请注意解决问题的过程

2658744907

刚刚接触matlab，多谢指点！！！

2658744907

ChaChing 发表于 2014-8-14 22:12
直觉是std中n与n-1之间的差异,所以好奇查了下kurtosis的算法(可以edit kurtosis查),不难发生其中是以mean的 ...

多谢指点