【光学】基于matlab涡旋光与球面波的干涉【含Matlab源码 597期】

本文介绍了如何使用matlab模拟涡旋光与球面波的干涉现象,提供了完整的源代码和运行结果。通过设置像素点、束腰半径和拓扑荷值,生成并显示干涉图像。代码适用于matlab 2014a版本,参考了相关书籍。
摘要由CSDN通过智能技术生成

⛄一、获取代码方式

获取代码方式1:
完整代码已上传我的资源:【光学】基于matlab涡旋光与球面波的干涉【含Matlab源码 597期】
点击上面蓝色字体,直接付费下载,即可。

获取代码方式2:
付费专栏Matlab物理应用(初级版)

备注:
点击上面蓝色字体付费专栏Matlab物理应用(初级版),扫描上面二维码,付费29.9元订阅海神之光博客付费专栏Matlab物理应用(初级版),凭支付凭证,私信博主,可免费获得1份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:1份本博客上传CSDN资源代码

⛄二、部分源代码

clear
clc
row=1080; %像素点设置成1080pix
g1=zeros(row); %生成1080阶零矩阵&#

  • 4
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
涡旋与平面波的干涉可以通过Matlab编写实现。下面是一个简单的代码实现: ```matlab % 设置参数 N = 200; % 图像大小 lambda = 632.8e-9; % 波长 k = 2*pi/lambda; % 波数 z = 0; % 传播距离 w0 = N/8; % 束半径 m = 5; % 涡旋的模式数 theta0 = pi/4; % 平面波的入射角度 % 生成网格 x = linspace(-1,1,N); y = linspace(-1,1,N); [X,Y] = meshgrid(x,y); % 生成高斯束 wz = w0*sqrt(1+(lambda*z/(pi*w0^2))^2); Rz = z*(1+(pi*w0^2)/(lambda*z)^2); w = w0*sqrt(1+(z/(pi*w0^2))^2); psi = atan(z/(pi*w0^2)); u = sqrt(2/pi)*wz/w0*exp(-(X.^2+Y.^2)/w^2).*exp(-1i*k*z).*exp(1i*psi); % 生成涡旋 theta = atan2(Y,X); phi = m*theta; v = exp(1i*phi).*u; % 生成平面波 p = exp(1i*k*(X*cos(theta0)+Y*sin(theta0))); % 干涉 w = u + v + p; % 显示图像 figure; subplot(2,2,1); imagesc(abs(u).^2); axis equal tight; colormap(jet); title('高斯束'); xlabel('x'); ylabel('y'); subplot(2,2,2); imagesc(abs(v).^2); axis equal tight; colormap(jet); title('涡旋'); xlabel('x'); ylabel('y'); subplot(2,2,3); imagesc(abs(p).^2); axis equal tight; colormap(jet); title('平面波'); xlabel('x'); ylabel('y'); subplot(2,2,4); imagesc(abs(w).^2); axis equal tight; colormap(jet); title('干涉图像'); xlabel('x'); ylabel('y'); ``` 这个代码生成一个大小为200x200的图像,其中包一个高斯束、一个模式数为5的涡旋和一个入射角度为45度的平面波。它们在相干条件下干涉,生成干涉图像。你可以调整参数来生成不同的干涉图像。
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

海神之光

有机会获得赠送范围1份代码

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值