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

xiannvjiejie

clear all
/ g  c/ ~7 G' l% D- J; |& ~' rI=10^45;
4 v/ Y( ?& W6 _% m2 Z4 H6 t. J3 RM=1.4;
$ E: ^, C! L: M3 J$ v& ~" |P0=0.01;%s
# p7 Z+ Y/ i! i7 s1 Y. ]0 Bu=1.7;
9 @, r& i8 ?8 g  ?4 D- `Om=2*pi/P0;
dt=630720000;%20year;
# t" n( m, D4 Y( r, f* |) Nc=3*10^10;%cm
- y% _- T0 o. o) s% @' ~for m=1:10000;
7 p0 D" u; V- a" E6 T" ?0 U5 n3 f    dOm=-(2*(u^2)*(Om^3)*dt)/(3*(c^3)*I);
    Om=Om+dOm;
; `3 h5 a- |% k9 g1 M$ i  m2 i) i    dotOm=dOm/dt;
    p(m)=2*pi/Om;
' T% X! A3 x  z7 l$ F    D(m)=(dotOm-dd)/dt;
    dd=dotOm;
3 Y- {% `0 f8 C( I9 @: |# |& g1 Y    d(m)=dotOm;
    dotp=-dotOm*2*pi/Om/Om;
    dp(m)=dotp;
/ E! W1 j$ t/ N5 `end
; J; }, N5 g! C; k& [% ?0 s) s" Thold on
. m( V! [% f/ O7 G" C, v2 l+ gplot(p,dp)


- d; h% L! v+ t% ~7 Y& p
我的dOm得到的值数量级在10的负49次方左右，所以循环过程中Om都是加0根本不会变化，如何提高精确度能够让Om在Om=Om+dOm中变化
) f8 T8 L2 T. z, C. S4 K/ ~1 A8 O  a

modengxian111

把你的单位改一下