方程数大于未知量个数,属于超定方程组,一般没有精确解,只能求最小二乘解.matlab中可以利用左除或伪逆法求解.
clear
c1=10496.108; c2=10500.275; c3=10504.442; c4=10508.609; c5=10512.776;
A=[1,c1,c1^2,c1^3; 1,c2,c2^2,c2^3; 1,c3,c3^2,c3^3; 1,c4,c4^2,c4^3; 1,c5,c5^2,c5^3;];
B=[-3043766.03 ;-3042809.53;-3041787.84;-3040701.05;-3039549.25];
xyzt=pinv(A)*B %伪逆法求解
%xyzt=A\B %左除法
xyzt =
1.0e+004 *
0.0005
1.8316
-0.0004
0.0000可是我把xyzt的数值带进去返算的结果不对啊？？？结果是对的，我既然放出答案，当然自己有验证过。我估计你是直接吧上面的结果带进去？上面是MTALAB默认的显示格式，只显示5位有效数字的。你再试试：xyzt =5.231005119 18316.49731 -3.59416421 0.0001735367694