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

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

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

    [LV.2]偶尔看看I

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

    EDA365欢迎您登录!

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

    x
    ( [: w  t. }% X& H' c; G: O
    计算一个方程里的参数,公式如图,为什么一直出错,运行不出来QAQ,麻烦大神指点一下,非常感谢!
    4 Q7 M7 `* a7 C( y# r. G3 G# w9 W1 O
    f1=@(x,a)(x.^(-a)-1)./a;
    " F0 |0 e# r  Q: M  {f2=@(x,a)diff(f1,x);4 ^. w  H8 X# P6 f8 Z+ C
    f3=@(x,a)f1(x,a)./f2(x,a);$ i. y) K# o) p0 ^. f
    f4=@(x,a)integral(@(x)f3(x,a),x,0,1);
    % w, k7 @9 ?9 I6 M. L& \* x; Ef5=@(a)1+4*f4(x,a)-6/2475;  Z3 c# E+ J1 a" c
    a=fsolve(f5);% H9 T7 a7 M$ f6 s% k. I& ~

    该用户从未签到

    2#
    发表于 2020-6-16 13:59 | 只看该作者
    clear; clc; close all;
    9 I8 i* _$ H  K: t& }9 s. Y& }syms t alpha real! d, Z3 y& g6 P+ H
    syms theta(t)- s- p; P. s; V
    theta = (t^(-alpha) - 1)/alpha;
      F6 G& J" G7 v. vDtheta = diff( theta, t );
    " ^% @7 h5 s# b  {2 Ztau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;
    $ _' _4 N( y3 |Alpha = solve([tau == 6/2457 ],[ alpha ])  s! G7 s1 o1 a4 W' W
    RHS_LHS = 4*int( ( (t^(-Alpha) - 1)/Alpha )/diff( ( (t^(-Alpha) - 1)/Alpha ), t ), t, 0, 1 ) + 1 - 6/2457
    # N6 ?8 K5 _, h9 m6 L" s
    8 |" T& p/ j" l' o  J. E4 FAlpha =- m) {9 E7 w, z8 ^0 E
    4/817
    ' G1 S3 |! x& @/ c" q7 ^
    % i: w* j" s: k7 G. V! Y" w6 JRHS_LHS =9 d2 `4 G5 U6 X% w6 S; E5 s
    0
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-7-25 10:54 , Processed in 0.109375 second(s), 26 queries , Gzip On.

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

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

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