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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
& ~/ m) X5 v( D- q惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
( E+ m- X9 b( L  L1 y2 r%目标函数：f
%不等式约束：g
%等式约束：h
. {  M+ q/ o+ |* v5 @, P; a%初始点：x0
%罚因子：r0
%缩小系数：c
! T# ?5 o5 t( Z; i) S* V, z, }%自变量向量：var
%精度：eps
( X/ @& @) O$ x: `%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf

) z8 W" \7 [/ o) a8 O" }% W! m9 }" P
目标函数：
function f=link_objfun1(x)
% m8 z8 h2 M2 b1 t8 B% wglobal L1 L5
f=0;
' r9 X, D* U' J' g6 L5 ?snti0=0;
for i=0:5
snti=snti0+pi/6*i/5;
: i* p& @* X" N/ G( g4 B& ifai=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2);
) l  `: q" ~4 X' Lpsai=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));
( g) O0 |: t8 O" [; O' aif snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
2 P- _  |! T. M. melseif snti>pi/2-x(5)
: D6 A9 X# u2 F% |* d   bati=-fai+psai+acos(x(4)/x(3));
4 L& r8 O! w/ Q: U4 m# cend
) m" o$ B) t' ~: Q; P5 df=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
* r9 |1 \% g8 t- J; \end

5 G" W9 r- z4 V& p6 Q
1 F. S1 F4 x8 p+ c+ Q请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
% _6 _! I% t0 d; b% v$ j  j运行如下：
syms t;
9 h' q( j2 g7 O8 r+ S  ^a=4;b=3;
8 E9 G4 Z$ L! `+ z( ^+ Ef=a*t;
g=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
/ O$ y, T6 L' n, e我将目标函数写成：
- Y$ R+ v; h, H: Xfunction f=link_objfun(s,t,a,m,n)
: V6 J) c+ S- S# U0 o# vglobal L1 L5
/ R  l- G1 I1 ^$ [0 _f=0;
5 h! T0 K3 G/ ~* ?snti0=0;
for i=0:5
  d- p: b) ?5 `  g  }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);
) T  S* V9 o! L7 U. r! [( Fpsai=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));
3 r- B0 W9 N) G+ ?) I3 Mif snti<=pi/2-n
   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
/ g. K+ L( Y9 Y$ |   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
/ ]3 X% E# `4 w7 b- B1 Y运行如下：
& y- p  Q: G0 A, w: ^7 Pclear all;
global L1 L5
( ~; N2 [! e7 i# AL1=3000;
% U! `( H+ b  F  k4 [: N# Z$ P/ R1 vL5=800;
syms s t a m n;
3 ]0 j. B% i9 of=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];
* T& A6 T5 s" J% n[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
% h6 Z. \+ F% T3 {/ P; j3 h??? Input argument "s" is undefined.
$ j8 i7 b8 y$ L" C: M
6 Q; H5 o$ J. @0 z7 w' zError in ==> link_objfun at 7
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);