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

Question

为什么我这个运行出来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

```

Answer 1

因为你的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计算，但是非常费时间