- 博客(17)
- 收藏
- 关注
RSS订阅原创 Matlab：拉盖尔-高斯光束光阑衍射
代码：clcclear allclose all%%N = 200; %取样点数lambda = 632e-6; %波长632nmk = 2*pi/lambda; %波数xo = linspace(-2,2,N); yo = linspace(-2,2,N);[x0,y0] = meshgrid(xo,yo);[theta,r] = cart2pol(x0,y0);w0 = 0.5; %束腰半径p = 0;z = 1000;
2020-08-14 10:53:11
10010
9
原创 Matlab：拉盖尔-高斯光束单缝衍射
代码：clcclear allclose all%% 拉盖尔-高斯光束N = 100; %取样点数lambda = 632e-9; %波长632nmk = 2*pi/lambda; %波数x0 = linspace(-2e-2,2e-2,N); y0 = linspace(-2e-2,2e-2,N);[X0,Y0] = meshgrid(x0,y0);[theta,r] = cart2pol(X0,Y0);w0 = 1e-3;
2020-08-13 19:35:56
2980
5
原创 Matlab：厄米-高斯光束合成拉盖尔-高斯光束
代码：clcclear allclose all%%N = 300; %采样点数lambda = 632e-9; %波长632nmk = 2*pi/lambda; %波数w0 = 3e-3; %束腰半径3mmZ_R = pi*w0^2/lambda; %瑞利长度row = linspace(-3*w0,3*w0,N); col = row;[x,y] = meshgrid(row,col);[theta,r] = ca
2020-05-28 20:59:27
12140
21
原创 Matlab：非局部非线性介质中两个同轴螺旋椭圆光束对光学图案的旋转
代码：clcclear allclose all%%N = 500; %采样点数A = 1; %振幅lambda = 632e-9; %波长632nmrow = linspace(-5,5,N); col = row;[x,y] = meshgrid(row,col);b = 2; c = 1;w_c = (b^2+c^2)/(2*b^2*c^2); %临界角速度THETA_c =
2020-05-27 15:59:34
1111
1
原创 Matlab：倍频过程中复合涡旋拓扑荷数的倍频效应
clcclear allclose all%%N = 300; %取样点数lambda = 632e-9; %波长632nmk = 2*pi/lambda; %波数x = linspace(-4,4,N);y = linspace(-4,4,N);[X,Y] = meshgrid(x,y);[theta,r] = cart2pol(X,Y);w0 = 3;p = 0; %径向量子数%% m_A = 0, m_B = 1
2020-05-16 09:45:37
1595
2
原创 Matlab：完美涡旋光束
代码：clcclear allclose all%% 振幅光学相位元件产生完美涡旋光场N = 300;row = linspace(-1.5e-6,1.5e-6,N); col = linspace(-1.5e-6,1.5e-6,N);[x,y] = meshgrid(row,col);[phi,rho] = cart2pol(x,y);lambda = 632e-9; %波长632nmk = 2*pi/lambda; %波数f = 0.5; %
2020-05-15 09:58:58
12858
7
原创 Matlab：拉盖尔-高斯光束杨氏双缝干涉
代码：clcclear allclose all%% L-G光束双缝干涉N = 300; %取样点数lambda = 632e-9; %波长632nmk = 2*pi/lambda; %波数x = linspace(-2,2,N);y = linspace(-2,2,N);[X,Y] = meshgrid(x,y);[theta,r] = ...
2020-05-14 16:25:28
2932
2
原创 Matlab：拉盖尔高斯光束的干涉
代码：clcclear allclose all%% 拉盖尔-高斯光束与平面光波干涉N = 300; %取样点数lambda = 632e-9; %波长632nmk = 2*pi/lambda; %波数x = linspace(-1e-2,1e-2,N);y = linspace(-1e-2,1e-2,N);[X,Y] = meshgrid(...
2020-05-05 08:52:13
8031
1
原创 Matlab：常见涡旋光束仿真
代码：clcclear allclose all%% 环形涡旋光束N = 200;lambda = 632e-9; %波长为632nmk = 2*pi/lambda; %波数w0 = 3; %束腰半径x = linspace(-10,10,N);y = linspace(-10,10,N);[X,Y] = meshgrid(x,y);[t...
2020-05-04 20:56:33
32491
43
原创 Matlab：贝塞尔高斯光束自由传输matlab仿真
程序部分：clcclear allclose all%%N = 100;lambda = 1064e-6; %波长1064nmk = 2*pi/lambda;x0 = linspace(-1,1,N); y0 = linspace(-1,1,N);[x1,y1] = meshgrid(x0,y0);[THETA,RHO] = cart2pol(x1,...
2020-05-04 09:56:40
22595
28
原创 Matlab：Fox_Li数值迭代法求激光器谐振腔本征模式
谐振腔的经典理论仅给出了部分简单腔型的模式解析解。对于平行平面腔以及在激光器的不断发展过程中所涌现的许多新型结构谐振腔通常是没有解析结果的 ,必须采用各种数值模拟方法进行求解 ,如 Fox-Li 迭代法、快速傅立叶变换法(FFT)、等效透镜波导法、特征向量法 、有限元法(FEM)和有限差分法(FDM)等 。...
2020-05-04 09:56:18
9039
原创 matlab:夫琅禾费衍射简单仿真
clcclear allclose all%% 光源部分N = 200;Diameter = 10;x = linspace(-5,5,N); y = linspace(-5,5,N);[x1,y1] = meshgrid(x,y);E1 = ones(N);E1(sqrt(x1.^2+y1.^2)>Diameter/2) = 0;I1 = E1.*conj(E1);I...
2020-05-04 09:55:45
8584
原创 matlab：杨氏双缝干涉
基本原理：杨氏干涉实验是两点光源干涉实验的典型代表。杨氏干涉实验以极简单的装置和巧妙构思实现了普通光源干涉。无论从经典光学还是从现代光学的角度来看，杨氏实验都具有十分重要的意义。杨氏双缝实验的装置如图2-18所示，按照惠更斯-菲涅耳原理，线光源S上的点将作为次波源向前发射次波（球面波），形成交叠的波场。在较远的地方放置一观察屏，屏上可以观测到一组几乎是平行的直线条纹。clcclose a...
2020-05-04 09:55:27
8489
2
原创 matlab：光场传输的各种算法
随着计算机的飞速发展，数值仿真在许多光学领域都有重要的应用。比如：1、在研究超分辨时需得到光场经特定结构后精确的聚焦情况2、研究激光相干合束时需计算阵列光束远场合束光斑的能量分布情况3、在设计衍射光学元件时需循环计算光场经衍射元件后的远场分布4、在设计激光谐振腔时需计算腔内光束经过多次反射后的稳定光场分布在这些研究中，普遍包含光场传输这一最基本的物理过程，因此研究光场传输快速、准确的实...
2020-05-04 09:54:20
17775
20
原创 Matlab：用空间光调制器产生不同拓扑荷的涡旋光束的相息图
clcclear allclose all%%c=-516:516;r=-516:516;l=1;[x,y]=meshgrid(c,r);[theta,r]=cart2pol(x,y);figure;for l=1:8 subplot(2,4,l) g=mod(l*theta,2*pi); imshow(g,[]) l=l+1;endfigu...
2020-05-04 09:52:42
8276
11
原创 Matlab：涡旋光束与平面光波、球面光波干涉
与平面光波干涉：clcclear allclose all%% 涡旋光束与平面光波干涉N = 200; %取样点数lambda = 632e-9; %波长632nmk = 2*pi/lambda; %波数x = linspace(-500,500,N);y = linspace(-500,500,N);[X,Y] = meshgrid(x,y...
2020-05-04 09:52:04
11975
12
原创 高斯光束的简单matlab仿真
clcclear allclose all%% 假设在z=0处有一振幅为A的高斯光束N = 100;lambda = 1064e-6; %波长为1064nmk = 2*pi/lambda; %波矢A = 1; %振幅w = 3; %高斯光束的束宽[x1,y1] = meshgrid(linspace(-10,10,N...
2020-05-04 09:50:52
23947
25
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人