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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
9 I& Q7 ~( R+ j( b6 Nfunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
; M" v" y: J% |9 l! A%目标函数：f
, N  W/ R0 E& q9 K$ i% }* L# ?. n%不等式约束：g
%等式约束：h
% o0 a4 w: b# ?1 [%初始点：x0
%罚因子：r0
%缩小系数：c
%自变量向量：var
%精度：eps
%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf
4 g# q( s/ W, I/ P3 y- g

$ Z% x" S7 D- n3 @( J3 Y目标函数：
function f=link_objfun1(x)
* X2 Y9 F( r8 |+ Y7 o( Oglobal L1 L5
" q" E& j- Q" Df=0;
snti0=0;
- a7 g7 M% [9 b* d; V2 Hfor 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)
: b" [( ?  C3 c7 m   bati=fai+psai+acos(x(4)/x(3));
9 `1 z$ {3 i2 U( Felseif snti>pi/2-x(5)
/ Z0 S7 e5 x- O' t6 U# D   bati=-fai+psai+acos(x(4)/x(3));
% i, l" O6 I+ M, g, qend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end

/ i* M8 r! A8 h0 x4 U
请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
$ P# a& Q" i; V9 K2 e' o% Qsyms t;
4 {3 R1 D1 @! N" m/ {a=4;b=3;
f=a*t;
: I8 c1 A' }9 X# n; lg=[t-b];
. T. ]8 v' @' G% ?( {/ ]% Y[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
( C, S3 k  J' ?3 w' fglobal L1 L5
f=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);
/ K, \8 h- O- m0 Gpsai=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));
9 w) t  c3 }0 U. n' s  Jif snti<=pi/2-n
   bati=fai+psai+acos(m/a);
0 B) U! x8 {) c9 s( celseif snti>pi/2-n
% f) W. G6 b$ x7 t   bati=-fai+psai+acos(m/a);
) i+ W  u9 a- F0 }+ bend
, k9 I: X- B7 H2 O" K2 qf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
" C) L! d; z) p; l0 z1 [end
6 v) P1 s: ]/ B0 |% i; K1 A, E运行如下：
clear all;
global L1 L5
3 a6 U! j( q9 v- ~/ \0 [( qL1=3000;
( U1 u, b" A. d! tL5=800;
syms s t a m n;
  I% y9 c+ Z+ j) G( Q9 [7 Kf=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.
* ]- f: `$ @$ R' J! s
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);