✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
更多Matlab仿真内容点击👇
⛄ 内容介绍
由于涡旋光束可以在不增加光谱带宽的前提下显著地提升信道的信息容量,而且能够更好地缓解湍流对光束传输时的负面影响.
⛄ 部分代码
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Gauss-Laguerre Acoustic Beam
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear; clc; close all;
% Model Parameters
l =-1; % topological charge;
n = abs(l)+0; % 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.