使用R solve(a,b)解方程 - 1

简介:
解普通方程组可以用函数solve(),solve()的基本用法是solve(A,b),其中,A为方程组的系数矩阵,b为方程组的右端。例如:

已知方程组 :
2x1+2x3=1
2x1+x2+2x3=2
2x1+x2=3


解法如下 : 
方程组等效于 :
2x1+0x2+2x3=1
2x1+1x2+2x3=2
2x1+1x2+0x3=3

> A
     [,1] [,2] [,3]
[1,]    2    0    2
[2,]    2    1    2
[3,]    2    1    0

> b=1:3
> b
[1] 1 2 3

> solve(A,b)
[1]  1.0  1.0 -0.5

即x1=1,x2=1,x3=-0.5。

这个解法用到公式 : 
A %*% solve(A,b)  ==  b
其中 : 
> A
     [,1] [,2] [,3]
[1,]    2    0    2
[2,]    2    1    2
[3,]    2    1    0

solve(A,b)
x1
x2
x3

b
1
2
3


[参考]
1. help(solve)
相关文章
|
2月前
|
机器学习/深度学习
PTA-方阵转置
方阵转置
18 1
|
11月前
|
C语言 C++
C++/C/PTA 找鞍点
一个矩阵元素的“鞍点”是指该位置上的元素值在该行上最大、在该列上最小。
87 0
欧拉计划Problem 5 最小公倍数
欧拉计划Problem 5 最小公倍数
【PTA】7-6 求最大公约数 (40point(s))
【PTA】7-6 求最大公约数 (40point(s))
166 0
|
机器学习/深度学习
【PTA】7-1 矩阵运算
【PTA】7-1 矩阵运算
2502 0
|
人工智能 算法 BI
三对角线性方程组(tridiagonal systems of equations)的求解
三对角线性方程组(tridiagonal systems of equations)   三对角线性方程组,对于熟悉数值分析的同学来说,并不陌生,它经常出现在微分方程的数值求解和三次样条函数的插值问题中。
1670 0