function [f,A]=gm(x)
%输入x,原始序列即x0(k);
%f-->初始序列x0(k)模拟后的随时间他t变化的序列
%A-->GM(1,1)模型的a,b参数组成的向量.
x1=ones(1,length(x));
z1=ones(1,length(x)-1);
for k=1:length(x)
x1(k)=sum(x(1:k));
end
for k=1:length(x)-1
z1(k)=mean(x1(k:k+1));
end
Y=x(2:length(x))';
B=[-z1',ones(length(z1),1)];
A=inv(B'**B'*Y;
beta=A(2)/(1+.5*A(1));
alpha=A(1)/(1+.5*A(1));
beta=num2str(beta);
alpha=num2str(alpha);
xx=num2str(x(1));
AA=num2str(A(1));
ft=strcat('(',beta,'-(',alpha,')*',xx,')*exp(-(',AA,'*(t-2)))');
f=inline(ft,'t');
ksjfkaj
sfakjldjf'sdjf
f
jsadofjkdsjfdsagf
dsagdsagdgsadgfd
sg
dsgadsag
sdgdsag
dsf
klajfdklsafk'gggggggggggggggaskjg
s
afkjkds
jgkagsag
dfaljgfd
sgnction [f,A]=gm(x)
%输入x,原始序列即x0(k);
%f-->初始序列x0(k)模拟后的随时间他t变化的序列
%A-->GM(1,1)模型的a,b参数组成的向量.
x1=ones(1,length(x));
z1=ones(1,length(x)-1);
for k=1:length(x)
x1(k)=sum(x(1:k));
end
for k=1:length(x)-1
z1(k)=mean(x1(k:k+1));
end
Y=x(2:length(x))';
B=[-z1',ones(length(z1),1)];
A=inv(B'**B'*Y;
beta=A(2)/(1+.5*A(1));
alpha=A(1)/(1+.5*A(1));
beta=num2str(beta);
alpha=num2str(alpha);
xx=num2str(x(1));
AA=num2str(A(1));
ft=strcat('(',beta,'-(',alpha,')*',xx,')*exp(-(',AA,'*(t-2)))');
f=inline(ft,'t');
ksjfkaj
sfakjldjf'sdjf
f
jsadofjkdsjfdsagf
dsagdsagdgsadgfd
sg
dsgadsag
sdgdsagnction [f,A]=gm(x)
%输入x,原始序列即x0(k);
%f-->初始序列x0(k)模拟后的随时间他t变化的序列
%A-->GM(1,1)模型的a,b参数组成的向量.
x1=ones(1,length(x));
z1=ones(1,length(x)-1);
for k=1:length(x)
x1(k)=sum(x(1:k));
end
for k=1:length(x)-1
z1(k)=mean(x1(k:k+1));
end
Y=x(2:length(x))';
B=[-z1',ones(length(z1),1)];
A=inv(B'**B'*Y;
beta=A(2)/(1+.5*A(1));
alpha=A(1)/(1+.5*A(1));
beta=num2str(beta);
alpha=num2str(alpha);
xx=num2str(x(1));
AA=num2str(A(1));
ft=strcat('(',beta,'-(',alpha,')*',xx,')*exp(-(',AA,'*(t-2)))');
f=inline(ft,'t');
ksjfkaj
sfakjldjf'sdjf
f
jsadofjkdsjfdsagf
dsagdsagdgsadgfd
sg
dsgadsag
sdgdsagnction [f,A]=gm(x)
%输入x,原始序列即x0(k);
%f-->初始序列x0(k)模拟后的随时间他t变化的序列
%A-->GM(1,1)模型的a,b参数组成的向量.
x1=ones(1,length(x));
z1=ones(1,length(x)-1);
for k=1:length(x)
x1(k)=sum(x(1:k));
end
for k=1:length(x)-1
z1(k)=mean(x1(k:k+1));
end
Y=x(2:length(x))';
B=[-z1',ones(length(z1),1)];
A=inv(B'**B'*Y;
beta=A(2)/(1+.5*A(1));
alpha=A(1)/(1+.5*A(1));
beta=num2str(beta);
alpha=num2str(alpha);
xx=num2str(x(1));
AA=num2str(A(1));
ft=strcat('(',beta,'-(',alpha,')*',xx,')*exp(-(',AA,'*(t-2)))');
f=inline(ft,'t');
ksjfkaj
sfakjldjf'sdjf
f
jsadofjkdsjfdsagf
dsagdsagdgsadgfd
sg
dsgadsag
sdgdsagnction [f,A]=gm(x)
%输入x,原始序列即x0(k);
%f-->初始序列x0(k)模拟后的随时间他t变化的序列
%A-->GM(1,1)模型的a,b参数组成的向量.
x1=ones(1,length(x));
z1=ones(1,length(x)-1);
for k=1:length(x)
x1(k)=sum(x(1:k));
end
for k=1:length(x)-1
z1(k)=mean(x1(k:k+1));
end
Y=x(2:length(x))';
B=[-z1',ones(length(z1),1)];
A=inv(B'**B'*Y;
beta=A(2)/(1+.5*A(1));
alpha=A(1)/(1+.5*A(1));
beta=num2str(beta);
alpha=num2str(alpha);
xx=num2str(x(1));
AA=num2str(A(1));
ft=strcat('(',beta,'-(',alpha,')*',xx,')*exp(-(',AA,'*(t-2)))');
f=inline(ft,'t');
ksjfkaj
sfakjldjf'sdjf
f
jsadofjkdsjfdsagf
dsagdsagdgsadgfd
sg
dsgadsag
sdgdsagnction [f,A]=gm(x)
%输入x,原始序列即x0(k);
%f-->初始序列x0(k)模拟后的随时间他t变化的序列
%A-->GM(1,1)模型的a,b参数组成的向量.
x1=ones(1,length(x));
z1=ones(1,length(x)-1);
for k=1:length(x)
x1(k)=sum(x(1:k));
end
for k=1:length(x)-1
z1(k)=mean(x1(k:k+1));
end
Y=x(2:length(x))';
B=[-z1',ones(length(z1),1)];
A=inv(B'**B'*Y;
beta=A(2)/(1+.5*A(1));
alpha=A(1)/(1+.5*A(1));
beta=num2str(beta);
alpha=num2str(alpha);
xx=num2str(x(1));
AA=num2str(A(1));
ft=strcat('(',beta,'-(',alpha,')*',xx,')*exp(-(',AA,'*(t-2)))');
f=inline(ft,'t');
ksjfkaj
sfakjldjf'sdjf
f
jsadofjkdsjfdsagf
dsagdsagdgsadgfd
sg
dsgadsag
sdgdsag