我正在尝试为 h（[MxM 矩阵]）求解以下方程：

我有 k，它是一个 2D [MxM] 对称矩阵，我有 x，它也是一个 2D [MxM] 对称 toeplitz（自相关）矩阵。

我知道这基本上是一个二维反卷积问题，但这不是我的领域，我不知道如何在 MATLAB 中做到这一点。

另外，如果可能的话，我更喜欢时域解决方案，但频域也可以！

尝试的解决方案#1：频域划分：

h_pred=ifft(fft2(k)./fft2(x));

我认为这应该可行，但由于某种原因，它一直给我一个想象的答案！

尝试的解决方案 #2：Toeplitz 矩阵求逆：

h_pred=inv(x)*k;

这效果更好，但是虽然恢复的 h_pred 更接近真实的 h，但它仍然很远......

请帮忙谢谢！