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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
9 v. t5 E1 }' G+ r+ z9 M7 k: l' ufunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
; [6 l- b6 ~* A2 e) C0 w# [' K* J, f%目标函数：f
" s3 \7 n  [8 n%不等式约束：g
%等式约束：h
/ [5 y) [  C% d% l8 G%初始点：x0
" u+ d1 U8 D/ ^0 e* K%罚因子：r0
%缩小系数：c
6 h, A- w/ E9 i%自变量向量：var
: A1 e8 V% U% N%精度：eps
0 P4 V' z( n" u. O6 p$ h0 i0 G! A%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf


/ \# g) p+ j" u8 s# o目标函数：
function f=link_objfun1(x)
global L1 L5
6 M5 m2 l& J( V7 Qf=0;
+ m5 u" V- L7 A: J7 v' h; G* G  Ksnti0=0;
for i=0:5
! s$ z  r4 @! Ssnti=snti0+pi/6*i/5;
5 U) \0 z8 I7 q4 r: nfai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
$ o% U9 `$ T, e2 Z* m& Cpsai=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));
9 D9 I. {0 w; w# q; Kelseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
end
/ c# F" V# g" u% }  L  V% }) B% j% t% Lf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end

/ C8 a  Q; n* D. i
6 w2 w+ T# E/ `. T请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
9 u+ ]) h7 z0 E/ W- [运行如下：
syms t;
, _2 _" K: V, ]8 U# Ja=4;b=3;
f=a*t;
- Q! p+ E. U) x, Jg=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
- H' w+ Q6 W6 h  }global L1 L5
f=0;
snti0=0;
8 }$ D: f9 n/ o" A* o$ Cfor 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
' W* I+ j& J. q2 X8 k* |9 ~   bati=fai+psai+acos(m/a);
: d- c6 ]! Z1 K2 p9 velseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
' e* q0 L# [1 l2 J) T运行如下：
clear all;
/ }1 Q1 c- k3 j. O+ `- C( }0 Uglobal L1 L5
' i. y) j+ h( q. M1 z4 a6 cL1=3000;
# [6 E; n- ^7 kL5=800;
# I: L5 w9 f5 E( D! v: }$ s( Jsyms 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])
/ P, ?# I2 `1 `??? Input argument "s" is undefined.
" Z3 t: E  r, T# ^- {6 }' K
Error in ==> link_objfun at 7
. H3 S4 P% R( `, @2 Jfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);