| 查看: 1653 | 回复: 1 | |||
[交流]
用MATLAB或者C语言生成一维随机粗糙面
|
| 现在急需要粗糙面程序用你模拟电磁散射，MATLAB或者C语言都可以，希望大虫们帮个忙，感激不尽，呵呵！ |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有4人回复
-
基金委咋了？2026年的指南还没有出来？
已经有9人回复
-
基金申报
已经有5人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有5人回复
-
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
推荐一本书
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用Matlab设计一种最佳配方需要用哪种语言可以实现
已经有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-函数实现积分分离PID控制器的仿真
已经有4人回复
- 求助matlab编程
已经有13人回复
- 编写快速的MATLAB代码
已经有230人回复
- 【求助】同样一个问题用matlab和fortran编程计算，结果相差比较大，不知道为什么
已经有4人回复
- 【求助】如何将matlab中的randn函数转变为C语言形式的啊？
已经有3人回复
- 【求助】帮我编个程啊~C语言，C++，Matlab，VB...都行~
已经有14人回复
- 【求助】Weickert的各向异性扩散方程的滤波方法的matlab程序代码
已经有8人回复
osd1888
铁杆木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 5460.8
- 红花: 1
- 帖子: 546
- 在线: 407小时
- 虫号: 3640095
- 注册: 2015-01-12
- 专业: 电磁场与电路
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+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