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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
8 Z# ~9 s6 e/ c惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
- Y$ V( X7 l4 X%目标函数：f
8 p' _. c1 |" k6 b, i+ w0 ]%不等式约束：g
%等式约束：h
' @" n: x" I1 E# T1 n%初始点：x0
& @7 {' b' x$ O" ]%罚因子：r0
%缩小系数：c
%自变量向量：var
%精度：eps
%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf

6 s& `+ p- {0 b3 r1 q
$ F* P  w+ B# Q/ I目标函数：
function f=link_objfun1(x)
global L1 L5
3 M3 ]3 F. H: M) q3 b  m1 t  Z3 nf=0;
) U; P) o+ [5 n: N- T) y: ksnti0=0;
3 S+ {/ s0 g) w: M5 f. h: Afor 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));
  j/ H# T* H/ x" Q% o7 _; j( I5 Kif snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
  C( N# M4 v8 y% T, l7 ^# d# `end
$ Q" M4 w% r- A) I5 e9 x. W$ T
0 P) E. R$ A6 o3 w) k, N9 V
. a; h' V( ?. ^# C2 W" |# i% S$ t# H4 X请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
7 h) Z& r" b! K( g9 u: Z& J运行如下：
& d8 w4 ~1 S7 m3 g. w; Bsyms t;
a=4;b=3;
8 p4 Q1 L" j% m) O) }. uf=a*t;
& t  X% X; i: R3 [g=[t-b];
- `7 B: g! g5 x+ [- f[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
. @6 h  L& a4 ?& K6 a; ~- n" Bf=0;
snti0=0;
3 l; y+ {) q0 o# \6 {( n9 efor i=0:5
/ ]3 G- ^0 e( ]3 r  E+ Y3 Fsnti=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);
% ?: b; W- l! s0 f: w; p* Nelseif snti>pi/2-n
. q4 @& T9 ~8 P# ]5 d2 J   bati=-fai+psai+acos(m/a);
0 s% H5 [" d! D4 e& r+ D2 l: O/ |( _* Uend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
9 h9 n! ~* L( `* @2 G% Aend
运行如下：
clear all;
global L1 L5
, }& V) R; z, g, }$ B* ]7 V# iL1=3000;
L5=800;
. Y# v' k5 B7 F/ U6 L1 Bsyms s t a m n;
f=link_objfun;
# e1 t- ]! U/ Y5 _$ Lg=[s;t;a;m;n;L1-s;t-s;a-s;a-s-t+L1;s+t+a-L1;s+t+a-m];
5 O. M! [- I/ d9 b[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
8 I: O$ D! z' M( W* D??? Input argument "s" is undefined.
! P4 `( E9 I6 P7 f
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);