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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
%不等式约束：g
%等式约束：h
%初始点：x0
%罚因子：r0
9 h0 v; V5 K8 y" Z+ K%缩小系数：c
%自变量向量：var
%精度：eps
# g+ K# E( d0 E* P) ]+ }: b, S%目标函数取最小值时的自变量的值：x
, W6 r% C- T. M3 q  C0 J%目标函数的最小值：minf
5 a6 q8 Y( f4 t. i
2 O. T5 a1 \1 R/ f7 [1 B
目标函数：
function f=link_objfun1(x)
global L1 L5
# l' P$ F0 y8 y6 e8 Xf=0;
& Z; _9 U# }% _! }9 I) {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);
3 X+ O3 e+ S4 C; L' ~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)
7 A) ^8 x6 i  J4 H0 u7 `0 s   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
# m5 Q1 [# ^% X3 S4 n  W) |end
+ R' w8 S% O* ~+ O$ D9 wf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
8 S, H6 F0 K0 A7 F

请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
& }4 M. ^, Y2 d/ m9 }0 }0 M) Isyms t;
" d' b) Q& A  i; B2 C6 l( D- Ba=4;b=3;
' h) d- O" n; N: Kf=a*t;
g=[t-b];
* |2 |( E1 R4 o[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
: u3 [+ [) w: X9 b4 e' pfunction f=link_objfun(s,t,a,m,n)
global L1 L5
& g" J* U5 |& a$ `6 T/ Z' _9 R1 `f=0;
: i$ x; }; ?  G- Ysnti0=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);
! S, a' M5 b- ]) u& c9 Fpsai=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));
* z' y, d# i9 H* ^if snti<=pi/2-n
   bati=fai+psai+acos(m/a);
; {2 @- u. b& i; G5 D! l8 Y( eelseif snti>pi/2-n
$ U% o! ]2 `' p0 t. U! f   bati=-fai+psai+acos(m/a);
end
) K: ^  n& n8 Y" w& k! S2 q5 Ff=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
6 B! P; \7 l: N  q* ^  ~& vend
1 V0 |3 L0 f8 L! L运行如下：
+ W0 p8 o3 V# oclear all;
4 `! {+ y  Z, f# B% s3 jglobal L1 L5
L1=3000;
L5=800;
) v: G: @2 n9 r" Q' rsyms s t a m n;
f=link_objfun;
0 ?9 S* I) X; c) dg=[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])
$ Y9 f4 T0 p: V5 R2 T7 H3 x# L) _??? Input argument "s" is undefined.

Error in ==> link_objfun at 7
9 O) e/ l/ j, Z0 L1 y; sfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);