数学建模社区-数学中国
标题:
四杆机构优化设计
[打印本页]
作者:
changfeng1937
时间:
2012-2-25 19:02
标题:
四杆机构优化设计
本人菜鸟一个，论文要用到优化设计，写了一个程序。运行错误提示：??? Error using ==> svd
$ A6 H! f9 c4 a" D! X
Input to SVD must not contain NaN or Inf.求大神们帮帮忙。很着急！！
+ M9 a8 [: H( Q
clc;
5 L$ }" F) |6 g! q T0 j
clear all;
4 @% G+ s& O+ N5 [
global L1 L5
) M7 J6 @/ o, l/ J8 r5 y* K }
L1=3000;
# q8 T$ j0 s, K1 m f1 z' q
L5=1800;
4 J( g/ x; f& H( Y4 K/ c" A5 Y, J
x0=[780;2000;800;800;pi/12];
6 G- K) A+ V# l# w9 J
options=optimset('largeScale','off','TolFun',1e-12);
; {( I" v2 O6 X( s, I0 Z
[x,fval,exitflag,output]=fmincon(@link_objfun1,x0,[],[],[],[],[],[],@link_confun1,options);
. a7 _+ n6 G% ]; O' c
2 X: S6 R( V* v/ N
6 Q6 z3 T4 [0 q) c% d7 y8 y
%目标函数
2 H/ X0 c6 A$ @ r
function f=link_objfun1(x)
; {+ w7 L! R- V; p) J$ }" F0 |- b
global L1 L5
' X% c/ D4 H# Y5 [3 j
f=0;
% t# |* N; e4 W" e# ^
snti0=0;
6 ?! z; o3 D1 |& c( q
for i=0:10
: s# W/ u, \1 ~) j- E4 x- l
snti=snti0+pi/6*i/10;
9 n. y( E* r B& F; G& q0 U
bati=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2)...
2 w+ t9 w* \9 S8 |2 [) X, k! \
+acos((L1^2+x(1)^2+x(3)^2-x(2)^2-2*L1*x(1)*sin(snti+x(5)))/2*x(3)*(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2)...
( }! u* g' }1 A/ q3 r; f9 }
+acos(x(4)/x(3));
$ {) I* r# O1 ~
f=f+(L1*cossnti+L5*cos(snti+bati)-L1)^2;
1 W/ C9 U9 F+ ^; e$ W
end
5 z9 F" x' d: |" T( `% K
3 \6 m: N `4 J+ w* S% Q
%约束函数
8 Q# V" W; h: ]) w) @! t
function [c ceq]=link_confun1(x)
; x9 Z' H2 Q9 N! R4 {1 X
global L1
( B% C: U B* Z, X0 {( t$ _9 k4 e
snti0=0;
% \/ M; K0 y. V8 w) m
bati0=acos((L1-x(1)*sin(snti0+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti0+x(5)))^1/2)...
( B& z3 _1 b3 b3 M: k
+acos((L1^2+x(1)^2+x(3)^2-x(2)^2-2*L1*x(1)*sin(snti0+x(5)))/2*x(3)*(L1^2+x(1)^2-2*L1*x(1)*sin(snti0+x(5)))^1/2)...
/ T( s3 }! |5 h! [$ E
+acos(x(4)/x(3));
# m1 f% H$ u) H8 A
c=[-x(1)
1 b1 y' ]9 M* G4 R8 [9 j
-x(2)
6 [+ L0 f+ T* |* E1 R1 i- V+ C
-x(3)
( j. b1 C, `6 ?; u
-x(4)
" R) d7 R# C f: d& U- x0 Y8 J
x(1)-L1
; d8 m3 k N: o' ?- f3 N% B1 K
x(1)-x(2)
( M6 U* [2 B7 n6 j
x(1)-x(3)
7 q, L7 N. W3 y$ ^6 g9 ~
x(1)+x(2)-x(3)-L1
0 N7 C& e( E- X* \9 N
L1-x(2)-x(3)-x(1)
: ^% u1 ?' q8 i( ^5 t$ K3 W! Z% [$ Z
x(3)-x(1)-x(2)-L1];
# c( _$ j }: Q7 Z3 R) t
ceq=bati0-pi/2;
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5