数学建模社区-数学中国
标题:
四杆机构优化设计
[打印本页]
作者:
changfeng1937
时间:
2012-2-25 19:02
标题:
四杆机构优化设计
本人菜鸟一个，论文要用到优化设计，写了一个程序。运行错误提示：??? Error using ==> svd
8 j' j4 k% c6 u, H, m0 W7 P
Input to SVD must not contain NaN or Inf.求大神们帮帮忙。很着急！！
. ?6 t' ]% u7 N( B! `
clc;
4 v' g7 W/ V/ `' ~7 \: ^
clear all;
0 a4 M. @$ k6 l% k7 {. |2 }0 z
global L1 L5
( M. ~0 `& X. ? k% I w4 `4 w
L1=3000;
$ e) m0 w2 g7 Z y: q% R: w% ]. k* R4 R
L5=1800;
: N% x! @+ t/ Z* z% A3 _+ O# Y
x0=[780;2000;800;800;pi/12];
/ Z; h. k1 q+ q" J8 ~+ w2 h
options=optimset('largeScale','off','TolFun',1e-12);
: @' m8 b2 C0 b3 @6 E* ]
[x,fval,exitflag,output]=fmincon(@link_objfun1,x0,[],[],[],[],[],[],@link_confun1,options);
. Q# |& F( {% q" h5 N
- r9 |/ W. ]9 n* _4 x. g4 r
7 t# L7 Q$ l: U2 i! G" _5 @* M+ p
%目标函数
) g' ~ I% d, F0 c& _
function f=link_objfun1(x)
* u$ n, Y, ^1 q% B6 \% j& O; L
global L1 L5
% m/ C2 z. y2 S6 O
f=0;
% {) I8 j4 A7 V: e
snti0=0;
% g$ b: [; _9 B5 u M# d+ s
for i=0:10
6 P% y, }0 U" P1 s: r' k
snti=snti0+pi/6*i/10;
y2 b" ^" g; a8 ~
bati=acos((L1-x(1)*sin(snti+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti+x(5)))^1/2)...
; `5 A4 |0 s, H5 @
+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)...
* e' E$ g; i: |) y& B; `6 R
+acos(x(4)/x(3));
$ {) _+ [) \1 f2 N- B/ e
f=f+(L1*cossnti+L5*cos(snti+bati)-L1)^2;
6 y/ z2 b) n0 U! p' V
end
. _- [# k( Q$ u! c( X2 X
7 W) ^3 j0 B0 q B
%约束函数
3 I4 q; l4 |7 S, p4 I t6 g1 C
function [c ceq]=link_confun1(x)
: O, E& s D/ V9 \2 G2 u8 d% y0 g
global L1
+ n+ w! P% e0 Q3 @: g/ q
snti0=0;
! G, q+ F* f. t* a
bati0=acos((L1-x(1)*sin(snti0+x(5)))/(L1^2+x(1)^2-2*L1*x(1)*sin(snti0+x(5)))^1/2)...
2 h ?& _+ b7 C! `. T* h# R
+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)...
. f4 o7 W3 {1 P3 m; G/ g+ {+ [
+acos(x(4)/x(3));
/ Q: M; e# O6 p! w& s9 V$ ^( r9 ]
c=[-x(1)
; d F4 ]6 f/ H# E* B* `5 w, `, Q
-x(2)
5 t1 z1 C7 T, a& s" X' v& ^
-x(3)
! l% I- o ~1 ?& U) ?
-x(4)
* y! H$ v; n+ f. g( s( M9 E& s! f
x(1)-L1
7 \) C1 R3 S2 y3 p) L+ [
x(1)-x(2)
) v8 @5 N( u' ]$ c) @" Q
x(1)-x(3)
M2 S% o* t3 u3 R# A7 b
x(1)+x(2)-x(3)-L1
4 J% o% z y& Q+ b o" s
L1-x(2)-x(3)-x(1)
! [4 R4 A% C) |
x(3)-x(1)-x(2)-L1];
$ U; Z$ @3 ^2 t% @! {# n* }3 E
ceq=bati0-pi/2;
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5