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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
5 r4 r( t$ e% i/ ]. d9 V" {% z惩罚函数调用格式如下：
7 Y, ~6 s; K/ Y* x6 [6 u- Nfunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
9 b1 [. W2 q/ J0 m%目标函数：f
%不等式约束：g
% Z8 K) Y, X6 f' w# N7 @. x: J1 T%等式约束：h
%初始点：x0
%罚因子：r0
: Z1 j: w3 W+ d+ v%缩小系数：c
%自变量向量：var
% l  R3 O, W( D" T2 f%精度：eps
/ C0 ~/ N8 p0 D- i: M* L%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf

* s! V7 ^* u$ }) D% _
! y' ~- p) x4 I% i+ Y* M目标函数：
, ?' p% h) V! ^% O% B# N2 mfunction f=link_objfun1(x)
global L1 L5
' L- H1 t+ E8 of=0;
! I& j' G! C( K6 l% w. Ysnti0=0;
, I  n! ]' G, Q  ^9 n% zfor i=0:5
$ L& }: k* L+ [% J1 Z' @snti=snti0+pi/6*i/5;
, h; @& m& x  R% a1 Pfai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
7 z9 \5 u/ E0 Cpsai=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)
. h: I9 }$ }, @( x   bati=fai+psai+acos(x(4)/x(3));
9 Y. U6 B7 S) R1 b4 lelseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
end
+ b& p( p# P% J. F, df=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end

  Q0 D+ o" z- p( G: d& T+ w- i
请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
6 w$ n2 q( n$ t& L& Z4 lsyms t;
: u" ^, y. I& x: K  u7 P) J% Xa=4;b=3;
9 t+ y( ^7 Q- n" _2 `2 S- B' d! z- ^f=a*t;
' G$ Z) w% n# d4 Vg=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
' k4 T# o' w% \1 g) b4 Uglobal L1 L5
f=0;
7 c' x' M  Z5 B* }snti0=0;
& v  B+ P& ?! hfor i=0:5
) C) {& ^. ^5 k* lsnti=snti0+pi/6*i/5;
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
% M0 U3 M6 k& wpsai=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));
' a9 o+ r7 }0 q% p5 a/ vif snti<=pi/2-n
2 y& Q, f- H' v) b   bati=fai+psai+acos(m/a);
2 }4 \" g+ E% relseif snti>pi/2-n
+ Z8 I0 I; f7 D( }   bati=-fai+psai+acos(m/a);
$ Y" z2 C( t2 X9 t0 j2 hend
( ^! m% L" }5 U4 Y* o) e0 d6 j( Xf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
运行如下：
clear all;
% b/ L+ d1 L6 t: Q+ B5 s8 Sglobal L1 L5
L1=3000;
L5=800;
  I1 i5 Y4 o$ g, Esyms s t a m n;
) o& m, L3 v# y9 {3 Q$ mf=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])
2 s0 p) P: f# p4 Y??? Input argument "s" is undefined.
" f; B7 n  f; N/ E. C
2 q) p$ |) h0 y# a6 JError in ==> link_objfun at 7
6 t9 A+ y: |& r- Vfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);