p81_polyfit.m

来自「数值分析, 同济大学教材<<现代数值数学和计算>>数值分」· M 代码 · 共 10 行

x = [0 0.25 0.50 0.75 1.00]';
y = [1.00 1.284 1.6487 2.1170 2.7183]';
x1 = ones(5, 1);
A = [x1 x x.^2]
Z = A \ y

p = polyfit(x, y, 2)