24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1835  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 春秋思良 的主题更新
信息提示
请填处理意见