matlab 分数阶0<m=p/q<1累加序列生成方法

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