24小时热门版块排行榜    

查看: 1351  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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 ]
5楼2010-10-09 18:14:33
已阅   回复此楼   关注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的回帖

jxguo86727

银虫 (小有名气)

引用回帖:
Originally posted by snoopyzhao at 2010-10-09 18:14:33:
Maxima 的结果。

[code]
(%i1) eq1: a*b/d=0.005754;
                                a b
(%o1)                           --- = 0.005754
                                 d
(%i2) eq2: a*c/(1-d)= ...

些许疑问,还请帮忙解决:
1、最后是不是做了三次数值逼近;2、数值逼近的方法是采用什么方法呢;3、设的收敛条件或是数值解的误差是什么。
多谢多谢!
6楼2010-10-09 23:24:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见