本实验要求开发一个2-D FFT程序包，使其可以应用于后续的几个实验。这个程序需要完成的功能是用因子（-1）^(x+y)乘以输入图像以实现滤波的中心变换，并还要求用一个实矩阵乘以一个复数矩阵通过调用两个图像的乘法程序来实现对应元素的相乘，同时计算反傅立叶变换，得到的结果乘以（-1）^(x+y)并取其实部最后计算频谱。实验中用到傅立叶变换的基本公式，通过实验我们可以更加深刻的理解频域滤波的基础。-The experiment calls for the development of a 2-D F

该实验的目的是开发一个 2-D FFT程序包。要求程序能完成下面的功能: (a) 用因子 (-1)x+y 乘以输入图像，以实现滤波的中心化变换； (b) 用一个实矩阵乘以一个复数矩阵，即用实矩阵中的元素同时乘以复数矩阵对应位置上的复数的实部与虚部。 可以通过调用两个图像的乘法程序来实现对应元素的相乘； (c) 计算反付立叶变换； (d) 结果乘以 (-1)x+y ，并取其实部； (e) 计算频谱。 -The purpose of this ex

Jd=diffusion(J,method,N,K) J :待扩散的原图像 method : lin : 线性扩散 (常数 c=1). pm1 : perona-malik, c=exp{-(|grad(J)|/K)^2} pm2 : perona-malik, c=1/{1+(|grad(J)|/K)^2} rmp : 复数 K :边缘门限参数 N :迭代次数 dt:时间步长 (0 < dt <= 0.25, 缺省值为0.2)