我在尝试对使用rfft2
创建的数组进行部分逆变换(沿一维合成)时遇到问题(当然,后来进行了过滤)。
问题是,我使用整个程序的真实输入,然后使用rfft2
从实域到傅立叶域来利用这种对称性是有意义的(我做了很多变换该程序)。这根据实数的(nx,ny)
矩阵生成(nx,ny//2-1)
(如果ny为奇数,则最后一个数字为+1,但您知道了要点)。我将这个过程命名为从(x,y)
域到(kx,ky)
域。我可以使用irfft2
还原此过程,并获取我的真实数据。
但是,在算法的一点上,我必须在中间域(kx,y)
中操作数据。在这个域中,我不知道系数(仍然很复杂)是否具有某种对称关系,但是我强烈怀疑它们必须验证(kx,y)
和(-kx,y)
处的系数是互共轭的。 / p>
我想到手动在(kx,ky)域中手动创建(nx,ny)复数数组(明确使用对称性),并在最后一个(ifft
后面简单地使用y
轴。但是,这样做很讨厌阅读,恕我直言,很容易出错,更不用说效率了,我可以为这一步而辞职。
实际的问题是,我可以有效地使用numpy
来计算该混合域的扩展系数吗?
干杯