Day6 Numerical simulation of optical wave propagation之通过随机介质(如大气湍流)的传播(二)
2.蒙特卡洛相位屏
大气折射率变化是一个随机的过程,通过大气的光程长度也同样是随机的。因此,湍流模型仅给出统计平均值,如折射率变量的结构函数和功率谱。
建立大气相位屏的问题就是产生随机过程独立表达式的问题,即相位屏的建立是通过把计算机产生的随机数变换为采样点网格上的两维相位值阵列来实现的,相位值阵列具有与湍流引起的相位变化相同的统计特性。
为了在有限的网格上生成相位屏,相位通常被写成各种基底函数的权重加和,常用于这一目的的基组为泽尔尼克多项式和傅里叶数列(FS),这两种基组各有优缺点。
(一)FT方法——最普遍的相位屏生成方法
湍流诱导光学相位为
改写成傅里叶级数
傅里叶系数Cn,m 服从具有零平均值和方差的环形复高斯统计,其中傅里叶系数Cn,m 的方差推导过程如下:
将9.78带入到9.77,可得如下结果
其中如下为改进的von-karman refractive-index PSD’s ,单位为cycles/m
总结程序过程:首先利用MATLAB的randn函数,生成零平均值和单位方差的高斯随机数,然后乘以(9.79)给出的方差均根值便得到了傅里叶系数Cn,m的随机曲线。进而利用二维DIFT得到相位屏。详细代码见ft_sh_phase_screen子函数。
(二)分谐波强化的FT方法——
上述FT方法中采用的二维DIFT不能生成准确的相位屏,最大偏差发生在大空间间隔处,即低空间频率处。——由(9.51)生成的相位PSD在较低空间频率具有很高的功率。实际上很多文献表明经常不能对足够低的空间频率进行采样来准确表征低阶模式,如倾斜。
基本原理:首先采用上述二维DIFT法生成相位屏,然后在相位屏中间的低频处这一小范围内采用分谐波法生成低频相位屏,即利用多个分谐波的和得到相位屏中间的低频处的相位,相应低频处相位的公式如下。详细代码见ft_phase_screen子函数。
注:分谐波是指频率等于一个周期性振荡基频的整分数的正弦分量。例如,频率等于基频二分之一的波称为二次分谐波,三分之一的波称为三次分谐波等。
最新文章
- cmder添加右键菜单
- UIImage 调整图片大小
- Laravel 5 基础(一)- Laravel入门和新建项目
- BZOJ_3670_[NOI2014]_动物园_(kmp)
- IT增值服务,客户案例(一)--山东青岛在职人士,2年.Net经验,转Java开发半年
- css3中display和box小结
- kubernetes入门之kube-proxy实现原理
- Python小问题汇总
- Linux安装配置VPN服务器
- openjudge8465:马走日 [搜索]
- HTML-移动开发技巧 响应式布局 弹性布局
- Android NDK Application.mk(中文翻译)
- JPA注解实体类,给表添加创建时间,更新时间,id的生成以及创建唯一约束
- PS抠图简单粗暴2种方法
- 洛谷P1402 酒店之王
- Android关于LinearLayout和RelativeLayout背景设置的区别
- Js_闭包跟作用域
- C# lambda表达式参数的正确使用姿势
- Maven 环境变量设置
- 哈希表(散列表)—Hash表解决地址冲突 C语言实现
热门文章
- 第52章 撤销端点(Revocation Endpoint) - Identity Server 4 中文文档(v1.0.0)
- JSON 与 JS 对象的区别与对比
- Android开发——使用intent传递对象
- noi.ac #289. 电梯(单调队列)
- 从.Net到Java学习第十二篇——SpringBoot+JPA提供跨域接口
- ionic3 Loading组件的用法
- jQuery省市区三级联动菜单
- angular部署到iis出现404解决方案
- shell脚本批量ssh登陆主机并执行命令
- 高德地图 Service 创建服务 USERKEY_PLAT_NOMATCH