| 5 | 1/1 | 返回列表 |
| 查看: 3499 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
23點(diǎn)雨新蟲 (小有名氣)
|
[求助]
求MATLAB生成二維高斯粗糙面的代碼 已有1人參與
|
|
| 各位大神 誰能提供生成二維高斯粗糙面的代碼,謝謝了! |

鐵桿木蟲 (正式寫手)
|
%********** 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' ) |
鐵桿木蟲 (正式寫手)
|
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)); |
新蟲 (小有名氣)

新蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +3 | kongweilin 2026-03-26 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +8 | Mr. Z 2026-03-25 | 8/400 |
|
|
[考研] 0856求調(diào)劑 +5 | zhn03 2026-03-25 | 6/300 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 7/350 |
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +4 | 13831862839 2026-03-24 | 5/250 |
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +7 | 722865 2026-03-23 | 7/350 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[基金申請(qǐng)] 請(qǐng)教下大家 2026年國家基金申請(qǐng)是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 材料學(xué)碩297已過四六級(jí)求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|