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

MATLAB数据非常小,输出的值为零怎么办

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
clear all, E+ }% ]& O' L/ ~
I=10^45;# y  ?9 Z1 S1 S) X
M=1.4;
- A. l& h, D2 _. u! I  jP0=0.01;%s
5 e2 V& H( v0 k6 p) Lu=1.7;
# @/ n3 E% X' B, eOm=2*pi/P0;  w; @5 d* t5 u8 b2 o1 S
dt=630720000;%20year;
2 R1 t' o2 O) l% V9 Sc=3*10^10;%cm2 z3 x. A3 j& v" R2 s( c
for m=1:10000;
2 l( M2 z( ?: P, T$ y    dOm=-(2*(u^2)*(Om^3)*dt)/(3*(c^3)*I);5 i; X- B3 e( X' o+ _: j1 b, {. E
    Om=Om+dOm;
$ T! Z1 [- N# w; A    dotOm=dOm/dt;8 Z! V9 U" I( i) P
    p(m)=2*pi/Om;
* n9 ^6 X; G" c. [+ O" g# L    D(m)=(dotOm-dd)/dt;# w! H) Y7 Z) k* j
    dd=dotOm;
2 E0 V" R$ r1 a9 j$ _: t, ^$ ~* c" {    d(m)=dotOm;
. ~* Z/ ]+ a; |& P    dotp=-dotOm*2*pi/Om/Om;. }# E) {1 z0 C. ~9 M
    dp(m)=dotp;! G+ H2 X( N! s* @. K) q- |$ b
end
5 J5 B/ Q  S- e, s3 `  hhold on
( T4 R. y' k" R0 l, B2 g+ q3 Z/ iplot(p,dp)$ ?* j* K/ ~# ~/ e. H9 T! G* H& ^0 s
7 ]" L6 M( @: k4 D; F  x

- f- o! f6 t- q0 k
7 T8 r# y% t0 r5 g; p( Z我的dOm得到的值数量级在10的负49次方左右,所以循环过程中Om都是加0根本不会变化,如何提高精确度能够让Om在Om=Om+dOm中变化* w. ~( j' [* T- }2 b1 U
1 p  m1 _3 h. |0 X2 j
  • TA的每日心情
    开心
    2023-1-3 15:10
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2022-7-19 15:21 | 只看该作者
    把你的单位改一下
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-6-19 10:59 , Processed in 0.078125 second(s), 23 queries , Gzip On.

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

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

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