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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
0 j, c) N; t6 i' d; ~2 d%目标函数：f
6 F4 _( [5 d: e$ O* i" C%不等式约束：g
%等式约束：h
/ _4 W+ q7 J. ?$ ]" ]%初始点：x0
+ U2 N. N- C# D' X* g) p  B# v" N%罚因子：r0
, h8 [+ ^, o0 S4 L%缩小系数：c
%自变量向量：var
%精度：eps
%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf


目标函数：
function f=link_objfun1(x)
global L1 L5
5 k! p% S. F( K+ ?/ Xf=0;
8 a# g0 Q& E) _/ c3 H/ i( Nsnti0=0;
for i=0:5
* k6 f% U" Q: h) {# E4 W* ysnti=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);
2 D& Z2 H' r' T" F/ bpsai=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)
/ o7 [' H$ b; M  P( b   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
. s* X' I- a" l% O0 y4 Z% \# [6 Lend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
9 _/ W. C. Y2 F- {, O, R/ b5 Cend

# X$ W1 r) d" B, i5 c- z
请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
* \9 R& L( ]# E2 y  F( b6 q2 Y运行如下：
; X: y  n% l& m- _. t" z! Xsyms t;
" M# ]9 Y$ o+ u* P! j" a4 oa=4;b=3;
# b6 Q. b+ t" l7 Vf=a*t;
g=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
1 a8 C: [5 q5 w2 m, l8 U; i我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
+ Z, M) R. E& Y4 oglobal L1 L5
9 X+ c, e5 }( H% M+ of=0;
snti0=0;
% X3 A& N0 ?* N4 dfor i=0:5
snti=snti0+pi/6*i/5;
fai=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
   bati=fai+psai+acos(m/a);
* o; M1 \( P$ @+ F: s; Z. aelseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
5 K* w" r3 g' ^" p8 c' Iend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
( l+ ]3 b  f! M9 o0 @end
' m, j  F  l+ f2 V% w运行如下：
clear all;
global L1 L5
L1=3000;
% O7 }" c, \* H# l/ yL5=800;
syms s t a m n;
- l) J$ L0 a3 M' _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.

* \0 s, D0 B1 ?/ |8 X6 kError in ==> link_objfun at 7
  f4 `; G/ K( M, s( Hfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);