【光波电子学】MATLAB仿真厄米特-高阶高斯光束以及绘制其光强分布图
目录
- 1 物理模型
- 2 MATLAB实现
- 2.1 光斑效果展示
- 2.2 高斯光束光强分布的三维模型
- 2.3 代码实现
- 3 知识拓展
- 3.1 高斯分布是什么?
- 3.2 单独实现高阶高斯分布
1 物理模型
在激光器的方形镜对称共焦腔中,厄米特—高斯分布描述了腔镜上自再现模的场分布,对其进行数值模拟,可以观察不同阶数下本征模的场分布情况
(x,y)表示镜面上的任意点坐标,腔长:L,光波波长:λ,原点选在镜面中心,则在近轴近似条件下 L,R>>α>>λ ,腔镜镜面上的电场复振幅可以表示为
umn≈CmnHm(X)Hn(Y)eX2+Y22u_{mn} \approx C_{mn}H_m(X)Hn(Y)e^{\frac{X^2+Y^2}{2}}umn≈CmnHm(X)Hn(Y)e2X2+Y2
其中,m=0, 1, 2, …和 n=0, 1, 2…,m、n 为横截面上 x、y 方向出现的节点;CmnC_{mn}Cmn为与 m、n 有关的常量Hm(X)、Hn(Y)为厄密多项式,可以表示为
Hm(X)=(−1)meX2dm(e−X2)dXmH_m(X) =(-1)^me^{X^2}\frac{d^m(e^{-X^2})}{dX^m}Hm(X)=(−1)meX2dXmdm(e−X2)
H0(X)=1,H1(X)=2X,H2(X)=4X2−2,H3(X)=−8x3+12xH_0(X) = 1,H_1(X) = 2X,H_2(X)=4X^2-2,H_3(X)=-8x^3+12xH0(X)=1,H1(X)=2X,H2(X)=4X2−2,H3(X)=−8x3+12x
也可以写为
umn≈CmnFm(X)F(Y)u_{mn}\approx C_{mn}F_m(X)F(Y)umn≈CmnFm(X)F(Y)
F(X)=Hm(X)e−X2/2F_(X) = H_m(X)e^{-X^2/2}F(X)=Hm(X)e−X2/2
F(Y)=Hn(X)e−Y2/2F_(Y) = H_n(X)e^{-Y^2/2}F(Y)=Hn(X)e−Y2/2
因此不同阶数的光场分布对应光斑分布可以表示为:
激光的模式常用微波中的表示方法,即TEMmnTEM_{mn}TEMmn ,其中TEM00TEM_{00}TEM00是基模,可以看出 m 和 n 的值越大,光场也越向外扩展。当 m=0,n=0 时,光场分布为基模高斯光束TEM00TEM_{00}TEM00,其传播场的振幅分布可写为
U00=exp(−21+ζ2x2+y2ws2)U_{00} = exp(-\frac{2}{1+\zeta ^2} \frac{x^2+y^2}{w_s^2})U00=exp(−1+ζ22ws2x2+y2)
ζ=2z/L,ws=λL/π\zeta = 2z/L,w_s = \sqrt{\lambda L/\pi}ζ=2z/L,ws=λL/π
z为传播距离,L为腔长,λ为波长z为传播距离,L为腔长,\lambda 为波长z为传播距离,L为腔长,λ为波长
其零点是对称共焦腔的中心。高斯光束的传播呈现双曲线型的形态,波长一定情况下,共焦腔内 z 处的截面半径大小仅取决于腔长 L。高斯光束的束腰半径可以写为
w0=λL2π=12wsw_0 = \sqrt{\frac{\lambda L}{2\pi}}= \frac{1}{\sqrt{2}}w_sw0=2πλL=21ws
2 MATLAB实现
2.1 光斑效果展示
2.2 高斯光束光强分布的三维模型
2.3 代码实现
clear all
close allx=[-5:0.1:5].*1e-3; % 坐标范围
lambda=632.8e-9; % 波长
L=10; % 腔长
X=x.*sqrt(2*pi./(lambda.*L));
[X1,Y1]=meshgrid(X,X); % X和Y的坐标矩阵%% 厄米特函数
% H0-H3
H0X=1;
H0Y=1;
H1X=2.*X1;
H1Y=2.*Y1;
H2X=4.*X1.^2-2;
H2Y=4.*Y1.^2-2;
H3X=-8.*X1.^3+12.*X1;
H3Y=-8.*Y1.^3+12.*Y1;%% 厄米特——高斯函数
% F0-F3
FX0=H0X.*exp(-X1.^2-2);
FY0=H0Y.*exp(-Y1.^2-2);
FX1=H1X.*exp(-X1.^2-2);
FY1=H1Y.*exp(-Y1.^2-2);
FX2=H2X.*exp(-X1.^2-2);
FY2=H2Y.*exp(-Y1.^2-2);
FX3=H3X.*exp(-X1.^2-2);
FY3=H3Y.*exp(-Y1.^2-2);%% 振幅分布u00=FX0.*FY0; % TEM00
subplot(3,4,1)
imagesc(abs(u00)) u10=FX1.*FY0; % TEM10
subplot(3,4,2)
imagesc(abs(u10)) u20=FX2.*FY0; % TEM20
subplot(3,4,3)
imagesc(abs(u20)) u30=FX3.*FY0; % TEM30
subplot(3,4,4)
imagesc(abs(u30)) u01=FX0.*FY1; % TEM01
subplot(3,4,5)
imagesc(abs(u01)) u11=FX1.*FY1; % TEM11
subplot(3,4,6)
imagesc(abs(u11)) u21=FX2.*FY1; % TEM21
subplot(3,4,7)
imagesc(abs(u21)) u31=FX3.*FY1; % TEM31
subplot(3,4,8)
imagesc(abs(u31)) u02=FX0.*FY2; % TEM02
subplot(3,4,9)
imagesc(abs(u02)) u12=FX1.*FY2; % TEM12
subplot(3,4,10)
imagesc(abs(u12)) u22=FX2.*FY2; % TEM22
subplot(3,4,11)
imagesc(abs(u22)) u33=FX3.*FY3; % TEM33
subplot(3,4,12)
imagesc(abs(u33)) colormap gray; %显示灰度图
3 知识拓展
3.1 高斯分布是什么?
说白了就是上图中右下角的山坡型的图形。这个山坡型就是下面三维高斯光束的横截面。所以高斯光束就和高斯分布结合起来了。
这时候,我将三种图对应起来,有没有发现共通点???
从左开始的第一张是高斯分布图,中间张是高斯光束光强三维分布图,最右边是高斯光束光强从三维图的顶部角度看上去的光斑。所以在下面要想实现光斑,只需要画出高斯分布的三维图形,从顶部角度视角展示图形就能得到光斑。
附上绘制这张图的代码
clear all
close allx=[-5:0.1:5].*1e-3; % 设定腔镜坐标
lambda=632.8e-9; % 设定激光器出光频率
L=10; % 设定激光腔腔长
X=x.*sqrt(2*pi./(lambda.*L));
Y= X
[X1,Y1]=meshgrid(X,X); % 获取X和Y的坐标矩阵% 设定厄米特函数,n和m只从0取到2 %
% 绘制三维的
H0X=1;
H1X=2.*X1;
H2X=4.*X1.^2-2;
H0Y=1;
H1Y=2.*Y1;
H2Y=4.*Y1.^2-2;% 三维公式
FX0=H0X.*exp(-X1.^2-2); % 生成厄米特——高斯函数
FX1=H1X.*exp(-X1.^2-2);
FX2=H2X.*exp(-X1.^2-2);FY0=H0Y.*exp(-Y1.^2-2);
FY1=H1Y.*exp(-Y1.^2-2);
FY2=H2Y.*exp(-Y1.^2-2);% 二维公式
H0X1=1;
H1X1=2.*X;
H2X1=4.*X.^2-2;
H0Y1=1;
H1Y1=2.*Y;
H2Y1=4.*Y.^2-2;F0X=H0X1.*exp(-X.^2-2); % 生成厄米特——高斯函数
F1X=H1X1.*exp(-X.^2-2);
F2X=H2X1.*exp(-X.^2-2);F0Y=H0Y1.*exp(-Y.^2-2);
F1Y=H1Y1.*exp(-Y.^2-2);
F2Y=H2Y1.*exp(-Y.^2-2);% 00
I2=FX0.^2.*FY0.^2; % 考虑X和Y两个方向分布
I1=F0X.^2.*F0Y.^2;
subplot(3,3,1)
plot(X,I1)
subplot(3,3,2)
surf(X1,Y1,I2)
subplot(3,3,3)
surf(X1,Y1,I2)
view(2)% 01
I2=FX0.^2.*FY1.^2; % 考虑X和Y两个方向分布
I1=F0X.^2.*F1Y.^2;
subplot(3,3,4)
plot(X,I1)
subplot(3,3,5)
surf(X1,Y1,I2)
subplot(3,3,6)
surf(X1,Y1,I2)
view(2)% 02
I2=FX0.^2.*FY2.^2; % 考虑X和Y两个方向分布
I1=F0X.^2.*F2Y.^2;
subplot(3,3,7)
plot(X,I1)
subplot(3,3,8)
surf(X1,Y1,I2)
subplot(3,3,9)
surf(X1,Y1,I2)
view(2)
3.2 单独实现高阶高斯分布
其中用到的高斯分布的公式
Hm(x)H_m(x)Hm(x)
F=e−x2+y2w2F = e^{-\frac{x^2+y^2}{w^2}}F=e−w2x2+y2
Hm2(x)F2H_m^{2}(x)F^2Hm2(x)F2
clear
clc%% 初始值
x = -2*pi:pi/100:2*pi;% 横坐标范围
w = 2 % F 函数中w取值
yy =x % 纵坐标范围%% 画光强分布
figure(2) % define figure
% 【1 1】个图
subplot(3,4,1);
plot([-10,10],[10,10],'linewidth',2)% 【1 2】个图
subplot(3,4,2);
y = fx(x,yy,w)
plot(x,y,'linewidth',2);
% 【1 3】个图
subplot(3,4,3); % subplot(x,y,n)x表示显示的行数,y表示列数,n表示第几幅图片
y = hf(x,yy,w,hx1(x),fx(x,yy,w))
plot(x,y,'linewidth',2);% 【1 4】个图
subplot(3,4,4);
y = hf2(x,yy,w,hx1(x),fx(x,yy,w))
plot(x,y,'linewidth',2);% 【2 1】个图
% 画45度斜线
subplot(3,4,5);
y = hx2(x)
plot(x,y,'linewidth',2);% 【2 2】个图
subplot(3,4,6);
y = fx(x,yy,w)
plot(x,y,'linewidth',2);% 【2 3】个图
subplot(3,4,7);
y = hf(x,yy,w,hx2(x),fx(x,yy,w))
plot(x,y,'linewidth',2);% 【2 4】个图
subplot(3,4,8);
y = hf2(x,yy,w,hx2(x),fx(x,yy,w))
plot(x,y,'linewidth',2);% 【4 1】个图
% 画y=x^2-10
subplot(3,4,9);
y = hx3(x)
plot(x,y,'linewidth',2);% 【4 2】个图
subplot(3,4,10);
y = fx(x,yy,w)
plot(x,y,'linewidth',2);% 【4 3】个图
subplot(3,4,11);
y = hf(x,yy,w,hx3(x),fx(x,yy,w))
plot(x,y,'linewidth',2);% 【4 4】个图
subplot(3,4,12);
y = hf2(x,yy,w,hx3(x),fx(x,yy,w))
plot(x,y,'linewidth',2);%% 函数
% 绘制水平线函数
function r=hx1(x)r =1
end
% 绘制45度直线函数
function r=hx2(x)r =2.*x
end% 绘制凹曲线函数
function r=hx3(x)r =4.*x.^2-2
end% 绘制F的函数
function r =fx(x,y,w)r = exp(-(x.*2+y.^2)./(w^2))
end% 绘制Hm(x)*F
function r =hf(x,y,w,source,fun)r = fun.*source
end
% 绘制Hm(x)^2*F^2的函数
function r =hf2(x,y,w,source,fun)r = (fun.^2).*(source.^2)
end
对于绘制出来的图形,并没有理想中的那么均匀,但是趋势是大致相同的,有两个原因
(1)原始教材中的图形,是理想化的图像,只是为了突出表述了每个图形的特点。
(2)本人绘制的图形,选择的参数只是随机选择的,比如F函数中w和y的取值仅仅是随机取的。如果想要实现最理想化的图形,需要慢慢的调参
【光波电子学】MATLAB仿真厄米特-高阶高斯光束以及绘制其光强分布图相关推荐
- 基于Matlab有限差分的高阶迭代最小二乘积分的波前重构算法
该算法根据x与y方向波前斜率,可对其波前进行重构. 一.算法验证 首先,分别得到x方向与y方向的波前斜率,分别如下图所示: 图1 X方向波前斜率 图2 Y方向波前斜率 基于有限差分的高阶迭代最小二乘积 ...
- matlab仿真高阶传递函数,Matlab/Simulink动力学系统建模与仿真(第2版)
绪论 第1章 系统建模与仿真基础 1.1 系统仿真模型框图表示法 1.1.1 基本仿真元件 1.1.2 简单仿真框图结构 1.2 拉普拉斯变换 1.2.1 拉普拉斯变换的定义及其性质 1.2.2 拉普 ...
- 使用matlab求高阶累积量
本文介绍如何使用matlab自带的高阶累积量函数求取一个随机过程的高阶累积量,运行demo之前确保matlab中已经安装了HOSA的工具箱. 安装过程中经常会碰到的两个问题: 1.info.xml出问 ...
- matlab凸透镜防仿真,西安邮电大学matlab仿真解剖.doc
西 安 邮 电 大 学 专业课程设计报告书 院系名称:电子工程学院学生姓名:李 群学号专业名称:光信息科学与技术班 级:光信1103实习时间:2014年4月8日至2014年4月18日 课程设计题目: ...
- 激光光束传输影响matlab仿真,matlab仿真光束的传输特性详解.doc
课程设计题目: 用matlab仿真光束的传输特性. 任务和要求 用matlab仿真光束通过光学元件的变换. 设透镜材料为k9玻璃,对1064nm波长的折射率为1.5062,镜片中心厚度为3mm,凸面曲 ...
- 定时同步 matlab仿真,高阶QAM定时同步的MATLAB仿真及其FPGA实现
目前,基于软件无线电的数字接收机,其定时同步主要采用异步采样恢复,即采样时钟独立工作,通过估算定时误差,控制内插滤波器内插出最佳采样时刻的值.常用的定时误差估计算法主要有:迟早门算法.米勒-穆雷算法. ...
- 【高阶累积量】基于高阶累积量的信噪比盲估计法的matlab仿真
1.软件版本 matlab2013b 2.系统描述 所谓高阶统计量,通常是指信号的高阶矩,高阶累积量等信号统计信息的统称.在这些高阶统计量中,高阶累积量具有十分重要的重用,不同阶数的累积量可以反映出信 ...
- 【调制识别】基于高阶累积量的调制方式识别算法matlab仿真,课题识别BPSK, QPSK, 16QAM, 64QAM
目录 一.理论基础 二.核心程序 三.仿真测试结果 一.理论基础 由于数字调制信号越来越多地应用于通信信号处理领域,因此对数字信号调制识别的研究也越来越多.传统的调制识别的判决方法有:决策判决法.高阶 ...
- 从光波叠加到条纹分布的matlab仿真,基于Matlab仿真算法的光源空间相干性研究
王帆, 张永安, 阳胜, 郭胤初, 张亚萍. . 基于Matlab仿真算法的光源空间相干性研究. 光学学报, 2017, 54(9): 092601-. Wang Fan, Zhang Yongan, ...
- matlab 实现高阶动态电路,高阶电路动态特性的仿真分析
作者/ 贺为婷 杨建华 西安工业大学 电子信息工程学院(陕西 西安 710032)本文引用地址:http://www.eepw.com.cn/article/201612/340860.htm 摘要: ...
最新文章
- 【2020】清华大学《高级机器学习》课件和专家特邀报告(附pdf下载)
- 多样性计算时代,怎样的技术生态才能满足发展需求
- 《我也能做CTO之程序员职业规划》写作过程重播之二
- js及jQuery实现checkbox的全选、反选和全不选
- 我们的2008。。。。。。
- Learning to Rank 中Listwise关于ListNet算法讲解及实现
- 八大基本数据类型对应的八大包装类(含对应面试题解析)
- linux 怎么管理文件夹,Linux 是如何管理目录文件?
- 如何选择正确的控制系统?PLC和DCS各有不同
- 基于tcp实现远程执行命令
- docker打包部署flask镜像
- 【RWM】基于matlab路由无线传感器网络模拟随机路点运动模型【含Matlab源码 1565期】
- CSDN资源下载积分规则
- Maven经验分享(一)安装部署
- 文字转语音怎么在线生成MP3格式的音频?
- 苏州大学计算机专业考研报录比,【图片】18年苏州大学计算机872考研经验分享【苏州大学研究生吧】_百度贴吧...
- 小菜鸡跟着小破站黑马程序员学习JAVA的真实经历、感受与收获
- 工程测量(建筑物施工、监测)
- 全球及中国智能X光安检机行业研究及十四五规划分析报告
- python 步数_用python如何修改微信和支付宝每天走路的步数
热门文章
- 网马(原理)及制作全过程
- 解决FTP服务器FileZilla server中文乱码问题
- 局域网ip扫描工具_IP Scanner Pro——ip局域网扫描工具
- vbs脚本学习整人Demo
- 透明背景flash html5,DW网页设计:添加透明背景的 Flash
- Tar的详细用法(转自Linux伊甸园)
- android媒体--stagefright概述
- Android车载应用开发与分析(11)- 车载Android应用开发入门指南
- 《Head First Android》读后感,电子书PDF下载
- plc c语言做报警,wincc声音报警C脚本版