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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 求解一个方程组，急！！！
52/1返回列表
查看: 1502  |  回复: 19
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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
要求转化为一个关于x（或y,z）的一元三次方程。
我自己手算过很多遍，不知道哪里出问题，都是一元四次方程。但是按照老师的说法是存在一元三次这种情况的，我一直没找到。课题进展就一直卡在这里下不去了。
希望有数学高手能够帮忙手算一下，不胜感激！
采用Mathematica计算出来的结果与手算结果一致，是不是表明只有一元四次这种情况？
回复此楼

» 猜你喜欢
1楼2014-07-25 10:03:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
14楼: Originally posted by 309065816 at 2014-07-28 08:49:54
问题确实没有解决，不过估计确实如你所说，后面也不会有人应助什么的，金币都给你了。...

楼主你是不是对mathematica的用法不很懂?
从前面的答案看是已经解决了你的疑问: 即结果是一元三次的多项式形式.
数学上没有任何特别之处. 多项式方程组的求解这种规模也没有什么难度,只是变量变化形式繁琐而已. 符号计算软件处理这种问题失手的机会几乎为零,除非软件有bug. 这里的符号计算的结果也很容易可视化成TeXForm,就跟用LaTeX写的公式看上去一样; 完全没有必要追求手写的结果. 你是希望别人手写了之后拍照贴出来吗? 如果有什么特别的要求,应该在提问或者补充说明的时候讲出来; 否则在别人看来很容易的事,可能恰恰就是你觉得很难没有解决的.
赞一下 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
15楼2014-07-28 11:45:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

Mr__Right

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
309065816: 金币+100 2014-07-26 09:44:00
最高次数的确是3
你用手算能算得出来才怪,所以,如果你手算了,很可能是有问题的.
得到的式子比较复杂, 比如, 其中一项是 x*Sqrt[a x^4+b x^3 + c x^2 + d]的形式, 也可以理解为不超过3次

用符号计算的时候, 为什么出现k1, k01, k2, k02并存? 还是本来相同? 不超过10的时候完全可以编号下去

解题思路:
假设三个多项式分别为 f1[x,y,z]==0, f2[x,y,z]==0, f3[x,y,z]==0
需要定义三个关于x,y,z函数f1,f2,f3

假设求关于x的多项式,则需要消去y, z
任取两个等式, 比如f2,f3解关于y,z方程组
Solve[{f2[x,y,z]==0,f3[x,y,z]==0},{y,z}]//FullSimplify
得到两组解,假设为Rule01, Rule01

然后把y,z的结果代入f1中化简,
Collect[f1[x,y,z]/.Rule01//FullSimplify//ExpandAll,x]

然后就得到想要的关于x的式子; 最高次数为3

c d k3 k5 + (-d k1 k3 + k02 k1 k4 + k1 k2 k4 - d k3 k4 - k1 k3 k4 +
    k01 k4^2 + k2 k4^2 + k1 k4 k5) x^2 + k1 k4^2 x^3 +
x (c d k1 k3 + k01 k02 k4 + c d k3 k4 - d k3 k5 + k01 k4 k5 +
    k2 k4 k5 +
    k4 \[Sqrt](k01^2 k02^2 + 2 k01^2 k02 k5 + 2 k01 k02 k2 k5 +
        k01^2 k5^2 + 2 k01 k2 k5^2 + k2^2 k5^2 + 2 k01 k02^2 k1 x +
        2 k01 k02 k1 k2 x - 2 k01 k02 k1 k3 x + 2 k01^2 k02 k4 x +
        2 k01 k02 k2 k4 x + 4 k01 k02 k1 k5 x + 2 k01 k1 k2 k5 x +
        2 k02 k1 k2 k5 x + 2 k1 k2^2 k5 x - 2 k01 k1 k3 k5 x +
        2 k1 k2 k3 k5 x + 2 k01^2 k4 k5 x + 4 k01 k2 k4 k5 x +
        2 k2^2 k4 k5 x + 2 k01 k1 k5^2 x + 2 k1 k2 k5^2 x +
        k02^2 k1^2 x^2 + 2 k02 k1^2 k2 x^2 + k1^2 k2^2 x^2 -
        2 k02 k1^2 k3 x^2 + 2 k1^2 k2 k3 x^2 + k1^2 k3^2 x^2 +
        4 k01 k02 k1 k4 x^2 + 2 k01 k1 k2 k4 x^2 +
        2 k02 k1 k2 k4 x^2 + 2 k1 k2^2 k4 x^2 - 2 k01 k1 k3 k4 x^2 +
        2 k1 k2 k3 k4 x^2 + k01^2 k4^2 x^2 + 2 k01 k2 k4^2 x^2 +
        k2^2 k4^2 x^2 + 2 k02 k1^2 k5 x^2 + 2 k1^2 k2 k5 x^2 -
        2 k1^2 k3 k5 x^2 + 4 k01 k1 k4 k5 x^2 + 4 k1 k2 k4 k5 x^2 +
        k1^2 k5^2 x^2 + 2 k02 k1^2 k4 x^3 + 2 k1^2 k2 k4 x^3 -
        2 k1^2 k3 k4 x^3 + 2 k01 k1 k4^2 x^3 + 2 k1 k2 k4^2 x^3 +
        2 k1^2 k4 k5 x^3 + k1^2 k4^2 x^4))
赞一下(3人) 回复此楼
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
2楼2014-07-26 08:00:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)
引用回帖:
2楼: Originally posted by Mr__Right at 2014-07-26 08:00:31
最高次数的确是3
你用手算能算得出来才怪,所以,如果你手算了,很可能是有问题的.
得到的式子比较复杂, 比如, 其中一项是 x*Sqrt的形式, 也可以理解为不超过3次

用符号计算的时候, 为什么出现k1, k01, k2, k02并 ...

先谢谢了。不过老师自己说过博士期间手算出来过，所以一直让我算。k1，k01这些是具有物理意义的，所以就是这么表示的。
你是用什么软件算的呢？matlab还是mathematics？
赞一下 回复此楼
3楼2014-07-26 09:43:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
309065816: 金币+10 2014-07-27 22:54:50
引用回帖:
3楼: Originally posted by 309065816 at 2014-07-26 09:43:47
先谢谢了。不过老师自己说过博士期间手算出来过，所以一直让我算。k1，k01这些是具有物理意义的，所以就是这么表示的。
你是用什么软件算的呢？matlab还是mathematics？...

1.很明显,你看我给出的代码是Mathematica表示的,你如果会用Mathematica应该看得出来;

2.手算除非每计算一步都用变量代换,把同类型系数用一个新的符号表示;否则容易出错; 手算计算效率低,而且碰上稍微有难度的非常吃力;

按照我说的方法应该很容易得到你导师说的结果.
你还有别的疑问吗?
赞一下 回复此楼
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
4楼2014-07-26 09:49:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭