24小时热门版块排行榜    

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

crazyiyy

木虫 (小有名气)

[求助] MATLAB解非线性方程组 已有1人参与

本人有个非线性方程组,已用1stOpt解出,但需要将结果导入到Matlab中,十分麻烦。因此想用Matlab将此方程组解出来,肯定路过大神帮忙!题目如下
1/((4.5+(8.5)*(x+y-x*y)+(33/2)*(x+y-x*y)^2+(33/2)*(x+y-x*y)^3+(33/2)*(x+y-x*y)^4)*(1-(x+y-x*y))/(1-(x+y-x*y)^5)+(1-(12/(48*2^3-4)))*(1-x)+6*(12/(48*2^3-4))+(1-(12/(48*2^3-4)))*(1-x)*(12*(1-y)+1))-z=0;
         (13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)=x;
         (1-(1-z)^8)/(2-(1-z)^9)=y;
注:其中x,y,z都是[0,1]数,一定需要用matlab解出来,感激不尽!
回复此楼

» 猜你喜欢

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

用心做好每一件事
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crazyiyy

木虫 (小有名气)

送红花一朵
引用回帖:
3楼: Originally posted by 月只蓝 at 2014-07-14 16:14:10
以下程序第对2楼中程序的改进,提高解的小数点后的位数。


function jie_fun
clear a;clc
format long
u0=;
u=fsolve(@myfun,u0)
f_val=myfun(u)



function f=myfun(u)
x=u(1);
y=u(2);
z=u(3);
...

您好! 您的答案的确跟1stOpt一致,还有个问题希望您能为我解惑。如果我需要对方程组中的某些系数进行变化,进而得到一组解,我应该如何操作。
例如 我需要对方程1中的48*2^3 变为48*2^x  其中x=1,2,3,4 这样我应该得到4组解,
用心做好每一件事
4楼2014-07-14 16:39:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crazyiyy

木虫 (小有名气)

送红花一朵
引用回帖:
5楼: Originally posted by 月只蓝 at 2014-07-14 16:45:00
我的处理办法是把x定义为全局变量。
以前已经有过相关代码了,参见:
http://muchong.com/bbs/viewthread.php?tid=6267100&authorid=1122189...

您好! 我还有最后一个问题,就是matlab种一定要给出初值, 请问初值的调试有没有什么方法或者经验啊,或者说如何判断当前初值得到的解是正解。
用心做好每一件事
7楼2014-07-14 17:06:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 crazyiyy 的主题更新
信息提示
请填处理意见