⛄一、获取代码方式
获取代码方式1:
完整代码已上传我的资源:【光学】基于matlab光波叠加【含Matlab源码 2371期】
点击上面蓝色字体,直接付费下载,即可。
获取代码方式2:
付费专栏Matlab物理应用(初级版)
备注:
点击上面蓝色字体付费专栏Matlab物理应用(初级版),扫描上面二维码,付费29.9元订阅海神之光博客付费专栏Matlab物理应用(初级版),凭支付凭证,私信博主,可免费获得1份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:1份本博客上传CSDN资源代码
⛄二、部分源代码
%% 光波在几种特殊条件下的叠加仿真
%
%% 驻波:两个频率、振幅大小及方向相同而传播方向相反的单色光波叠加
%设叠加的两个波表达式为:E1=acos(kz+wt), E2=acos(kz-wt+delte)
format short
lambda = 7e-4; %入射光波长(mm)
delte = pi/3; %两个波的相位差弧度制
a = 3; %两个波的振幅
k = 2pi/lambda; %波数
z = linspace(-10e-4, 10e-4, 1000); %波的位移分量
figure(‘NumberTitle’, ‘off’, ‘Name’, ‘驻波仿真’);
plot(z, 2acos(kz+delte/2), ‘r’, ‘Linewidth’, 1.2);
grid on
xstr={[‘位移(mm)’]; [‘入射光波长:’ num2str(lambda*1e6) ‘nm’ ’ 相位差:’ num2str(delte) ’ 振幅:’ num2str(a)]};
xlabel(xstr, ‘Fontsize’, 11)
ylabel(‘振幅(负数代表相位相反)’, ‘Fontsize’, 11)
title(‘驻波(Stationary Wave)仿真示意图’,‘FontSize’,14);
%% 椭圆偏振光:两个频率相同且正交的单色光波叠加
%设两个叠加的光波表达式为 Ex=a1cos(kz-wt) Ey=a2cos(kz-wt+delte)
delte = pi/3; %两个波的相位差弧度制
a1 = 4;
a2 = 4; %两个波的振幅
figure(‘NumberTitle’, ‘off’, ‘Name’, ‘圆偏振光合矢量末端运动轨迹’);
thelte = linspace(0, 2pi, 500);
plot(a1cos(thelte), a2cos(thelte+delte),‘g’, ‘Linewidth’,1.2);
title(‘椭圆偏振光(Elliptically Polarized Light)仿真示意图’,‘FontSize’,14);
phi = 0.5atand(2a1a2/(a12-a22)*cos(delte)); %椭圆长轴与x轴夹角(°)
xstr = {[‘相位差:’ num2str(delte) ’ 振幅:’ num2str(a1) ‘,’ num2str(a2)] ;['椭圆长轴与x轴夹角{\phi} = ’ num2str(phi) ‘°’]};
xlabel(xstr, ‘Fontsize’, 11);
axis equal
grid on
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 门云阁.MATLAB物理计算与可视化[M].清华大学出版社,2013.
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除