clc
load('E:\hua_1\程序\網絡延時測試程序客戶端--UDP協議\TestSock_Cnt\平移gamma分布\TCPtime.txt');
[num1,num]=size(TCPtime) ;   %num 數據量
num
MeanTCPtime=mean(TCPtime)    %平均值
VarianceTCPtime = std(TCPtime,1)^2  %方差


SUM=0;
for i = 1:num
    tempx = (TCPtime(i) - MeanTCPtime)^3; 
    SUM = SUM+tempx;
end
SkenessTCPtime = (SUM/num) / (VarianceTCPtime^(3/2))  %扭量
syms x y z
[x,y,z] = solve('x*y+z-MeanTCPtime','x*(y^2)-VarianceTCPtime','2/(x^(1/2))-SkenessTCPtime','x','y','z');   %矩估計
aph=eval(x)              %參數輸出
beta=eval(y)
gamma=eval(z)