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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
0 h4 k& e2 m7 C6 F. H3 Cfunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
% A- B+ _& \/ w7 N%不等式约束：g
2 ^3 y4 N$ N! e/ @%等式约束：h
& w( j- P" Q5 d% N- D% ~( F%初始点：x0
2 U3 u. y* n: O( w8 @%罚因子：r0
7 X- R8 C# P9 H% Q& S# V* R6 @%缩小系数：c
1 g8 ^2 X; W) S# ?%自变量向量：var
%精度：eps
$ m% k4 `+ S8 j! G%目标函数取最小值时的自变量的值：x
( [0 V5 `9 P4 t3 n9 ^1 z" w1 W%目标函数的最小值：minf
" [2 Q6 a+ d* L' N& ]

目标函数：
function f=link_objfun1(x)
global L1 L5
f=0;
snti0=0;
- o( |0 j9 w" c: A+ |+ I' dfor 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);
4 w+ n& q, Z% G8 R5 I' ~/ apsai=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)
   bati=-fai+psai+acos(x(4)/x(3));
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
' `5 H6 Q' Z( G; Qend
/ I3 E- _" P& E7 N8 [

7 J1 X0 r5 v! M请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
syms t;
a=4;b=3;
f=a*t;
8 n8 H7 [" T6 v9 o4 z( |' n0 o. xg=[t-b];
5 W) `& p# q" ]9 D) T! \6 A7 B& I& t, o[x,minf]=minNF(f,[5],g,10,0.5,[t])
  `1 B, K/ U3 E& \- K7 @- ]我将目标函数写成：
6 P( I5 Y6 L: V! pfunction f=link_objfun(s,t,a,m,n)
global L1 L5
2 S. V1 s) Y1 tf=0;
: i" A' R8 L( z1 T& ssnti0=0;
- |9 W5 P, p# [6 j# ?2 w" L- vfor i=0:5
snti=snti0+pi/6*i/5;
: y$ D4 }* B% Rfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
& \& P7 l$ e% z6 `1 ^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
   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
* E  {  Q/ ~( q/ b2 |: x% b, j   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
7 ~# q& q' ?% B/ |4 C5 aend
运行如下：
5 g5 \! e* \  X" s9 `clear all;
global L1 L5
L1=3000;
& c4 N0 A) V  a" E& q) `! V+ SL5=800;
( F/ P7 e5 X* zsyms s t a m n;
f=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])
??? Input argument "s" is undefined.

Error in ==> link_objfun at 7
6 ?4 _& C0 y7 T! Sfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);