应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 【求助】帮忙求解四元方程组
51/1返回列表
查看: 1352  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jxguo86727

银虫 (小有名气)

[交流] 【求助】帮忙求解四元方程组 已有4人参与

ab/d=0.005754
ac/(1-d)=4.169×〖10〗^(-6)
(0.2-a)/(0.25-b)=4/d
(0.2-a)/(0.03-c)=24/7d
上述四式构成的方程组，求解未知数a、b、c、d的值，a的取值范围0~0.2，b的取值范围0~0.25，c的取值范围0~0.03，d的取值范围0~1。

PS:貌似应该用数值法求解，如果哪位大虾能帮忙计算出来，还望把数值计算方法大体上一说，还有收敛条件、设定的初始值、计算次数等。

不胜感激，必有重谢哦。。。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-10-09 15:13:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wayn1bz

铁虫 (小有名气)

mathematica求解

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
dongdong3881(金币+2): 多谢指导 2011-02-22 13:31:37

In[27]:= Solve[{a b/d == 0.005754,
  a c/(1 - d) == 4.169*10^-6, (0.2 - a)/(0.25 - b) ==
   4/d, (0.2 - a)/(0.03 - c) == 24/7 d}, {a, b, c, d}]

Out[27]=
{{c -> -0.0000465351, b -> 0.19178, a -> 0.0451108,
  d -> 1.50353}, {c -> 0.0000333019, b -> -0.0353554, a -> 0.542452,
  d -> -3.33309}, {c -> 0.0000616971, b -> -0.115831, a -> -0.187562,
  d -> 3.77572}, {c -> 0.0298603, b -> -41.2829, a -> -0.0820862,
  d -> 588.939}, {c -> 0.0299882, b -> -41.3689, a -> 0.282085,
  d -> -2028.08}}
赞一下(2人) 回复此楼
工作是一种谋生，科研同样是一种谋生手段。
10楼2011-02-22 08:44:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

jxguo86727

银虫 (小有名气)
引用回帖:
Originally posted by zxzj05 at 2010-10-09 15:39:31:
试试maple

能麻烦大虾帮忙算一下吗？
赞一下 回复此楼
3楼2010-10-09 16:26:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)
我好久没用了
机器上也没安装
不好意思
赞一下 回复此楼
储氢家族欢迎储氢研究者！
4楼2010-10-09 17:12:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★ ★ ★ ★ ★
jxguo86727(金币+20): 2010-10-09 23:17:03
mingdong(金币+5):感谢交流！ 2010-10-12 10:19:50
Maxima 的结果。
CODE:
(%i1) eq1: a*b/d=0.005754;
                                a b
(%o1)                           --- = 0.005754
                                 d
(%i2) eq2: a*c/(1-d)=4.169*10^(-6);
                          a c
(%o2)                    ----- = 4.1689999999999994E-6
                         1 - d
(%i3) eq3: (0.2-a)/(0.25-b)=4/d;
                                 0.2 - a    4
(%o3)                            -------- = -
                                 0.25 - b   d
(%i4) eq4: (0.2-a)/(0.03-c)=24/(7*d);
                                0.2 - a    24
(%o4)                           -------- = ---
                                0.03 - c   7 d
(%i5) solve([eq1,eq2,eq3,eq4],[a,b,c,d]);

rat: replaced -0.005754 by -433/75252 = -0.005753999893691

rat: replaced -4.168999999999999E-6 by -37/8875030 = -4.168999992112703E-6

rat: replaced 0.2 by 1/5 = 0.2

rat: replaced 0.25 by 1/4 = 0.25

rat: replaced 0.2 by 1/5 = 0.2

rat: replaced 0.03 by 3/100 = 0.03
(%o5) [[a = 0.18571751923947, b = 0.22416573626101,
c = - 1.3996884948583785E-4, d = 7.235228539576366],
[a = 0.01412678085481, b = 0.22439931350114, c = 1.3252708943556937E-4,
d = 0.55092802172929], [a = 1.3977106364338493E-4, b = 0.24969694088565,
c = 0.029646433041302, d = 0.0060654157927555]]

[ Last edited by snoopyzhao on 2010-10-9 at 23:22 ]
赞一下(1人) 回复此楼
5楼2010-10-09 18:14:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭