系数矩阵 A 的行列式为范德蒙 (Vandermonde) 行列式,
且五个坐标点的横坐标各不相同,则该行列式不等于零,
所以方程组有唯一解。
syms a0 a1 a2 a3 a4
eq1=sym(‘a0=-27’);
eq2=sym(‘a0+a1+a2+a3+a4=0');
eq3=sym(' a0+2*a1+2^2*a2+2^3*a3+2^4*a4=21 ');
eq4=sym(' a0+3*a1+3^2*a2+3^3*a3+3^4*a4=0 ');
eq5=sym(' a0+4*a1+4^2*a2+4^3*a3+4^4*a4=-75 ');
[a0 a1 a2 a3 a4]=solve(eq1,eq2,eq3,eq4,eq5)
MATLAB 代码为:
A=[1,0,0,0,0;1,1,1,1,1;1,2,2^2,2^3,2^4;1,3,3^2,3^3,3^4;1,4,4^2,4^3,4^4]
y=[-27;0;21;0;-75]
U=rref([A,y])
或