[发明专利]一种基于快速傅里叶变换的二维离散数据平滑方法

说明书

本发明具体涉及一种基于快速傅里叶变换的二维离散数据平滑方法，包括以下步骤：将待平滑的二维离散数据通过二维傅里叶变换由空间域变换到波数域，以得到原数据角谱；以补零的方式对原数据角谱进行扩展，得到扩展角谱；将扩展角谱通过二维傅里叶逆变换由波数域变换到空间域，以得到平滑后的二维离散数据。本发明通过快速傅里叶变换将二维离散数据由空间域变换到波数域，并在波数域中进行补零扩展，最后逆变换回空间域，实现了二维离散数据的快速平滑。

技术领域

本发明属于二维离散数据平滑领域，具体涉及一种基于快速傅里叶变换的二维离散数据平滑方法。

背景技术

在工程实践中，往往受测试条件的影响，只能获取少量的测试数据，如何根据少量测试数据获得更加平滑的结果具有重要的实际意义。

目前已有学者提出采用插值方法平滑离散采样数据，如拉格朗日拟合法，分段三次埃米尔特拟合法，最小二乘法，三次样条曲线法等。这些曲线平滑的方法各有其优点，同时它们也都存在一定的局限性，只能在特定的条件下才能使用，不具有普适性，但是它们都极大的提高了曲线的光滑性与精确性。三次样条函数计算较简单，能通过所有采样点，且有二阶连续导数，是较理想的插值工具。武婕等人在“基于三次样条插值的算数傅里叶谐波分析法”中，对三次样条插值法进行了改进，解决了动态曲线的三次样条插值问题，使得测得的离散点可以形成一条平滑的曲线。

卜兆林在“关于FFT补零充位问题的进一步讨论”中，对一维傅里叶变换补零问题进行了探讨，然其重点为对一维测试数据进行补零，进而获得频谱值。钱辉等为了提高图像处理的运算速度，提出了一种结合频谱移位的二维傅里叶变换的FPGA实现方法，通过改变二维傅里叶变换的列变换起始位置，从矩阵的中间列开始进行列变换，使得二维傅里叶变换频谱的零频从矩阵的四个顶角移到中间位置，从而实现在完成二维傅里叶变换的同时完成一半的频谱移位，但其频谱移位是简单进行行、列变换，并未改变频谱幅值。

发明内容

本发明针对现有技术中的不足，提供一种基于快速傅里叶变换的二维离散数据平滑方法，旨在解决二维离散数据的快速平滑问题，采用的技术方案如下：

一种基于快速傅里叶变换的二维离散数据平滑方法，包括以下步骤：

步骤1：将待平滑的二维离散数据通过二维傅里叶变换由空间域变换到波数域，以得到原数据角谱；

步骤2：以补零的方式对步骤1中得到的原数据角谱进行扩展，得到扩展角谱；

步骤3：将步骤2中得到的扩展角谱通过二维傅里叶逆变换由波数域变换到空间域，以得到平滑后的二维离散数据。

进一步地，步骤1中原数据的角谱根据以下公式计算获取：