MATLAB 中这个for循环这样有什么问题吗

为什么我这个运行出来G2_si只有61个值而且全是0


```c#
T=8; %ns t_rt=2*pi/w_FSR
r_s=1;%MHz
t0=3.7;%ns
dT=0.165; %ns探测器分辨率

for t=-100:1:100 %ns
    if t >= t0/2
        G2_si(:,1)=exp(-1*r_s*T-4*(1*T-t).^2/(dT)^2);
        for j=1:1:60
            G2_si(:,j+1)=G2_si(:,j)+exp(-j*r_s*T-4*(j*T-t).^2/(dT)^2);
        end
    else
        G2_si(:,1)=exp(-1*r_s*T-4*(1*T-t).^2/(dT)^2); 
        for j=1:1:60
            G2_si(:,j+1)= G2_si(:,j)+exp(-j*r_s*T-4*(-j*T-t).^2/(dT)^2);
        end
    end
end

```

因为你的exp函数里面的值为-1*r_s*T-4*(1*T-t).^2/(dT)^2,我看第一步这个值就是-5000多,相当于是exp(-5000),这个值在matlab双精度表示就是0,因为matlab最多表示到10的-308次方.exp(-5000)早就是0了。建议仔细检查是不是数据的问题,比如纳秒应该跟GHz才能相乘才对,把量级搞对!!
然后,如果必须要这么小的计算结果,那么你可以用sym计算,但是非常费时间