部分代码（内含3个例子，1分绝对划算）： function yc0=GM1(x0,N) %x0 is the original data , %N is the number you want to predict. n=length(x0); for i=1:n x1(i)=sum(x0(1:i)); end % for k=2:n % 紧邻均生成 z % z(k)=0.5*x1(k)+0.5*x1(k-1);
function [yc0] gm11 x0 number n length x0 ; global N for i 1:n x1 i sum x0 1:i ; end % for k 2:n % 紧邻均生成 z % z k 0 5 x1 k +0 5 x1 k 1 ; % end z x1 1:n 1 +x1 2:n 2; % for i 1:n 1 % b i 1 z i+1 ; % y i x0 i+1 ; % end b : 1 z"; b : 2 1; y x0 2:n ;
