pianduqiaodu.m

自己编写的比较全面的故障诊断matlab函数程序

%**************************************************************************************************************************************************************
%   本程序为函数功能，目的是为求取一列信号数据的偏度系数和峭度系数
%      
%   中国北方发动机研究所试验测试中心  苏志霄   2006年10月
%**************************************************************************************************************************************************************
%
function R=pianduqiaodu(data)
n=size(data,1);
junzhi=mean(data);                               %计算均值
fangcha=mean((data-junzhi).^2);                  %计算方差
bianliang=(data-junzhi)/fangcha;
piandu=sum(bianliang.^3)*sqrt(1/(6*n));          %计算偏度系数
qiaodu=(sum(bianliang.^4)/n-3)*sqrt(n/24);       %计算峭度系数
R(1)=piandu;
R(2)=qiaodu;