Motivated by the discrete dipole approximation (DDA) for the scattering of electromagnetic waves by a dielectric obstacle that can be considered as a simple discretization of a Lippmann-Schwinger style volume integral equation for time-harmonic Maxwell equations, we analyze an analogous discretization of convolution operators with strongly singular kernels. For a class of kernel functions that includes the finite Hilbert transformation in 1D and the principal part of the Maxwell volume integral operator used for DDA in dimensions 2 and 3, we show that the method, which does not fit into known frameworks of projection methods, can nevertheless be considered as a finite section method for an infinite block Toeplitz matrix. The symbol of this matrix is given by a Fourier series that does not converge absolutely. We use Ewald's method to obtain an exponentially fast convergent series representation of this symbol and show that it is a bounded function, thereby allowing to describe the spectrum and the numerical range of the matrix. It turns out that this numerical range includes the numerical range of the integral operator, but that it is in some cases strictly larger. In these cases the discretization method does not provide a spectrally correct approximation, and while it is stable for a large range of the spectral parameter $\lambda$, there are values of $\lambda$ for which the singular integral equation is well posed, but the discretization method is unstable.
翻译:由离散的 dipole 近似值(DDAD) 驱动, 将电磁波通过电磁波散布为电磁波的电磁屏障碍, 它可以被视为时间- 调和 Maxwell 方程式中Lipmann- Schwinger 风格体积整体方程式的简单离散化, 我们分析的是具有强烈单核内核的聚合操作器的类似离散操作器。 对于包含 1D 中的有限 Hilbert 转换和 DAD 2 和 3 的 Maxwell 体积整体操作器的主要部分的内核功能, 我们显示, 与已知的投影法不相容, 但是, 仍然可以被视为一个无限块块块块的离析方块 Toeplitz 矩阵中的有限部分方法。 这个矩阵的符号由不绝对趋同的 Fourier系列提供。 我们使用 Ewald 的方法来获得这个符号的指数性快速趋同序列表示, 显示它是一个捆绑的功能, 从而可以描述矩阵的频谱和数字范围。 它包括集操作的数值范围, 但是, 它在某情况下, 最不稳定的光化法度范围 。</s>