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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
  m* [. C# N$ i- a# `% N3 e惩罚函数调用格式如下：
/ L: M, m0 Y3 @- T$ Z3 }% _2 F5 Hfunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
) C; d2 P  E7 t) b5 |1 k%目标函数：f
%不等式约束：g
( b* C4 E6 l# u0 z8 ?/ D%等式约束：h
( I- S* I- }7 F. b& K9 f%初始点：x0
4 W; N( z- e! U0 d! y; d%罚因子：r0
* Y+ X' G8 i4 F% T  ?# n0 B5 }0 g%缩小系数：c
%自变量向量：var
%精度：eps
%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf
8 b# r( a8 J1 C

2 K" B9 K8 O3 o( B0 C目标函数：
function f=link_objfun1(x)
; h2 D. x( t" a* U7 Rglobal L1 L5
9 h, k. P- Y' [7 b' @f=0;
0 V$ C( n0 X' W2 j# |) G5 ]: Isnti0=0;
  B# C9 q0 _3 e& }for 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);
psai=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));
) [/ ^. u: ]; R9 q! ]+ i: qif snti<=pi/2-x(5)
( F# A  \* @1 F: m: Y( C7 Y$ r   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
" O+ ~2 R4 f, K6 s   bati=-fai+psai+acos(x(4)/x(3));
  U7 p8 P) x7 uend
! n" l, V' w( N, gf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end

5 w/ ]* h; [5 i7 D  G6 H+ d8 h/ V
# r) P+ e( ~9 h# d; }请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
3 h: q/ L# S4 q# ~' w% P运行如下：
syms t;
a=4;b=3;
( B+ J% @& |. y+ J+ W; bf=a*t;
g=[t-b];
& }8 R/ Z3 E4 u) j[x,minf]=minNF(f,[5],g,10,0.5,[t])
9 c$ X3 T* R" P0 D我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
$ X* E2 d. ]# Y5 ]" P( O6 E; \; Vglobal L1 L5
1 C1 ^0 t6 b7 s# H/ o' }f=0;
snti0=0;
for i=0:5
$ M; n' N, E; l# wsnti=snti0+pi/6*i/5;
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
6 D  k3 A9 c7 i6 S+ Y# 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));
if snti<=pi/2-n
- \. q+ }: o* }  K/ J   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
- t; G: a/ r; ~* g8 nend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
/ n  Z& _7 Q+ p  Uend
, |; S6 G$ ^+ b2 P运行如下：
: ]; x4 ^* K4 _  g) Vclear all;
global L1 L5
( ]" ?% [7 J; i  o8 {/ |L1=3000;
L5=800;
syms s t a m n;
/ {7 G, e6 \; h/ I# T- Bf=link_objfun;
' R. \# i' U5 g! n0 y6 J( P- E; X$ ^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])
??? Input argument "s" is undefined.

& b" M; o, A9 S! Z# e2 |Error in ==> link_objfun at 7
  \0 S+ `- h4 H" M3 ~& g; Bfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);