clear all;clc;close all
%〈灰色系统理论及其应用〉第7版,刘思峰 P193例9.2.1
%对X00数据序列,生成m(m=p/q<1)阶累加序列X01。
X00=[247.839,273.021,289.014,285.208,288.818,297.078];
% k=6;
X01=[];
for k = 1:length(X00)
s=0;
for i=1:k
m=0.1; %m-代表分数阶
product=1;
if i<k
for j=0:(k-i-1)
product=product*(m+j);
end
s=s+ product/factorial(k-i)*X00(i);
else
s=s+X00(k);
end
end
X01=[X01,s];
end
X01
%〈灰色系统理论及其应用〉第7版,刘思峰 P193例9.2.1
%对X00数据序列,生成m(m=p/q<1)阶累加序列X01。
X00=[247.839,273.021,289.014,285.208,288.818,297.078];
% k=6;
X01=[];
for k = 1:length(X00)
s=0;
for i=1:k
m=0.1; %m-代表分数阶
product=1;
if i<k
for j=0:(k-i-1)
product=product*(m+j);
end
s=s+ product/factorial(k-i)*X00(i);
else
s=s+X00(k);
end
end
X01=[X01,s];
end
X01