为什么我这个运行出来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计算,但是非常费时间