数学建模社区-数学中国
标题:
matlabd的fsolve函数求解方程组
[打印本页]
作者:
追梦人1101
时间:
2013-1-26 19:23
标题:
matlabd的fsolve函数求解方程组
function fun=Untitled6(x)
* J9 ~7 |; y$ A' Y, L/ U
theta=pi/18;
/ h" _' L; \1 e8 H4 C
k=0.97; %k为万向节的传输效率
: ^3 s. u0 ]6 J' K
R=13; %R为轮子的半径
; V5 q6 r5 o: a& A" I
r=4;
& }+ y5 }3 O# a- e: [: a% b* c
D=250;
" \+ n9 P4 `+ {
fr=0.2; %fr为轮子的滚动摩阻
& V% n7 A$ X1 k; S2 y6 h
mu=0.4; %mu为轮与有机玻璃的摩擦系数
7 j# f& k+ _% k
Wf=5; %Wf为电线所产生的负载
% f- H4 U) o' G" b
G1=5;
7 G% g1 M) Z1 ]
G2=80; %G2为电机，导向装置和加工头的重力和N
* H' I5 z# L! {% b( K0 J) ]# F
alpha=11.5*pi/180;
4 ^5 m3 z! x1 K1 }
beta=33.5*pi/180;
+ ~7 D I; E$ ]. y( r
gamma=45*pi/180;
& d) ?) N+ ?) W; l) j
rho=1.205;%rho为气体密度kg/m3
0 o2 b2 y$ [# i7 z9 x6 ?
V=300;%V为气体速度m/s
8 q7 K$ D" T! n% a, K. I+ m* [( O
d=10;
" f+ X- {! }- m7 A
Q=V*pi*d^2/4;%Q为气体流量
! ?4 D% e% P- |* ^* ]6 ]4 u
1 y' ]- C+ n1 E9 m$ u/ y" {
f1=x(1);
8 W. y0 I. F, p, @+ N$ ]
%f1为轮与管壁之间的滚动摩擦
% f% {+ t7 }8 b( A4 J& p) Q6 Z
f2=x(2);
( l" n; e% w+ O0 T9 L
%f2为轮与轮轴之间的摩擦
6 ]" m5 A3 I* c0 J# s
Fp=x(3);
{% r2 n; q, P! X2 z$ C7 k! ?
% Hx为由于周向力引起的滑动趋势的摩擦力
5 |$ p, L- k* ]; P- W3 E
Fx=x(4);
* q: M8 r* B) x* M5 }
%Hy为由于G和W引起的沿管道轴线方向滑动趋势的摩擦力
' e. }0 N& T+ z. z. A3 M# M6 x
Fy=x(5);
5 B6 B6 | S7 ]
%Fx为轴对轮的周向分力
# S" T( U2 Z7 \7 W1 X
Fz=x(6);
% {; h/ i5 B3 y1 B- v, S. g: e: P
%Fy为轴对轮沿管道轴线的分力
; M5 ~) I! q, v
Hx=x(7);
2 `7 y2 U9 O( ]' ?$ J
%Fz为轴对轮沿管道径向的分力
& G9 s$ E4 s- E6 v8 L* d
Hy=x(8);
) ~# a( X+ Y! T# v4 F
%Fp为牵引力分布在每个导向轮上的力
9 }7 y2 j+ y: X2 {$ c8 l
Nz=x(9);
D) @4 k$ m' P" p3 Q
Tm=x(10);
' h2 J' N) y3 q2 ^" |; v: ?
! ?/ Y( }" L; Y1 M7 I6 a
: W) h$ l1 `' X8 x; ?+ E5 L
%Nx=x(9);
" O2 M3 K8 I, T: Q" Z
%Nx为轴对轮的周向分力
( i! w/ d( l, G8 s0 U7 D [' q
%Ny=x(10);
+ |+ G8 |6 n( N- y# l7 s9 Y/ R
%Ny为W和G在管道轴向的分力与滚动摩擦的和
8 O( k' C0 x! s# k1 k% K4 m
%f=x(11);
( E8 k7 u2 ^8 r9 L) c; i1 c
%f为由电机扭转和气体冲量引起的滚动趋势的摩擦力
3 ^$ E X) B9 H& W; D7 j1 q
2 C3 D: g1 h5 f1 R
fun(1)=Hy+(f1+f2)*sin(theta)-Fy;
* x3 |' x/ E: G
fun(2)=(f1+f2)*cos(theta)+Fx-Hx;
: w! v/ c+ P' }; p9 Q; d2 F+ t
fun(3)=Hy*sin(theta)*R+f1*R-Hx*R*cos(theta)+f2*r;
- y- s( t# d3 K* n: D
fun(4)=Hx^2+Hy^2-(mu*Fz)^2;
, ?5 S. d5 q: Z/ Z2 l, c
fun(5)=k*Tm-3*Hx*D/2+3*f1*cos(theta)*D/2;
3 r k' l( g& p- R7 t& C }
fun(6)=k*Tm-3*Fx*D/2;
# c% ~. u! J+ \, u# {& I0 c
fun(7)=Fp-(G2*cos(gamma)/6+Nz*fr/R+Wf/6);
7 n* n7 d4 b) c: b# w$ n
fun(8)=(Tm+3*rho*Q*V*D/2)/(D*6/2)-mu*Nz;
1 Y! F( t5 P+ s
fun(9)=2*Fp/cos(alpha)+G1*cos(beta)/3-Fy;
0 l8 N3 F/ l" i
fun(10)=3/2*D*Fx/k-Tm;
% _5 I* _8 P$ }7 K" C! h
8 a# }: U0 \7 \. d, D
( r* R a/ O; M! I4 L* Y4 V
在求解时，matlab显示:
' ^$ C& ?. c2 j
"Solver stopped prematurely.
# [# d4 O& D6 I8 |
fsolve stopped because it exceeded the function evaluation limit,
. M- b+ v. K1 d0 z
options.MaxFunEvals = 1000 (the default value)."
* j8 e% `7 S( ^( k5 n5 f8 ]% {6 z
即超过了fsolve函数的迭代次数限制（1000），该怎么解决呢？
作者:
o晴o
时间:
2017-8-26 10:08
& _( F, t/ [5 ~7 N$ t0 M
作者:
2249039642
时间:
2017-11-12 16:49
这个算的是啥问题，我觉得matlab解非线性真心不怎么行
# v2 p# C* v9 \$ y {1 e6 ]
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
