Page 319 - 数学建模算法与应用
P. 319

第十一章  时间序列分析研究


                   end


                   % 显示一次移动平均的结果

                   disp(‘ 一次移动平均值 yhat1:’);
                   disp(yhat1);


                   % 计算一次移动平均的长度

                   m2 = length(yhat1);


                   % 计算二次移动平均
                   yhat2 = zeros(1, m2 - n + 1);

                   for i = 1:m2 - n + 1
                       yhat2(i) = sum(yhat1(i:i + n - 1)) / n;
                   end



                   % 显示二次移动平均的结果
                   disp(‘ 二次移动平均值 yhat2:’);
                   disp(yhat2);



                   % 绘制原始数据
                   figure;
                   plot(1:m1, y, ‘*’);
                   title(‘ 原始数据 ’);

                   xlabel(‘ 时间 ’);
                   ylabel(‘ 值 ’);
                   grid on;



                   % 计算趋势平滑系数
                   a21 = 2 * yhat1(end) - yhat2(end);
                   b21 = 2 * (yhat1(end) - yhat2(end)) / (n - 1);



                                                                                      309
   314   315   316   317   318   319   320   321   322   323   324