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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
%目标函数：f
* y1 x: C2 F% W/ Z; L/ S%不等式约束：g
%等式约束：h
%初始点：x0
9 `; e+ y8 A8 R. U7 _%罚因子：r0
' f2 |* s5 j. I3 s+ o" r& k# l+ H%缩小系数：c
%自变量向量：var
$ a. Q2 f, i# v( ~%精度：eps
3 j/ T& N. f" x0 g& M6 o%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf


. E% m8 m' d4 q: c& T" t目标函数：
: L6 @9 n2 {4 s/ h( ?% ?function f=link_objfun1(x)
global L1 L5
f=0;
9 q) `0 D8 t( v! B& r5 f1 w, ysnti0=0;
+ q- W, I0 c; S( mfor i=0:5
6 W- l0 u0 q2 A5 M$ Wsnti=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);
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));
! J( E) s  j7 I+ s7 uif snti<=pi/2-x(5)
# ]( E/ j# y5 k- E0 Z3 ^   bati=fai+psai+acos(x(4)/x(3));
! C1 W/ h& P& G" X* Oelseif snti>pi/2-x(5)
   bati=-fai+psai+acos(x(4)/x(3));
2 a# A$ m0 K* N  e" U8 ~5 L( t% Eend
( J2 A$ O; G1 J! j" P9 b# a4 \f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
5 x' ]- O. T: k, Eend
6 T, C( ~4 @) P) _
* ^/ m* i, H& g  s9 h
& A) m: P2 I) H3 H% L8 n$ F. x请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
! ], s5 b' f3 {7 p6 b2 v运行如下：
syms t;
a=4;b=3;
& y0 v- p+ r+ s2 H9 c4 f8 ^& _f=a*t;
g=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
% r, u! E! o% ?' b" b我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
f=0;
snti0=0;
: d1 c' W4 D7 i) }: K( yfor i=0:5
+ d+ P4 H/ \! u6 E1 O: v8 Tsnti=snti0+pi/6*i/5;
* K( g( n! c4 v, v+ c* }- ?0 V% Nfai=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));
3 I& \% y0 U0 X3 m3 A! {0 Dif snti<=pi/2-n
   bati=fai+psai+acos(m/a);
; a4 a8 @( b' Z. E' s+ O' @2 uelseif snti>pi/2-n
( L- @' D0 O3 p0 q( x' A% a   bati=-fai+psai+acos(m/a);
: U1 d2 N7 X" y$ Xend
  \9 E/ |# |2 P# J9 ]0 uf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
! ?2 \( T8 N% S5 F4 A" }end
% D& O' F) E/ f) b0 i运行如下：
clear all;
. A! R! ~, ?& B8 hglobal L1 L5
L1=3000;
L5=800;
syms s t a m n;
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 B* }7 ?7 n0 g2 `9 y[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
??? Input argument "s" is undefined.
! C( t( t+ k' |" T( a+ }" p
Error in ==> link_objfun at 7
5 ^0 {! H- f$ S/ \' Xfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);