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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
8 B+ j8 ~$ G- a# [7 U惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
% Z! \9 Z( `- l, ]/ Z; A%目标函数：f
, q% W% K2 P$ i%不等式约束：g
' ?5 e4 y% N$ ]4 D%等式约束：h
0 [3 _! R9 H3 U; }- F8 ?; x%初始点：x0
%罚因子：r0
%缩小系数：c
%自变量向量：var
  |" M, @" ~% S5 |' g%精度：eps
/ E' m( ~) r3 y# ]%目标函数取最小值时的自变量的值：x
%目标函数的最小值：minf
! B5 ?6 ~$ y0 f% z, m6 X

目标函数：
- ~# T; b# }/ pfunction f=link_objfun1(x)
" L5 v3 x" u( }/ l; [# ]global L1 L5
f=0;
snti0=0;
for i=0:5
" _' m" a7 [! B. V* psnti=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));
. v6 J6 z: Z6 s8 G+ \2 p6 mif snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
! A7 P. z- K: X9 Felseif snti>pi/2-x(5)
/ E8 R2 l+ g2 E6 B$ U" D1 O   bati=-fai+psai+acos(x(4)/x(3));
end
! n' c: C6 }* L4 D) \f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
& a& T1 J  E6 Q+ c% _6 u' Bend


/ ]. q/ y9 }' i8 c. v8 J; r请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
运行如下：
* N" P# o( B: b' z8 C2 Zsyms t;
a=4;b=3;
( q: A/ b, b) j+ B) [- C+ W+ Ff=a*t;
g=[t-b];
* X, L- `; [6 u' @  ~2 Q[x,minf]=minNF(f,[5],g,10,0.5,[t])
) i  B7 @# w) t1 {7 `我将目标函数写成：
function f=link_objfun(s,t,a,m,n)
global L1 L5
f=0;
9 J, z& D: n& n8 Usnti0=0;
2 \* c6 _5 ]: R+ u1 F' [5 e( N4 |! @for i=0:5
" @  D+ X6 S' R5 ?' Qsnti=snti0+pi/6*i/5;
4 L' R$ E* e: W8 E6 n- P/ N: W- gfai=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));
& N4 R4 M% R  D: eif snti<=pi/2-n
! r0 ], j* j( r7 e   bati=fai+psai+acos(m/a);
elseif snti>pi/2-n
6 m, r+ J  g8 k' {7 l   bati=-fai+psai+acos(m/a);
end
, s9 p2 C+ r0 }6 g% R7 Vf=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
5 b" B1 Q- B# D' e7 R6 }end
3 _  I, c* X+ c4 o6 `+ u- M运行如下：
" D% @3 _" e0 _clear all;
& Y. O1 b5 f5 I, U/ Lglobal L1 L5
& S# e0 }* C' p& s: PL1=3000;
' ?. _+ ?" C. ~, T* e; QL5=800;
5 x# ]9 Q3 Z& Z6 W. R$ Jsyms s t a m n;
7 N! ?( R. f1 b+ a" S6 m5 A6 }f=link_objfun;
* X8 v2 z$ Z/ I) q+ O  R- v3 U( I; C9 ?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.

# G3 H# W0 V: I% k- {8 @. }2 dError in ==> link_objfun at 7
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);