使用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)
目录
相关文章
|
算法 安全 Java
z3-solver求解器
一个非常高级的工具,SMT求解器。应用领域非常广,解各类方程,解各类编程问题(例如解数独),解逻辑题等都不在话下。
|
9月前
|
Python
NumPy 差分、最小公倍数、最大公约数、三角函数详解
NumPy 差分 离散差分意味着相邻元素之间的减法。 例如,对于 [1, 2, 3, 4],离散差分将是 [2-1, 3-2, 4-3] = [1, 1, 1] 要找到离散差分,使用 diff() 函数。
|
10月前
|
Java
HDU-2199-Can you solve this equation
HDU-2199-Can you solve this equation
47 0
|
10月前
|
Java
HDU-2199-Can you solve this equation?
HDU-2199-Can you solve this equation?
54 0
欧拉计划Problem 5 最小公倍数
欧拉计划Problem 5 最小公倍数
129 0
|
人工智能 算法 BI
三对角线性方程组(tridiagonal systems of equations)的求解
三对角线性方程组(tridiagonal systems of equations)   三对角线性方程组,对于熟悉数值分析的同学来说,并不陌生,它经常出现在微分方程的数值求解和三次样条函数的插值问题中。
1760 0