数学建模社区-数学中国
标题:
matlabd的fsolve函数求解方程组
[打印本页]
作者:
追梦人1101
时间:
2013-1-26 19:23
标题:
matlabd的fsolve函数求解方程组
function fun=Untitled6(x)
8 ~- b/ q$ m8 x. p. L
theta=pi/18;
" h3 t" O- B5 Z, A
k=0.97; %k为万向节的传输效率
& |) U% b: j. @* n l; s& z
R=13; %R为轮子的半径
& `1 z% a$ p8 Y6 h
r=4;
3 Q6 E/ C. P! s8 W) h( d; A# A
D=250;
5 a& o) d& M' K; r
fr=0.2; %fr为轮子的滚动摩阻
0 c- o/ k* S& A0 W
mu=0.4; %mu为轮与有机玻璃的摩擦系数
/ T# r3 q( }# D7 T* h
Wf=5; %Wf为电线所产生的负载
% e, P) {# }/ K
G1=5;
: ~! C" C7 d, A/ A0 Z+ U
G2=80; %G2为电机，导向装置和加工头的重力和N
, G$ B, n8 ?: U% q: v
alpha=11.5*pi/180;
) W% J* n& j5 n1 c( f, q3 u
beta=33.5*pi/180;
, k7 ]' K4 @7 r/ q1 X3 w
gamma=45*pi/180;
4 D g# W' ~( f$ D5 s# t
rho=1.205;%rho为气体密度kg/m3
! M9 _1 ~( F6 d I) K: p& G% s( h* d
V=300;%V为气体速度m/s
# Q% |) N9 F' n6 \+ z# C% W
d=10;
0 v$ v% A# B" F3 V$ L, c
Q=V*pi*d^2/4;%Q为气体流量
4 X0 w" F. {5 q7 U- L" N/ J
; a& C9 H( E1 U2 G/ I" g
f1=x(1);
3 b4 e& H8 u; H# S' T) ^/ R1 `
%f1为轮与管壁之间的滚动摩擦
* _2 N3 z9 _' r3 f4 z, L- G& P
f2=x(2);
7 ?- l+ V4 {# v- k
%f2为轮与轮轴之间的摩擦
$ H* X$ h3 y. j1 S2 W+ T
Fp=x(3);
5 b2 F* f( \' A# i9 F4 R% D/ h9 Q: W2 l
% Hx为由于周向力引起的滑动趋势的摩擦力
8 A. Z4 @. t r$ ^# E& [# G
Fx=x(4);
+ \" ~' ~" h1 R d
%Hy为由于G和W引起的沿管道轴线方向滑动趋势的摩擦力
$ x' t" a) ?9 l8 e
Fy=x(5);
8 i _; C+ o9 L: |' k# [$ E1 H
%Fx为轴对轮的周向分力
/ S( D; C2 p1 @# j
Fz=x(6);
* d) t2 P% v& a+ u) {
%Fy为轴对轮沿管道轴线的分力
" M- N% B( Z: R7 n) l6 o1 c
Hx=x(7);
; G: }9 X9 s8 h/ x3 z! r
%Fz为轴对轮沿管道径向的分力
S- A, P. }: T1 c6 S
Hy=x(8);
/ U+ B) b0 b) f
%Fp为牵引力分布在每个导向轮上的力
4 M, @" h; K0 M3 L5 G7 F4 d. w+ J% g
Nz=x(9);
! o) s+ e; G% M9 l
Tm=x(10);
1 }& g/ |. }7 q! W9 ]7 o
8 j! }1 ?) r4 S0 m1 }$ L
9 w% Z Y& D" q% \% t$ U
%Nx=x(9);
( ^* L$ O- a: G7 ^
%Nx为轴对轮的周向分力
2 X$ x! M! \2 x/ Z" e1 y$ T
%Ny=x(10);
. p* D: \: H- p R) t& E8 n; \
%Ny为W和G在管道轴向的分力与滚动摩擦的和
. t6 |9 t) G& P" Z( O& |
%f=x(11);
& n. e7 q. X1 w" j; h
%f为由电机扭转和气体冲量引起的滚动趋势的摩擦力
; f+ E7 P: `6 s- Q
V& |. V9 X& Z0 [
fun(1)=Hy+(f1+f2)*sin(theta)-Fy;
2 A: s9 X$ t6 |8 m% C, P
fun(2)=(f1+f2)*cos(theta)+Fx-Hx;
3 ?6 p; a7 l# z
fun(3)=Hy*sin(theta)*R+f1*R-Hx*R*cos(theta)+f2*r;
- V- H" `) d. _# ~7 I. U
fun(4)=Hx^2+Hy^2-(mu*Fz)^2;
8 m5 N! i7 ?$ S
fun(5)=k*Tm-3*Hx*D/2+3*f1*cos(theta)*D/2;
9 z. T* d- k8 y7 j) _8 m( _
fun(6)=k*Tm-3*Fx*D/2;
# |' a" |$ z; @, b8 U7 M$ C
fun(7)=Fp-(G2*cos(gamma)/6+Nz*fr/R+Wf/6);
# \' g# U8 F4 B5 Z3 z
fun(8)=(Tm+3*rho*Q*V*D/2)/(D*6/2)-mu*Nz;
8 K6 D; i! ]/ m. i
fun(9)=2*Fp/cos(alpha)+G1*cos(beta)/3-Fy;
, s0 {1 d9 G K) d: ?1 C, G
fun(10)=3/2*D*Fx/k-Tm;
: o. D- ^ {3 h9 l/ i* _
2 g9 z. O, g$ ~/ R
$ V' s/ e( k4 S7 H0 h; x3 K
在求解时，matlab显示:
* z0 @$ |/ J: Z/ V. E9 [
"Solver stopped prematurely.
: b* i' s) s1 y# p
fsolve stopped because it exceeded the function evaluation limit,
A, D6 h/ Y1 Q& y
options.MaxFunEvals = 1000 (the default value)."
, g4 O5 c7 x* }8 y! w8 T" S
即超过了fsolve函数的迭代次数限制（1000），该怎么解决呢？
作者:
o晴o
时间:
2017-8-26 10:08
# ]. ~1 I' z, B) L/ S, G
作者:
2249039642
时间:
2017-11-12 16:49
这个算的是啥问题，我觉得matlab解非线性真心不怎么行
+ S4 o Q& u$ T. ^$ \5 y% G
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
# ]. ~1 I' z, B) L/ S, G