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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
, z# [' D4 G7 @/ I: U. F, j# \function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
$ ?3 a4 S. b+ Q8 i) o+ R0 T%目标函数：f
%不等式约束：g
%等式约束：h
%初始点：x0
& L+ Q1 i" L) ?0 p%罚因子：r0
%缩小系数：c
%自变量向量：var
%精度：eps
4 J$ C6 Q; S* o( c+ {%目标函数取最小值时的自变量的值：x
6 o2 V/ O) p, ]! c, l%目标函数的最小值：minf


目标函数：
1 M8 q. I0 S0 f3 G% W) t) m2 R4 Tfunction f=link_objfun1(x)
global L1 L5
f=0;
snti0=0;
for i=0:5
snti=snti0+pi/6*i/5;
/ b- m( f6 n/ V" i$ u  ofai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
8 _' V9 L" a, o" M/ }% _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)
4 s' E  F4 k5 |, y: a5 D   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
; Q' d) d8 ?# y7 @$ Rend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
$ f/ \& ^' j/ d" a9 y

* T7 e, ~" O) e' a% f9 S: E请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
syms t;
( a9 R! o) V& X# i/ i, U' oa=4;b=3;
f=a*t;
$ m7 p! v/ f3 K8 J* Ig=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
6 n/ R; ^3 w' X- u. q9 [5 w我将目标函数写成：
9 ^8 ]5 W2 {. k7 {/ u( [function f=link_objfun(s,t,a,m,n)
; ]  E7 Q! I! Lglobal L1 L5
! N; v% j2 M' W+ p% T4 ~f=0;
snti0=0;
/ y- k6 d* i+ {for i=0:5
. k, Q8 X7 k0 b% _& V# msnti=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
% A% e4 u( `  U   bati=fai+psai+acos(m/a);
' K1 y& Y  E2 ^; G' u) ^' Z' Z0 h) }elseif snti>pi/2-n
/ h& J- f9 W3 e  d- Y   bati=-fai+psai+acos(m/a);
end
6 t/ v& I2 B7 V$ e5 mf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
运行如下：
9 }' E" h0 l+ x& Xclear all;
global L1 L5
1 U! O% I  q! p7 YL1=3000;
L5=800;
syms s t a m n;
. W5 n' E$ r* l+ }0 A! @4 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];
& ^% ?4 E% L. I8 y/ l[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
??? Input argument "s" is undefined.

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);