没有合适的资源?快使用搜索试试~ 我知道了~
首页MATLAB 光学仿真入门.pdf
MATLAB 光学仿真入门.pdf
5星 · 超过95%的资源 需积分: 23 146 下载量 172 浏览量
更新于2023-04-28
评论 19
收藏 4.55MB PDF 举报
matlab光学仿真入门,里面包含了光波通过不同介质时反射率和透射率的变化、高斯光束性质相关的仿真。
资源详情
资源评论
资源推荐
MATLAB 光学仿真入门
光的电磁理论基础
光波由光疏介质进入光密介质
分析光波由光疏介质进入光密介质时反射率和透射率的变化。光疏介质空气 ,光密介质石英玻
璃 ,作出 、 分量的振幅反射率和振幅透射率以及他们的绝对值随入射角度的变化曲
线。
n1=1, n2=1.45;
theta=0:0.1:90;
a=theta*pi/180;
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-
(n1/n2*sin(a)).^2));
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-
(n1/n2*sin(a)).^2));
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
figure(1);
subplot(1,2,1);
plot(theta,rp,'-',theta,rs,'--
',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
legend('r_p','r_s','|r_p|','|r_s|')
xlabel(' heta_i')
ylabel('Amplitude')
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
axis([0 90 -1 1])
grid on
subplot(1,2,2);
plot(theta,tp,'-',theta,ts,'--
',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
legend('t_p','t_s','|t_p|','|t_s|')
xlabel(' heta_i')
ylabel('Amplitude')
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
axis([0 90 0 1])
grid on
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
由图可知:
1. 当入射角 ,即垂直入射时, 、 和 、 都不为 ,表示存在反射波和折射波。
2. 当入射角 ,即掠入射时, , ,即没有折射光波。
3. 、 随 的增大而减小, 随 的增大而增大,直到等于 。
4. 先随 的增大而减小,到达一特定的值 时,有 ,即反射波中此时没有 分量,
只有 分量,产生全偏振现象,然后随着 的增大, 不断增大,直到等于 。
作出 、 分量的能流反射率和能流透射率以及他们的平均值随入射角度的变化曲线。
n1=1, n2=1.45;
theta=0:0.1:90;
a=theta*pi/180;
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-
(n1/n2*sin(a)).^2));
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-
(n1/n2*sin(a)).^2));
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Rp=abs(rp).^2;
Rs=abs(rs).^2;
Rn=(Rp+Rs)/2;
Tp=n2*sqrt(1-(n1/n2*sin(a)).^2)./(n1*cos(a)).*abs(tp).^2;
Ts=n2*sqrt(1-(n1/n2*sin(a)).^2)./(n1*cos(a)).*abs(ts).^2;
Tn=(Tp+Ts)/2;
figure(1);
subplot(1,2,1);
plot(theta,Rp,'-',theta,Rs,'-.',theta,Rn,'--','LineWidth',2)
legend('R_p','R_s','R_n')
xlabel('\theta_i')
ylabel('Amplitude')
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
axis([0 90 0 1])
grid on
subplot(1,2,2);
plot(theta,Tp,'-',theta,Ts,'-.',theta,Tn,'--','LineWidth',2)
legend('T_p','T_s','T_n')
xlabel('\theta_i')
ylabel('Amplitude')
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
axis([0 90 0 1])
grid on
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
由图可知:
1. 当入射角 时,垂直入射时能流反射率 、 和 、 都不为 ,此时存在反射光
波。
2. 随着 的增大, 不断增大至 , 不断减小至 ,但始终有 。
3. 随着 的增大, 先减小,直至一特定的值 时变为 ,而后随着 的增大不断增大到 。
的过程正好相反,在入射角为 时为 ,且始终有 。
光波由光密介质进入光疏介质
分析光波由光密介质进入光疏介质时反射率和透射率的变化。光疏介质空气 ,光密介质石英玻
璃 ,作出 、 分量的振幅反射率和振幅透射率以及他们的绝对值随入射角度的变化曲线。
与上述过程相同,只需要将折射率互换。此处的光波变化分析略去。
在图中, 后, 和 以及 和 产生了很大的差异。因为 , 如果
则 ,计算得到的 和 将变为复数。但作图时只取了
实部。
[theta' rp' rs']
ans =
0.0000 + 0.0000i -0.1837 + 0.0000i 0.1837 + 0.0000i
0.1000 + 0.0000i -0.1837 + 0.0000i 0.1837 + 0.0000i
0.2000 + 0.0000i -0.1837 + 0.0000i 0.1837 + 0.0000i
0.3000 + 0.0000i -0.1837 + 0.0000i 0.1837 + 0.0000i
0.4000 + 0.0000i -0.1837 + 0.0000i 0.1837 + 0.0000i
0.5000 + 0.0000i -0.1837 + 0.0000i 0.1837 + 0.0000i
......
43.0000 + 0.0000i 0.5449 + 0.0000i 0.7542 + 0.0000i
43.1000 + 0.0000i 0.5754 + 0.0000i 0.7727 + 0.0000i
43.2000 + 0.0000i 0.6108 + 0.0000i 0.7938 + 0.0000i
43.3000 + 0.0000i 0.6531 + 0.0000i 0.8185 + 0.0000i
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
剩余12页未读,继续阅读
StriverLite
- 粉丝: 114
- 资源: 22
上传资源 快速赚钱
- 我的内容管理 收起
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
会员权益专享
最新资源
- zigbee-cluster-library-specification
- JSBSim Reference Manual
- c++校园超市商品信息管理系统课程设计说明书(含源代码) (2).pdf
- 建筑供配电系统相关课件.pptx
- 企业管理规章制度及管理模式.doc
- vb打开摄像头.doc
- 云计算-可信计算中认证协议改进方案.pdf
- [详细完整版]单片机编程4.ppt
- c语言常用算法.pdf
- c++经典程序代码大全.pdf
- 单片机数字时钟资料.doc
- 11项目管理前沿1.0.pptx
- 基于ssm的“魅力”繁峙宣传网站的设计与实现论文.doc
- 智慧交通综合解决方案.pptx
- 建筑防潮设计-PowerPointPresentati.pptx
- SPC统计过程控制程序.pptx
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功
评论5