 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3491 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
23點(diǎn)雨新蟲 (小有名氣)
|
[求助]
求MATLAB生成二維高斯粗糙面的代碼 已有1人參與
|
|
| 各位大神 誰(shuí)能提供生成二維高斯粗糙面的代碼，謝謝了！ |
回復(fù)此樓» 猜你喜歡
08工學(xué)調(diào)劑
已經(jīng)有17人回復(fù)
-
340求調(diào)劑
已經(jīng)有4人回復(fù)
-
311求調(diào)劑
已經(jīng)有3人回復(fù)
-
0805 316求調(diào)劑
已經(jīng)有5人回復(fù)
-
食品專碩 一志愿雙一流 328
已經(jīng)有4人回復(fù)
-
材料調(diào)劑
已經(jīng)有6人回復(fù)
-
300分，材料，求調(diào)劑，英一數(shù)二
已經(jīng)有5人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有13人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 急：求用Fortran或者其他語(yǔ)言編寫一個(gè)程序，處理下面的數(shù)據(jù)�。。。。。。。。。�
已經(jīng)有37人回復(fù)
- 高斯方程積分編程matlab或C語(yǔ)言
已經(jīng)有15人回復(fù)
- matlab中如何產(chǎn)生逆高斯分布的隨機(jī)數(shù)？
已經(jīng)有9人回復(fù)
- 【Matlab編程】高斯光束經(jīng)透鏡聚焦，焦點(diǎn)處的三維光場(chǎng)分布
已經(jīng)有111人回復(fù)
- 圖像處理高斯模糊化的MATLAB代碼
已經(jīng)有5人回復(fù)
- 現(xiàn)在做到高斯函數(shù)的誤差函數(shù)一步，又是一個(gè)方程不知道用matlab怎么寫
已經(jīng)有14人回復(fù)
- 請(qǐng)教向量Q滿足高斯分布，用matlab如何表示？
已經(jīng)有11人回復(fù)
- MATLAB中產(chǎn)生高斯白噪聲，涉及到awgn和wgn函數(shù)
已經(jīng)有3人回復(fù)
- 【分享】立方晶系標(biāo)準(zhǔn)投影圖+投影圖自動(dòng)生成軟件【已搜索無(wú)重復(fù)】】
已經(jīng)有220人回復(fù)
- 【求助】如何用Matlab產(chǎn)生高斯白噪聲
已經(jīng)有6人回復(fù)
23點(diǎn)雨
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2183.8
- 散金: 23
- 紅花: 1
- 帖子: 76
- 在線: 34.6小時(shí)
- 蟲號(hào): 2655666
- 注冊(cè): 2013-09-14
- 性別: MM
- 專業(yè): 中藥抗腫瘤藥理
送紅花一朵 
3樓2015-03-05 18:17:50
osd1888
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 5460.8
- 紅花: 1
- 帖子: 546
- 在線: 407小時(shí)
- 蟲號(hào): 3640095
- 注冊(cè): 2015-01-12
- 專業(yè): 電磁場(chǎng)與電路
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +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
osd1888
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 5460.8
- 紅花: 1
- 帖子: 546
- 在線: 407小時(shí)
- 蟲號(hào): 3640095
- 注冊(cè): 2015-01-12
- 專業(yè): 電磁場(chǎng)與電路
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
23點(diǎn)雨: 金幣+20, ★★★★★最佳答案, 太謝謝啦！很有幫助！ O(∩_∩)O~ 2015-03-10 11:05:11
23點(diǎn)雨: 金幣+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' ) |
» 本帖已獲得的紅花（最新10朵）
4樓2015-03-07 09:34:21
23點(diǎn)雨
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2183.8
- 散金: 23
- 紅花: 1
- 帖子: 76
- 在線: 34.6小時(shí)
- 蟲號(hào): 2655666
- 注冊(cè): 2013-09-14
- 性別: MM
- 專業(yè): 中藥抗腫瘤藥理
5樓2015-03-10 15:02:53
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +5 | @taotao 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0854電子信息求調(diào)劑 +6 | α____ 2026-03-22 | 7/350 |
|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
|
[考研] 316求調(diào)劑 +7 | 梁茜雯 2026-03-19 | 7/350 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|