版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）

返回列表

23点雨

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 2183.8
散金: 23
红花: 1
帖子: 76
在线: 34.6小时
虫号: 2655666
注册: 2013-09-14
性别: MM
专业: 中药抗肿瘤药理

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

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

回复此楼

» 猜你喜欢

三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复

» 本主题相关价值贴推荐，对您同样有帮助:

急：求用Fortran或者其他语言编写一个程序，处理下面的数据！！！！！！！！！！已经有37人回复
高斯方程积分编程matlab或C语言已经有15人回复
matlab中如何产生逆高斯分布的随机数？已经有9人回复
【Matlab编程】高斯光束经透镜聚焦，焦点处的三维光场分布已经有111人回复
图像处理高斯模糊化的MATLAB代码已经有5人回复
现在做到高斯函数的误差函数一步，又是一个方程不知道用matlab怎么写已经有14人回复
请教向量Q满足高斯分布，用matlab如何表示？已经有11人回复
MATLAB中产生高斯白噪声，涉及到awgn和wgn函数已经有3人回复
【分享】立方晶系标准投影图+投影图自动生成软件【已搜索无重复】】已经有220人回复
【求助】如何用Matlab产生高斯白噪声已经有6人回复

O(∩_∩)O~

1楼2015-02-03 11:05:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

osd1888

铁杆木虫 (正式写手)

应助: 3 (幼儿园)
金币: 5460.8
红花: 1
帖子: 546
在线: 407小时
虫号: 3640095
注册: 2015-01-12
专业: 电磁场与电路

【答案】应助回帖

感谢参与，应助指数 +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=1

N/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的回帖

23点雨

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 2183.8
散金: 23
红花: 1
帖子: 76
在线: 34.6小时
虫号: 2655666
注册: 2013-09-14
性别: MM
专业: 中药抗肿瘤药理

送红花一朵

引用回帖:

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的回帖

osd1888

铁杆木虫 (正式写手)

应助: 3 (幼儿园)
金币: 5460.8
红花: 1
帖子: 546
在线: 407小时
虫号: 3640095
注册: 2015-01-12
专业: 电磁场与电路

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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点雨

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 2183.8
散金: 23
红花: 1
帖子: 76
在线: 34.6小时
虫号: 2655666
注册: 2013-09-14
性别: MM
专业: 中药抗肿瘤药理

能不能按照你发的第一个代码的方法，改成二维的呢？？

赞一下

回复此楼

O(∩_∩)O~

5楼2015-03-10 15:02:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

23点雨

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 2183.8
散金: 23
红花: 1
帖子: 76
在线: 34.6小时
虫号: 2655666
注册: 2013-09-14
性别: MM
专业: 中药抗肿瘤药理

送红花一朵

引用回帖:

4楼: Originally posted by osd1888 at 2015-03-07 09:34:21
%********** 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,
...

能不能按照你发的第一个代码的方法，改成二维的呢？？

赞一下

回复此楼

O(∩_∩)O~

6楼2015-03-10 16:13:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 23点雨的主题更新

返回列表