函數gmcal=gm1(x)
sizexd2 = size(x，2)；
%找到數組長度
k = 0；
對於y1=x
k = k+1；
如果k & gt1
x 1(k)= x 1(k-1)+x(k)；
累計發電量百分比
z 1(k-1)=-0.5 *(x 1(k)+x 1(k-1))；
%z1維減去1用於計算b。
yn 1(k-1)= x(k)；
其他
x 1(k)= x(k)；
結束
結束
%x1，z1，k，yn1
sizez1=size(z1，2)；
% size(yn 1)；
z2 = z 1 '；
z3 = ones(1，sizez 1)'；
YN = yn 1 '；%換位
%YN
b =；gm1(x)
%其中5999，5903，5848，5700，7884可以用已知的歷史數據代替，無論多少。