24小时热门版块排行榜    

查看: 3084  |  回复: 13

309065816

木虫 (正式写手)

[求助] 请教高手,如何解非线性方程组!!! 已有1人参与

(1)D*(c-x)-k1*x*(1-y-z)+k01*y-k4*x*z+k5*z=0
(2)k1*x*(1-y-z)-k01*y-k2*y+k02*z-k3*y*z=0
(3)k2*y-k02*z+k3*y*z-k4*x*z-k5*z=0
上述方程组中,D,k1~k5均为一些变量参数,想得到x,y,z关于这些参数的表达式。
想通过代入消元或者其他方法将上述方程化为一元多次方程的形式。例如,将y,z用x表示,最终将方程化为关于x的一元多次方程。
我自己计算过,由于方程中含有一些相同项,通过将(1)(2)(3)式相加得到,z关于x的表达式;
将(1)(2)相加,再将z关于x的表达式代入,可得到y关于x的表达式;
再将上述y.z关于x的表达式代入到(3)中,结果为0,没有得到关于x的表达式。
上述过程应该是有问题的,但是具体问题不知道在哪?正确的做法应该是怎么样的?
想问下具体的过程该是怎么样的?
我最终想要的就是将上述3个方程化为一个一元多次方程,请高手指点一下!
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
309065816: 金币+10, ★有帮助, 好的,谢谢了,我先试试,对的话会再给金币的。 2014-07-01 21:15:31
309065816: 金币+10 2014-07-08 12:54:26
先将x看作常量,通过(1)、(2)解得y=y(x)、z=z(x),再代入(3),这样就应该不会出现楼主述及的问题了。如果这样做仍会出现楼主遇到的问题,那只能说明这三个方程不是相互独立的,至少一个方程会有另外两个推出。这样一来,方程组的解如果存在的话就会有无穷多个。
3楼2014-07-01 21:12:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

309065816

木虫 (正式写手)

根据要求,希望能够手动推导得出一元多次方程,matlab中的solve函数求解精度不能满足要求,且存在问题
2楼2014-07-01 20:47:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hank612

至尊木虫 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
309065816(feixiaolin代发): 金币+10 2014-07-09 19:58:01
用 Mathematica 的Groebner Basis 功能吧, 一步到位
自己google它的原理, 反正到处都是
请教高手,如何解非线性方程组!!!
Emuch018.png

We_must_know. We_will_know.
4楼2014-07-02 10:27:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)

引用回帖:
4楼: Originally posted by hank612 at 2014-07-02 10:27:51
用 Mathematica 的Groebner Basis 功能吧, 一步到位
自己google它的原理, 反正到处都是

Emuch018.png

嗯,十分感谢。不过老师是要求我自己手算推导。所以还是要花时间动动笔,不过你这个结果可用来作为参考。
5楼2014-07-02 12:51:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daul_0124

银虫 (初入文坛)

k01,k02何意?
学习
6楼2014-07-02 21:15:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)

引用回帖:
6楼: Originally posted by daul_0124 at 2014-07-02 21:15:31
k01,k02何意?

本来应该是k-1,k-2.表示逆反应速率常数,为了简便起见,就用了k01,k02
7楼2014-07-03 00:19:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)

引用回帖:
4楼: Originally posted by hank612 at 2014-07-02 10:27:51
用 Mathematica 的Groebner Basis 功能吧, 一步到位
自己google它的原理, 反正到处都是

Emuch018.png

想给你发金币,但是你不是应助回帖,不知道怎么才能把金币发给你
8楼2014-07-08 12:55:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)

引用回帖:
4楼: Originally posted by hank612 at 2014-07-02 10:27:51
用 Mathematica 的Groebner Basis 功能吧, 一步到位
自己google它的原理, 反正到处都是

Emuch018.png

我手算推导的结果跟你计算的一样。这个结果应该只有一个吧?你能不能再帮我在试试分别用y,z表示的一元多次方程呢?一元四次方程求解不是很方便。我基本都是手算。
9楼2014-07-08 12:57:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)

引用回帖:
4楼: Originally posted by hank612 at 2014-07-02 10:27:51
用 Mathematica 的Groebner Basis 功能吧, 一步到位
自己google它的原理, 反正到处都是

Emuch018.png

我安装了一个mathematics9.0,点文件-新建-程序包,然后按照你照片里面的输入,然后运行程序包,没有输出结果,
=GroebnerBasis[{D*(c-x)-k1*x*(1-y-z)+k01*y-k4*x*z+k5*z,
k1*x*(1-y-z)-k01*y-k2*y+k02*z-k3*y*z,
k2*y-k02*z+k3*y*z-k4*x*z-k5*z},{x,D,c,k1,k01,k4,k5,k2,k02,k3},{y,z}]
提示:
Syntax::tsntxi: "=GroebnerBasis[{D*(c-x)-k1*x*(1-y-z)+k01*y-k4*x*z+k5*z,k1*x*(1-y-z)-k01*y-k2*y+k02*z-k3*y*z,k2*y-k02*z+k3*y*z-k4*x*z-k5*z},{x,D,c,k1,k01,k4,k5,k2,k02,k3},{y,z}]" 不完整;需要更多输入.
Syntax::sntxi: 不完全表达式;需要更多输入 .
这是怎么回事呢?因为之前一直用matlab,对这个不了解,所以用起来比较费劲。
10楼2014-07-08 15:49:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 309065816 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗? 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[职场人生] 40+中年人失业后,还能找到好工作吗? +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见