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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
4 R+ m1 x7 ~* I3 r6 i, A1 o* Z惩罚函数调用格式如下：
/ i# z0 V4 d; L) g4 ufunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
8 ^7 r5 r/ Z, x%目标函数：f
%不等式约束：g
%等式约束：h
7 Z) g7 E: k6 {5 M4 l%初始点：x0
2 Q0 U% d! r; m' J0 n3 h%罚因子：r0
%缩小系数：c
%自变量向量：var
%精度：eps
%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf

/ C, t! k9 q7 }$ x: m
' \, q/ \0 L6 c3 j目标函数：
7 _$ F4 n) F2 k& W$ g1 x$ a4 pfunction f=link_objfun1(x)
! k2 N- G3 Y8 i* j1 Oglobal L1 L5
; i4 F- H; |# p3 T6 Af=0;
; F& U+ q5 |- e3 w% e4 Gsnti0=0;
for i=0:5
snti=snti0+pi/6*i/5;
2 f! B6 f' w% G& cfai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
5 F3 A! y3 W/ u" w& j! r$ V6 A7 hpsai=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));
2 _/ N7 D* X. |$ S/ Uif snti<=pi/2-x(5)
7 ~4 @, W/ A! G4 P( X   bati=fai+psai+acos(x(4)/x(3));
) H9 m( U- L4 F+ A( {' e- S' \elseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
" X( v0 G! I% F/ ~/ Tend
$ `3 ~! k6 {2 K! A# M% C& Q5 J2 hf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
, z+ I8 S+ }$ a# d: q  S; w9 H& eend

* {2 k' {* W4 l) F7 t1 R/ H/ {8 [; u
4 O  u" i( |/ {" |1 |( S请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
( M. F5 N% W7 `$ `5 E8 |& _; }运行如下：
syms t;
$ t; v0 e5 D9 _8 S; E0 _a=4;b=3;
f=a*t;
* ]! i: U* v9 j* d7 j  }# @/ Lg=[t-b];
7 |$ w* W8 |* E0 A[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
+ O% z. X$ Q, b$ \; Pglobal L1 L5
3 v9 i( U3 L( o1 j$ ]8 ~f=0;
snti0=0;
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);
4 O# R% @; i& ?0 n8 K  R, h$ Spsai=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
4 p3 s4 l/ h0 n6 A: |   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
4 b1 N! C2 Y$ O9 A/ |! G) G3 ]end
4 p4 b, L3 Q% B# R; u. a1 K+ `f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
" B0 \$ N! l5 R. i1 N8 |+ L# Fend
; H6 U4 Z/ t- c运行如下：
' F% M+ J9 w4 h" o. \" n% \1 Iclear all;
4 U( m' A0 K' p7 |1 E- i. \global L1 L5
L1=3000;
; V4 |8 N2 ]# [6 KL5=800;
syms s t a m n;
f=link_objfun;
. f& P) h$ T9 U8 N2 Rg=[s;t;a;m;n;L1-s;t-s;a-s;a-s-t+L1;s+t+a-L1;s+t+a-m];
% w: b- Y: P: `" v# S[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
??? Input argument "s" is undefined.
* i4 A( m' R6 i! h; }' |; K
; a: t) H9 w$ |5 g/ o$ xError in ==> link_objfun at 7
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);