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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
/ P" z; c) U: B! f; I0 q% ]7 m惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
%不等式约束：g
%等式约束：h
%初始点：x0
" N5 u# o" u; P  _%罚因子：r0
  I% N/ D* A" W/ q- M& F0 v' A%缩小系数：c
2 Q1 x# p" r: f' H9 P, S6 n%自变量向量：var
: d: o$ f' m+ ]1 n/ s%精度：eps
%目标函数取最小值时的自变量的值：x
1 b5 {" s0 ?" b! P( e5 F! Y* ?%目标函数的最小值：minf
# e) ^  a2 f. K$ X  d9 q  r

& w1 w& u' R3 b- k1 e: c目标函数：
; ^! u9 D% j1 J/ [; G4 H+ w9 M' vfunction f=link_objfun1(x)
' X" a+ k$ G$ }) ?# Aglobal L1 L5
f=0;
# u8 l$ u2 `, N, _snti0=0;
9 i5 y2 a- w- r5 wfor i=0:5
1 }9 E6 \. t% tsnti=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);
& W6 O& z! o5 A9 ?$ apsai=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)
   bati=fai+psai+acos(x(4)/x(3));
1 {+ t' y) n6 ~. Z" S, Relseif 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;
end

0 z. M+ }2 M3 u; ^
- }( p- l/ {+ J! j6 [请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
% M" l7 Q" A, E6 e7 A$ S/ t" Z运行如下：
syms t;
a=4;b=3;
f=a*t;
g=[t-b];
2 O8 o. R! h" f9 ?6 r4 {/ O' @[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
0 a) D) d2 p/ s. Qglobal L1 L5
- z+ l( t1 g  o3 l) d# D# gf=0;
  I7 a  h: J) W& D- |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);
+ a: Z  [: `' [9 ]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));
+ Y% f! E2 u( F4 g4 K* r$ z+ T1 ~if snti<=pi/2-n
/ d6 v; `- P3 J8 [+ N" X   bati=fai+psai+acos(m/a);
; _. W  t5 k$ selseif snti>pi/2-n
   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
运行如下：
/ l8 m' c% R) J: ~0 m2 Hclear all;
3 k6 U# O+ w) ], e! Z8 {global L1 L5
4 U' m; e& ]5 B- r, w8 LL1=3000;
! ?5 u8 ]# ?% W$ D/ P* \L5=800;
) \- t1 i; L1 z/ isyms s t a m n;
* H9 s! X9 W+ J, rf=link_objfun;
g=[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])
??? Input argument "s" is undefined.

Error in ==> link_objfun at 7
! u5 j- [: x3 ~8 }1 h# hfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);