clc?%清屏，以使結果獨立顯示
format?long;?%設置計算精度
x=[10.9011.5711.3611.4911.82];
if?length（x（:1））==1?%對輸入矩陣進行判斷，如不是一維列矩陣，進行轉置變換
????x=x‘;
end
n=length（x）;?%取輸入數據的樣本量
z=0;
for?i=1:n?%計算累加值，并將值賦予矩陣be
????z=z+x（i:）;
????be（i:）=z;
end
for?i=2:n?%對原始數列平行移位
????y（i-1:）=x（i:）;
end
for?i=1:n-1?%計算數據矩陣B的第一列數據
????c（i:）=0.5*（be（i:）+be（i+1:））;
end
for?j=1:n-1?%計算數據矩陣B的第二列數據
????e（j:）=1;
end
for?i=1:n-1?%構造數據矩陣B
????