在fmincon中出现错误

pTDbn25

在fmincon中出现错误
' |! u# j8 m3 U完整代码：clear all;
: I1 C; ^8 j8 p3 S" C! |close all;
clc;
x0=[480;100;120];
6 h; v5 }. R- A5 P$ mlb=[300;60;90];
/ f( V& z0 S. Z( Yub=[650;140;150];
4 U# v& Y5 H: v# u' p3 j0 K4 J" P[x,fval]=fmincon('fun',x0,[],[],[],[],lb,ub,'con1');

* s) R! g& X- E# q' G7 l! W7 f+ A, o8 [错误代码：

错误使用 fmincon (line 619)
/ Z: D" A+ T0 k7 a2 A7 OSupplied objective function must return a scalar value.
) K' [  `6 [3 |
出错 Untitled (line 7)
8 q) P3 a* V7 [3 e* S; X[x,fval]=fmincon('fun',x0,[],[],[],[],lb,ub,'con1');

9 D  X" D$ Q4 p4 z; n  G$ S

nkkopd

来学习一下

kekek

clear all;
close all;
/ p" i* m2 M  y9 j$ e% Q3 o) J+ }clc;
x0=[480;100;120];
5 z1 {/ Z/ t' Q; \lb=[300;60;90];
ub=[650;140;150];
4 x! W0 E3 R& z  N5 _7 n" }, R[x,fval]=fmincon(@fun,x0,[],[],[],[],lb,ub,@con1);
) {( ~7 v* Z5 \2 K5 q2 rfunction f=fun(x)
. h5 M  N; c/ q: S+ Wf=0.25*pi*7.8*10^(-6)*(x(1)+x(3))*(x(2)^2-30^2);
1 L1 z/ {3 ^2 ]6 p- R9 a: [2 u5 fend
1 X8 F# R" s5 p) L8 _function [c,ceq]=con1(x)
$ ^+ G. \+ H8 W' B4 {! cc=(64*15000*x(3)^2*(x(1)+x(3))/(3*pi*2.1*10^5*x(2)^4-30^4)-0.05);
ceq=[];
end

pTDbn25

function f=fun(x)
$ ?4 U5 ^1 D0 |: k" m/ s5 q$ \/ |f(x)=0.25*pi*7.8*10^(-6)*(x(1)+x(3))*(x(2)^2-30^2);
5 s( x; C* N0 x9 z
function [c,ceq]=con1(x)
( Y  i6 {: i: T# Fc=(64*15000*x(3)^2*(x(1)+x(3))/(3*pi*2.1*10^5*x(2)^4-30^4)-0.05);
ceq=[];