matlab模拟高斯随机粗糙表面

[manbuyunduan89]

http://www.cnki.com.cn/Article/CJFDTotal-RHMF200610015.htm，参考文献是陈辉、胡元中的《粗擦表面计算机模拟》，，这个是我模拟的图像结果，，文献中的结果
clear
clc

N=512;
L=0.05;
delta=0.05;%表面均方根粗糙度
betax=10;%x方向的相关长度
betay=10;
C=1;%功率谱密度

dx=L/N;dy=dx;
NN=-N/2:N/2-1;
[Nx,Ny]=meshgrid(NN,NN);
taux=dx.*Nx;tauy=dy.*Ny;

%%生成具有指定自相关函数的粗糙表面
eta=randn(N,N);%高斯分布白噪声
A=fft2(eta);%傅里叶变换
R=zeros(N,N);
R=delta^2*exp(-2.3*((taux/betax).^2+(tauy/betay).^2).^0.5);%自相关函数
Gz=1/(2*pi^2).*fft2(R);%功率谱密度函数
H=(Gz/C).^0.5;%传递函数
Z=H.*A;%表面高度的傅里叶变换
z=ifft2(Z);%表面高度分布

figure(1)
meshz(abs(z));
title('rough surface');
axis square

有没有哪位大神能帮忙看看，问题出在哪儿

[chengyi0923]

我也遇到同样的问题啊，求高手指点

[qq_14970105]

楼主，请问你的问题解决了没有，求程序啊。