光波在介质中界面上的反射及透射特性的仿真
光课程设计——光波在介质中界面上的反射及透射特性的仿真教学提纲
光课程设计——光波在介质中界面上的反射及透射特性的仿真西安邮电大学光学报告学院：电子工程学生姓名：专业名称：光信息科学与技术班级：光信1103班设计名称：光波在介质中界面上的反射及透射特性的仿真一、课程设计目的1．掌握反射系数及透射系数的概念；2．掌握反射光与透射光振幅和相位的变化规律；3．掌握布儒斯特角和全反射临界角的概念。
二、任务与要求对n1=1、n2=1.52及n1=1.52、n2=1的两种情况下，分别计算反射光与透射光振幅和相位的变化，绘出变化曲线并总结规律。
三、课程设计原理光在介质界面上的反射和折射特性与电矢量的振动方向密切相关。
由于平面光波的横波特性，电矢量可在垂直传播方向的平面内的任意方向上振动，而它总可以p s m E Et E E r imtm m im rm m ,,,0000===分解成垂直于入射面振动的分量和平行于入射面振动的分量，一旦这两个分量的反射、折射特性确定，则任意方向上的振动的光的反射、折射特性也即确定。
菲涅耳公式就是确定这两个振动分量反射、折射特性的定量关系式。
（1）s 分量和p 分量垂直入射面的振动分量- -s 分量平行入射面的振动分量- -p 分量定义：s 分量、p 分量的反射系数、透射系数分别为（2）反射系数和透射系数定义：s 分量、p 分量的反射系数、透射系数分别为p s m E Et E E r imtm m im rm m ,,,0000===（3）菲涅耳公式已知界面两侧的折射率21n n 、和入射角1θ，就可由折射定律确定折射角2θ；进而可由菲涅耳公式求出反射系数和透射系数。
绘出如下按光学玻璃（n=1.5）和空气界面计算，在21n n <（光由光疏介质射向光密介质）和21n n >（光由光密介质射向光疏介质）两种情况下，反射系数、透射系数随入射角1θ的变化曲线。
(a)光由光疏介质射向光密介质 (b)光由光密介质射向光疏介反射光与入射光中s,p 分量的相位关系： （1）n1＜n2时，光疏入射光密 s 分量的反射系数s r ：反射光中的s 分量与入射光中的s 分量相位相反；反射光中的s 分量相对入射光中的s 分量存在一个π相位突变(rs ϕ=π)； p 分量的反射系数p r ：在1θ<B θ范围内,p r ＞0，反射光中的p 分p 量与入射光中的分量相位相同(rp ϕ=0)；在1θ>B θ范围内,p r <0，反射光中的p 分量相对入射光中的p 分量有π相位突变(rp ϕ=π)；（2）n1>n2时，光密入射光疏 s 分量的反射系数s r ：入射角1θ在0到C θ(临界角，12/sin n n C =θ）的范围内，s 分量的反射系数s r ＞0。
基础光学第4章光波在界面的反射和折射规律课件
无论 n1 n2 或 n1 n2
透射光1’和2’振动方
向相同。即无半波损失。
只要光线2存在，光线1
和2的振动方向总是相
反的，即1和2的光程之
间存在半波损失。
光在多层透明介质界面的反射和折射
n1 < n2< n3 或n3 < n2< n1
时，光线1和2之间的光程
没有半波损失。
当折射率不按顺序排列时，
p
s
n2 n1
t p ts
2n1
n2 n1
入射
反射
约定
n1 < n2 n1 > n2
rp
rs
tp
ts
+
+
+
+
+
约定
实际
实际
+
反射
入射
n1
n2
约定
实际
透射
(a) n1 < n2
n1 < n2 时反射光与入射光振动方向相反
n1 > n2 时反射光与入射光振动方向相同
在任何情况下，透射光的方向和入射光相同
中的多次反射，分别求光从空气（折射率为1.0）正入射到玻
璃上表面，以及光从玻璃下表面射出时的振幅反射率、光强
反射率、振幅透射率和光强透射率。
【解】 正入射： i1 i2 0
n n
2n1
rp 2 1 rs ， t p ts
n2 n1
n2 n1
2
n2 n1
在反射和折射过程中，p, s两个分量的振动是相互独立的
4.2 菲涅尔反射和折射公式
n cos i1 n1 cos i2
光波场在2种介质界面反射时的半波损失
光波场在2种介质界面反射时的半波损失
庹有康;陈希明
【期刊名称】《重庆邮电大学学报（自然科学版）》
【年(卷),期】2004(016)006
【摘要】在波动光学中常常会遇到光波场在2种介质分界面上反射时的半波损失问题,而大多数教科书并未能从理论上阐明半波损失产生的条件和原因.从电磁场的基本理论出发,分析和讨论了光波场在两种界面上反射时的半波损失问题,从本质上全面、合理地解释了在实验中产生半波损失的实验结果.
【总页数】3页(P36-38)
【作者】庹有康;陈希明
【作者单位】重庆邮电学院,重庆,400065;重庆邮电学院,重庆,400065
【正文语种】中文
【中图分类】O431.1
【相关文献】
1.光从晶体到各向同性介质界面的反射和透射 [J], 裴芳芳;陈西园
2.论光在介质表面反射时的半波损失问题--兼与张静江先生商榷 [J], 陶志伟;刘涵哲
3.半波损失——平面偏振光在界面反射时的相位变化 [J], 龙海燕
4.半波损失——平面偏振光在界面反射时的相位变化 [J], 龙海燕;
5.光波入射到光疏—光密界面时反射光与折射光电矢量分析 [J], 陶志伟;高飞
因版权原因，仅展示原文概要，查看原文内容请购买。
光波在介质中界面的反射及透射的特性仿真实验题
1. 如何确定入射面?答：入射光与反射光以及法线共同构成的平面即入射面2.什么是临界角?临界角是光疏到光密，还是光密到光疏时发生?答：临界角就是全反射角，他指的是光线由光密介质入射到光疏介质时正好发生全反射时的入射角。
3.利用全反射现象能否产生圆偏振光?答；利用全反射现象可以产生圆偏振光，一个偏振光在一定角度上经过两次全反射可以产生圆偏振光，菲涅耳棱镜就是利用这个原理所制成的。
4.解释反射系数及透射系数的概念。
答：当电磁波由一个磁导率为μ1、介电常数为ε1的均匀介质，进入另一个具有磁导率为μ2、介电常数为ε2的均匀介质时，一部分电磁波在界面上被反射回来，另一分电磁波则透射过去。
反射波与透射波的振幅同入射波振幅之比，分别称之为反射系数与透射系数。
5.根据仿真曲线解释反射及透射光的相位变化规律。
答：图中反应了他们的相位的变化规律，例如图三所示在布儒斯特角处它的相位发生了π的跃变，而根据一个确定的波的表达式来看它是由余弦函数的的变化来确定的，而rp在菲涅耳表达式中是两个确定的余弦函数之比，所以rp由正变为负的时候，其中有一个余弦函数肯定相位发生了变化（奇变偶不变，符号看象限），且在布儒斯特角处，而在全反射角处也会发生变化，而且是逐渐变化的，这是因为当入射角逐渐增大的时候，它满足一个公式tan(fai/2)=-√((sin θ)^2-n^2)/cosθ),从公式可以看出相位会随着入射角的变化而渐变，当θ=π/2时，tan(fai/2)为无穷，所以fai=π。
6.试说明布儒斯特角的概念。
答：布儒斯特角，又称偏振角，是自然光经电介质界面反射后，反射光为线偏振光所应满足的条件。
7.试分析布儒斯特角与临界角哪个大。
答：临界角大于布儒斯特角，我们从它们的公式可以简单的推导出来，布儒斯特角为arctan(n2/n1),全反射角为arcsin(n2/n2), 假设n2/n1=x,因为有光密入射到光疏，所以n2>n1，因此x>1，此时布儒斯特角为arctan(x),全反射角为arcsin(x),我们对它两个同时求导得到：（arctan(x)）’=1/(1+x^2),而（arcsin(x)）’=1/√(1+x^2),由此我们可以得出全反射角公式的倒数大，也就是说，在相同变量的情况下它的数值大，从而我们也就说明了临界角大于布儒斯特角。
第二章-理想平板介质光波导中的光传播特性及仿真
第2章介质光波导分析方法2.1 平板介质光波导一般概念2.1 平板介质光波导一般概念波动理论法则是把平板介质光波导中的光波看作是满足波导边界条件的麦克斯韦方程组的解。
2.2 平板光波导分析的射线法振幅反射率和附加相移振幅反射率和附加相移S 波（TE 波——电矢量平行于界面）振幅反射率：光传播过程相位变化：光波不仅在介质中传播过程中相位会发生改变，在界面上反射时相位也会变化。
对于θ1 ＜ θ1c ，界面上发生全反射，此时上式的分子和分母中第二个平方根内为负数，因此得到的振幅反射率r 为复数。
1.106分子分母同乘k ）振幅反射率和附加相移振幅反射率和附加相移s 波（ TE波——电矢量平行于界面） 附加相移为：p 波（ TM波——磁矢量平行于界面）在界面发生全反射时引起的附加相移为：（1.145）（1.144）界面：n1、n2、n3的界面，不是入射面平板波导中的其他光场均可视为TEM 模：模式只有横向分量，而无纵导模特征方程导模特征方程入射光线两次反射后与入射光线同方向传输特征方程特征方程A、B 两点的距离为:C、D两点的距离为:光线CD 还经历了两个附加相移:分别是介质1、3 界面处全反射的附加相移ϕ3 和介质1、2 界面处全反射的附加相移ϕ2.平板光波导的特征方程：特征方程特征方程 界面处的附加相移会因入射光偏振方向的不同而有所差异，因此就能够得到两个不同模式下的特征方程电矢量平行于界面的导波式中：特征方程特征方程同样地，磁矢量平行于界面的导波TM 模的特征方程（代入ΦM2和ΦM3） ：这里采用的是简单光线传播的射线理论。
实际上，从麦克斯韦方程出发，结合介质界面处的边界条件也可以推导出以上特征方程。
引入的几个重要参数——都是θ的函数,得到一个光波模式的波矢就可以求解其他引入的几个重要参数——都是θ的函数,得到一个光波模式的波矢就可以求解其他纵向波矢横向波矢衰减系数< n 1kn 2k << n 1k2.3 平板光波导中的TE模TE模的电磁理论求解TE模的电磁理论求解平板光波导中的TE模仅有E y由麦克斯韦方程：(2.30)TE模的电磁理论求解TE模的电磁理论求解的式子因此可以将H的分量表示为Ey代入式(2.30)，可以得到关于Ey的波动方程，j = 1;2;3 表示分别是在芯层、衬底和覆盖层。
光波在介质界面上的反射和折射 PPT课件
(2) 大小
i t / 2 n1 sini n2 sint
tan B
n2 n1
n21
(3) 应用
3、全反射
•
设光波从光密介质射向光疏介质(n1>n2)，
折射角θ2大于入射角θ1。当sinθ1＝n2/n1时，θ2
为90o，这时折射角沿界面掠过。若入射角再增
大，使sinθ1>n2/n1 ，这时不能定义实数的折射 角。使θ2=90o的入射角θ1称为临界角，记作θc 即
E0ip cosi E0rp cosr＝Et0 p cost
2、反射系数和透射系数
rp
E0rp E0ip
n2 cos1 n1 cos2 n2 cos1 n1 cos2
tan(1 2 ) tan(1 2 )
tp
E0tp E0ip
2n1 cos1 n2 cos1 n1 cos2
2 cos1 sin 2 sin(1 2) cos(1 2)
RT 1
四、反射率和透射率
3、反射率随入射角变化关系
R随入射角θ1的变化关系
11日出生于苏格兰杰德伯勒，1800年毕业于爱
丁堡大学，曾任“爱丁堡杂志”、“苏格兰杂 志”、“爱丁堡百科全书”编辑，爱丁堡大学
教授、校长等。1815年被选为皇家学会会员， 1819年获冉福德奖章。
•
布儒斯特主要从事光学方面的研究。1812
年发现当入射角的正切等于媒质的相对折射率 时，反射光线将为线偏振光（现称为布儒斯特
Ets Eis Ers
n H2 H1 0
n Htp cost Hip cosi Hrp cosr 0
Hip cosi H rp cosr Htp cost
7
光波在各向同性介质界面的反射和折射 ppt课件
ppt课件
17
(2)大角度入射(掠射)的反射特性
由图1-24(a)，有
n1＜n2，光疏到光密。θ 1≈900的掠射情况。
rs 0, rp 0
在入射点处，反射光矢量Er与入射光矢量Ei方向近似相 反，将产生半波损失。 n1>n2，光密到光疏。掠射θ 1≈900>θ c。全反射。 在入射点处，反射光产生半波损失的条件：
ki sin i kr sin r , ki sin i kt sin t n1 sin i n1 sin r , n1 sin i n2 sin t
反射定律
T 1-21
折射定律
描述光在介质面上的传播方向
ppt课件 3
1.2.2 菲涅耳公式
描述入射光、反射光和折射光 之间的振幅、相位关系。 1.s分量和p分量 垂直入射面的振动分量- -s分量 T 1-23 平行入射面的振动分量- -p分量 规定分量和分量的正方向如图所示 2.反射系数和透射系数 定义：s分量、p分量的反射系数、透射系数分别为
① n1＜n2，光疏到光密。先考察θ 1=00的正入射情况。 由图1-24(a)，有
rs 0, rp 0
考虑P30 T1-23，有关光场振动正方向的规定，则有
可见：在入射点处，合成的反射光矢量Er相对入射光场Ei反 向，相位发生π突变，或半波损失。 对于θ 1非零、小角度入射时，都将近似产生π相位突变，或 半波损失。
入射光中s分量和p分量的透射率(不相同)为
n2 cos 2 2 sin 21 sin 2 2 Ts ts n1 cos1 sin 2 (1 2 )
n2 cos 2 2 sin 21 sin 2 2 Tp tp 2 n1 cos1 sin (1 2 ) cos2 (1 2 )
平面光波在电介质表面的反射特性仿真研究
１ Байду номын сангаас— １
《 电子设 计 工程 ） ０ ２年 第 １ ２１ １期
透射波 ： Ｅ＝ ，ｘ ｋ（ｉ —ｏ ）【 ， ｏ ｐｉ ２ｓ ｃｓ Ｅｅ ［ ｎ 一ｏ 胡 （） ４
１２ ８ ３年 ） 。其 中 ， （ １ 和 式 （ ３ 是 反 射 公 式 ， （ ２ 和 式 式 １） １） 式 １）
Ｅ ｅｐｉ １ｓ ０ — Ｏ０ｙ一 ｔ ； ｘ ｋ（ｉ ｒ ＣＳｒ）ｔ】 ［ ｎｘ Ｏ
（） ３
基 金 项 目： 军队 重 点科 研 项 目（ Ｊ ５ ３ ） Ｋ ０ １８ 作 者 简 介 ： 昊 鹏 （９ ７ ） 男 ， 宁 沈 阳人 ， 士 研 究 生 。研 究 方 向 ： 确 制 导 武 器 的作 战 使 用与 仿 真 。 王 １８ 一 ， 辽 硕 精
ｃ ｓ 矿ｃ ｓ ｏ ０悱 ｃｓ ０ ｏ （ ６）
失问题 。
２ 平面 光 波在 电 介质 表 面 的 反射 和 折 射
２１ 电矢 量 平行 入 射 面 ．
平 面 光 波 的 电 矢 量 平 行 于入 射 面 ， 此 其 电场 只有 Ｐ分 因 量 ， 磁 场 垂 直 于 入 射 面 ， 此 只 有 Ｓ分 量 ， 面 光 波 传 输 其 因 平 方 向矢 量 ｜ 在 入 射 面 内 ， 与 ： 平 行 。 以 、 和 曰所 确定 ｊ ｝ 轴 Ｅ
ｃｓｉ ｅｐｉ ｉｉ０－ ｃｓ Ｅ０ ｐｉ ｉｉ０－ｔ ｏＯ ￣ｘ ｋｓ ｎｒ Ｉ一 ｏ ，ｘ ｋ ｓ ｃ 】 Ｅ ［ ｓ （ ） ｅ ［ｓ ｎ ｏ ＝ ｏ０ ￣ｘ ｋｓｉＯ－ ｔ ｃｓ， ｅｐｉ ２ ｎｒｏ］ Ｅ ［ ｓ ｔ 再 结 合 ｊ＝ ￣ ｌ 和 ｋ＝ ￣ Ｊ 代 人 （ ） 简 可 得 ： ｝ ２ｒ １ ｎ ２２ ｎ Ａ， ５化
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
西安邮电大学光学报告学院：电子工程学生姓名：专业名称：光信息科学与技术班级：光信1103班设计名称：光波在介质中界面上的反射及透射特性的仿真一、课程设计目的1．掌握反射系数及透射系数的概念；2．掌握反射光与透射光振幅和相位的变化规律；3．掌握布儒斯特角和全反射临界角的概念。
二、任务与要求对n1=1、n2=1.52及n1=1.52、n2=1的两种情况下，分别计算反射光与透射光振幅和相位的变化，绘出变化曲线并总结规律。
三、课程设计原理光在介质界面上的反射和折射特性与电矢量的振动方向密切相关。
由于平面光波的横波特性，电矢量可在垂直传播方向的平面内的任意方向上振动，而它总可以分解成垂直于入射面振动的分量和平行于入射面振动的分量，一旦这两个分量的反射、折射特性确定，则任意方向上的振动的光的反射、折射特性也即确定。
菲涅耳公式就是确定这两个振动分量反射、折射特性的定量关系式。
（1）s分量和p分量p s m E Et E E r imtm m im rm m ,,,0000===垂直入射面的振动分量- -s 分量平行入射面的振动分量- -p 分量定义：s 分量、p 分量的反射系数、透射系数分别为（2）反射系数和透射系数定义：s 分量、p 分量的反射系数、透射系数分别为p s m E E t E E r imtm m im rm m ,,,0000===（3）菲涅耳公式已知界面两侧的折射率21n n 、和入射角1θ，就可由折射定律确定折射角2θ；进而可由菲涅耳公式求出反射系数和透射系数。
绘出如下按光学玻璃（n=1.5）和空气界面计算，在21n n <（光由光疏介质射向光密介质）和21n n >（光由光密介质射向光疏介质）两种情况下，反射系数、透射系数随入射角1θ的变化曲线。
(a)光由光疏介质射向光密介质 (b)光由光密介质射向光疏介反射光与入射光中s,p 分量的相位关系： （1）n1＜n2时，光疏入射光密 s 分量的反射系数s r ：反射光中的s 分量与入射光中的s 分量相位相反；反射光中的s 分量相对入射光中的s 分量存在一个π相位突变(rs ϕ=π)； p 分量的反射系数p r ：在1θ<B θ范围内,p r ＞0，反射光中的p 分p 量与入射光中的分量相位相同(rp ϕ=0)；在1θ>B θ范围内,p r <0，反射光中的p 分量相对入射光中的p 分量有π相位突变(rp ϕ=π)；（2）n1>n2时，光密入射光疏 s 分量的反射系数s r ：入射角1θ在0到C θ(临界角，12/sin n n C =θ）的范围内，s 分量的反射系数s r ＞0。
反射光中的s 分量与入射光中的s 分量同相位，rs ϕ=0；入射角1θ>C θ时，发生全反射，1212cos sin 2tan θθϕn rs--=（21/n n n =）； p 分量的反射系数p r ：在1θ＜B θ范围内，p r <0，反射光中的p 分量相对入射光中的p 分量有π相位突变（rp ϕ=π)；在B θ<1θ<C θ范围内，p r >0，反射光中的p 分量与入射光中的p 分量相位相同(rp ϕ=0)；入射角1θ>C θ时，发生全反射，12122cos sin 2tan θθϕn n rp--=； 四、课程设计步骤（流程图）（1）定义变量n1,n2,f1.（2）给变量赋值,其中n1=1,n2=1.52,还有一种情况其中n1=1.52,n2=1（3）设计for 循环，使f1每循环一次加 /1000，实现在f1每变化一次下，得出相应的反射系数，透射系数的值，从而得出程序的循环（4）根据程序仿真结果五、仿真结果分析102030405060708090-1-0.500.51Fn1<n2s/p 分量与相位的关系010203040506070809001234Ff r s010203040506070809001234Ff r p0102030405060708090-11234n1>n2s/p 分量与相位的关系F01020304050607080901234Ff r s01020304050607080901234Ff r p结论：光在介质面上的反射、透射特性有三个因素决定：入射光的偏振态，入射角，界面两侧介质的折射率。
(1)光波由光疏介质射向光密介质（n1<n2）a.n1<n2时，反射系数rs<0,说明反射光中的s 分量与入射光中的s 分量相位相反。
(即frs=π)b.而p分量的反射系数rp在f1<fb范围内，rp>0,说明反射光中的p 分量与入射光中的p分量相位相同。
（即frp=0）c.在f1>fb范围内，rp<0,说明反射光中的p分量与入射光中的p分量π相位突变。
（即frp=π）(2)光波由光密介质射向光疏介质（n1>n2）a.入射角f1在0到fc的范围内，s分量的反射系数rs>0,说明反射光中s分量与入射光中的s分量同相位。
（即frs=0）b.P分量的反射系数rp在f1<fb范围内，rp<0,说明反射光中的p分量相对入射光中的p分量有π相位突变。
（即frp=π）c.在fb<f1<fc范围内，rp>0,说明反射光中的p分量与入射光中的p分量相位相同。
六、仿真小结光在介质界面上的反射、透射特性由三个因素决定：（1）入射光的偏振态；（2）入射角;（3）界面两侧介质的折射率。
由rs、rp、ts、tp随入射角的变化曲线可知，在入射角从0度到90度的变化范围内，不论光波以什么角度入射至界面，也不论界面两侧折射率大小如何，s分量和p分量的透射系数t总是取正值，因此，折射光总是与入射光同相位。
通过本次实验，掌握了反射系数及透射系数的概念，反射光与透射光振幅和相位的变化规律，布儒斯特角和全反射临界角的概念。
七、程序clear all;%n1=1;%n2=1.52;n1=1.52;n2=1;n=n2./n1;if n1<n2subplot(1,3,1)qa=0:pi/100:pi/2;qb=asin(n1.*sin(qa)./n2);rs=-sin(qa-qb)./sin(qa+qb);rp=tan(qa-qb)./tan(qa+qb);ts=2.*cos(qa).*sin(qb)./sin(qa+qb);tp=2.*cos(qa).*sin(qb)./sin(qa+qb)./cos(qa-qb);plot(qa*180./pi,rs,'r',qa*180./pi,rp,'c',qa*180./pi,ts,'b',qa*180./pi ,tp,'g')legend('rs','rp','ts','tp')%rssubplot(1,3,2)for qa=0:pi/1000:pi/2qb=asin(n1.*sin(qa)./n2);rs=-sin(qa-qb)./sin(qa+qb);if rs<=0Frs=pi;elseFrs=0;endplot(qa*180./pi,Frs,'r')hold onendlegend('Frs')%rpsubplot(1,3,3)for qa=0:pi/1000:pi/2qb=asin(n1.*sin(qa)./n2);rp=tan(qa-qb)./tan(qa+qb);if rp<=0Frp=pi;elseFrp=0;endplot(qa*180./pi,Frp,'b')hold onendlegend('Frp')elsesubplot(1,3,1)qc=asin(n2./n1);qa=0:0.0001:qc;qb=asin(n1.*sin(qa)./n2);rs=-sin(qa-qb)./sin(qa+qb);rp=tan(qa-qb)./tan(qa+qb);ts=2.*cos(qa).*sin(qb)./sin(qa+qb);tp=2.*cos(qa).*sin(qb)./sin(qa+qb)./cos(qa-qb);plot(qa*180./pi,rs,'r',qa*180./pi,rp,'c',qa*180./pi,ts,'b',qa*180./pi ,tp,'g')hold onqa=qc:0.0001:pi/2;tp=0;ts=0;rs=1;rp=1;plot(qa*180./pi,rs,'r',qa*180./pi,rp,'c',qa*180./pi,ts,'b',qa*180./pi ,tp,'g')hold onlegend('rs','rp','ts','tp')%rsqc=asin(n2./n1);subplot(1,3,2)for qa=0:pi/1000:qcqb=asin(n1.*sin(qa)./n2);rs=-sin(qa-qb)./sin(qa+qb);if rs<=0Frs=pi;elseFrs=0;endplot(qa*180./pi,Frs,'r')hold onendqa=qc:pi/1000:pi/2;Frs= 2.*atan(sqrt(sin(qa).^2-(n.^2))./cos(qa));plot(qa*180./pi,Frs,'r')hold onlegend('Frs')%rpsubplot(1,3,3)for qa=0:pi/1000:qc;qb=asin(n1.*sin(qa)./n2);rp=tan(qa-qb)./tan(qa+qb);if rp<=0Frp=pi;elseFrp=0;endplot(qa*180./pi,Frp,'b')hold onendqa=qc:pi/1000:pi/2;Frp= 2.*atan(sqrt(sin(qa).^2-(n.^2))./cos(qa)./n.^2); plot(qa*180./pi,Frp,'b')hold onlegend('Frp')end。