找回密码
 注册
关于网站域名变更的通知
查看: 374|回复: 1
打印 上一主题 下一主题

麻烦大神指点一下,非常感谢!

[复制链接]
  • TA的每日心情
    开心
    2019-11-29 15:39
  • 签到天数: 4 天

    [LV.2]偶尔看看I

    跳转到指定楼层
    1#
    发表于 2020-6-16 11:33 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

    EDA365欢迎您登录!

    您需要 登录 才可以下载或查看,没有帐号?注册

    x
    0 J# ^' y9 U: }& q# }" b8 Y. V
    计算一个方程里的参数,公式如图,为什么一直出错,运行不出来QAQ,麻烦大神指点一下,非常感谢!
    ) }! L- f- ]9 \- G- L
    9 ~* a' F- H' N2 |) m; p- B1 [f1=@(x,a)(x.^(-a)-1)./a;; K/ a3 z% b! r; Z
    f2=@(x,a)diff(f1,x);& U; J0 @% X) @) x0 u
    f3=@(x,a)f1(x,a)./f2(x,a);
    6 A7 D: F! s: e( E9 Zf4=@(x,a)integral(@(x)f3(x,a),x,0,1);
    2 _* v9 J2 o. }( [f5=@(a)1+4*f4(x,a)-6/2475;
    " v% w4 H8 ^% @+ |$ _+ s/ @a=fsolve(f5);
      Q# B* z+ v+ J! r

    该用户从未签到

    2#
    发表于 2020-6-16 13:59 | 只看该作者
    clear; clc; close all;- E  q' S: F# e
    syms t alpha real
      ^7 }2 [; t5 ~$ E2 P- W; }3 usyms theta(t)% H5 D8 j1 r  N5 g1 ^) ^' ?/ a
    theta = (t^(-alpha) - 1)/alpha;5 n6 F0 }8 O4 T2 d, z! u4 a
    Dtheta = diff( theta, t );7 f: Z* I, }7 E( g
    tau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;/ t  d' v4 v% e8 v6 T
    Alpha = solve([tau == 6/2457 ],[ alpha ])
    + i# f7 m& p" Y& [% \3 N3 vRHS_LHS = 4*int( ( (t^(-Alpha) - 1)/Alpha )/diff( ( (t^(-Alpha) - 1)/Alpha ), t ), t, 0, 1 ) + 1 - 6/2457  y' I' }4 F8 ~8 G7 ]. s9 ^

    $ v2 n" U6 h6 Q, W4 I0 sAlpha =
    3 R% b: U/ K' A1 ]& P2 j4/8172 b; f( [+ J3 \- `/ R$ ]: |$ d

    1 P9 f* M2 A& a, xRHS_LHS =* b. e$ [( G! I" R+ M( S
    0
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-8-24 08:29 , Processed in 0.125000 second(s), 26 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表