应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 求MATLAB生成二维高斯粗糙面的代码
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
51/1返回列表
查看: 3397  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

23点雨

新虫 (小有名气)

[求助] 求MATLAB生成二维高斯粗糙面的代码已有1人参与

各位大神 谁能提供生成二维高斯粗糙面的代码，谢谢了！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
O(∩_∩)O~
1楼2015-02-03 11:05:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

23点雨

新虫 (小有名气)
送红花一朵
引用回帖:
2楼: Originally posted by osd1888 at 2015-02-04 11:26:42
function =rsgeng(N,rL,h,lc,seed);
% RSGENG generates 1D Gaussian random rough surfaces with Gaussian Spectrum.
% INPUT:
% N=total number of sample points
% rL=rough surface length
% h=rms height ...

谢谢！这个我早就有了，想要一个二维高斯面的，也就是在X,Y方向都变化的。你能提供么？
赞一下(1人) 回复此楼
O(∩_∩)O~
3楼2015-03-05 18:17:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

osd1888

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
function [f,df,x]=rsgeng(N,rL,h,lc,seed);
% RSGENG generates 1D Gaussian random rough surfaces with Gaussian Spectrum.
% INPUT:
% N=total number of sample points
% rL=rough surface length
% h=rms height
% lc=correlation length
% seed=seed of random number generator
%
% OUTPUT:
% f=rough surface profile
% df=df/dx
% x=sample points on the surface
%

randn('seed',seed);
y=randn(N,1);
for n=1N/2-1);
  bh(n)=(y(2*n-1)+i*y(2*n))/sqrt(2);
end;

bhc=conj(bh);
bhf=fliplr(bhc);
bi=[bh y(N-1) bhf y(N)];
kx=2*pi*[-N/2+1:1:N/2]/rL;
y1=sqrt(wk(kx,h,lc));
y=y1*sqrt(2*pi*rL);
b=y.*bi;
xs=[b(N/2+1:1:N) b(1:1:N/2)];
xt=[xs(N),xs(1:1:N-1)];
ft=ifft(xt,N);
ft=ft*N/rL;
fs=[ft(2:1:N),ft(1)];
f=[fs(N/2+1:1:N) fs(1:1:N/2)];
f=real(f);
dx=rL/N;
x=[-N/2+1:1:N/2]*dx;
n=2:N-1;
df1=(f(n+1)-f(n-1))/(2*dx);
df=[(f(2)-f(N))/(2*dx),df1,(f(1)-f(N-1))/(2*dx)];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Gaussian spectral density %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y=wk(kx,h,lc)

y=h^2*lc*exp(-(kx*lc*0.5).^2)/(2*sqrt(pi));
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

23点雨
2楼2015-02-04 11:26:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

osd1888

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
23点雨: 金币+20, ★★★★★最佳答案, 太谢谢啦！很有帮助！ O(∩_∩)O~ 2015-03-10 11:05:11
%********** INPUTS ***********************************
%SELECT GRID SIZE: N=NUMBER OF GRID NODES
N=128;   %For efficiency (FFTs) N should be a power of 2
%SELECT GRID SPACING, CORRELATION LENGTH,
%AND SURFACE HEIGHT VARIANCE, ALL IN METERS
%(NOTE: delx should be smaller than lc for )
%(      good surface representation)
lambda=10./3.;             %radiation wavelength (m)
lc  =1.000*lambda;         %transverse correlation length (m)
zdev=0.100*lambda;         %surface height variance (m)
delx=lambda/7.;            %cell size (grid spacing) (m)
%********** END INPUTS ********************************
%CREATE 2-D GAUSSIAN FILTER
%normalize correlation length to cel1 size
lc=lc/delx;
%x and y are dimensionless for now
x=linspace(-(N-1)/2,(N-1)/2,N);y=x;
%x and y coords.  for each grid pt.
[X,Y]=meshgrid(x,y);
%2-D filter
G=exp((-X.^2 -Y.^2)/(lc.^2/2.))/(sqrt(pi)*lc/2.);
%scale x and y axes to be meters for future use
x=x*delx;y=y*delx;
%CREATE GAUSSIAN DISTRIBUTED, UNCORRELATED RANDOM SURFACE
R=randn(size(G)) ;
%USE 2-D CIRCULAR CONVOLUTION TO MAKE GAUSSIAN CORRELATED SURFACE
%WITH CORRELATION LENGTH lc (celIs) aud STANDARD DEVIATION 1.
Z=ifft2(fft2(R).*fft2(G));%Circular convolution of R ald G via FFT
Z=real(Z);               %Convolution of two real arrays is real
Z=Z*zdev;                 %scale to obtain desired std.dev.(m)
%PLOT THE SURFACE
figure
mesh(x,y,Z)                  %random surface
axis ([min(x),max(x),min(y),max(y),min(min(Z))*2,max(max(Z))*2]);
xlabel('x (m)'); ylabel('y (m)'); zlabel('z (m)');
title ('Random Surface' )
赞一下(3人) 回复此楼

» 本帖已获得的红花（最新10朵）

23点雨
4楼2015-03-07 09:34:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

23点雨

新虫 (小有名气)
能不能按照你发的第一个代码的方法，改成二维的呢？？
赞一下 回复此楼
O(∩_∩)O~
5楼2015-03-10 15:02:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭