回 帖 发 新 帖 刷新版面

主题:[讨论]求助matlab解超越方程组问题

以下是非线性方程组,其中x1min和x2min是要求解的变量,Rou为已知量,其余的量如Proz4、a12、Matrix1等都是解一次方程就要变一次的量:

eq1=((cos(x2min)+1/2*sin(Proz4+x1min)^2*(1-cos(x2min)))*cos(x1min)+(-1/2*sin(Proz4+x1min)*cos(Proz4+x1min)*(1-cos(x2min))-1/2*2^(1/2)*sin(x2min))*sin(x1min))*Matrix1+(-(cos(x2min)+1/2*sin(Proz4+x1min)^2*(1-cos(x2min)))*sin(x1min)+(-1/2*sin(Proz4+x1min)*cos(Proz4+x1min)*(1-cos(x2min))-1/2*2^(1/2)*sin(x2min))*cos(x1min))*Matrix2+(1/2*sin(Proz4+x1min)*(1-cos(x2min))-1/2*2^(1/2)*cos(Proz4+x1min)*sin(x2min))*Matrix3-sin(Rou)*sin(Pro2)*cos(a12)-sin(Rou)*cos(Pro2)*cos(a22)+cos(Rou)*cos(a32);

eq2=((-1/2*sin(Proz4+x1min)*cos(Proz4+x1min)*(1-cos(x2min))+1/2*2^(1/2)*sin(x2min))*cos(x1min)+(cos(x2min)+1/2*cos(Proz4+x1min)^2*(1-cos(x2min)))*sin(x1min))*Matrix1+(-(-1/2*sin(Proz4+x1min)*cos(Proz4+x1min)*(1-cos(x2min))+1/2*2^(1/2)*sin(x2min))*sin(x1min)+(cos(x2min)+1/2*cos(Proz4+x1min)^2*(1-cos(x2min)))*cos(x1min))*Matrix2+(-1/2*cos(Proz4+x1min)*(1-cos(x2min))-1/2*2^(1/2)*sin(Proz4+x1min)*sin(x2min))*Matrix3-sin(Rou)*sin(Pro2)*cos(b12)-sin(Rou)*cos(Pro2)*cos(b22)+cos(Rou)*cos(b32);

之前有用fsolve求解过,但是得出来的解反带回原方程组只有前几个点是正确的,不知道是否与初值的设置有关?如果是,要怎么解决?matlab里面还有其他解超越方程组的函数吗?请高手指导,谢谢啦!

回复列表 (共1个回复)

沙发

http://www.zhonghuabz168.com本文来自:http://www.dream8.com.cn/article.php?pd=762&id=21329

我来回复

您尚未登录,请登录后再回复。点此登录或注册