24小時熱門版塊排行榜

返回列表

23點(diǎn)雨

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2183.8
散金: 23
紅花: 1
帖子: 76
在線: 34.6小時
蟲號: 2655666
注冊: 2013-09-14
性別: MM
專業(yè): 中藥抗腫瘤藥理

[求助] 求MATLAB生成二維高斯粗糙面的代碼已有1人參與

各位大神誰能提供生成二維高斯粗糙面的代碼，謝謝了！

回復(fù)此樓

» 猜你喜歡

271材料工程求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

急：求用Fortran或者其他語言編寫一個程序，處理下面的數(shù)據(jù)�。。。。。。。。。� 已經(jīng)有37人回復(fù)
高斯方程積分編程matlab或C語言已經(jīng)有15人回復(fù)
matlab中如何產(chǎn)生逆高斯分布的隨機(jī)數(shù)？已經(jīng)有9人回復(fù)
【Matlab編程】高斯光束經(jīng)透鏡聚焦，焦點(diǎn)處的三維光場分布已經(jīng)有111人回復(fù)
圖像處理高斯模糊化的MATLAB代碼已經(jīng)有5人回復(fù)
現(xiàn)在做到高斯函數(shù)的誤差函數(shù)一步，又是一個方程不知道用matlab怎么寫已經(jīng)有14人回復(fù)
請教向量Q滿足高斯分布，用matlab如何表示？已經(jīng)有11人回復(fù)
MATLAB中產(chǎn)生高斯白噪聲，涉及到awgn和wgn函數(shù) 已經(jīng)有3人回復(fù)
【分享】立方晶系標(biāo)準(zhǔn)投影圖+投影圖自動生成軟件【已搜索無重復(fù)】】已經(jīng)有220人回復(fù)
【求助】如何用Matlab產(chǎn)生高斯白噪聲已經(jīng)有6人回復(fù)

O(∩_∩)O~

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

osd1888

鐵桿木蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 5460.8
紅花: 1
帖子: 546
在線: 407小時
蟲號: 3640095
注冊: 2015-01-12
專業(yè): 電磁場與電路

【答案】應(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));

贊一下

回復(fù)此樓