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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
%不等式约束：g
%等式约束：h
%初始点：x0
+ c2 G7 L5 E$ X5 a- N; l%罚因子：r0
& e. Q: i9 B& e5 ~. N' a%缩小系数：c
2 ~( Y' }" n4 g%自变量向量：var
%精度：eps
6 v# F9 g" ~$ j- x+ W* W%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf
- Y6 f- L' d% x

目标函数：
function f=link_objfun1(x)
& s" e& p2 G( H8 l! E4 ~9 |global L1 L5
f=0;
  ?) \  j: H1 Tsnti0=0;
* S  ]" [1 k. \9 Z: u, b4 }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));
, G6 \* o2 k9 B5 _3 A9 f; Y" w- J) j# gif snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
+ b- q% L3 f4 C, Felseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
end
2 W4 \/ L# _' b6 E; R6 nf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
4 d0 d0 U* \" s% R" `5 a1 nend
/ U5 }# H& l& n" v# o7 g) m

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

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
2 \) Z# [6 J' |" p% n- A% ]1 N运行如下：
3 ^) v0 k4 w( @. P+ }' m. l" n) ^0 Csyms t;
5 F1 s  s& W  o" ^: {/ F, Ma=4;b=3;
f=a*t;
g=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
, p- ^) ?0 K1 f9 Lfunction f=link_objfun(s,t,a,m,n)
global L1 L5
( P" x4 }. t, df=0;
. K/ ^' o! \: m5 r1 g6 gsnti0=0;
for i=0:5
9 A8 Z/ c& k: ]# J" w4 nsnti=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);
  V- S2 I% }3 a8 j4 l  {elseif snti>pi/2-n
! d/ d6 r* P; ^/ I4 ~: I5 U   bati=-fai+psai+acos(m/a);
end
: N7 R  c& a4 cf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
6 c2 K' n( W- v' Send
0 n% r5 E6 ]9 G4 W运行如下：
clear all;
- V& {) Q( _0 {' E8 x2 L4 `global L1 L5
L1=3000;
L5=800;
. y4 ^" n: z, U9 ]$ `$ _4 wsyms s t a m n;
/ P' f* r! v! A% [- ?$ ^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];
  Y* e7 i1 T5 i/ T* c. W[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
7 \. a# H7 R+ }2 h" h??? Input argument "s" is undefined.
+ m6 z: Q% Y+ f, z
Error in ==> link_objfun at 7
) w. ]! c4 N0 V3 j9 J5 ffai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);