9楼: Originally posted by duyoudu at 2013-06-17 10:40:40
应该是没问题的。fsolve的算法是个通用算法，可以解答线性和非线性的，虽然效率不如专门的求解函数高。但这个局部最优的问题是存在的，可以考虑函数工具箱中提供的全局求解的函数，具体是哪个你可以看看帮助。...

9楼: Originally posted by duyoudu at 2013-06-17 10:40:40
应该是没问题的。fsolve的算法是个通用算法，可以解答线性和非线性的，虽然效率不如专门的求解函数高。但这个局部最优的问题是存在的，可以考虑函数工具箱中提供的全局求解的函数，具体是哪个你可以看看帮助。...

Parameter x(6);
Constant d=0.045,D1=0.05,k=1.4,ff=0.004,L=100,P1=1500000,T1=300,Pa=101325,M=29,R=8.314;
Function ((k+1)/2)*ln((x2^2*x3)/(x1^2*x4))-(1/(x1^2)-1/(x2^2))+4*ff*L/D1;
        (D1^2/D1^2)*x5*sqrt((2*k*M)*((Pa/x5)^(2/k)-(Pa/x5)^((k+1)/k))/(R*x6*(k-1)))-x2*x5*sqrt((k*M)/(R*x6));
        (D1^2/D1^2)*x5*sqrt((2*k*M)*((Pa/x5)^(2/k)-(Pa/x5)^((k+1)/k))/(R*x6*(k-1)))-x1*P1*sqrt((k*M)/(R*T1));
        1+((k-1)*x1^2)/2-x3;
        1+((k-1)*x2^2)/2-x4;
        x6-(x3*T1)/x4;

13楼: Originally posted by dingd at 2013-06-17 14:25:43
用1stOpt吧，其它软件得到正解的概率几乎为0：

Parameter x(6);
Constant d=0.045,D1=0.05,k=1.4,ff=0.004,L=100,P1=1500000,T1=300,Pa=101325,M=29,R=8.314;
Function ((k+1)/2)*ln((x2^2*x3)/(x1^2*x4))-(1/ ...

11楼: Originally posted by xiangzigood at 2013-06-17 11:22:57
全局求解？您能给推荐一个不，刚解除matlab不久。。...