请教一下，有关路面功率谱估计？

大鹏之举

为什么我做出得来功率谱估计是这样的？下面是我的程序？感觉和标准谱差别太大了？
A=importdata('');%导如高程数据
y=mean(A);%求取左轮高程的平均值
x=A-y;%去均值
C=importdata('‘);%导入含有采样车速度的数据块
z=averagespeed(C);%求取采样车的平均车速
N=20000; %有20000个点.
t=1000/z %求取采样时间
fs=20000/t;%求取采样频率
nfft=2^nextpow2(N)
P=abs(fft(x,nfft).^2/N);
d=0:round(nfft/2-1)
f=d*fs/nfft
s=f/z
figure(2)
loglog(s,P(d+1));
hold on
for i= 1 :8  ;
    x=16*4^(i-1);
            %参考空间频率n0下的路面功率谱密度
n0=0.1;             %参考空间频率n0
L=409.6;
l=0.1;
N=L/l;              %采样点数
n1=0.01;            %空间频率范围n1--nu
nu=3;
w=2;                %频率指数
n=linspace(0.01,3,N);
Gx=x*(n/n0).^(-w);  %路面功率谱密度拟合表达式
loglog(n,Gx)
hold on
end

[ 本帖最后由 大鹏之举 于 2009-3-21 21:09 编辑 ]