matlab 条件方程组的解,solve 时解方程组的限制条件问题

本帖最后由 oldlybaby 于 2017-5-28 14:43 编辑

简单来说，需要求解a1,a2,a3,但只有两个关于a1,a2,a3的方程f1,f2，附加条件是a1+a2+a3最小，请问怎么求解方程组，我的程序(方程有点长)如下

syms a1 a2 a3 ;复制代码

f1=cos(a3)*(10*sin(a1)*(cos(a2) - 1) + 10/((10*sin(a1)*sin(a2) - 10*cos(a1)*(cos(a2) - 1) + cos(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - sin(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)))^2/(10*sin(a1)*(cos(a2) - 1) - sin(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - cos(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + 10*cos(a1)*sin(a2))^2 + 1)^(1/2) - sin(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - cos(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + 10*cos(a1)*sin(a2)) - sin(a3)*(10*sin(a1)*sin(a2) - 10*cos(a1)*(cos(a2) - 1) + cos(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - sin(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + (10*(10*sin(a1)*sin(a2) - 10*cos(a1)*(cos(a2) - 1) + cos(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - sin(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2))))/(((10*sin(a1)*sin(a2) - 10*cos(a1)*(cos(a2) - 1) + cos(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - sin(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)))^2/(10*sin(a1)*(cos(a2) - 1) - sin(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - cos(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + 10*cos(a1)*sin(a2))^2 + 1)^(1/2)*(10*sin(a1)*(cos(a2) - 1) - sin(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - cos(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + 10*cos(a1)*sin(a2)))) - (cos(a3) - 1)*(10*sin(a1)*(cos(a2) - 1) - sin(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - cos(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + 10*cos(a1)*sin(a2)) + sin(a3)*(10*sin(a1)*sin(a2) - 10*cos(a1)*(cos(a2) - 1) + cos(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - sin(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)))-T(1);复制代码

f2=cos(a3)*(10*sin(a1)*sin(a2) - 10*cos(a1)*(cos(a2) - 1) + cos(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - sin(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + (10*(10*sin(a1)*sin(a2) - 10*cos(a1)*(cos(a2) - 1) + cos(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - sin(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2))))/(((10*sin(a1)*sin(a2) - 10*cos(a1)*(cos(a2) - 1) + cos(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - sin(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)))^2/(10*sin(a1)*(cos(a2) - 1) - sin(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - cos(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + 10*cos(a1)*sin(a2))^2 + 1)^(1/2)*(10*sin(a1)*(cos(a2) - 1) - sin(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - cos(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + 10*cos(a1)*sin(a2)))) - (cos(a3) - 1)*(10*sin(a1)*sin(a2) - 10*cos(a1)*(cos(a2) - 1) + cos(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - sin(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) + sin(a3)*(10*sin(a1)*(cos(a2) - 1) + 10/((10*sin(a1)*sin(a2) - 10*cos(a1)*(cos(a2) - 1) + cos(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - sin(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)))^2/(10*sin(a1)*(cos(a2) - 1) - sin(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - cos(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + 10*cos(a1)*sin(a2))^2 + 1)^(1/2) - sin(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - cos(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + 10*cos(a1)*sin(a2)) - sin(a3)*(10*sin(a1)*(cos(a2) - 1) - sin(a2)*(10*cos(a1) - (8*cos(a1))/(sin(a1)*(cos(a1)^2/sin(a1)^2 + 1)^(1/2))) - cos(a2)*(10*sin(a1) - 8/(cos(a1)^2/sin(a1)^2 + 1)^(1/2)) + 10*cos(a1)*sin(a2))-T(2);复制代码

[a1, a2,a3]=solve(f1,f2);复制代码

QQ截图20170528144323.png (18.77 KB, 下载次数: 0)

2017-5-28 14:37 上传