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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
, w* v/ @# F7 M! @  R- o; D%目标函数：f
0 v/ E: I4 v5 F1 S# L%不等式约束：g
%等式约束：h
%初始点：x0
0 }& _: \* M2 N, b7 Y%罚因子：r0
. J" C' V; C8 O1 E& ]( J$ H8 B%缩小系数：c
%自变量向量：var
7 O5 L& v6 Q( H: M%精度：eps
%目标函数取最小值时的自变量的值：x
3 r# H$ P* Q/ V% D0 ~%目标函数的最小值：minf

  d( o3 Q# ~4 m7 e5 S: }+ E2 u
9 e. D5 v; c  ]( Z  F5 B3 q目标函数：
( p- B2 V" F, W1 x! a1 i  efunction f=link_objfun1(x)
0 a  G$ i& N% Q" t" nglobal L1 L5
* t; L& L& d3 j, qf=0;
! ]" U& s" n* Q7 Z  msnti0=0;
8 g5 E# J6 L% a! n( k2 Vfor i=0:5
snti=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);
7 G. g, d) t/ P1 V3 i# K2 ^5 Vpsai=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)
# s' m8 ^2 {% Y! c  A   bati=fai+psai+acos(x(4)/x(3));
* B( a; O; N& k; p* ?elseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
end
1 p- F8 w! |7 d7 v: y; Sf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
7 Z/ R4 w) x3 X% n8 n

6 U9 s! r! ]/ O( d6 R0 Q1 G( K请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
' I$ Y' b$ V1 {/ W* esyms t;
6 \% y* l% Z0 ka=4;b=3;
f=a*t;
g=[t-b];
- o, z6 N3 S* P+ t) [/ K[x,minf]=minNF(f,[5],g,10,0.5,[t])
  x8 a! s* o* N3 f: D我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
$ ]( E; a8 D1 If=0;
snti0=0;
. z+ f- k" Z- q; v2 ufor i=0:5
snti=snti0+pi/6*i/5;
; y5 `. ?/ e5 V, p7 Hfai=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);
) d, p; s* }  e" melseif snti>pi/2-n
4 o4 X" S; a/ H   bati=-fai+psai+acos(m/a);
end
1 P2 s* ^1 y& K% X# O$ j+ R4 `f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
! ~, {; i2 V( }+ P) f* a8 |& Jend
* e- d! W! G1 ?% c9 L9 ~- U运行如下：
  i0 s' W, A, \# `6 Iclear all;
global L1 L5
L1=3000;
' k2 {; f5 R# m. S! r" EL5=800;
+ T4 P# Q! W: i9 {+ [: g5 Osyms s t a m n;
! I5 c5 Q8 g. h% {6 I- ~! {f=link_objfun;
: s: M# l" V4 N: U2 X% l* E" ag=[s;t;a;m;n;L1-s;t-s;a-s;a-s-t+L1;s+t+a-L1;s+t+a-m];
6 Q1 y" k& T2 X[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
??? Input argument "s" is undefined.
' V) a( m7 x. @% {2 v* b- x
1 @  |* U3 `" x# C% zError in ==> link_objfun at 7
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);