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

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

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

    [LV.2]偶尔看看I

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

    EDA365欢迎您登录!

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

    x

    1 e- h5 w9 b3 `% F计算一个方程里的参数,公式如图,为什么一直出错,运行不出来QAQ,麻烦大神指点一下,非常感谢!
    9 |; V& {) S2 {3 ?0 b8 j) [$ h' h) o- ]% ]% A1 u" s6 a' B5 Q9 p
    f1=@(x,a)(x.^(-a)-1)./a;
    $ v; T/ C5 Q2 x' W/ i$ A1 x+ }f2=@(x,a)diff(f1,x);
    . V1 @- s) m% c$ [# I' W; ~f3=@(x,a)f1(x,a)./f2(x,a);' Q+ q9 b: x( d+ N
    f4=@(x,a)integral(@(x)f3(x,a),x,0,1);% m% W3 [- I. _1 F  [! v4 }
    f5=@(a)1+4*f4(x,a)-6/2475;
    / D# J2 e0 o8 b) ^7 u9 ma=fsolve(f5);; B' J1 L3 o3 @& {% }2 z# |# {% n

    该用户从未签到

    2#
    发表于 2020-6-16 13:59 | 只看该作者
    clear; clc; close all;) b% K, f$ Z9 W" R+ _+ P7 ?
    syms t alpha real
    9 e) a' p$ B6 Q4 v* q" v1 jsyms theta(t)
    & G7 l/ p6 I; Z8 g  d+ G4 F. Atheta = (t^(-alpha) - 1)/alpha;
    + m3 x! r1 ]! m( y; ZDtheta = diff( theta, t );
    % u1 }0 K( W, u* vtau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;
    / W) b8 w3 A+ l2 \; E7 PAlpha = solve([tau == 6/2457 ],[ alpha ])/ e* W- Z, }5 G! v+ w
    RHS_LHS = 4*int( ( (t^(-Alpha) - 1)/Alpha )/diff( ( (t^(-Alpha) - 1)/Alpha ), t ), t, 0, 1 ) + 1 - 6/2457
    - y+ O/ |, W1 Q6 u0 b
    ( y/ k6 \5 _1 pAlpha =5 {% o( C, _. R: e
    4/817) v: t' T0 A0 U6 p. i2 L

    $ l4 M$ o4 j: \: w/ \" aRHS_LHS =
    , O9 p5 j! P7 C% W0
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-24 11:10 , Processed in 0.171875 second(s), 26 queries , Gzip On.

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

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

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