# w=1800
# a+b=w
# a>b
# a:w=b:a约等于0.618
# 求a的值
#导入sbmpb
from sympy import *
#定义变量
a = Symbol('a')
b = Symbol('b')
print(solve([a+b-1800,a/1800-0.618,b/a-0.618],[a,b]))
怎么解方程···········?
- 写回答
- 好问题 0 提建议
- 追加酬金
- 关注问题
- 邀请回答
-
4条回答 默认 最新
- 琐琐闻 2021-03-11 20:59关注
由于浮点数计算的精度问题会导致方程无解
二元一次方程式的求解只需要两条方程,例如a:w(1800)=0.618,b:a=0.618可以确定a与b的值
我们将其代入代码进行计算
>>> from sympy import * >>> a = Symbol('a') >>> b = Symbol('b') >>> print(solve([a/1800-618/1000,b/a-618/1000],[a,b])) {a: 1112.40000000000, b: 687.463200000000}
此时已经得出唯一解a=1112.4,b=687.4632,并且a+b=1799.863,与a+b=1800有一定的误差;故将此唯一解代回a+b=1800方程时会否定该解,造成无解的现象。
所以说只需要删除三条方程中的任意一条即可。
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 无用评论 打赏 举报
悬赏问题
- ¥15 关于#matlab#的问题:期望的系统闭环传递函数为G(s)=wn^2/s^2+2¢wn+wn^2阻尼系数¢=0.707,使系统具有较小的超调量
- ¥15 FLUENT如何实现在堆积颗粒的上表面加载高斯热源
- ¥30 截图中的mathematics程序转换成matlab
- ¥15 动力学代码报错,维度不匹配
- ¥15 Power query添加列问题
- ¥50 Kubernetes&Fission&Eleasticsearch
- ¥15 報錯:Person is not mapped,如何解決?
- ¥15 c++头文件不能识别CDialog
- ¥15 Excel发现不可读取的内容
- ¥15 关于#stm32#的问题:CANOpen的PDO同步传输问题