代码如下：
SystemOpen@
Export["a.dat",
First@Cases[
Normal@Plot[-NIntegrate[-w Cos[
w t] Exp[-w^2/4]/2/(Log[2]/0.2^2) (BesselJ[0, 2 w] +
Re[HankelH1[0, 2 w]/Pi/w/(-HankelH1[1, w/2])]*
NIntegrate[(0.5 - 4 Cos[y]) BesselJ[1,
Sqrt[0.25 + 16 - 4 Cos[y]] w]/
Sqrt[0.25 + 16 - 4 Cos[y]], {y, 0, Pi}] +
Sum[Re[2 HankelH1[k, 2 w]/Pi/
w/(HankelH1[k, w/2] 2 k/w -
HankelH1[k + 1, w/2])] NIntegrate[(0.5 -
4 Cos[y]) BesselJ[1,
Sqrt[0.25 + 16 - 4 Cos[y]] w] Cos[k y]/
Sqrt[0.25 + 16 - 4 Cos[y]], {y, 0, Pi}], {k,
1, \[Infinity]}]), {w, 0, \[Infinity]}], {t, 0, 10},
PlotRange -> {-0.05, 0.12}] Line[t_] -> t, Infinity]]