%**************************************************************************************************************************************************************
%   本程序為函數功能，目的是為求取一列信號數據的偏度系數和峭度系數
%      
%   中國北方發動機研究所試驗測試中心  蘇志霄   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;