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

为什么Matlab计算结果和单片机相差非常大?

[复制链接]
  • TA的每日心情
    开心
    2022-1-29 15:03
  • 签到天数: 2 天

    [LV.1]初来乍到

    跳转到指定楼层
    1#
    发表于 2022-8-24 13:56 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

    EDA365欢迎您登录!

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

    x

    遇到一个很奇怪的问题,下边这个公式使用Matlab和EXCEL计算结果相同且正确

    使用单片机和csharp测试又是另一个结果

    -(2*L^4*b - L^4*c + 8*L^2*b^3 + 2*L^2*c^3 + 8*c^3*h^2 - 4*L^2*b*c^2 - 4*L^2*b^2*c - 8*L^2*b*h^2 - 4*L^2*c*h^2 + 2*L*h*((L + c)*(L - c)*(L^2 - 4*b^2 + 4*h^2)*(L^2 - 4*b^2 + 8*b*c - 4*c^2 + 4*h^2))^(1/2))/(4*(4*L^2*b^2 - 4*L^2*b*c + L^2*c^2 - 4*L^2*h^2 + 4*c^2*h^2))

    单片机和csharp计算a=-9.58444976;

    double L =    270;
    ! y. g0 s3 h1 i1 Ddouble b =    14.33;9 X1 v! F6 r' w( R: [/ h
    double c =    20.25;' C1 r$ L! j, }$ S" Z
    double h =    120.21;

    a =  -(2*pow(L,4.00)*b - pow(L,4.00)*c + 8.00*pow(L,2)*pow(b,3.00) + 2*pow(L,2.00)*pow(c,3.00) + 8*pow(c,3)*pow(h,2) - 4*pow(L,2)*b*pow(c,2) - 4*pow(L,2)*pow(b,2)*c - 8*pow(L,2)*b*pow(h,2) - 4*pow(L,2)*c*pow(h,2) + 2*L*h*pow(((L + c)*(L - c)*(pow(L,2) - 4*pow(b,2) + 4*pow(h,2))*(pow(L,2) - 4*pow(b,2) + 8*b*c - 4*pow(c,2) + 4*pow(h,2))),(1/2)))/(4*(4*pow(L,2)*pow(b,2) - 4*pow(L,2)*b*c + pow(L,2)*pow(c,2) - 4*pow(L,2)*pow(h,2) + 4*pow(c,2)*pow(h,2)))1 u8 g/ z5 n) s3 g
    ;" k: ?: c  F; v- U5 E; P7 ?  i* V% E

    EXCEL和Matlab计算a=126.3747;

    270    L
    ' k" ]/ }% `: e0 G- ]( _2 S5 [14.33    b
    0 v1 x9 J3 R! e! F0 k# v20.25    c6 }! K, W# {* ?; [
    120.21    h

    = -(2*A1^4*A2 - A1^4*A3 + 8*A1^2*A2^3 + 2*A1^2*A3^3 + 8*A3^3*A4^2 - 4*A1^2*A2*A3^2 - 4*A1^2*A2^2*A3 - 8*A1^2*A2*A4^2 - 4*A1^2*A3*A4^2 + 2*A1*A4*((A1 + A3)*(A1 - A3)*(A1^2 - 4*A2^2 + 4*A4^2)*(A1^2 - 4*A2^2 + 8*A2*A3 - 4*A3^2 + 4*A4^2))^(1/2))/(4*(4*A1^2*A2^2 - 4*A1^2*A2*A3 + A1^2*A3^2 - 4*A1^2*A4^2 + 4*A3^2*A4^2))


    * z7 z: X& `: E6 s1 |8 Z& G

    不知道这个问题是怎么产生的,正在找原因。

    , B% U2 c; C. m
  • TA的每日心情
    开心
    2022-1-29 15:03
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
     楼主| 发表于 2022-8-24 14:33 | 只看该作者
    问题找到, 里边有一个(1/2)改成(1.00/2.00)后解决
    ) p! E! T3 l9 K& s, {解决方法是把公式分解,一段一段对比

    该用户从未签到

    3#
    发表于 2022-8-24 14:40 | 只看该作者
    会不会溢出了
  • TA的每日心情
    开心
    2022-1-29 15:03
  • 签到天数: 2 天

    [LV.1]初来乍到

    4#
     楼主| 发表于 2022-8-24 14:59 | 只看该作者
    问题找到,是2个整数相除结果的小数部分被舍掉了
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-8-24 01:42 , Processed in 0.125000 second(s), 23 queries , Gzip On.

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

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

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