24小时热门版块排行榜    

查看: 2014  |  回复: 3

春秋思良

金虫 (正式写手)

[求助] 重金求助高手解决二元高次方程组MatlAB求解问题

em=1;f1=0.05;f2=0.1;
a=70;
lt=1/2/(a^2-1)^1.5*(a^2*atan(sqrt(a^2-1))-sqrt(a^-1));
ll=1-2*lt;
ef1=1.8-0.4i;ef2=2.9-0.8i;
1+1/(3/(2*f1*(x(1)-em)/(em+ll*(x(1)-em))+f1*(x(2)-em)/(em+lt*(x(2)-em)))-1/3)-ef1/em=0; %方程(1);
1+1/(3/(2*f2*(x(1)-em)/(em+ll*(x(1)-em))+f2*(x(2)-em)/(em+lt*(x(2)-em)))-1/3)-ef2/em=0; %方程(2);
联立(1)(2)求方程组的解,式中想想x(1)和x(2)为待求量。其他为已知量。
我用matlab中solve求解时没有显式解,后来应用[x,fval]=fsolve(@myfun,x0) 求解时,不同的初始值x0得到不同的解。上述方程应该为二元四次方程组,可能会有4组解,如何使用MatlAB把各个解都计算出来呢?请高手帮忙啊!

原始公式.JPG



文本.JPG

» 猜你喜欢

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

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

1225214

金虫 (正式写手)

【答案】应助回帖


06233053: 金币-1, 求助帖,请勿灌水~ 2012-12-01 15:36:46
唉。。。。难啊。。你去问问叫兽吧。。。。。

[ 发自手机版 https://muchong.com/3g ]
2楼2012-11-21 22:12:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

我觉得你的方程有点问题啊,按你大括号里的方程(2Gx+Gy)/3的值不同啊,这是不可能的吧,不管你上面的u是什么参数,一组u只能得到一个相同的Gx和Gy。
showmethemoney
3楼2012-11-22 11:47:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

春秋思良

金虫 (正式写手)

引用回帖:
3楼: Originally posted by csgt0 at 2012-11-22 11:47:37
我觉得你的方程有点问题啊,按你大括号里的方程(2Gx+Gy)/3的值不同啊,这是不可能的吧,不管你上面的u是什么参数,一组u只能得到一个相同的Gx和Gy。

谢谢你的提醒,图上方程列的有点问题,那个f是浓度,应该为方程组中的两个方程对应两个不同浓度的G。问题已经解决了,用1stopt算的。

[ 发自手机版 https://muchong.com/3g ]
4楼2012-11-23 00:14:10
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 春秋思良 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] fileCode有新解读? +9 Tide man 2026-08-08 17/850 2026-08-08 14:56 by archvillain
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +5 a089 2026-08-07 5/250 2026-08-08 14:04 by alian_214
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 10:39 by 3OOjAIS77qg2
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见