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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
! x; y( n3 t) K2 w; S' r8 Ifunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
0 T, X8 O) C! o2 l" H%不等式约束：g
) L6 l3 A; s* E3 H, E" g9 ~3 [/ J%等式约束：h
%初始点：x0
: ]! r0 o8 K' D! W1 g9 Z%罚因子：r0
%缩小系数：c
4 A, D% x  |! F%自变量向量：var
. D) J2 z0 A: e! H6 p1 }" U%精度：eps
%目标函数取最小值时的自变量的值：x
# U7 n$ M2 G5 U% b. U# `. o+ v7 K2 Q%目标函数的最小值：minf

& U. Z# d% ~: E+ N" r
目标函数：
function f=link_objfun1(x)
7 z$ o" f! e1 G& [* Vglobal L1 L5
5 }5 M  `5 k3 |9 l8 }$ k& j- m: Bf=0;
, {: c% W& ^1 o+ l! `snti0=0;
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));
0 v6 ?- B; W& W7 p4 E; Dif snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
1 |2 s( ^6 @( ^/ Y6 K. y% Qend

0 Q  w6 J5 G! U, O
( V8 x0 k& v' z% L) X请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
$ i0 z. r, {: Y2 M# ]/ C: X运行如下：
syms t;
1 J6 ]* @7 a: Z! f3 z3 W: ?- Za=4;b=3;
6 C2 |% M! ?. z. i1 T0 r; Mf=a*t;
) U- e: `( g0 |% s1 Qg=[t-b];
6 o# c' J# e# [" o[x,minf]=minNF(f,[5],g,10,0.5,[t])
: h# q2 r' q% B: i$ b2 Q+ l/ e8 i我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
4 X! p1 D* T2 ^  l" i* U( y- n9 Bf=0;
snti0=0;
for i=0:5
snti=snti0+pi/6*i/5;
% T* G1 C/ e+ L- Ofai=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
" C# W- q. m' H5 L( U  x/ ]   bati=fai+psai+acos(m/a);
& u6 x; j) F3 s9 \elseif snti>pi/2-n
5 I9 l) ?5 |. N2 c   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
& A: H* S+ o3 G3 ]  ]- y! H7 t运行如下：
clear all;
6 K: R/ ~. o2 f* j0 I, N6 Tglobal L1 L5
L1=3000;
L5=800;
syms s t a m n;
5 W2 E' J2 |+ [$ ]2 @$ E& v, Qf=link_objfun;
1 x# R9 a- I+ C! c& Sg=[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])
9 M0 |. |" t% N??? Input argument "s" is undefined.

. x! K! F3 M/ Q# R7 ~' tError in ==> link_objfun at 7
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);