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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
4 p& u; b' \4 C) a1 S" [7 y惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
7 o7 b) F& r& W0 i/ q; [( N7 n9 @* i9 D%目标函数：f
7 ?( z1 D) X) u5 L0 {2 d: J# C%不等式约束：g
%等式约束：h
# n" Q5 v* y6 y+ z%初始点：x0
%罚因子：r0
%缩小系数：c
%自变量向量：var
%精度：eps
+ T$ y8 {: ^4 g% t%目标函数取最小值时的自变量的值：x
9 T6 _2 X6 v5 |% w; R%目标函数的最小值：minf
  b: B+ I+ e2 v: F$ ~; @
0 I8 A: V1 u; ]4 W. v' [+ o9 `* U
目标函数：
function f=link_objfun1(x)
# @: p2 r) z6 dglobal L1 L5
f=0;
snti0=0;
; t& O' w! _. i8 dfor i=0:5
: Z# {$ v, P& Z5 psnti=snti0+pi/6*i/5;
" I6 ~1 }% v) X  V  Pfai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
# }+ w6 h2 X% ~, o4 @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)
7 s' P; W; R7 V0 c. C/ F$ V   bati=fai+psai+acos(x(4)/x(3));
, Y+ T5 k* A) \" ?# E5 _$ Zelseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
" R/ \/ R  D* Hend
* [1 B' y" u* `! A. uf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
; N9 l1 x3 p$ K  L; `7 o1 j/ V9 Aend
# O5 E% V0 u8 {9 F0 |" j+ h+ q
. ]" O' h  q( x& f
% L2 l( A1 R9 B3 o3 g% e$ r4 X请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
. k: v3 U$ x1 E- \, A7 z) y! c运行如下：
syms t;
. c2 R7 x$ K2 M: y6 Z' Xa=4;b=3;
f=a*t;
) D  e; C+ T6 S8 S/ ^* Pg=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
& ?# U, H! X+ Y6 _" cfunction f=link_objfun(s,t,a,m,n)
7 L* U" a& n: x/ a" l! a, Mglobal L1 L5
2 k& B/ T( z0 k0 C) i9 {) c. |/ af=0;
7 R5 f' s6 [2 ksnti0=0;
for i=0:5
( N7 V5 R7 F7 e* vsnti=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);
elseif snti>pi/2-n
8 e3 u) f. e4 B7 i% `' y. K. N   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
& q; r" a/ M; N  q* u& K8 Z. ~( Jend
运行如下：
$ ]3 a# ^% b8 f5 t) X8 ~, lclear all;
global L1 L5
) o& ~+ H# T6 E9 o9 B. r- iL1=3000;
  h, p; a0 A8 U( {5 y; {; ZL5=800;
syms s t a m n;
& K3 |: Y; s% o, ^# ]0 ]f=link_objfun;
1 G" {; Y- i$ ~( K+ bg=[s;t;a;m;n;L1-s;t-s;a-s;a-s-t+L1;s+t+a-L1;s+t+a-m];
: m% D" }- y: W% J5 [" ][x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
4 x- h( P7 d. V??? Input argument "s" is undefined.

# @) T6 i- w1 _( u& U0 lError in ==> link_objfun at 7
/ B& D3 \, [6 D- C) d" y5 xfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);