数学建模社区-数学中国
标题:
matlabd的fsolve函数求解方程组
[打印本页]
作者:
追梦人1101
时间:
2013-1-26 19:23
标题:
matlabd的fsolve函数求解方程组
function fun=Untitled6(x)
' m) a5 W3 E- p+ x+ \) }+ G
theta=pi/18;
1 @, ~0 ]3 T$ M* l; n; X
k=0.97; %k为万向节的传输效率
+ v) {# E. ], B7 U
R=13; %R为轮子的半径
. i" L! Z" \' M
r=4;
( F* F& N* G1 t4 X. O$ s
D=250;
# ~5 {6 ?, Y. H# q
fr=0.2; %fr为轮子的滚动摩阻
1 h) M- ]2 }4 W ~4 H8 Z3 @2 q
mu=0.4; %mu为轮与有机玻璃的摩擦系数
: A2 u& l* N4 n- @& d* Q
Wf=5; %Wf为电线所产生的负载
4 Q- }. V) A# s
G1=5;
9 M! s3 c0 E R2 O* H' }5 n: G0 E7 h
G2=80; %G2为电机，导向装置和加工头的重力和N
0 C/ E( h4 X+ \0 f) j
alpha=11.5*pi/180;
; x0 z. a4 O* O0 |( j! Q9 s
beta=33.5*pi/180;
0 E: ^. C$ v+ D, Q7 N7 R% d" E
gamma=45*pi/180;
. [: F% ^' I, A; S) V6 K- q
rho=1.205;%rho为气体密度kg/m3
: d2 I. G, W. h5 Y# G
V=300;%V为气体速度m/s
2 Q4 ^/ N4 P- @0 d- Z6 F
d=10;
) ~8 e: S3 o9 i/ F5 `
Q=V*pi*d^2/4;%Q为气体流量
, O4 ]' y0 ^. Y6 E R' l" a
2 [+ s- |) r) W; \
f1=x(1);
0 A3 M8 f% q' z; X
%f1为轮与管壁之间的滚动摩擦
9 `3 x6 d& x" g" s
f2=x(2);
$ k' A' Y" t2 p' E( M1 F! b
%f2为轮与轮轴之间的摩擦
7 y: d/ d2 v( N/ D% ]% f
Fp=x(3);
, x/ s0 v( m& B' I3 K
% Hx为由于周向力引起的滑动趋势的摩擦力
; L% E6 j. ^/ f& Y# ~7 z
Fx=x(4);
" f* V: b+ D S- v$ e( f7 G$ H
%Hy为由于G和W引起的沿管道轴线方向滑动趋势的摩擦力
& U2 ^5 ^5 A9 T' P4 @
Fy=x(5);
' [9 V8 K" ?2 `; b: k, D
%Fx为轴对轮的周向分力
; ~2 L, F# o y0 o& R! C
Fz=x(6);
4 j' q6 t% o0 K$ m
%Fy为轴对轮沿管道轴线的分力
5 I0 v- t& `1 J4 c6 M
Hx=x(7);
$ R Y" [2 k, }7 t' ~
%Fz为轴对轮沿管道径向的分力
7 L( E* X1 `3 {* C
Hy=x(8);
' q8 I2 n' V) g# L5 }
%Fp为牵引力分布在每个导向轮上的力
3 ]% K* ?" r" T7 {4 {- b
Nz=x(9);
* J0 x3 u. P; f% ~' Q# W
Tm=x(10);
8 ~% P2 W8 V; `3 \8 |$ W+ ?; e7 a
( J2 r# z. h6 W) L# |: Z+ c
. K/ U) ]$ S6 d( O% ~* l& T( o/ h
%Nx=x(9);
$ z" q! T6 F0 ]2 x: Q7 Q$ A
%Nx为轴对轮的周向分力
# M% L. d) Y# Q6 s7 y; {! U. N% V
%Ny=x(10);
, U" _3 e7 J: l2 C" S% K
%Ny为W和G在管道轴向的分力与滚动摩擦的和
% p6 \. U+ Y1 G+ q& L
%f=x(11);
2 G( ^' t* r3 C+ @1 J3 m& `
%f为由电机扭转和气体冲量引起的滚动趋势的摩擦力
: Q; F! J" n( G x
b7 P. D+ L' e5 e$ w
fun(1)=Hy+(f1+f2)*sin(theta)-Fy;
3 b$ p2 n( j! j' R8 b
fun(2)=(f1+f2)*cos(theta)+Fx-Hx;
* f/ s, f3 x6 _+ P0 n$ K' [! j
fun(3)=Hy*sin(theta)*R+f1*R-Hx*R*cos(theta)+f2*r;
: X8 A/ c0 i. a% T# t
fun(4)=Hx^2+Hy^2-(mu*Fz)^2;
5 y8 l+ s" N+ T
fun(5)=k*Tm-3*Hx*D/2+3*f1*cos(theta)*D/2;
1 N' Y/ \+ r8 ^ ~6 S
fun(6)=k*Tm-3*Fx*D/2;
" P6 u2 L- q O: m- U4 a. J
fun(7)=Fp-(G2*cos(gamma)/6+Nz*fr/R+Wf/6);
9 d2 W8 B, U8 k# g" p
fun(8)=(Tm+3*rho*Q*V*D/2)/(D*6/2)-mu*Nz;
# D! d1 W% P6 N. E- d
fun(9)=2*Fp/cos(alpha)+G1*cos(beta)/3-Fy;
% m6 Q: Y+ o3 j. G
fun(10)=3/2*D*Fx/k-Tm;
( g/ e1 s$ f- m' d2 Y: ]
1 H3 Q/ j# ] v/ T; o
2 H$ `3 H) {& a0 o$ m# B
在求解时，matlab显示:
8 I' D" Z/ _8 R3 W3 |3 x* I. g
"Solver stopped prematurely.
1 v- k5 `9 O6 h# a4 g; T
fsolve stopped because it exceeded the function evaluation limit,
! E" X$ b9 j3 P
options.MaxFunEvals = 1000 (the default value)."
- m' _6 H6 c1 P
即超过了fsolve函数的迭代次数限制（1000），该怎么解决呢？
作者:
o晴o
时间:
2017-8-26 10:08
1 d5 z4 ]* I) B2 f% j" E
作者:
2249039642
时间:
2017-11-12 16:49
这个算的是啥问题，我觉得matlab解非线性真心不怎么行
+ Q% k+ z. G0 f1 L
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
1 d5 z4 ]* I) B2 f% j" E