EDA365电子论坛网

标题: 在fmincon中出现错误 [打印本页]

作者: pTDbn25    时间: 2020-12-10 10:17
标题: 在fmincon中出现错误
在fmincon中出现错误0 L( ~9 {: [: p! i$ A! C
完整代码:clear all;6 |* X: |. B0 \7 U( `
close all;
* t* j! T4 x5 m0 C8 x: ?clc;
" e' \' f, V" a$ ^; ^- G5 T8 z+ f/ Jx0=[480;100;120];
$ m  h6 n% T3 C' Q# c( vlb=[300;60;90];
# W( p1 p0 m- _* ^- \4 ~ub=[650;140;150];
8 ^" Q' @8 {) Y$ g* ][x,fval]=fmincon('fun',x0,[],[],[],[],lb,ub,'con1');
. g9 N9 @: c$ P4 B; {+ |) L! i/ i; m
错误代码:+ ?& \5 j1 G3 M& B. Z0 K# v

& R: |1 e# [7 D错误使用 fmincon (line 619); e$ \  j1 D% M, C
Supplied objective function must return a scalar value.3 u- }3 A4 \# Y3 j, `$ G
% l! k0 s8 B9 y$ R
出错 Untitled (line 7)% Q% b0 v) \& w1 u* X
[x,fval]=fmincon('fun',x0,[],[],[],[],lb,ub,'con1');
9 q7 {! R+ T4 V, x  V, ~' R* y7 n8 A. [

8 m' ]# v, U* ?. ^% x' ~$ }1 |
作者: pTDbn25    时间: 2020-12-10 10:54
function f=fun(x)
( [, c: Y, T+ b  Bf(x)=0.25*pi*7.8*10^(-6)*(x(1)+x(3))*(x(2)^2-30^2);
  Z! N& G8 G4 V' l; Y
2 m; U! x7 k) cfunction [c,ceq]=con1(x)3 f1 v* y5 N3 F
c=(64*15000*x(3)^2*(x(1)+x(3))/(3*pi*2.1*10^5*x(2)^4-30^4)-0.05);
- k/ N4 [/ T, Y! bceq=[];
作者: kekek    时间: 2020-12-10 14:31
clear all;! W; Q! c4 s& p. p. N1 g
close all;
! M2 i& u! X  O/ P8 Fclc;& z8 J+ |. p% l9 M8 l5 m
x0=[480;100;120];; \& D7 ?' W- ^2 |9 F
lb=[300;60;90];
. R- r6 m) h$ r3 R, Rub=[650;140;150];
' H/ \; z( M% r" T9 c. V  X[x,fval]=fmincon(@fun,x0,[],[],[],[],lb,ub,@con1);0 j3 e3 n. N/ {; C8 G
function f=fun(x)0 a; z1 X, `. s6 |$ T8 r2 j9 n" @, G
f=0.25*pi*7.8*10^(-6)*(x(1)+x(3))*(x(2)^2-30^2);* j0 X3 A% c. \' g
end' J( `* @' C' S+ e/ q
function [c,ceq]=con1(x)4 H* q' p  e$ N- z3 c" d
c=(64*15000*x(3)^2*(x(1)+x(3))/(3*pi*2.1*10^5*x(2)^4-30^4)-0.05);
: Y- Y# H. `# Z- q3 Z: y. @/ ~! Eceq=[];- |$ b, J" t1 k
end
作者: nkkopd    时间: 2020-12-10 14:37
来学习一下




欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2