24小时热门版块排行榜    

查看: 341  |  回复: 0

dage16wo

木虫 (小有名气)

[求助] fmincon优化的问题

function f=objtctfun(x)
global G1 L1 G4 L4 Gg Gd Gw G05 L05 G04 L04 G03 L03 G02 L02 G01 L01 k l a c
f=(G1*L1+G4*L4+(Gg+Gd)*L1+Gw*L4+G05*L05+G04*L04+G03*L03+G02*L02+G01*L01)*cos(k*pi/180)*(l+x(3))/(sqrt(a^2*c^2+x(2)^2*x(1)^2-2*a*c*x(1)*x(2)));

以上为目标函数



function [c,ceq]=conobject(x)
global a c1 k l h1 k1 k2 pi
c=[(l+x(3))-sqrt(a^2+x(1)^2)-sqrt(x(2)^2+c1^2);
sqrt(x(2)^2+c1^2)-sqrt(a^2+x(1)^2)-(l+x(3));
sqrt(x(2)^2+c1^2)-sqrt(a^2+x(1)^2)-(l+x(3));
(l+2*x(3))-(sqrt(a^2+x(1)^2)+sqrt(x(2)^2+c1^2));
(sqrt(x(2)^2+c1^2))-(sqrt(a^2+x(1)^2)+(l+2*x(3)));
sqrt(x(2)^2+c1^2)-(sqrt(a^2+x(1)^2)+(l+2*x(3)));
x(1)-h1;
c1-x(1);
x(2)-k2;
k1-x(2);
a^2+x(1)^2+c1^2+x(2)^2-2*(a*x(2)+x(1)*c1)*cos(k*pi/180)+2*sqrt(a^2*c1^2+x(2)^2*x(1)^2-2*a*c1*x(1)*x(2))*sin(k*pi/180) -(l+2*x(3))^2;
];
ceq=a^2+x(1)^2+c1^2+x(2)^2-2*(a*x(2)+x(1)*c1)-(l+x(3))^2;



以上为约束函数



G1=410  L1=21550  G4=200   L4=20730   Gw=49.79    Gg=199.16   Gd=100  G05=0    L05=0    G04=0  L04=0  G03=830.72  

L03=17707   G02=1214.5  L02=11517   G01= 1659.95   L01=4445

a=100   c1=390   k=30   l=500   h1=4000   k1=1000  k2=2500  

初始值x0=[500:500:500]  非线性优化为题


出现问题

??? Error using ==> vertcat
CAT arguments dimensions are not consistent.

Error in ==> conobject at 3
c=[(l+x(3))-sqrt(a^2+x(1)^2)-sqrt(x(2)^2+c1^2);

Error in ==> fmincon at 654
        [ctmp,ceqtmp] = feval(confcn{3},X,varargin{:});

Caused by:
    Failure in initial user-supplied nonlinear constraint function evaluation. FMINCON cannot continue.
回复此楼
loveisalwayshere!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dage16wo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 11:18 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 11:09 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 6/300 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 14/700 2026-11-01 10:43 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文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 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,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 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 起风了的快乐
信息提示
请填处理意见