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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
9 B% d3 p- Z" {5 @- V6 W, o惩罚函数调用格式如下：
) x, ]4 E3 h+ b2 `1 M% h! lfunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
/ V5 y; C6 m& n$ g* b2 S%目标函数：f
%不等式约束：g
%等式约束：h
) S: z0 Q1 V  S$ ]6 `5 q, ^%初始点：x0
7 w4 q6 i$ l& f2 c%罚因子：r0
%缩小系数：c
0 H  S( Q: X+ B! @- p5 ?/ J7 f%自变量向量：var
7 M+ k) X+ {) n0 _; B4 C) A%精度：eps
%目标函数取最小值时的自变量的值：x
- L8 j0 B3 q" N9 k8 c%目标函数的最小值：minf


目标函数：
function f=link_objfun1(x)
global L1 L5
f=0;
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));
if snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
- @* ?6 Q& ?2 |) Lelseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
/ m" d7 D; @# _  N  {0 x3 oend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
  W, y3 X$ n* w7 V2 w/ M  s* R4 Mend
5 w/ Y- Y/ j" ]; G6 i
0 ~  @# S2 w- A7 }0 b
: b0 ~4 l. ]$ v. X9 `4 R9 Q. Y请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
( d6 q' F* \6 V. X) c) n5 ^syms t;
a=4;b=3;
7 P( |6 S1 `7 Q! R: p0 Hf=a*t;
g=[t-b];
% m" O+ N) B( H1 Z! \7 ~! `: k[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
3 a! m* p' m, x! y5 N( Z' Jfunction f=link_objfun(s,t,a,m,n)
, C8 K( d- d! ~: Wglobal L1 L5
, l# O# ^$ x- }0 i* i% K7 x& Z# If=0;
* l/ V  T9 P5 ?" J2 S; p/ \snti0=0;
for i=0:5
snti=snti0+pi/6*i/5;
- H/ X: K% w0 q; R/ J! s; t: F) Vfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
7 E" L6 b) v8 N6 ]9 L; c; C  _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));
! q& P# Y3 D; j9 M+ `if snti<=pi/2-n
   bati=fai+psai+acos(m/a);
; Z5 y( ?& d7 {0 ~8 Q, }elseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
! p! h3 D' i! mend
" K) z' j4 O" k+ n9 W, I+ f9 Df=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
' ^0 _9 P7 D- n3 \" t运行如下：
clear all;
global L1 L5
L1=3000;
' [* a! h, n- C' pL5=800;
syms s t a m n;
f=link_objfun;
! `3 J+ z! J/ y$ @+ Cg=[s;t;a;m;n;L1-s;t-s;a-s;a-s-t+L1;s+t+a-L1;s+t+a-m];
/ t7 h$ B8 p2 a[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
$ n+ L- \  u# L: t& ???? Input argument "s" is undefined.
/ F* o' i* @- K7 z! l
- e- M* c) B  T/ V4 Q8 P) ]Error in ==> link_objfun at 7
) Z6 c) a' b. I7 o1 o+ `+ E3 H7 ^fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);