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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
/ G5 z0 T. C+ y; M3 S$ Ufunction [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
1 }/ D- [- H" r2 C) z7 G/ ?# {%目标函数：f
/ f6 n/ s7 `' _5 a& U%不等式约束：g
0 d9 K/ ^1 u! [" y4 i%等式约束：h
%初始点：x0
, U, D! ]. w. d$ D- g: x%罚因子：r0
%缩小系数：c
%自变量向量：var
%精度：eps
* w/ P. ^5 |' S5 x% n; S%目标函数取最小值时的自变量的值：x
( l, o& B* p! o- L8 x- j# M%目标函数的最小值：minf
7 ?3 o  q$ {6 p$ @4 e- M# s

  t; z/ r" Z- I* ~目标函数：
0 I. l6 L; `! K8 zfunction f=link_objfun1(x)
global L1 L5
f=0;
) B8 |+ j- l& d! V6 i0 {$ E8 Csnti0=0;
for i=0:5
snti=snti0+pi/6*i/5;
6 p5 ^: c5 D2 _# Y: n$ Gfai=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));
if snti<=pi/2-x(5)
2 v9 h8 z) A7 n0 n) E* x5 }   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
+ q: T7 y% B* _; A  o- U8 H   bati=-fai+psai+acos(x(4)/x(3));
3 u3 z' r- s4 f# i5 A: y9 M: H5 }8 cend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end


请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
syms t;
a=4;b=3;
8 i8 t& L( k5 a) hf=a*t;
) J2 |9 |. _1 U$ S7 ig=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
, t4 J( `( y8 T3 D9 C9 tfunction f=link_objfun(s,t,a,m,n)
global L1 L5
f=0;
snti0=0;
2 q# |  Z9 N: ]# P! Rfor i=0:5
snti=snti0+pi/6*i/5;
# f6 w: ?1 {7 Z" `, Gfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
2 w& o% L( c3 J5 B- N5 G/ W' ypsai=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));
6 ]5 M  H0 @) {% `9 s! _2 rif snti<=pi/2-n
% V6 ]% B- j# b- w- a8 y4 O   bati=fai+psai+acos(m/a);
' ^) J5 a/ L# I) helseif snti>pi/2-n
- Z2 B  {7 h  \, h% {4 g0 K- H" ?   bati=-fai+psai+acos(m/a);
) Q# ?" M+ g$ N( G( {  s8 Bend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
  m$ n& M7 s7 C4 R, H  B* s2 E/ Zend
运行如下：
clear all;
global L1 L5
8 J  M; t8 w5 d' R6 D9 A* _L1=3000;
8 c& e. E; D9 QL5=800;
syms s t a m n;
* u. g7 X( g9 w3 }8 kf=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.

2 _5 ]7 p% o9 Q3 a+ U2 `2 GError in ==> link_objfun at 7
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);