%------------------------检验预处理后的数据是否符合AR建模要求,计算自相关和偏相关系数---------------% %---------------------------------------计算自相关系数-----------------------------------% R0=0; for i=1:88 R0=Y(i)^2/88+R0; %标准化处理后的数据的方差 end for k=1:20 %R 协方差 R(k)=0; for i=k+1:88 R(k)=Y(i)*Y(i-k)/88+R(k); end end x=R/R0; %自相关系数x = 协方差/方差 %画图 figure; plot(x) title('DC14->DC10的自相关系数分析图');
时间: 2023-08-15 15:59:10 浏览: 190
这段代码是用来计算预处理后数据的自相关系数,并绘制自相关系数分析图。
首先,代码计算了标准化处理后数据的方差R0。通过循环遍历每个数据,将数据的平方除以数据总数并求和,得到方差R0。
然后,代码使用另一个循环计算了自相关系数R。循环从1到20,对于每个延迟k,计算了Y(i)和Y(i-k)的乘积,并将结果除以数据总数求和,得到协方差R(k)。
接下来,代码将协方差R除以方差R0,得到自相关系数x。
最后,代码绘制了自相关系数分析图,以x为纵轴。标题指明了分析的对象是"DC14->DC10"。
你可以运行这段代码来查看自相关系数分析图。
阅读全文