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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
7 o# [# d( F2 Y/ l: N+ [3 ?function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
9 Y, V# A& M5 K5 M%目标函数：f
9 v" m% J& a( Q4 `7 N$ K%不等式约束：g
%等式约束：h
%初始点：x0
# a4 J+ g" ?1 k) n$ I%罚因子：r0
%缩小系数：c
% P, G  u" B- p$ o6 h1 e6 Y%自变量向量：var
%精度：eps
! l' B" p  o5 p* V: l%目标函数取最小值时的自变量的值：x
1 Y. u' x! e, C! d' z) ~: w%目标函数的最小值：minf


( [5 {- \( L  ~4 [! i8 Z5 V! {目标函数：
+ T# u. A9 {5 Z2 v/ X* u( B) gfunction f=link_objfun1(x)
global L1 L5
f=0;
snti0=0;
" u$ b1 b9 l  {0 a  M# ufor i=0:5
7 v% A$ R9 N+ |  z7 x: `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));
if snti<=pi/2-x(5)
2 D) h7 R% W3 [2 D7 W   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
3 Q) }% O5 e- [, E- b- Gend
4 w; X/ s8 ]0 N( F! Gf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end


& ?* i7 j: U! K7 N请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
. U4 T6 F6 k9 d3 n) a) Ssyms t;
: A7 f! P. e  sa=4;b=3;
$ R$ O. s; o" O+ H! l+ F/ o( S1 Zf=a*t;
. w" `4 s+ V$ Q0 b9 y6 g% d% vg=[t-b];
4 n+ E) A& z" X7 @# \[x,minf]=minNF(f,[5],g,10,0.5,[t])
! h# _4 @& Z3 l2 J, |  X& Z我将目标函数写成：
! R% P, J5 I  e" S$ s; |6 Jfunction f=link_objfun(s,t,a,m,n)
global L1 L5
% H" S' |) W0 p$ w2 O( N( Kf=0;
% N1 B# l0 y, K" X4 Msnti0=0;
1 _7 B4 n; C! }0 x3 l1 wfor 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);
5 @. I2 N* s! r# w3 N. F& opsai=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
# i! u& q5 O' R4 }1 B; R! X0 J2 Y   bati=fai+psai+acos(m/a);
9 ~0 W0 C) n) ~8 [3 x  ]elseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
end
1 Z  a. T5 }  P6 r5 d( Cf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
! J# V) [- |  B0 G运行如下：
clear all;
global L1 L5
L1=3000;
L5=800;
5 \. o' u$ W' g6 r  _! z% psyms s t a m n;
f=link_objfun;
% l! Y- g, g4 [, L9 n- r( @+ Eg=[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])
6 }9 ?6 s6 Y- Q5 E  Y* m??? Input argument "s" is undefined.

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);