程序比较长,数据也很多,截其中一部分吧,具体数据就不写出来了
就是很简单的矩阵运算,求解一个四元方程,得到的解回带回去就不等于0.
b(1,1)=-U1(1,5);
b(2,1)=-U12(1,5);
b(3,1)=-U123(1,5);
b(4,1)=-U123(3,5);
A(1,1)=U1(1,2);
A(2,1)=U12(1,2);
A(3,1)=U123(1,2);
A(4,1)=U123(3,2);
A(1,2)=U1(1,4);
A(2,2)=U12(1,4);
A(3,2)=U123(1,4);
A(4,2)=U123(3,4);
A(2,3)=U2(1,4);
A(3,3)=U3(1,1)*U2(1,4)+U3(1,2)*U2(2,4)+U3(1,3)*U2(3,4)+U3(1,4)*U2(4,4);
A(4,3)=U3(3,1)*U2(1,4)+U3(3,2)*U2(2,4)+U3(3,3)*U2(3,4)+U3(3,4)*U2(4,4);
A(3,4)=U3(1,4);
A(4,4)=U3(3,4);
X1=A\b;
XL(1,1)=A(3,1)*X1(1,1)+A(3,2)*X1(2,1)+A(3,3)*X1(3,1)+A(3,4)*X1(4,1)-b(3,1)
涉及到的就是这一小部分,都是很简单的运算。
XL(1,1)应该等于0,但结果却不是。
单独计算了A(3,1)*X1(1,1)+A(3,2)*X1(2,1)+A(3,3)*X1(3,1)+A(3,4)*X1(4,1)和b(3,1),两者结果相等,但是相减后就不等于0了。
情况就是这样的,非常简单的情况,但是却一直得不到正确结果。谢谢解答!