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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
9 e/ O0 s+ G4 i. f0 O1 w$ h! A* ^惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
4 [2 V  x! N/ V( o, B, E%目标函数：f
%不等式约束：g
7 J3 v# T& _+ Z%等式约束：h
) [. C3 l- S0 ~5 e) @%初始点：x0
9 T4 j8 O5 S& t; ?$ h%罚因子：r0
%缩小系数：c
%自变量向量：var
: q2 Q6 b: t  a%精度：eps
! v% f8 d- R8 D6 T! H( Y%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf
8 Y. n% D" L' `0 c, h1 U  f
# v) ~/ A! w9 Q$ R+ Q! a+ R
目标函数：
function f=link_objfun1(x)
global L1 L5
8 D; f) m& j& B% l( ^' Y8 |  Sf=0;
snti0=0;
for i=0:5
( z* j+ M/ F/ w) Asnti=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));
if snti<=pi/2-x(5)
' h- a$ R6 ^9 Y5 u   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
5 \2 m" t2 D# [3 Q6 p% K   bati=-fai+psai+acos(x(4)/x(3));
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end

: T* Q) X1 U' C! D; n
# O3 b4 x# G3 D' ]# k请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
- G' r" d1 ]4 v0 {6 E6 V; M运行如下：
syms t;
a=4;b=3;
7 r' U+ j2 r; k$ R' r+ Xf=a*t;
g=[t-b];
" P$ t0 s) b. p) m8 g2 E2 e[x,minf]=minNF(f,[5],g,10,0.5,[t])
- ]/ ]* P( Z- c我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
" D4 E8 z& v$ q$ t+ nf=0;
snti0=0;
for i=0:5
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);
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
   bati=fai+psai+acos(m/a);
7 x% d5 y7 W+ V3 Celseif snti>pi/2-n
  Q) W' J+ _  z3 v' z   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
$ j/ s8 |) T+ g* v) rend
运行如下：
clear all;
global L1 L5
L1=3000;
L5=800;
syms s t a m n;
0 c: j1 {1 m9 g6 b8 _. [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];
. F4 G2 o/ l; k1 F$ R1 _( }[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
??? Input argument "s" is undefined.
) a/ i2 a+ m5 U8 p) ]
5 p6 k1 L8 S( Y: t$ h% v, Z, I/ vError in ==> link_objfun at 7
# Q" @. T& ]$ _  n' Nfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);