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

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

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

    [LV.2]偶尔看看I

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

    EDA365欢迎您登录!

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

    x
    8 t% U) _) J  Z0 F( n2 q
    计算一个方程里的参数,公式如图,为什么一直出错,运行不出来QAQ,麻烦大神指点一下,非常感谢!
    / B1 U4 f4 G- E- r8 p6 l
    / [' o6 v( W9 W: i  X8 lf1=@(x,a)(x.^(-a)-1)./a;
    ) P3 X) W6 v( z5 F0 t  ~f2=@(x,a)diff(f1,x);
    " H9 `; u9 L& w1 n  bf3=@(x,a)f1(x,a)./f2(x,a);7 A' N- n/ V' ~3 [( y7 H
    f4=@(x,a)integral(@(x)f3(x,a),x,0,1);
    : N, P' K5 y& P  P$ e1 hf5=@(a)1+4*f4(x,a)-6/2475;
    & a) q) b! Q. [; d0 Ha=fsolve(f5);5 p. F  d+ }9 Y

    该用户从未签到

    2#
    发表于 2020-6-16 13:59 | 只看该作者
    clear; clc; close all;
    1 B# H8 g9 c9 Y% Gsyms t alpha real; E) w% P9 M) A+ V5 h3 [% A
    syms theta(t)  l% i4 k: |; l- E
    theta = (t^(-alpha) - 1)/alpha;' g- {" M; T) G6 X+ Z
    Dtheta = diff( theta, t );7 a( S) y* G2 |9 ]
    tau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;1 R- Y" @& G! n4 B2 J) U6 W. [
    Alpha = solve([tau == 6/2457 ],[ alpha ])
    2 O2 y3 ~  q2 g2 MRHS_LHS = 4*int( ( (t^(-Alpha) - 1)/Alpha )/diff( ( (t^(-Alpha) - 1)/Alpha ), t ), t, 0, 1 ) + 1 - 6/2457
    : i8 |6 {" P: u
    8 j1 `2 X5 G: ^0 r! Q3 c9 JAlpha =
    ( I4 [0 `4 J' v7 o4/817& y. y7 w9 Z0 K+ |  K) {
    ' f) n; Z! U4 x4 `  Q! X+ c+ B$ i
    RHS_LHS =
    1 ~4 g( J: X- t8 z5 A: x) B& J0
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

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

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

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

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