function jie_fun
u0=[0.5 0.5 0.5];
u=fsolve(@myfun,u0)
f_val=myfun(u)



function f=myfun(u)
x=u(1);
y=u(2);
z=u(3);
f(1)=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;
f(2)= (13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)-x;
f(3)=(1-(1-z)^8)/(2-(1-z)^9)-y;

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



function f=myfun(u)
x=u(1);
y=u(2);
z=u(3);
f(1)=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;
f(2)= (13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)-x;
f(3)=(1-(1-z)^8)/(2-(1-z)^9)-y;

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);
...

4楼: Originally posted by crazyiyy at 2014-07-14 16:39:24
您好！ 您的答案的确跟1stOpt一致，还有个问题希望您能为我解惑。如果我需要对方程组中的某些系数进行变化，进而得到一组解，我应该如何操作。
例如 我需要对方程1中的48*2^3 变为48*2^x  其中x=1,2,3,4 这样我应 ...

4楼: Originally posted by crazyiyy at 2014-07-14 16:39:24
您好！ 您的答案的确跟1stOpt一致，还有个问题希望您能为我解惑。如果我需要对方程组中的某些系数进行变化，进而得到一组解，我应该如何操作。
例如 我需要对方程1中的48*2^3 变为48*2^x  其中x=1,2,3,4 这样我应 ...

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

7楼: Originally posted by crazyiyy at 2014-07-14 17:06:19
您好！ 我还有最后一个问题，就是matlab种一定要给出初值， 请问初值的调试有没有什么方法或者经验啊，或者说如何判断当前初值得到的解是正解。...

8楼: Originally posted by 月只蓝 at 2014-07-14 18:37:07
把算出来的解代回原方程，看能有多接近于0，3楼正是给出了这些数值。
此外，对于多解的情况，最终解的如果有物理意义要求，举例说，比如未知数是开尔文温度T，是不可能小于0的，在解里面就应当取满足实际物理意义 ...

9楼: Originally posted by mjjackey at 2015-03-02 16:40:35
用fsolve，貌似也只能设置三种算法，没办法设置更多的算法“调用具有全局优化能力的算法，比如说遗传算法，”...