版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

☆皓月當空☆

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 73.7
散金: 10
紅花: 1
帖子: 38
在線: 6.4小時
蟲號: 2173071
注冊: 2012-12-04
專業(yè): 激光

[交流] 用MATLAB或者C語言生成一維隨機粗糙面

現(xiàn)在急需要粗糙面程序用你模擬電磁散射，MATLAB或者C語言都可以，希望大蟲們幫個忙，感激不盡，呵呵！

回復此樓

» 猜你喜歡

287求調劑已經有5人回復
22408 344分求調劑一志愿華電計算機技術已經有3人回復
289求調劑已經有3人回復
317求調劑已經有10人回復
廣西大學家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調劑）已經有3人回復
學校已經提交到NSFC，還能修改嗎？已經有8人回復
08工學調劑已經有5人回復
281求調劑（0805）已經有25人回復
085600材料與化工已經有6人回復
265求調劑已經有8人回復

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

用Matlab設計一種最佳配方需要用哪種語言可以實現(xiàn) 已經有19人回復
這個方程組該怎么求解呢？用matlab or c語言？已經有3人回復
算法程序用matlab，C語言，軟件平臺，硬件平臺運行耗時之間的關系已經有9人回復
高斯方程積分編程matlab或C語言已經有15人回復
怎么在一個0矩陣中產生隨機分布的1小矩陣已經有4人回復
matlab 編程教程《C語言與matlab接口編程與實例》已經有7人回復
大家一般都到哪個網站下載matlab 和C語言程序？已經有7人回復
求助高手將以下C語言翻譯為Matlab語言，急急急！已經有12人回復
怎樣用Matlab調用已有的C語言程序呢？已經有3人回復
利用matlab中的S-函數(shù)實現(xiàn)積分分離PID控制器的仿真已經有4人回復
求助matlab編程已經有13人回復
編寫快速的MATLAB代碼已經有230人回復
【求助】同樣一個問題用matlab和fortran編程計算，結果相差比較大，不知道為什么已經有4人回復
【求助】如何將matlab中的randn函數(shù)轉變?yōu)镃語言形式的��？已經有3人回復
【求助】幫我編個程啊~C語言，C++，Matlab，VB...都行~ 已經有14人回復
【求助】Weickert的各向異性擴散方程的濾波方法的matlab程序代碼已經有8人回復

1樓 2014-11-26 17:15:07

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

osd1888

鐵桿木蟲 (正式寫手)

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

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

function [f,df,x]=rsgeng(N,rL,h,lc,seed);
% 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));

贊一下

回復此樓

2樓2015-02-04 21:13:16

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

相關版塊跳轉我要訂閱樓主 ☆皓月當空☆ 的主題更新