dataprocess.m

在mathematica环境下运行

Sk4=Skewness[a4]
Ku1=Kurtosis[a1]
Ku2=Kurtosis[a2]
Ku3=Kurtosis[a3]
Ku4=Kurtosis[a4]
\[Phi]=2\[Pi]/n;
s1=\[Phi] a1^2/2;
s2=\[Phi] a2^2/2;
s3=\[Phi] a3^2/2;
s4=\[Phi] a4^2/2;
ss1=Apply[Plus,s1];
ss2=Apply[Plus,s2];
ss3=Apply[Plus,s3];
ss4=Apply[Plus,s4];
d1=ss1/(10000 \[Pi])
d2=ss2/(10000 \[Pi])
d3=ss3/(10000 \[Pi])
d4=ss4/(10000 \[Pi])