求助地形谱分析

s7udY

请教一个大家一个问题。

假如我有一个n*1的单波数地形数据h，水平分辨率（r）是n*1的矩阵（各点分辨率不一样，单位是km）。求地形谱的话，先对h做快速傅里叶变换（h1=fft（h）），h1的模代表振幅（但是跟真实的振幅有差别，要对h1*2/n，其中h1（1）只除以n就可以了，即：h2=2*abs（h1）/n，h2（1）=h2（1）/2，h3=h2*h2），h3就是纵坐标（单位是m*m/km），我的一个疑问是，是不是还要对h3除以2pi，才能化成m*m*cycle/km？

接下来求波数，波数（f）是（0：n-1）.*r/n，我想问这里的f是不是也要乘以2pi？

我的matlab程序是：
h1=fft(h);
h2=2*abs(h1)/n;h2(1)=h2(1)/2;h3=h2*h2/2/pi;
f=2*pi*(0:n-1).*r/n;
loglog(f,h3)

s7udY

再加一句，我现在觉得不需要对h2除以n再乘以2了，新的程序如下：
h1=fft(h);
h2=abs(h1);
h3=h2*h2/2/pi;
f=2*pi*(0:n-1).*r/n;
loglog(f,h3)