三元一次方程大家应该是不陌生的,形如
aX + bY + cZ = d
的就是,其中X,Y,Z是未知的变量,a,b,c,d 都是已知的常量,通常呢,需要至少3个没有线性关系的已知等式才能求唯一解。例如,已知:
2X + 3Y - 4Z = 345
5X - 6Y + 7Z = 361
8X + 9Y - 10Z = 235
就可以求出唯一解:X= -984.7667 Y= -61.2 Z= 327.5667
看起来确实有点难度哦!
不知道多少人还记得数学解法,主要是利用消元思想使三元变二元,再变一元。我搜索了一下,是如下3个步骤:
①利用代入法或加减法,消去一个未知数,得出一个二元一次方程组;
②解这个二元一次方程组,求得两个未知数的值;
③将这两个未知数的值代入原方程中较简单的一个方程,求出第三个未知数的值,把这三个数写在一起的就是所求的三元一次方程组的解。
那么,如果是要在R里面操作呢,这个时候矩阵的思想就可以帮上忙了。我先公布代码:
b=c(