找回密码
 注册
Simdroid-非首页
查看: 65|回复: 0

求助matlab解超越方程组

[复制链接]
发表于 2012-4-17 21:02:37 | 显示全部楼层 |阅读模式 来自 福建厦门
以下是超越方程组,其中x1min和x2min是要求解的变量,Rou为已知量,其余的量如Proz4、a12、Matrix1等都是循环一次就要变一次的量,每个循环中Pro2从0~360变化,求出每个Pro值对应的x1min和x2min:
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还有其他解超越方程组的函数吗?
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-10-1 21:36 , Processed in 0.030565 second(s), 12 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表