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