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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
9 z  `% E$ _+ }, t) U" e%不等式约束：g
%等式约束：h
. [" Z) R: C6 W) s%初始点：x0
3 g9 L  U$ E/ d% T7 ^# B%罚因子：r0
%缩小系数：c
2 e2 V/ `" F2 Y2 X5 V! d% N%自变量向量：var
%精度：eps
- h! D5 @; a- k9 J% l% z) s, Q%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf

3 t: x3 p" A+ ?
4 W5 i5 D3 U) R+ Q! j* M' M目标函数：
function f=link_objfun1(x)
global L1 L5
f=0;
snti0=0;
2 B9 I! I( L% C3 nfor i=0:5
4 G( n' d& \) d3 Psnti=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 n3 @* e" j; Q- Y/ wpsai=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)
   bati=fai+psai+acos(x(4)/x(3));
6 Z* C+ J/ ?& |( c0 T4 melseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
end
( w; f; h9 i  d' r/ u0 n% wf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
% H; z( H+ P  U" r

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

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
: W7 Q4 t% J) G* t运行如下：
" q. q3 |( ]7 j9 q: H0 L: n3 n& T0 Esyms t;
: ?, G2 \1 R6 W0 Ha=4;b=3;
0 J9 N* \: q! P8 T/ Zf=a*t;
# r/ }7 L2 Q( r( N5 n( Lg=[t-b];
7 Y4 f1 [. ~7 U7 o* X' v[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
5 u. k) P) L3 Qfunction f=link_objfun(s,t,a,m,n)
9 E0 u9 i! A# _9 g3 tglobal L1 L5
f=0;
* }, C7 C4 q! D. z( T  ?snti0=0;
for i=0:5
# _( P- S: X  S' y$ i/ L5 Msnti=snti0+pi/6*i/5;
# ~! G3 U, K' b4 \+ z5 ?3 M* Vfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
6 p8 h2 w" n9 t9 ~) K% }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
5 S8 X, Q, }4 x2 b   bati=fai+psai+acos(m/a);
1 a6 @6 y, H' ^elseif snti>pi/2-n
- A; z* k  K* e7 m0 Y% p   bati=-fai+psai+acos(m/a);
end
0 Q8 P# N' K, o5 ]7 K, I- f; qf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
运行如下：
clear all;
8 T! C6 E& u, R5 _. A' ~global L1 L5
L1=3000;
L5=800;
; v3 W; z6 t8 a( F3 `; i& ~; usyms s t a m n;
9 `  _2 l7 \: e; sf=link_objfun;
. |  n# i% g4 F  S- [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])
7 T+ ]0 ?1 K, `% {9 Z, S/ x% e??? Input argument "s" is undefined.

2 ]& G) C8 r0 q4 {! }  jError in ==> link_objfun at 7
8 K+ ^- g/ K' v7 @) yfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);