[交流] MATALB生成具有指定自相关函数的粗糙表面 已有5人参与

参考文献是陈辉、胡元中的《粗擦表面计算机模拟》，，这个是我模拟的图像结果，，文献中的结果见附件 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 有没有哪位大神能帮忙看看，问题出在哪儿 更多 0 1390271227_63743.jpg 1390271252_165409.jpg

回复此楼