某条河流的一个水文站从1915年到1973年记录的最大径流量见下，共59个数据。

15600 8960 10400 10600 10800 9880 9850 10900 8810 9960 12200 7510 8640 6380 6810 8820 14400 7440 7240 6430 11000 7340 9260 5290 9130 7480 6980 9650 7260 8750 9900 7310 9040 7310 8850 7840 10700 6190 9610 7580 9990 6150 8250 6030 8980 6180 9630 9490 2340 11100 5090 10900 6490 12600 6640 7430 6760 10000 9300

求自相关函数和偏相关函数(k=1,2,3,.....,15)，保留3位有效数字。

用matlab编写程序，计算出，，并做出其图像。

在Matlab编辑器中编写以下程序：

clear ;
z = [15600 8960 10400 10600 10800 9880 9850 10900 8810 9960 ...12200 7510 8640 6380 6810 8820 14400 7440 7240 6430 11000 ...7340 9260 5290 9130 7480 6980 9650 7260 8750 9900 7310 9040 ...7310 8850 7840 10700 6190 9610 7580 9990 6150 8250 6030 8980 ...6180 9630 9490 2340 11100 5090 10900 6490 12600 6640 7430 6760 ...10000 9300] ;
%计算z的平均值
m=mean(z);
%z的标准化
for k = 1:59 w(k) = z(k) - m;
endx = 0:1:15;
figure;
subplot(2,1,1);
PartialACF=parcorr(z,15);                %用这个函数，直接可以求得偏相关函数的值
plot(x,PartialACF,'-b.');                %画曲线
xlabel('k取值');
ylabel('偏相关函数值');
title('偏相关函数');subplot(2,1,2);
ACF = autocorr(z,15);                     %用这个函数直接可以求得自相关函数的值
plot(x,ACF,'-b.');                        %画曲线
xlabel('k取值');
ylabel('自相关函数值');
title('自相关函数');

 若要求得得自相关函数和偏相关函数得值取有效数字3位，可以使用vpa(a,b)这个函数。

在这个函数语句中，a 是所求得的这一组数,也就是所求对象；b是保留几位有效数字。

在本题目中，可以在命令窗口输入：

 vpa([PartialACF],3)

 

 vpa([ACF],3)

即可得出最后结果。