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

