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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
# D6 v$ V2 I3 l, V9 H8 o%不等式约束：g
4 a8 J; }% w- e1 k9 W%等式约束：h
3 c; j% J0 d; M  [%初始点：x0
2 R$ Z: V$ m, r%罚因子：r0
, y4 q* n0 Y9 N; R7 F%缩小系数：c
%自变量向量：var
3 [( F0 Q: |1 r( `%精度：eps
%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf

# h4 e8 d6 k, v
/ T4 y- F+ Q9 F' O4 x( J6 {$ p* P# x目标函数：
) ?2 O0 `' A) `9 y  y; ]& Sfunction f=link_objfun1(x)
' t9 [% j8 O; R/ \& I0 xglobal L1 L5
f=0;
snti0=0;
1 D! ^/ Q& |( D, l  qfor i=0:5
& p( }) U( T& @snti=snti0+pi/6*i/5;
& _- b$ w& i% v% L1 V+ H& cfai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
! i: e/ G. a# Q! y8 [; _2 w" H0 I5 Dpsai=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));
if snti<=pi/2-x(5)
5 y6 o, x) |6 K6 c* Z6 b; w   bati=fai+psai+acos(x(4)/x(3));
" P; u1 B8 ?8 G9 {6 i- aelseif snti>pi/2-x(5)
' I: n' J2 N+ {/ D& P& Y. G& ^1 q   bati=-fai+psai+acos(x(4)/x(3));
end
6 z/ q$ N! \4 b' O  U: l2 if=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end


, M: B  D2 t: r3 v: ?: s3 m2 z. R请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
+ V# _7 v# H# Z8 D0 g, n! n1 Z运行如下：
/ r5 p, K3 f, _: s3 V9 D( `0 E; }) Rsyms t;
a=4;b=3;
f=a*t;
: N/ D/ c* ]" Z: Mg=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
4 ~" U2 n9 {2 @0 z2 A0 Bfunction f=link_objfun(s,t,a,m,n)
global L1 L5
+ s" ?0 v& w5 z' If=0;
) H% Z8 Z6 Y6 p, Ksnti0=0;
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);
" i& @* `" [4 j) W( P  Rpsai=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));
' h+ p4 z6 `' Y3 R% _9 dif snti<=pi/2-n
" C4 U/ Y. I- m0 u' ^  O. k6 }9 y   bati=fai+psai+acos(m/a);
5 g- i+ S  `1 X/ Y+ zelseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
. t1 r; t% F, E0 g) v1 H" F. ~end
运行如下：
clear all;
+ K: D- c$ D: X+ _, b! C  P  Tglobal L1 L5
4 t1 S5 U# K1 @L1=3000;
1 x+ P/ U; A) M7 V9 bL5=800;
7 f4 O( V% N+ v) Y& q; M- usyms s t a m n;
) p  P: U, W4 Yf=link_objfun;
4 I! O5 X  r4 }3 rg=[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.
4 V! s  l: z. e0 F# O0 K) q
8 t; ^* U$ J  BError in ==> link_objfun at 7
+ U7 e' Y+ u8 @3 ?6 I. wfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);