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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
' C; c; _& F% I& n+ y2 K惩罚函数调用格式如下：
6 m: P! z+ z2 r- Cfunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
. A: O6 {# G5 H, ?  `7 {5 T0 k* J%不等式约束：g
%等式约束：h
%初始点：x0
% U- _. F$ D- h' b%罚因子：r0
%缩小系数：c
( V( `$ P9 V* y# t: @* ^9 R1 q8 f%自变量向量：var
7 _3 g) D6 G9 z1 n' c9 U%精度：eps
%目标函数取最小值时的自变量的值：x
' v1 F5 ~( d% ?9 B) F  Y* g%目标函数的最小值：minf

  g+ t5 d) r. V- z
目标函数：
function f=link_objfun1(x)
global L1 L5
f=0;
snti0=0;
) D. r# z; s9 y9 B, [for i=0:5
snti=snti0+pi/6*i/5;
% v  p7 M+ }) x  w: [4 kfai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
8 p5 w$ x4 T1 S( F8 L/ v9 t% d+ qpsai=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));
% E& N' j! x: v$ R- n7 \% @if snti<=pi/2-x(5)
# k! t- D  l" {1 ^( ]   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
: N8 E  h4 A. M4 Y4 Iend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end


/ c: w& ~1 j' f. E# K3 U0 p请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
syms t;
' L- j4 x" c# W4 f  A, u5 b: X8 H; ca=4;b=3;
f=a*t;
g=[t-b];
9 g& j# I, O' ~% `1 `[x,minf]=minNF(f,[5],g,10,0.5,[t])
: g) z) _. G4 h. G我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
f=0;
snti0=0;
7 Q$ ?9 O+ o: p0 ?. M5 ~. Y+ A' I6 [for i=0:5
snti=snti0+pi/6*i/5;
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));
if snti<=pi/2-n
   bati=fai+psai+acos(m/a);
( c, j4 ?6 L- J* l* U0 ^/ Q" qelseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
7 ]- S: ~/ @/ ]( Q3 aend
! ]: g4 \; J0 w0 jf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
# A3 U4 G; l3 w4 e" |& ]运行如下：
, c2 d' |7 L% q" H3 mclear all;
global L1 L5
L1=3000;
L5=800;
syms s t a m n;
8 p% C9 l- f0 g' ?- Jf=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];
1 u, Z5 M2 u! V5 m' o$ F[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
??? Input argument "s" is undefined.

Error in ==> link_objfun at 7
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);