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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
! o/ L2 ~: U0 k7 R; _1 t+ |3 X2 ~$ |惩罚函数调用格式如下：
6 `( o8 F( X3 y6 [3 `& b  ~, efunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
2 m( ~/ E% k; F" @8 F5 G0 m/ H%目标函数：f
% B3 Z# {+ D" w6 d6 e%不等式约束：g
%等式约束：h
%初始点：x0
  i3 P/ p1 m0 F) N: Q3 m%罚因子：r0
%缩小系数：c
- k( \) d7 |5 L) Q% R2 a; x: Q%自变量向量：var
1 E* R# f! u% e2 T%精度：eps
, ]( `; Y9 R& z; q4 z%目标函数取最小值时的自变量的值：x
) b% `, }( P0 q. m9 R2 n( e%目标函数的最小值：minf
5 D. w! e  K$ l* [( {7 m/ L0 Q

, H6 N0 d; M8 E8 ?目标函数：
7 }- `4 h( p- dfunction f=link_objfun1(x)
; }* K  }5 P: A- D3 {global L1 L5
, s6 u# c" R) x# |; ?% c/ mf=0;
snti0=0;
/ \# F3 X- J) A4 z, s  J, }' zfor i=0:5
snti=snti0+pi/6*i/5;
! q  E6 o( G4 ^$ o; M$ mfai=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));
if snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
% z! S" s, f( L; u( w6 Zelseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
  ?( j& i; t' ?5 ^3 M1 v( h
/ K/ D8 a3 i- d/ V2 M
1 G- B( _7 l5 p" {2 U& u" \请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
8 q3 G2 I9 C' l' ]$ f运行如下：
syms t;
a=4;b=3;
f=a*t;
g=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
% R3 |' J  w% }  F' Y6 g: I0 mfunction f=link_objfun(s,t,a,m,n)
2 @* u2 _5 A8 Q3 C8 a/ k7 Vglobal L1 L5
f=0;
snti0=0;
for i=0:5
snti=snti0+pi/6*i/5;
; Q' d; Q5 j# x* z6 v7 Bfai=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
0 G+ N( M" ]9 g( g! s! K   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
3 {, h4 T$ e8 I3 Pend
1 t9 E- e2 e3 J- v# N, L5 qf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
运行如下：
clear all;
global L1 L5
2 \; ]# S4 K9 R+ k' |0 WL1=3000;
" X# G1 {7 ~; c- T4 \L5=800;
syms s t a m n;
* U, V* \. i9 L: @; g; [  Xf=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];
0 s9 ~5 x2 O& @+ @+ K' f. o: M' \[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
??? Input argument "s" is undefined.
) H  b5 Z8 V& ?8 D
1 O3 i) k+ h& k- Y5 H  R( OError in ==> link_objfun at 7
2 U1 |. V, s0 \, \6 zfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);