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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
; k; {% \0 }# c" `( F4 Vfunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
; j" }- J, F6 h; ?7 Q! w%不等式约束：g
7 }: A* J) a+ y%等式约束：h
%初始点：x0
2 Z' N# K& T. G; X% Y%罚因子：r0
" K. V  U8 g4 z; p/ ~' @4 B%缩小系数：c
%自变量向量：var
7 p; F, U" u# I5 W%精度：eps
%目标函数取最小值时的自变量的值：x
/ a  U+ X! @/ R9 O# A, {2 U2 N%目标函数的最小值：minf
1 z8 J% f( G5 \% N5 t# @

  E/ G1 \9 z' U: J0 }目标函数：
4 ^- X( u+ M- }5 `$ x3 Ufunction f=link_objfun1(x)
global L1 L5
f=0;
snti0=0;
for i=0:5
( E1 f6 m3 B8 c' N2 \  [" \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));
7 E- o, a/ M. _* N% ~% n7 V& o% ~if snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
& n1 Z" b1 j* o* g, h   bati=-fai+psai+acos(x(4)/x(3));
+ t$ x9 k& z7 M" J+ G: send
4 _+ ^4 E1 s( \- A0 df=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
+ ^! e: K( @$ j! H$ F( Bend
- n. M; |! Y: z. C8 F: U0 L1 s
! I9 o' n: v, j; k. N5 x
请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
4 f8 o3 I, a, m+ ?' g0 `$ s运行如下：
syms t;
a=4;b=3;
; F# V9 T: Y3 [f=a*t;
g=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
( P1 c2 j% C9 Q, `* F8 m: s0 X我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
8 e9 {5 K4 K' J; d# P& Kglobal L1 L5
f=0;
snti0=0;
for i=0:5
" B; O. o9 c! [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);
0 V/ w: X9 O: X" wpsai=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
% b7 [* O# f" e. I+ S' z; z6 N   bati=fai+psai+acos(m/a);
( p: \4 n. M: @, i8 ielseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
end
9 v! t( [/ S' f% }f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
2 N6 a' r5 T5 d. a) Tend
# M: b* k5 U  o. n6 S6 c运行如下：
. z4 |0 q) D. ]4 mclear all;
' `7 A" O* c! s( O; @global L1 L5
% B# z7 @1 b! K; |8 C* FL1=3000;
L5=800;
syms 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])
??? Input argument "s" is undefined.

Error in ==> link_objfun at 7
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);