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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
2 b8 T# F& E, x2 [惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
) s8 K' I4 D& T. t%不等式约束：g
( X  R) \: i  W%等式约束：h
- I/ i1 ]) k4 C" d%初始点：x0
%罚因子：r0
%缩小系数：c
& _2 x8 Z% @: p4 P% y%自变量向量：var
. H) l9 y7 j5 Y! Z0 D%精度：eps
%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf


  K8 x0 l1 C2 {目标函数：
function f=link_objfun1(x)
global L1 L5
f=0;
2 H# S8 [& `7 S) f' gsnti0=0;
for i=0:5
snti=snti0+pi/6*i/5;
( Q5 y2 t4 X, B! ]fai=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));
4 S0 _- e: n  y9 N) c! C! v" R- R1 J) }if snti<=pi/2-x(5)
  k) d: k3 B. J$ x! r5 @   bati=fai+psai+acos(x(4)/x(3));
: }/ l* d, q5 @% Relseif snti>pi/2-x(5)
7 b: o, f$ A( j7 [! |   bati=-fai+psai+acos(x(4)/x(3));
% A0 p. B6 `, ?: a2 Tend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end

+ s7 c* D  Q# x9 a0 g
) u# b- e0 x$ h5 D" L- q请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
( R2 m6 I5 N7 Q3 ~/ _7 tsyms t;
- [2 S! Z5 \  ^8 `a=4;b=3;
. m$ w# j7 E; g: r9 cf=a*t;
1 g0 X7 y5 q; Y* ag=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
, P: v! A7 F' N$ {! d& t5 A我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
( B7 J. u/ p% b$ w4 P8 Hf=0;
  i9 V0 S# K  }4 ?# ?snti0=0;
for i=0:5
1 q5 }% L  S+ I8 k( t" z: `% Isnti=snti0+pi/6*i/5;
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
/ G) \( S6 W/ C/ S$ x+ npsai=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
   bati=-fai+psai+acos(m/a);
end
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
$ V& v: s! }& o+ ^end
) J: ~4 ?8 W9 u运行如下：
clear all;
global L1 L5
L1=3000;
L5=800;
( P+ G2 ~3 o# U; g+ M, fsyms s t a m n;
/ |) [" @! A) f9 O) D6 m& P' g* |f=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];
1 R$ p/ ?4 O6 z( ~6 ~1 v. ^: G& T[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
2 i0 ]! T: |3 ???? Input argument "s" is undefined.

8 T, B( @1 x8 e3 G8 M7 t, W0 f( `Error in ==> link_objfun at 7
- D# U$ q5 o3 k3 yfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);