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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
: m% ~, c& B2 h$ C  j" Q惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
& U* d& ]# o2 c%不等式约束：g
%等式约束：h
; K) V. R7 W2 y. S% H; |%初始点：x0
%罚因子：r0
2 O0 P  H: a/ u' N0 u' X- S%缩小系数：c
  m* k7 h* M8 `" l) ^' c3 Z; z7 ?2 Q%自变量向量：var
%精度：eps
%目标函数取最小值时的自变量的值：x
9 n3 L$ D4 @, O+ k( r4 H; `3 T) g%目标函数的最小值：minf

/ j- q# S3 g' o+ _* K
: r" r$ y, l+ H; d目标函数：
function f=link_objfun1(x)
global L1 L5
f=0;
snti0=0;
" ~+ f: j6 U/ Kfor i=0:5
: t% L. ^0 c4 w) Z: _snti=snti0+pi/6*i/5;
1 U% k) h# `; {2 a- e7 C* bfai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
7 x; h! J1 G4 T- l( i5 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));
. L. S  e- N' }$ }% oif snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
- O* y. Q; m( {" R7 Celseif snti>pi/2-x(5)
$ {8 {+ u) q* }   bati=-fai+psai+acos(x(4)/x(3));
9 X0 r# Y' w. j" _3 }& |end
- v: W- e: _) Q, u5 X! [f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
$ {7 I  I2 _& q

请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
" W+ H. U" U4 hsyms t;
a=4;b=3;
f=a*t;
; j) K' n8 m/ d" yg=[t-b];
8 _$ ]  n/ s* g9 [9 E" l[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
f=0;
snti0=0;
0 }) r" Y. \$ [" Yfor i=0:5
/ R+ Y) I7 }" W- t8 B# m$ ksnti=snti0+pi/6*i/5;
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
- y# G2 [, F/ spsai=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));
" b; u: ^. m( }4 Y/ f% o5 A0 ?  dif snti<=pi/2-n
   bati=fai+psai+acos(m/a);
% _1 Q: z, U. @1 jelseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
end
5 M7 y. v6 s5 L/ [* Y0 k; |f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
3 G1 P" N! \- T/ hend
1 g3 k  I* u2 b: p运行如下：
clear all;
global L1 L5
L1=3000;
L5=800;
* t% W* N# d( j- b' b/ U- d$ X5 Jsyms 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];
[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
. P) i0 a! G# @6 P& E- P* E??? Input argument "s" is undefined.
# _% w3 S% U2 G" D# n
Error in ==> link_objfun at 7
4 p$ K7 {" h+ ?; Y: R: xfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);