matlab | 生成低通滤波器

引言

生成的低通滤波器如下：(在菲涅尔衍射积分中，我们可认为它是一个衍射面上的衍射孔)

算法

以下算法可直接生成如上图所示的低通滤波器。

r=512,c=r;                             
a=zeros(r,c);                          
a(r/2-r/4:r/2+r/4,c/2-c/4:c/2+c/4)=1; 

分析

上述算法难点在于对矩阵的索引，下面给出详细的分析：
r/2-r/4:r/2+r/4等效为r/4:3r4，也就是提取矩阵的第r/4到3r4中的元素。
c/2-c/4:c/2+c/4表示对列进行同样的操作。
提取出来一个小矩阵，让小矩阵所有元素都为1。
便得到了低通滤波器，