matlab 分數階0

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