????sum=zeros（131）;
????z=zeros（131）;
????A=xlsread（‘E:\C.xlsx‘）;
????x=A（:1）;%時間序
????y=A（:2）;%徑流數據列
????n=length（x）;
????m=length（y）;
????s=0;%設置和為0
????for?i=1:n
????????for?j=1:n