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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
& t( U- R" S% o& `" e7 Q% I惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
5 o! n7 @$ s1 m" M$ C%目标函数：f
! E1 T7 P5 A, e! ?; ^! H3 A%不等式约束：g
%等式约束：h
5 P( X& D7 l0 e5 w8 Q$ ~2 O. C* S%初始点：x0
%罚因子：r0
/ P5 q/ ]7 L+ L%缩小系数：c
( I6 {& F+ G3 M. g; t%自变量向量：var
1 M5 _8 t/ R2 y7 f3 C%精度：eps
: j0 x/ K3 v( X9 R, i+ [%目标函数取最小值时的自变量的值：x
( u6 }( r/ n8 Q%目标函数的最小值：minf


; G; k5 {# R$ o( n$ u" _目标函数：
! F8 Y- S" ~. M6 a9 `3 hfunction f=link_objfun1(x)
) v# s) y/ |2 x# fglobal L1 L5
0 ~* P+ |  u- d% Z! S2 [f=0;
snti0=0;
8 n$ R- J& S6 R5 o; Pfor i=0:5
- Q+ X/ K- J$ m* x; I& ?6 I+ vsnti=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);
/ y$ _2 N8 h3 ipsai=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));
elseif snti>pi/2-x(5)
$ s4 f  g- {" k9 Q' ?; z1 q   bati=-fai+psai+acos(x(4)/x(3));
end
/ E( g* ^6 ~8 ~1 |f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
* @/ Q, h1 B9 g. @1 `9 }end
( q# T+ O: m3 Z/ l/ t, [$ T

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

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
syms t;
  {# W0 l" a* B7 `- Ka=4;b=3;
f=a*t;
g=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
7 Y/ m# [' X4 x1 u; Gf=0;
( h  p* Y2 |8 s5 z' A  hsnti0=0;
3 \- F9 I' e- o& y) [. {for i=0:5
snti=snti0+pi/6*i/5;
: M+ \2 H5 U" c3 v2 gfai=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
4 ]9 X; N9 x" t# ~: g  U: G4 ^# w  c- p   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
0 }4 K! m' ]  k& Z   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
  y  Y% t0 U, h0 jend
运行如下：
clear all;
global L1 L5
* f- t+ M8 X6 N4 iL1=3000;
L5=800;
syms s t a m n;
7 R- ], K; z1 r9 M* S9 Zf=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];
( l* ?4 z5 L& E3 s. l[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
7 f! Y. c) x% V: m??? Input argument "s" is undefined.

. U  S9 r5 T; K2 s, C5 vError in ==> link_objfun at 7
* w$ t3 @+ u2 z& J) x6 J, `  T6 o0 Q1 Hfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);