《量化投资:以MATLAB为工具》

MATLAB技术论坛

 找回密码
 注册账号
查看: 139|回复: 0
收起左侧

[提问] 谁有高版本的1stopt帮我解个超定方程

[复制链接]

该用户从未签到

发表于 2017-8-12 09:45:58 | 显示全部楼层 |阅读模式
我想用高版本的1stopt解超定方程,12个方程,6个未知数,有人可以帮我解吗??1.5版本的只能显示一个结果,我想看到所有的结果。谢谢了!!

Parameter  theta1[-3.2289,3.2289],theta2[-2.7053,0.6109],theta3[-2.2689,2.6878],theta4[-6.1087,6.1087],theta5[-2.2689,2.2689],theta6[-6.1087,6.1087];
function     cos(theta2 - theta3)*cos(theta1)*cos(theta5) - sin(theta1)*sin(theta4)*sin(theta5) - cos(theta1)*cos(theta2)*cos(theta4)*sin(theta3)*sin(theta5) + cos(theta1)*cos(theta3)*cos(theta4)*sin(theta2)*sin(theta5)=0;
            sin(theta6)*(cos(theta5)*sin(theta1)*sin(theta4) + cos(theta2 - theta3)*cos(theta1)*sin(theta5) + cos(theta1)*cos(theta2)*cos(theta4)*cos(theta5)*sin(theta3) - cos(theta1)*cos(theta3)*cos(theta4)*cos(theta5)*sin(theta2)) - cos(theta6)*(cos(theta4)*sin(theta1) - cos(theta1)*cos(theta2)*sin(theta3)*sin(theta4) + cos(theta1)*cos(theta3)*sin(theta2)*sin(theta4))=1;
             cos(theta6)*(cos(theta5)*sin(theta1)*sin(theta4) + cos(theta2 - theta3)*cos(theta1)*sin(theta5) + cos(theta1)*cos(theta2)*cos(theta4)*cos(theta5)*sin(theta3) - cos(theta1)*cos(theta3)*cos(theta4)*cos(theta5)*sin(theta2)) + sin(theta6)*(cos(theta4)*sin(theta1) - cos(theta1)*cos(theta2)*sin(theta3)*sin(theta4) + cos(theta1)*cos(theta3)*sin(theta2)*sin(theta4))=0;
             (cos(theta1)*(136*cos(theta2 - theta3) + 240*sin(theta2 - theta3) - 129*sin(theta2) + 52))/200=-0.295;
              cos(theta1)*sin(theta4)*sin(theta5) + cos(theta2 - theta3)*cos(theta5)*sin(theta1) - cos(theta2)*cos(theta4)*sin(theta1)*sin(theta3)*sin(theta5) + cos(theta3)*cos(theta4)*sin(theta1)*sin(theta2)*sin(theta5)=0;
            sin(theta6)*(cos(theta2 - theta3)*sin(theta1)*sin(theta5) - cos(theta1)*cos(theta5)*sin(theta4) + cos(theta2)*cos(theta4)*cos(theta5)*sin(theta1)*sin(theta3) - cos(theta3)*cos(theta4)*cos(theta5)*sin(theta1)*sin(theta2)) + cos(theta6)*(cos(theta1)*cos(theta4) + cos(theta2)*sin(theta1)*sin(theta3)*sin(theta4) - cos(theta3)*sin(theta1)*sin(theta2)*sin(theta4))=0;
             cos(theta6)*(cos(theta2 - theta3)*sin(theta1)*sin(theta5) - cos(theta1)*cos(theta5)*sin(theta4) + cos(theta2)*cos(theta4)*cos(theta5)*sin(theta1)*sin(theta3) - cos(theta3)*cos(theta4)*cos(theta5)*sin(theta1)*sin(theta2)) - sin(theta6)*(cos(theta1)*cos(theta4) + cos(theta2)*sin(theta1)*sin(theta3)*sin(theta4) - cos(theta3)*sin(theta1)*sin(theta2)*sin(theta4))=-1;
             (sin(theta1)*(136*cos(theta2 - theta3) + 240*sin(theta2 - theta3) - 129*sin(theta2) + 52))/200=0;
              sin(theta2 - theta3)*cos(theta5) - cos(theta2 - theta3)*cos(theta4)*sin(theta5)=-1;
              cos(theta2 - theta3)*cos(theta6)*sin(theta4) + sin(theta2 - theta3)*sin(theta5)*sin(theta6) + cos(theta2 - theta3)*cos(theta4)*cos(theta5)*sin(theta6)=0;
               sin(theta2 - theta3)*cos(theta6)*sin(theta5) - cos(theta2 - theta3)*sin(theta4)*sin(theta6) + cos(theta2 - theta3)*cos(theta4)*cos(theta5)*cos(theta6)=0;
             (17*sin(theta2 - theta3))/25 - (6*cos(theta2 - theta3))/5 + (129*cos(theta2))/200 + 27/40=-0.005;
您需要登录后才可以回帖 登录 | 注册账号

本版积分规则

QQ|网站地图|MATLAB技术论坛|Simulink仿真论坛 ( 陕ICP备08102094号  

GMT+8, 2017-10-20 07:35 , Processed in 0.087653 second(s), 26 queries , Gzip On.

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

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