数学建模社区-数学中国
标题:
四杆机构优化设计
[打印本页]
作者:
changfeng1937
时间:
2012-2-25 19:02
标题:
四杆机构优化设计
本人菜鸟一个，论文要用到优化设计，写了一个程序。运行错误提示：??? Error using ==> svd
4 H, l- X5 u( f& z
Input to SVD must not contain NaN or Inf.求大神们帮帮忙。很着急！！
6 M T" F% W4 s+ |1 m* B5 F4 V; a6 F
clc;
; v" ]' w# p6 d5 H2 |7 ^( k+ Z
clear all;
: s" Y4 }- H# r8 Q
global L1 L5
* |# @ W: R* ^* m' k7 v u6 A7 e% E
L1=3000;
% G. i; J( ?: i% y
L5=1800;
. j+ n4 T0 ^, ]# J) ?
x0=[780;2000;800;800;pi/12];
: R8 E; K% l4 p
options=optimset('largeScale','off','TolFun',1e-12);
/ j0 ^# ] Z: e9 M: N9 k1 Q
[x,fval,exitflag,output]=fmincon(@link_objfun1,x0,[],[],[],[],[],[],@link_confun1,options);
p9 b2 n' q5 ]
, {, n& f& U2 n8 Y
; m5 H0 c/ X6 z' J
%目标函数
6 t! |. A2 I& B/ U* l2 O$ U/ e
function f=link_objfun1(x)
. B* I5 Q& a; |% V8 A$ ^% n
global L1 L5
( U' I: A% M) d7 O, K
f=0;
2 p! c( H3 m. P5 U# U" h& x6 q
snti0=0;
; W: _6 N( {; ]: e7 q) b9 x8 r
for i=0:10
. F2 [' A& [) Z5 O8 J" |
snti=snti0+pi/6*i/10;
# w6 R- e* s2 s& 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)...
6 b& A0 W7 A7 H% C9 V5 N- s
+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)...
' w; O( D, ~8 W# i" f# q* ^( z
+acos(x(4)/x(3));
+ V) K; q7 F( ?' h; p: k% ^
f=f+(L1*cossnti+L5*cos(snti+bati)-L1)^2;
; ~! x6 y' A# G2 H1 r6 Z* m: s' P
end
1 G. r; M: z6 O
, d5 a. C1 r0 V+ |% `
%约束函数
& a& c6 c7 C& ]# J* e
function [c ceq]=link_confun1(x)
5 w6 U) K7 C7 M) t* T* j; W; b
global L1
) ~/ A1 ^5 ^, ~+ s g6 x3 W
snti0=0;
- [9 f o1 z2 q8 K3 Z
bati0=acos((L1-x(1)*sin(snti0+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti0+x(5)))^1/2)...
' v3 J3 d) h7 d. h9 X8 }
+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)...
{; { J- l2 f2 U" Z8 V
+acos(x(4)/x(3));
+ @2 k" x- K% U5 q0 F
c=[-x(1)
$ [; G$ l G0 x5 b5 f6 ~$ e
-x(2)
1 C) ]* B% b% x4 H% }
-x(3)
, I* W) \7 g3 ?
-x(4)
- a" M y0 W T, U- O9 Z0 ?
x(1)-L1
! @7 Q5 H. @) b
x(1)-x(2)
# ^2 R9 ^* @$ Y
x(1)-x(3)
x4 V2 y( D9 A7 E$ l, d% d# x
x(1)+x(2)-x(3)-L1
e" G, d! s' H5 |% {* m R8 _- Y
L1-x(2)-x(3)-x(1)
) g1 q8 |: f, |! @1 h7 v' H K* d
x(3)-x(1)-x(2)-L1];
0 c4 f# p7 t& e. c; Z$ a4 q0 V
ceq=bati0-pi/2;
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5