QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 10340|回复: 2
打印 上一主题 下一主题

matlabd的fsolve函数求解方程组

[复制链接]
字体大小: 正常 放大

2

主题

6

听众

296

积分

升级  98%

  • TA的每日心情
    奋斗
    2017-4-29 14:31
  • 签到天数: 34 天

    [LV.5]常住居民I

    2013挑战赛参赛者

    群组2013认证赛C题讨论群组

    群组第四届cumcm国赛实训

    跳转到指定楼层
    1#
    发表于 2013-1-26 19:23 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    2体力
    function fun=Untitled6(x)
    0 E. Z% }2 S# p8 R0 n/ K theta=pi/18;* s8 S) O$ Q4 R9 B( E2 L
    k=0.97; %k为万向节的传输效率
    ' j8 m$ [9 i& T4 W9 U5 e+ G* p1 i" ]7 F R=13; %R为轮子的半径  h% f8 Z- f6 \8 o
    r=4;
    5 H4 k9 I# Q) ? D=250;
    5 N6 W- {. E7 h, S" A5 t* r9 g8 c) B" i fr=0.2; %fr为轮子的滚动摩阻
    1 ?' |* L( r/ @& o mu=0.4; %mu为轮与有机玻璃的摩擦系数
    * W8 j, P' q5 w! I" S8 i$ q3 t Wf=5; %Wf为电线所产生的负载
    0 T# W$ Z  u  O7 A3 m8 R G1=5;
    1 C# O( Z: d- C G2=80; %G2为电机,导向装置和加工头的重力和N
    ' g& L5 P& e0 B; ` alpha=11.5*pi/180;& V3 J' q  V& y
    beta=33.5*pi/180;4 d' I& `0 E# x
    gamma=45*pi/180;
    4 E6 i, L* Q" o3 W2 \ rho=1.205;%rho为气体密度kg/m3- L& \6 L- [! m/ A0 Z* I, _9 o
    V=300;%V为气体速度m/s
    / t, E2 ?( O8 }5 O9 a! { d=10;6 j( X% O. E6 N9 z
    Q=V*pi*d^2/4;%Q为气体流量
    3 C2 ?* q# R  Q- k3 X
    : R/ r* h' |$ d9 I5 u3 [( ]+ Lf1=x(1);8 ~3 j. ~7 M) J
    %f1为轮与管壁之间的滚动摩擦8 x: ~! g# ~  ]/ y+ I5 |+ V
    f2=x(2);( t, `8 U' V; n0 t8 d
    %f2为轮与轮轴之间的摩擦7 K5 k* `! S3 s' W3 Y) Q. K
    Fp=x(3);+ ]/ A/ O' d: g! d8 a6 s9 T0 Z" L
    % Hx为由于周向力引起的滑动趋势的摩擦力
    - H9 D, n3 H$ T7 ? Fx=x(4);$ M# F9 ^' R% ?( O- y1 P1 j2 [
    %Hy为由于G和W引起的沿管道轴线方向滑动趋势的摩擦力. ]1 x7 J7 D7 ~/ m/ i1 ]7 F3 ^2 z
    Fy=x(5);
    0 e; Z" t1 I2 i8 F( w. b0 K( \! k& } %Fx为轴对轮的周向分力
    1 ^( ?  {6 f" q) o- E% h Fz=x(6);; t! }2 X1 l! J1 Z# X* P& a0 w4 Z
    %Fy为轴对轮沿管道轴线的分力; Y& ]/ x0 K6 ?2 n5 w8 N
    Hx=x(7);- Q+ ^* R! T/ ]: W, F, \- r. ^
    %Fz为轴对轮沿管道径向的分力
    ) `7 C( A" ?  g$ e- o( J& K* J% X" D1 j Hy=x(8);
    , P1 x+ K0 g( v  C %Fp为牵引力分布在每个导向轮上的力6 K5 }( w' _! B  L0 N: [
    Nz=x(9);
    * f  o1 G6 S9 @; ?4 _) D Tm=x(10);* x/ ^5 u. }7 d. n1 Y
    $ m0 D5 ?( n: k+ N: ^! g
    . B" L' O7 y0 B, x) |
    %Nx=x(9);: p2 z& ^$ N$ T7 i3 s8 o
    %Nx为轴对轮的周向分力
    * J' Q$ F' |) _$ R' t9 P: b %Ny=x(10);
    9 ]' I: c, @! ]4 G( c %Ny为W和G在管道轴向的分力与滚动摩擦的和9 `* `; x: }3 a7 m) U% @1 c
    %f=x(11);9 n0 ~3 [$ k" Z! @
    %f为由电机扭转和气体冲量引起的滚动趋势的摩擦力
    $ ]2 U1 j" N4 U: ?; t 5 C, e* T  \8 A! n, a
    fun(1)=Hy+(f1+f2)*sin(theta)-Fy;; M! h6 H5 t% J! l5 k$ ]6 d. w% w
    fun(2)=(f1+f2)*cos(theta)+Fx-Hx;  P& X9 Y( d/ C8 P; {
    fun(3)=Hy*sin(theta)*R+f1*R-Hx*R*cos(theta)+f2*r;
    % g/ R& a5 a) {, t1 Y fun(4)=Hx^2+Hy^2-(mu*Fz)^2;2 o. r0 [1 n$ ^2 f+ J  F
    fun(5)=k*Tm-3*Hx*D/2+3*f1*cos(theta)*D/2;; i8 l* V& I1 u2 _
    fun(6)=k*Tm-3*Fx*D/2;
    5 A7 f4 r+ S) R4 q7 ^; D fun(7)=Fp-(G2*cos(gamma)/6+Nz*fr/R+Wf/6);
    & G7 H6 W* R+ C% @5 C4 [ fun(8)=(Tm+3*rho*Q*V*D/2)/(D*6/2)-mu*Nz;
    ; e. h& ^9 i( G) G- s. @ fun(9)=2*Fp/cos(alpha)+G1*cos(beta)/3-Fy;
    $ N8 p  c1 X. x- c' O fun(10)=3/2*D*Fx/k-Tm;" j! {7 ?2 `& u/ e1 x$ K) m# V  f

    , t1 `0 _+ B/ K! _
    " i3 ]8 H* h9 F. a/ \( x在求解时,matlab显示:7 K$ n1 e2 Q% F& x4 |
    "Solver stopped prematurely.
    6 r' B$ [+ ?8 {. U) o! z6 l: @& Z fsolve stopped because it exceeded the function evaluation limit,
    " z4 r; G3 v; T+ p& W5 {; o0 y; V" { options.MaxFunEvals = 1000 (the default value)."
    " X: N* u# e" Q7 b' r9 H( ~ 即超过了fsolve函数的迭代次数限制(1000),该怎么解决呢?

    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    o晴o        

    1

    主题

    5

    听众

    32

    积分

    升级  28.42%

  • TA的每日心情
    慵懒
    2017-8-26 10:05
  • 签到天数: 7 天

    [LV.3]偶尔看看II

    自我介绍
    黑黑的小白一枚

    邮箱绑定达人

    回复

    使用道具 举报

    5

    主题

    3

    听众

    105

    积分

    升级  2.5%

  • TA的每日心情
    擦汗
    2018-2-9 16:38
  • 签到天数: 22 天

    [LV.4]偶尔看看III

    国际赛参赛者

    群组Matlab讨论组

    群组F题

    群组E题

    群组D题

    群组C题

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-7-19 11:07 , Processed in 0.694654 second(s), 66 queries .

    回顶部