下载了《精通MATLAB最优化计算》书上的优化程序，用它来做杆机构的优化，我的目标函..

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
9 e' `/ f* R+ ]) P5 m) B惩罚函数调用格式如下：
" F: ?- @4 g2 }  g6 X% xfunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
- \, q( x% t+ t4 \. q. @%目标函数：f
0 M% C7 e) s( s- R' k%不等式约束：g
%等式约束：h
%初始点：x0
3 f! G( f- [& X: v; N%罚因子：r0
, w, q# ^9 [) s. u%缩小系数：c
%自变量向量：var
* R6 L8 h- l. t$ O) i! E% G%精度：eps
; t, |8 P: f- w1 a2 }* ~, K%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf
. V8 ?* A/ ]/ B: x% ?6 d) G

) ]) c& S! E' o5 _& x目标函数：
function f=link_objfun1(x)
global L1 L5
f=0;
- Z8 u0 C' ~9 {snti0=0;
for i=0:5
& D! c" e+ B0 S! Ksnti=snti0+pi/6*i/5;
fai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
* Z5 l: j" _/ i( S4 L0 k6 L8 |5 `psai=acos(2*(L1^2+x(1)^2+x(3)^2-x(2)^2-2*L1*x(1)*sin(snti+x(5)))/(x(3)*(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2));
2 t6 Q( K2 }) _4 S+ t4 h1 Y3 S( W; `$ {" [if snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
% a4 V& h* s8 s, @" H9 Y   bati=-fai+psai+acos(x(4)/x(3));
4 X. E; R$ E( iend
  r7 [9 y& ^% h6 b9 w$ U% |3 Wf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end


$ Z: m% w8 t" ?/ V9 J请诸位帮帮忙，我的体力不多。

厚积薄发

没看出来您的问题是什么？

changfeng1937

我懂的不多，生搬硬套。书上的例题
4 }. d& ?* b! h& a. m, p$ c运行如下：
syms t;
5 [2 A- `. b8 n+ za=4;b=3;
f=a*t;
6 V+ @7 b( S3 l* P/ B5 [g=[t-b];
0 L3 c  W9 `$ w" x+ z! c* }  T[x,minf]=minNF(f,[5],g,10,0.5,[t])
+ k9 l% \7 j8 h7 L9 L' ?) D我将目标函数写成：
0 `% ?8 T7 }' A5 M7 Wfunction f=link_objfun(s,t,a,m,n)
global L1 L5
f=0;
9 s. ^: r3 P/ p+ X  R3 csnti0=0;
0 p; i4 W, D2 M4 W4 R; I, rfor i=0:5
2 p3 [4 B2 |" v( c( G: m) bsnti=snti0+pi/6*i/5;
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
( X& J" Q2 r6 C7 O6 Npsai=acos(2*(L1^2+s^2+a^2-t^2-2*L1*s*sin(snti+n))/(a*(L1^2+s^2-2*L1*s*sin(snti+n))^1/2));
if snti<=pi/2-n
  b# ]) d7 i/ C! j5 k   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
! {8 P  X$ R8 a( n6 c   bati=-fai+psai+acos(m/a);
end
0 [- c" c; s& a" Z( Gf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
; C3 V" T8 e* t3 S& ~! O+ [end
3 V. e. Z  h: _7 y9 B- L运行如下：
( j3 G. _+ J' V( Q$ v' zclear all;
global L1 L5
L1=3000;
L5=800;
syms s t a m n;
f=link_objfun;
g=[s;t;a;m;n;L1-s;t-s;a-s;a-s-t+L1;s+t+a-L1;s+t+a-m];
2 w$ n0 K. p6 ]' I& ~5 E: ~[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
0 I6 I$ S; J8 E5 D/ e??? Input argument "s" is undefined.
, w- |% }) `3 ~/ `0 W! v+ u' B4 L
& R, \2 X5 ?, [  ]; V( e, w  jError in ==> link_objfun at 7
$ }/ d4 m# Y9 j* m) q% Kfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);