| 查看: 352 | 回复: 0 | ||
[求助]
函数最值问题！！！！
|
|
小虫，刚接触MATLAB,需求一函数的最值，但约束条件有些复杂，顿时凌乱，望高手指导！ 以下是我的文件内容： function f = fun(x); x0=[1 1 1 1 10*pi/10 10*pi/180 10*pi/180 4*pi/180 8*pi/180 80*pi/180 175*pi/180]; T=1300; M=11500; R=10; U=0.16; fun=@(x)1/1.6900e-007*abs(abs(((T/R)*(x(1)*tan(x(5))-x(2)*tan(x(6))+x(3)*tan(x(7))*cos(x(9))-x(4)*tan(x(8)))-(2*M/R^2)*x(3)*sin(x(9)))*(sin(x(11))... -U*cos(x(11)))+((2*M/R^2)*(x(1)-x(2)+cos(x(9))-x(4))+(T/R)*x(3)*tan(x(7))*sin(x(9)))*(cos(x(10))+U*sin(x(11)))./(1.0256*sin(x(11)-x(10))))... -abs(((2*M/R^2)*(x(1)-x(2)+cos(x(9))-x(4))+(T/R)*x(3)*tan(x(7))*sin(x(9)))))+1/2.0331e-006*(abs(abs(((T/R)*(x(1)*tan(x(5))-x(2)*tan(x(6))+x(3)*tan(x(7))*cos(x(9))-x(4)*tan(x(8)))-(2*M/R^2)*x(3)*sin(x(9)))*(sin(x(10))... -U*cos(x(10)))+((2*M/R^2)*(x(1)-x(2)+cos(x(9))-x(4))+(T/R)*x(3)*tan(x(7))*sin(x(9)))*(cos(x(10))+U*sin(x(10)))./((1+U^2)*sin(x(10)-x(11))))... -abs((T/R)*(x(1)*x(5)-x(2)*x(6)+x(3)*tan(x(7))*cos(x(9))-x(4)*tan(x(8)))-(2*M/R^2)*x(3)*sin(x(9))))) A=[]; B=[]; Aeq=[1 1 1 1 0 0 0 0 0 0 0]; Beq=[10]; lb=[1,1,1,1,5*pi/180,8*pi/180,8*pi/180,10*pi/180,8*pi/180,85*pi/180,185*pi/180]; ub=[5,5,5,5,16*pi/180,14*pi/180,14*pi/180,14*pi/180,12*pi/180,95*pi/180,200*pi/180]; [x,fval]= fmincon(fun,x0,A,B,Aeq,Beq,lb,ub) c=(T/R)*(x(1)*tan(x(5))-x(2)*tan(x(6))+x(3)*tan(x(7))*cos(x(9))-x(4)*tan(x(8)))-(2*M/R^2)*x(3)*sin(x(9)) d=(2*M/R^2)*(x(1)-x(2)+cos(x(9))-x(4))+(T/R)*x(3)*tan(x(7))*sin(x(9)) p=((T/R)*(x(1)*tan(x(5))-x(2)*tan(x(6))+x(3)*tan(x(7))*cos(x(9))-x(4)*tan(x(8)))-(2*M/R^2)*x(3)*sin(x(9)))*(sin(x(11))... -U*cos(x(11)))+((2*M/R^2)*(x(1)-x(2)+cos(x(9))-x(4))+(T/R)*x(3)*tan(x(7))*sin(x(9)))*(cos(x(10))+U*sin(x(11)))./((1+U^2)*sin(x(11)-x(10))) w=((T/R)*(x(1)*tan(x(5))-x(2)*tan(x(6))+x(3)*tan(x(7))*cos(x(9))-x(4)*tan(x(8)))-(2*M/R^2)*x(3)*sin(x(9)))*(sin(x(10))-U*cos(x(10)))+((2*M/R^2)*(x(1)-x(2)+cos(x(9))-x(4))+(T/R)*x(3)*tan(x(7))*sin(x(9)))*(cos(x(10))+U*sin(x(10)))./((1+U^2)*sin(x(10)-x(11))) 其中要使c,d大于0;p,w小于0，不知道如何添加约束，还望高手指导》》》 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : laomo.m
2012-05-28 10:09:00, 2.61 K
» 猜你喜欢
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
听说今天filecode变了
已经有45人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
-
filecode
已经有19人回复
-
基金中了
已经有11人回复
-
固定端突然变了，今天
已经有8人回复
找到一些相关的精华帖子，希望有用哦~
毕设，数学相关专业，求指导！
已经有18人回复
- [求助]高等数学--最值问题
已经有7人回复
- 【求助】在圆柱内分别堆积立方体，八面体，球形的问题
已经有5人回复
- 【求助】函数最小值问题？
已经有9人回复
- 【求助】请高人帮忙三元二次方程
已经有11人回复
- 【求助】计算两段空间圆弧的最短距离---(已解决)
已经有11人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
100