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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
+ y3 u+ `8 d; ^' u/ u! u%目标函数：f
%不等式约束：g
%等式约束：h
* {( U) K- d5 R% {2 L7 u( x0 s%初始点：x0
- l9 z& m; G2 j+ [%罚因子：r0
, S8 W- J- m( U2 F0 [%缩小系数：c
%自变量向量：var
%精度：eps
%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf
; q: @% D4 S7 E

目标函数：
function f=link_objfun1(x)
global L1 L5
1 X: ~" L0 y( kf=0;
0 H- Z, I0 b7 i- |! o) s3 @4 @1 dsnti0=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);
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));
$ }2 k, a- Y. ?5 h- x' h/ Y: tif snti<=pi/2-x(5)
7 k" H5 C' [- g& V   bati=fai+psai+acos(x(4)/x(3));
! z% L5 j% Z# }9 b, O4 ~elseif snti>pi/2-x(5)
  q2 Q8 w' _& i" g3 t   bati=-fai+psai+acos(x(4)/x(3));
end
5 Y/ H# k# K" O2 S1 ff=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
8 O9 r& d) G' d' x8 R' W, Fend


) j3 ]( U; n- }3 q% f请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
! {& a3 ]2 _1 t- N! \. O: }syms t;
a=4;b=3;
f=a*t;
g=[t-b];
' S1 B( y' g: M1 I% Y+ e[x,minf]=minNF(f,[5],g,10,0.5,[t])
0 l! ~  _/ Y$ w5 I" b+ z/ }' l我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
! _( I2 J! O: p8 {  xf=0;
+ a9 O& ]8 @* P1 y! e/ _snti0=0;
for i=0:5
snti=snti0+pi/6*i/5;
+ O' T  {& t" f" k2 ]fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
* [  a: S9 t8 `; Q9 Lpsai=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
- ?* K4 d  j0 @) {1 @   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
! N! c2 i: m+ u, B& K% qend
运行如下：
clear all;
2 W7 U' R5 i7 h) w0 o6 F! [global L1 L5
- X- {) Y8 b) m: k& f6 ~L1=3000;
L5=800;
  d9 Y: K! v( r& c6 j" M0 tsyms s t a m n;
( d# @9 B, w+ ^7 Z$ \4 ]5 K% a$ ef=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])
" G3 F7 c. }( Q+ a5 t8 Q/ ^( i% s??? Input argument "s" is undefined.
8 w( _' R8 Q5 \
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);