z=[16?12?15?10?9?17?11?16?10?14];?%任給一組數(shù)據(jù)，通過下面函數(shù)求解樣本均值z1、樣本自協(xié)方差函數(shù)r、樣本自相關(guān)函數(shù)p（SACF）和樣本偏自相關(guān)函數(shù)q（SPACF）

N=size（z）;???%size（A1）返回矩陣A的行數(shù)??size（A2）返回矩陣A的列數(shù)
n=N（12）;
z1=sum（z）/n;?%樣本均值z1

s0=0;
for?i=1:n
????s0=（z（i）-z1）^2+s0;
end
r0=s0/n;

%%樣本自協(xié)方差函數(shù)r
for?k=1:（n-1）
????nr（k）=0;
????for?t=1:（n-k）
????????nr（k）=（z（t）-z1）*（z（t+k）-z1）+nr（k）;
????end?
end???
r=nr./n;

%%樣本自相關(guān)函數(shù)p
p=（r./r0）‘;%p=vpa（（r./r0）‘3）;?%vpa（nm）對n保留小數(shù)點(diǎn)后m位

%%偏自相關(guān)函數(shù)q
for?m=1:（n-2）
????q（11）=p（1）;
????D=0;T=0;
????for?h=1:m
????????D=p（m+1-h）*q（mh）+D;
????????T=p（h）*q（mh）+T;
????end
????????q（m+1m+1）=（p（m+1）-D）/（1-T）;
????for?j=1:m
????????q（m+1j）=q（mj）-q（m+1m+1）*q（mm+1-j）;
????end?????
end
q

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-10-30?20:42??程序源代碼\
?????文件????????6148??2017-10-29?20:56??程序源代碼\.DS_Store
?????文件??????220389??2017-10-29?20:56??程序源代碼\2017最新MATLAB從入門到精通直播課程.pdf
?????目錄???????????0??2017-10-30?20:41??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\
?????文件????????1636??2017-10-30?04:41??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM96A.M
?????文件?????????495??2017-10-30?04:41??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM96B.M
?????文件?????????290??2017-10-30?04:41??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM96Bfun.m
?????文件????????1067??2017-10-30?04:41??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM97A.M
?????文件?????????314??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM97Afun.m
?????文件?????????631??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM97Aoptim.m
?????文件????????4353??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM97B.M
?????文件????????4353??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM97Bcount.m
?????文件?????????393??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM97Brule.m
?????文件????????1483??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM98A1.M
?????文件????????1375??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM98A2.M
?????文件?????????862??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM98A3.M
?????文件????????8108??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM98B.C
?????文件?????????751??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM98B.M
?????文件?????????272??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM99A.M
?????文件?????????284??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM99A1.M
?????文件?????????645??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\JM99AFUN.M
?????文件?????????376??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\jm98a3fun.m
?????文件?????????412??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\jm99a1simu.m
?????文件?????????508??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\jm99asmfun.m
?????文件?????????592??2017-10-30?04:42??程序源代碼\中國大學(xué)生數(shù)學(xué)建模競賽題解\新建?文本文檔.txt
?????目錄???????????0??2017-10-30?20:42??程序源代碼\圖形\
?????文件????????5528??2017-10-30?04:34??程序源代碼\圖形\DRAW.M
?????文件????????3626??2017-10-30?04:34??程序源代碼\圖形\MESH.M
?????文件????????1567??2017-10-30?04:34??程序源代碼\圖形\MESHF.M
?????文件????????4385??2017-10-30?04:34??程序源代碼\圖形\PLOT.M
?????文件????????1937??2017-10-30?04:34??程序源代碼\圖形\PLOT3.M
............此處省略81個(gè)文件信息