\[ScriptCapitalN] = 1024; \[ScriptCapitalM] = 20;
uncorrelated =
Parallelize[
Table[RandomVariate[
NormalDistribution[0,
25*\[Pi]]], {\[ScriptCapitalN]}, {\[ScriptCapitalN]}]];
kernel = (1/(\[Pi]*(\[ScriptCapitalM]/2)^2))*
Table[circ[
Sqrt[(i - \[ScriptCapitalM]/2 +
1)^2 + (j - \[ScriptCapitalM]/2 + 1)^2]/(\[ScriptCapitalM]/
2)], {i, 1, \[ScriptCapitalM]}, {j, 1, \[ScriptCapitalM]}];
smooth = ListConvolve[kernel, uncorrelated, 1];
ListPlot[smooth[[All, \[ScriptCapitalN]/2]],
AxesLabel -> {"Pixel Number", "Smoothed Phase"}]

书上的图

因为运行一次时间太长了，这是把N改小后的结果。表上都没有数据

你circ的定义呢？