基于Matlab实现LG涡旋光束

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

⛄ 内容介绍

由于涡旋光束可以在不增加光谱带宽的前提下显著地提升信道的信息容量,而且能够更好地缓解湍流对光束传输时的负面影响.

⛄ 部分代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%

%                       Gauss-Laguerre Acoustic Beam

%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear; clc; close all;

% Model Parameters

 l =-1;         % topological charge;

% radial index; n=|l|,|l|+2,|l|+4 ...

 D = sqrt(2);   % is a constant for normalization;

% Discrete domain

xc=-6.2:0.2:6.2; %[-]

yc=-6.2:0.2:6.2; %[-]

Z=0.6; %[-] a XY-slice in the z-direction

[X,Y] = meshgrid(xc,yc);

[TH,R] = cart2pol(X,Y);

% Analytical functions

G = @(r,z) D./sqrt(1+z.^2).*exp(-r.^2./(1+z.^2)).*exp(-1i/4*(z.*r.^2)./(1+z.^2));

A = @(r,z) (sqrt(2)*r./sqrt(1+z.^2)).^abs(l).*LaguerreL((n-abs(l))/2,abs(l),2*r.^2./(1+z.^2));

PHI = @(th) exp(1i*l*th);

PSI = @(z) exp(-1i*(n+1)*atan(z));

P = @(th,r,z,t) G(r,z).*A(r,z).*PHI(th).*PSI(z).*exp(-1i*t);

% Compute profile for a seleted time 't':

p1=P(TH,R,Z,0);

% Plot a single slice of the presure profile

figure(1); fontsize=12;

set(gcf,'position',[100,100,600,200])

subplot(121), imagesc(xc,yc,flipud(abs(p1))); colorbar;

xlabel('x','interpreter','latex','fontsize',fontsize);

ylabel('y','interpreter','latex','fontsize',fontsize);

title('magnitude','interpreter','latex','fontsize',fontsize);

subplot(122), imagesc(xc,yc,flipud(angle(p1))); colorbar;

xlabel('x','interpreter','latex','fontsize',fontsize);

ylabel('y','interpreter','latex','fontsize',fontsize);

title('phase at $t_0$','interpreter','latex','fontsize',fontsize);

⛄ 运行结果

⛄ 参考文献

[1]朝鲁门. 涡旋光的产生与传播特性研究[D]. 内蒙古大学, 2014.

⛄ 完整代码

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料