应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 函数最值问题！！！！
11/1返回列表
查看: 371  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

myl5865115

新虫 (初入文坛)

[求助] 函数最值问题！！！！

小虫，刚接触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

» 猜你喜欢
1楼 2012-05-28 10:10:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 myl5865115 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 12:19 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:56 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 10:48 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 10:27 by BCdvBAEpt3zm
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 10:27 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 04:04 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 3/150 2026-11-01 04:04 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
关闭
请填处理意见
关闭