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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
& y; }# d/ `4 w4 F% U%目标函数：f
%不等式约束：g
%等式约束：h
) C2 Y9 h, `% B3 H0 J%初始点：x0
' ^& A/ W% ]" ]" D0 P' [0 k% _%罚因子：r0
%缩小系数：c
$ M3 G6 j  L  O# k- P, i%自变量向量：var
6 |8 T4 t9 P7 m, |" g, k3 o%精度：eps
! {. l2 E! w( j4 S. ]0 K8 ?) K4 ?%目标函数取最小值时的自变量的值：x
9 s- i& F5 c# E- @- v5 N%目标函数的最小值：minf

2 x) J! f9 y7 [& x$ I
/ H' s/ @( ~# {' T1 D9 ]) M目标函数：
function f=link_objfun1(x)
) U5 a% v! a2 I* W5 [global L1 L5
* X  Z/ F# u" l8 x7 P; nf=0;
! o: N9 M' K6 |  C0 D! E; {& T1 a; Xsnti0=0;
( Q7 ]  D7 _6 S/ D. \for i=0:5
$ ]! x" [/ ^6 M- S: m8 vsnti=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));
, F3 n9 _- g- }% J* k9 Zif snti<=pi/2-x(5)
  I! C: j5 p* S0 {( _4 P   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
  u& j" ^- i5 }, x- jend
' y& J& _2 D4 S( y6 R3 a) yf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end

7 b$ Y: Q' Z, `4 p+ R
请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
: b* J) ^$ x1 ], }# S& qsyms t;
a=4;b=3;
f=a*t;
8 h0 e- A5 @* ~% z4 ^3 g4 _g=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
7 d  k7 a/ h% ^8 df=0;
snti0=0;
for i=0:5
: z0 e& }! V& d; @5 gsnti=snti0+pi/6*i/5;
  s7 H4 S& t1 d/ Z, R8 wfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
! Q: ]/ ^9 [( i/ Q9 Ypsai=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  }; q' y' }# Yif snti<=pi/2-n
   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
# r3 J9 q; u+ {   bati=-fai+psai+acos(m/a);
end
. l, V- z" V" O7 H# J$ q& \* Yf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
$ u0 k! y0 Q* B: N9 Hend
5 C0 g4 F- Q) W5 N5 D# @运行如下：
clear all;
9 b0 _% P8 O/ H5 Wglobal L1 L5
L1=3000;
L5=800;
9 Q- l; G. ], g- Usyms s t a m n;
- o: a  g4 I, V7 ?% gf=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.

  ?( e1 w0 P' C3 D# z+ wError in ==> link_objfun at 7
5 h- [( O2 T, o7 rfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);