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

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

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

    [LV.2]偶尔看看I

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

    EDA365欢迎您登录!

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

    x

    ; P' V- O1 i, A! M计算一个方程里的参数,公式如图,为什么一直出错,运行不出来QAQ,麻烦大神指点一下,非常感谢!  i6 ^5 e3 Q/ d0 b% h0 f2 A5 z: ~

    . E; L/ s1 A# j& i0 _6 n' u: }f1=@(x,a)(x.^(-a)-1)./a;
    : j) S3 g$ Q1 Zf2=@(x,a)diff(f1,x);, b/ h5 V0 t- ]% X6 W7 ]
    f3=@(x,a)f1(x,a)./f2(x,a);2 m6 A' u. F3 M: C# O
    f4=@(x,a)integral(@(x)f3(x,a),x,0,1);
    ( e4 i2 q5 L, n/ x- W' _1 if5=@(a)1+4*f4(x,a)-6/2475;6 v( D: N. d6 M% t3 W% V/ X
    a=fsolve(f5);
    ) G4 a7 r3 \. m0 D' }$ o6 I

    该用户从未签到

    2#
    发表于 2020-6-16 13:59 | 只看该作者
    clear; clc; close all;
    5 A5 G4 B* `; [syms t alpha real2 I) N3 b2 l. V" N- n( O
    syms theta(t)
    + M' [& b: _/ Q6 w! q; |6 t, dtheta = (t^(-alpha) - 1)/alpha;3 t3 F) K% \: v* Y+ K9 |- r- o& F
    Dtheta = diff( theta, t );
    ( t1 X( e+ i! x* l8 atau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;
    4 l' f% C5 e' M( EAlpha = solve([tau == 6/2457 ],[ alpha ])& e& {+ g$ n8 P( x, U/ r0 ~7 U
    RHS_LHS = 4*int( ( (t^(-Alpha) - 1)/Alpha )/diff( ( (t^(-Alpha) - 1)/Alpha ), t ), t, 0, 1 ) + 1 - 6/2457
    7 S4 O6 {$ L. x, S% w1 l3 D+ V
    Alpha =
    6 j0 t0 m3 {8 K- q4/817
    8 [0 p/ J2 R: O! a' H6 K* o
    + R6 F% J/ p& TRHS_LHS =: Q2 d& x. k# N. u: ]" A# `, j
    0
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-5 22:53 , Processed in 0.156250 second(s), 26 queries , Gzip On.

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

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

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