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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
* E2 T( @1 b* r%目标函数：f
%不等式约束：g
%等式约束：h
%初始点：x0
%罚因子：r0
%缩小系数：c
%自变量向量：var
%精度：eps
%目标函数取最小值时的自变量的值：x
- {3 \2 F3 X$ ^9 x! C3 B%目标函数的最小值：minf
7 G0 a5 t3 e7 z1 H4 @, ?
2 a7 f  B- _8 _& n
/ o3 b& x! O$ O) Q目标函数：
function f=link_objfun1(x)
0 N* v1 ^9 U* e& r) ^! G' Oglobal L1 L5
6 I9 [# ]! E, [' G1 ef=0;
9 k1 a% P0 D$ Lsnti0=0;
/ E" w) G3 y( d# }! U- gfor 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);
& b3 E; r/ c/ I+ W5 spsai=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)
8 n; c0 h& G4 ?% k8 `   bati=fai+psai+acos(x(4)/x(3));
& E9 A* I( a# y- T5 yelseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
end
! m: Y" I' V: K9 T: k: gf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end


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

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
, M7 f1 T; U' ]' h0 {+ i  R运行如下：
syms t;
a=4;b=3;
4 ?' d) @( x. `3 C2 Df=a*t;
g=[t-b];
1 R2 |5 I7 D1 t4 N4 r2 Z[x,minf]=minNF(f,[5],g,10,0.5,[t])
& l- m& S; b  P6 u# E; H. n4 A我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
& i3 O' [6 U" T3 N9 w1 uglobal L1 L5
. @4 d! \4 I0 W  l. N( mf=0;
snti0=0;
/ E; d' Z" {+ `for 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));
+ T7 ?% z- B; d# M( Dif snti<=pi/2-n
" C$ N7 H) x4 i1 t# T% {! ?& K2 T   bati=fai+psai+acos(m/a);
! N5 p' i. F& ?; Zelseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
end
2 G. y9 F" u- Y  s$ \' o6 ^f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
运行如下：
clear all;
& Y6 X; z7 h; \) K: T9 X, ]global L1 L5
L1=3000;
L5=800;
syms s t a m n;
" @/ I. T! L( s; H" {6 D8 Jf=link_objfun;
0 p' C% O( z( \+ g# Pg=[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])
( a5 a% b  Z3 ]1 d" |??? Input argument "s" is undefined.

" I" |4 c; Y. p; `Error in ==> link_objfun at 7
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);