| 查看: 3482 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
23点雨新虫 (小有名气)
|
[求助]
求MATLAB生成二维高斯粗糙面的代码 已有1人参与
|
|
| 各位大神 谁能提供生成二维高斯粗糙面的代码，谢谢了！ |
回复此楼» 猜你喜欢
289求调剂
已经有3人回复
-
317求调剂
已经有10人回复
-
广西大学家禽遗传育种课题组2026年硕士招生（接收计算机专业调剂）
已经有3人回复
-
学校已经提交到NSFC，还能修改吗？
已经有8人回复
-
08工学调剂
已经有5人回复
-
281求调剂（0805）
已经有25人回复
-
085600材料与化工
已经有6人回复
-
265求调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有21人回复
-
295材料求调剂，一志愿武汉理工085601专硕
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急：求用Fortran或者其他语言编写一个程序，处理下面的数据！！！！！！！！！！
已经有37人回复
- 高斯方程积分编程matlab或C语言
已经有15人回复
- matlab中如何产生逆高斯分布的随机数？
已经有9人回复
- 【Matlab编程】高斯光束经透镜聚焦，焦点处的三维光场分布
已经有111人回复
- 图像处理高斯模糊化的MATLAB代码
已经有5人回复
- 现在做到高斯函数的误差函数一步，又是一个方程不知道用matlab怎么写
已经有14人回复
- 请教向量Q满足高斯分布，用matlab如何表示？
已经有11人回复
- MATLAB中产生高斯白噪声，涉及到awgn和wgn函数
已经有3人回复
- 【分享】立方晶系标准投影图+投影图自动生成软件【已搜索无重复】】
已经有220人回复
- 【求助】如何用Matlab产生高斯白噪声
已经有6人回复
osd1888
铁杆木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 5460.8
- 红花: 1
- 帖子: 546
- 在线: 407小时
- 虫号: 3640095
- 注册: 2015-01-12
- 专业: 电磁场与电路
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
23点雨: 金币+20, ★★★★★最佳答案, 太谢谢啦！很有帮助！ O(∩_∩)O~ 2015-03-10 11:05:11
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' ) |
23点雨4楼2015-03-07 09:34:21
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朵）
2楼2015-02-04 11:26:42
23点雨
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2183.8
- 散金: 23
- 红花: 1
- 帖子: 76
- 在线: 34.6小时
- 虫号: 2655666
- 注册: 2013-09-14
- 性别: MM
- 专业: 中药抗肿瘤药理
3楼2015-03-05 18:17:50
23点雨
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2183.8
- 散金: 23
- 红花: 1
- 帖子: 76
- 在线: 34.6小时
- 虫号: 2655666
- 注册: 2013-09-14
- 性别: MM
- 专业: 中药抗肿瘤药理
5楼2015-03-10 15:02:53
