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

changfeng1937

想用里边的惩罚函数程序来算，我的目标函数是用for循环写的。不会导入。
3 E/ k9 Q( c& l惩罚函数调用格式如下：
function [x,minf] = minMixFun(f,g,h,x0,r0,c,var,eps)
; f& o2 s3 k1 y$ H& {%目标函数：f
%不等式约束：g
( j, F6 X+ m( R( l  y1 F%等式约束：h
%初始点：x0
%罚因子：r0
%缩小系数：c
%自变量向量：var
/ r6 F  J9 `3 r. D%精度：eps
%目标函数取最小值时的自变量的值：x
( k2 T9 M/ A' G- a$ k. K" [%目标函数的最小值：minf
3 l- ~  S( `5 t7 h" `
7 ]6 @+ n+ G2 R) _' C) W$ i
* F. m2 r: I% e0 a# V" P* s: L( x目标函数：
, R+ Q' M0 z# I1 e9 v, a8 Wfunction f=link_objfun1(x)
global L1 L5
3 [3 B: X% X; e. Z+ L' v# rf=0;
; \4 Z9 Y; _3 S% R: Jsnti0=0;
for i=0:5
' p4 h/ L% K' Q' O$ {$ g1 F% [% csnti=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);
' V5 U( w* ~4 b5 Y. Z) W/ y5 T$ r: j+ Kpsai=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));
2 K) C- d/ Q& tif snti<=pi/2-x(5)
   bati=fai+psai+acos(x(4)/x(3));
elseif snti>pi/2-x(5)
" K, V, O$ u! y" ~& ~' S& H   bati=-fai+psai+acos(x(4)/x(3));
, T( a+ j: `+ L# |1 nend
- X2 F* {$ O3 N4 d& R& i9 If=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
! y. F6 J7 _) U& w* T8 T: O* _end


6 ?: r) E# Q0 F  Z$ M4 R请诸位帮帮忙，我的体力不多。

厚积薄发

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

changfeng1937

我懂的不多，生搬硬套。书上的例题
2 ?9 p( D4 h, v/ Z2 S+ k. H1 Q) U运行如下：
syms t;
a=4;b=3;
, ]; _& e3 ]; M( Af=a*t;
3 h# N% ~! f/ F, Wg=[t-b];
[x,minf]=minNF(f,[5],g,10,0.5,[t])
0 N" y1 S  M+ R我将目标函数写成：
, y7 q. B" `% E- J: Ofunction f=link_objfun(s,t,a,m,n)
! }- D! z& B8 _3 ]" Xglobal L1 L5
f=0;
/ X; M1 j& y* B8 d# r' Y- F# d$ Osnti0=0;
8 q- G, E, G% t9 ]for i=0:5
snti=snti0+pi/6*i/5;
fai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);
( X5 F+ J0 }1 [( G, w+ i: \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));
6 G; }5 d* n+ Q2 x( qif snti<=pi/2-n
   bati=fai+psai+acos(m/a);
) J2 P- }% s0 {3 c- _elseif snti>pi/2-n
8 ]* T% W/ y2 L   bati=-fai+psai+acos(m/a);
" T; o( g4 Z+ e, ~8 F( n/ send
f=f+(L1*cos(snti)+L5*cos(snti-bati)-L1)^2;
0 g' Z5 J) E4 ]& X) Z( Jend
运行如下：
clear all;
global L1 L5
! e" W; X4 E% {$ k0 I! N; [L1=3000;
; G7 x/ G% t4 ]L5=800;
% m: I4 }8 t& Q3 S1 ^( X3 Csyms 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];
/ X4 t, D' i; b4 \: P# f3 E[x,minf]=minNF(f,[780 2000 800 800 pi/12],g,15,0.6,[s t a m n])
. G4 \8 m- r; O  G??? Input argument "s" is undefined.

/ b+ _& a1 X, T$ `' B" CError in ==> link_objfun at 7
: a4 X% G; g! Sfai=acos((L1-s*sin(snti+n))/(L1^2+s^2-2*L1*s*sin(snti+n))^1/2);