矩阵的每一行代表一个方程,m行代表m个线性联立方程。 超定方程组:方程个数大于未知量个数的方程组。对于方程组Ra=y,R 为nm矩阵,如果R 超定方程一般是不存在解的矛盾方程。例如,如果给定的三点不在一条直线上, 我们将无法得到这样一条直线,使得这条直 线同时经过给定这三个点。 也就是说给定的条件(限制)过于严格, 导致解不存在。在实 验数据处理和曲线拟合问题中,求解超定方程组非常普遍。比较常用的方法是最小二乘法。 形象的说,就是在无法完全满足给定的这些条件的情况下,求一个最接近的解。 曲线拟合的最小二乘法要解决的问题,实际上就是求以上超定方程组的最小二乘解的 问题。 欠定方程组: 方程个数小于未知量个数的方程组。 对于方程组Ra=y,R 为nm矩阵,且n1e-8 %设置定位精度 j=1:8Dj=xyz(j,:)-P'; dj=norm(Dj); %计算接收机到第j 颗卫星向量及其长度 G(j,1:3)=Dj/dj;G(j,4)=-1; L(j,:)=dj-ro(j); %计算几何矩阵第j 行及方程右端值 end %解几何方程error=abs(b(1:3)); %计算位值修正量绝对值 P=P+b(1:3);k=k+1; %修正定位数据 end toc %程序运行计时结束 format long %显示满足精度的定位数据Format 计算结果 elapsed_time 9.613338291042342e+005-5.674076369901314e+006 2.740537661301808e+006 Column -3.600000000248986e+004