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

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

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

    [LV.2]偶尔看看I

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

    EDA365欢迎您登录!

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

    x

    5 Q3 z: p3 U6 h1 M% U# J计算一个方程里的参数,公式如图,为什么一直出错,运行不出来QAQ,麻烦大神指点一下,非常感谢!( k2 B, _+ ]4 @6 S: t

    ) I' Q- a* a! [% F& ^8 ^f1=@(x,a)(x.^(-a)-1)./a;% [8 `6 E& q0 q# N5 L+ W
    f2=@(x,a)diff(f1,x);6 I& l' }8 z/ l, p1 l
    f3=@(x,a)f1(x,a)./f2(x,a);
    , P" L; R" P. if4=@(x,a)integral(@(x)f3(x,a),x,0,1);. ~1 R6 K3 I0 t9 O2 O) r5 I
    f5=@(a)1+4*f4(x,a)-6/2475;0 P* A* u& R8 P. ], o' ?& I
    a=fsolve(f5);( L7 c/ B# z: H5 J3 @

    该用户从未签到

    2#
    发表于 2020-6-16 13:59 | 只看该作者
    clear; clc; close all;3 @( a$ ~+ ?' C& I0 p& L
    syms t alpha real
    + c9 `# ?  {% ^, G$ N# gsyms theta(t)( `* X0 S1 R  e2 |! U5 e
    theta = (t^(-alpha) - 1)/alpha;
    & _6 }  z  M; ~8 Z) TDtheta = diff( theta, t );
    4 J' r' t" B# L( ~7 o( A2 |tau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;8 }# b0 W2 t7 j9 m4 I1 s
    Alpha = solve([tau == 6/2457 ],[ alpha ])
    % X8 R! ^3 ]! o  m% `; F; |$ E; oRHS_LHS = 4*int( ( (t^(-Alpha) - 1)/Alpha )/diff( ( (t^(-Alpha) - 1)/Alpha ), t ), t, 0, 1 ) + 1 - 6/2457
    9 i! a, J9 |' g8 n# X; Z
    + Q+ u; M3 S, T2 ]) ^! B/ IAlpha =) U. O6 r# B, H$ D5 u
    4/8178 x) m2 W; l" Z/ n7 k( b

    6 q' n3 L' F* E- V# J9 J5 lRHS_LHS =0 @; W4 ~1 S2 q, c' z% i; H' R3 c
    0
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

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

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

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

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