做好备份ARMA GM

clc;clear;close all;
A=[0.02,0.04,0.09,0.12,0.15,0.12,0.08,0.05,0.09,0.11,0.19,0.24,0.35,0.33,0.38,0.42,0.55,0.68,0.45,0.32,0.33,0.24,0.27,0.32,0.36,0.44,0.46,0.57,0.62,0.75];


plot(1:30,A,'-ro'); %画图
hold on
B=[0.4,0.42,0.435,0.454,0.475,0.50,0.53,0.58,0.62,0.65]
plot(21:30,B,'-*');


hold on
C=[0.33,0.24,0.27,0.32,0.36,0.44,0.46,0.57,0.62,0.75]


D=C-B
plot(21:30,D,'-gd');
% B=[0.33,0.34,0.25,0.31,0.35,0.41,0.48,0.48,0.62,0.63]
% plot(21:30,B,'-*');
% % 
%  hold on
% C=[0.33,0.24,0.27,0.32,0.36,0.44,0.46,0.57,0.62,0.75]
%  
%  D=C-B
 plot(21:30,D,'-gd');
 legend('实际值','GM')
 xlabel('时间刻度')
 ylabel('冲击地压态势')