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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
, Z2 x  W/ y+ @: K% Z5 _, N3 O' i惩罚函数调用格式如下：
5 w1 h: z+ g- T, \function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
/ n; R( X$ G% @& H# Y* J%不等式约束：g
$ M' F  k: b" G: U; Y" Q- }7 R7 c%等式约束：h
%初始点：x0
%罚因子：r0
) \# h# s" \/ ~- d1 R%缩小系数：c
) C7 u& J5 x6 g; ~  ?%自变量向量：var
! m2 ]2 u6 X( _/ I%精度：eps
%目标函数取最小值时的自变量的值：x
& @% y# Y/ D# q5 S%目标函数的最小值：minf

8 ^; k8 [0 j% ^5 A8 @& e6 H% L
目标函数：
function f=link_objfun1(x)
global L1 L5
f=0;
snti0=0;
" F$ j7 ?( M. s) q; @4 L- dfor i=0:5
snti=snti0+pi/6*i/5;
) R* D) m, \( n5 x# s( s% ^fai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
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));
7 n. n9 {; `! Z  W1 G1 Hif snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
5 `# a/ V- D# pelseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
+ \, }( ]( t9 |' G/ y! @$ kend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
7 N: t: k; ~4 L! hend
/ ~6 R! Z: X! w5 P7 u6 s  v
7 T3 v- Z, Y' j- m) O; Z: t
/ D+ k1 w7 N, N2 q: O( ~+ y请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
8 {0 e$ V) T; l5 M+ w' V2 h运行如下：
syms t;
6 k1 D; X8 ~! q  R* i9 K, Ea=4;b=3;
' B8 D2 R# {2 }7 V& l( f5 Ff=a*t;
g=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
7 x& ]* f* B  _1 ^$ b0 f我将目标函数写成：
9 x1 z8 p+ o( Y+ k2 N; Ofunction f=link_objfun(s,t,a,m,n)
global L1 L5
f=0;
: L; X. M" B  ~# D" K" s. `snti0=0;
for i=0:5
9 x) r; m  P# Q" |4 R% [snti=snti0+pi/6*i/5;
2 i4 {- \4 E5 `( v8 B4 ^fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
psai=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));
7 r8 S: D. _; X. b6 u6 Uif snti<=pi/2-n
# ?9 W+ J  a! q   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
; q2 m: C8 Q% g   bati=-fai+psai+acos(m/a);
) o- |, X3 V4 h6 e% t5 a" |end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
运行如下：
clear all;
global L1 L5
$ q0 _4 {" j8 A$ M/ h( T' cL1=3000;
  I. _# m* S2 Y) P% \* r$ B: P2 mL5=800;
syms s t a m n;
f=link_objfun;
6 O+ F8 ?5 k: O. t8 t3 @g=[s;t;a;m;n;L1-s;t-s;a-s;a-s-t+L1;s+t+a-L1;s+t+a-m];
[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
??? Input argument "s" is undefined.
( d* I. k1 H  Q. S: ^
( y, q2 H& z8 b! l  h& ]! C7 r/ {Error in ==> link_objfun at 7
1 x7 s# @. C- ]fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);