www.gusucode.com > 灰色GM(1，3)神经网络matlab源码程序 > GM1_3.m

    function [ y ] = GM1_3( x0,x1,x2,x3 )
%UNTITLED12 Summary of this function goes here
%   Detailed explanation goes here

x10=x0;
x20=x1;
x30=x2;
x40=x3;

x11=leijia(x10);
x21=leijia(x20);
x31=leijia(x30);
x41=leijia(x40);

z1=jinlinxulie(x10,0.5);

B=[-z1(2:end)' x21(2:end)' x31(2:end)' x41(2:end)'];
Y=[x10(2:end)]';

D=inv(B'*B)*B'*Y;

n=length(x10);
x0(1)=x10(1);
for k=2:n
    x0(k)=D(2:end)'*[x21(k) x31(k) x41(k)]'-D(1)*z1(k);
end
y=x0;
end