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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
  \) V2 C0 [5 t# l惩罚函数调用格式如下：
6 s7 V/ Y+ F5 b- t- zfunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
, n% U: |3 D2 B1 M; x%目标函数：f
/ q& G$ f. F! r%不等式约束：g
! `9 i% E1 w5 \4 X%等式约束：h
%初始点：x0
%罚因子：r0
%缩小系数：c
% [6 ]. A* B) ~% A%自变量向量：var
, p, O9 G2 l6 j1 E' D/ X%精度：eps
%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf
) M; h# X- Q+ x: x0 o
# ?; I* A7 o  z7 r1 ^; j
% z4 _4 J* o' l; [( \5 X& j& g5 p目标函数：
function f=link_objfun1(x)
global L1 L5
! _+ `7 j2 D9 }8 F3 Hf=0;
snti0=0;
for i=0:5
  r, p1 o& ]8 s; _, k# W8 Gsnti=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);
* V8 v0 B: K; q* ~  f+ z" N& ypsai=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));
+ `. B" U' e: [. i- o6 W' i$ Oif snti<=pi/2-x(5)
! p- {1 v$ g2 S1 A9 v   bati=fai+psai+acos(x(4)/x(3));
" E$ Z  a: b% i! [& N* n8 Gelseif snti>pi/2-x(5)
6 E: o$ F5 M  _5 J/ x9 ^0 A   bati=-fai+psai+acos(x(4)/x(3));
end
8 x# h( A, N# v/ j4 s$ A4 R' {f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
9 M* `/ f# \5 y. C2 [7 Xend


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

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
: i8 X0 b! g4 d" t( T1 Isyms t;
a=4;b=3;
/ i9 k8 A& o4 A4 ~f=a*t;
# d* c4 N6 T0 @) _2 \% G4 k: ]$ Pg=[t-b];
# k( T8 j* A; ]2 V% Z: d* z[x,minf]=minNF(f,[5],g,10,0.5,[t])
, ~: H$ K+ |2 Q- C+ [" H; x0 B$ t' l我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
5 W+ L. S9 e: b$ Fglobal L1 L5
& a# x. H  I9 Lf=0;
snti0=0;
1 F8 }/ M1 L! y% ]% y  Afor i=0:5
snti=snti0+pi/6*i/5;
; j1 E9 j/ c8 w# zfai=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));
/ e/ C, p8 r. ]/ ]) _if snti<=pi/2-n
7 j+ Y; J: k5 x! S# X   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
0 Z1 f( D9 ^& }" C( Z, N   bati=-fai+psai+acos(m/a);
  W5 u+ O, w$ n0 x  h" X5 K4 J- Kend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
( q) o# ^% n4 K% l: u% k运行如下：
* Z, r9 X3 m2 `* l. K& B. A3 Bclear all;
- P, d( l- |, Z2 H; }9 q, w3 rglobal L1 L5
" ?' l; P" E; c/ i! a& @' {. fL1=3000;
L5=800;
syms s t a m n;
! t3 K  Q5 m% ~1 if=link_objfun;
! z4 N  X- c- h5 U# f' D3 D1 G) Xg=[s;t;a;m;n;L1-s;t-s;a-s;a-s-t+L1;s+t+a-L1;s+t+a-m];
- g2 r6 s' U/ {4 \[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
+ Q: V- R& H" }! tfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);