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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
5 e0 n* `6 L7 P* @2 C1 f! m% e惩罚函数调用格式如下：
, q/ `! k; @+ p, t* @, X3 {function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
2 q( f9 g: ^0 n1 j* H1 {" @%目标函数：f
%不等式约束：g
%等式约束：h
( @0 B% A/ R8 r3 {3 a% H%初始点：x0
%罚因子：r0
" f2 G; x: j# I) Q' r4 g8 q; v4 @%缩小系数：c
" F3 R7 A5 @( @! e/ e! B* d! b1 W2 M%自变量向量：var
%精度：eps
%目标函数取最小值时的自变量的值：x
: t* ?, {4 o8 r) J: e%目标函数的最小值：minf


8 }( F5 U1 @2 X& r! n目标函数：
function f=link_objfun1(x)
global L1 L5
7 ]; @& j+ _; Gf=0;
snti0=0;
, D1 N6 f2 V! i: d/ ]9 n. Sfor i=0:5
snti=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));
6 O4 T* h* S0 `- I' J" t& Oif snti<=pi/2-x(5)
* f# e, K6 [) @0 }   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
9 u' }0 ^$ E& }0 O! i4 r   bati=-fai+psai+acos(x(4)/x(3));
end
& Y, G. A. S# i/ Ff=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
! T, J5 ^" N+ D! k
! y! P; i* ^3 y8 C/ o
/ _& t* [0 L7 b# E请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
2 @: i7 x0 G8 d3 j5 V! v, d, `  Y# V运行如下：
syms t;
/ R( W# Q5 y  P- N- pa=4;b=3;
% O1 Q  D( D9 Y, Hf=a*t;
# Q+ m; w7 |/ j8 Gg=[t-b];
" f. G. g, U+ e* [: Y[x,minf]=minNF(f,[5],g,10,0.5,[t])
我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
- ~# c5 B6 t) n$ K/ s8 d+ a0 bf=0;
snti0=0;
3 F, x( s$ F$ h8 _for i=0:5
4 t5 t, m# U- T- K! Esnti=snti0+pi/6*i/5;
9 U# t+ z1 h! ~% zfai=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));
if snti<=pi/2-n
   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
, j$ e* a: q0 D2 [* l   bati=-fai+psai+acos(m/a);
9 X( ?: `) A8 _7 a1 Q. {% |8 l) Tend
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
end
运行如下：
1 D% X# C6 I/ T  P& v7 {9 lclear all;
global L1 L5
L1=3000;
L5=800;
( z! ~! G2 W  e$ c. ?# @7 \8 psyms s t a m n;
/ O1 Z+ U$ h0 Of=link_objfun;
. t; c# ?7 L: [5 eg=[s;t;a;m;n;L1-s;t-s;a-s;a-s-t+L1;s+t+a-L1;s+t+a-m];
% V$ Y$ ]6 T- l8 B+ L[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
??? Input argument "s" is undefined.
6 n1 Y; }1 C; N9 Z4 M7 |. l- R
Error in ==> link_objfun at 7
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);