2光波在介质中界面上的反射及透射特性的仿真
光波在介质界面上的反射和折射 菲涅耳公式
在讨论过程中，不计吸收、散射等能量损耗，因此， 入射光能量在反射光和折射光中重新分配，而总能 量保持不变。
2. 3 反射率和透射率 (Reflectivity and transmissivity) 若有一个平面光波以入射角1 斜入射介质分界面， 平面光波的强度为 Ii，则每秒入射到界面上单位面积 的能量为
sin (1 2 ) rs =sin (1 2 )
(134)
(Eis Ers )n1cos1 Ets n2 cos 2 (133)
3. 菲涅耳公式 利用类似方法，可以推出 p 分量的反射系数和透射系 数表示式， 这就是著名的菲涅耳公式：
sin(1 2 ) n1 cos 1 n2 cos 2 rs = sin(1 2 ) n1 cos 1 n2 cos 2 2 cos 1 sin 2 2n1 cos 1 ts = sin(1 2 ) n1 cos 1 n2 cos 2
Ers sin ( 2 1 ) Eis sin ( 2 1 )
(ki kr ) r 0 (121)
(134)
m s, p
（ 128）
E0 rm rm E0 im
(129)
3. 菲涅耳公式 由 （134）式和（133）式消去 Ers，经运算整理得
2n1 cos 1 ts =n1 cos 1 n2 cos 2 (135)
(145) (146)
n2 cos 2 2 sin 21 sin 2 2 Ts ts n1cos1 sin 2 (1 2 )
(147)
n2 cos 2 2 sin 21 sin 2 2 Tp tp 2 n1cos1 sin (1 2 ) cos 2 (1 2）
光波在分界面上的反射和折射
❖ 利用全反射时的能量特性，改变光的传播方向、传递能量。 ❖利用倏逝波特性产生的受抑全反射效应能制成光调制器或光输出耦 合器。 ❖ 利用全反射时的位相变化，选取适当的折射率和入射角，可改变 入射光的偏振状态。
(b)
玻璃菱体(n=1.51)，当 1=54037 (或48037) 时，有Δ =450。
因此，垂直菱体入射的线偏振光，若其振 动方向与入射面的法线成450角，则在菱体内上 下两个界而进行两次全反射后，s 分量和 p 分 量的相位差为900，因而输出光为圆偏振光。
See TIR from a fingerprint valley and FTIR from a ridge.
利用倏逝波特性产生的受抑全反射效应能制成光调制器或 光输出耦合器。
在斜面间的空气隙内的隐失波场的耦合作用下，光波可以从 一块棱镜透射到另一块棱镜，透射量的多少与间隔有关
棱镜波导耦合器：可以用来将光信号方便有效地耦合进薄膜波导中，或者将 在薄膜波导中传播的光信号引出波导。
利用全反射时的位相变化，选取适当的折射率和入射角，可改变 入射光的偏振状态。
k2xxt
沿x方向传播
等位相面（ 为x;： t） t0
k2xxk2xsint
c
倏逝波波长 ： 2 k2x sini
(2n)
光波的等相位面为一系列平行于z轴的等距直线
k1
’
④倏逝波等振幅面
等振 A 2幅 ex p z 面 ] [C ： zC 等相面 y
光波的等振幅面为一系列平行于x轴的等距直线
x 等 幅 面
透射波是一个沿x方向传播，振幅在z方向按指数规律衰减的复杂波； 其振幅在介质表面法向方向作指数迅速衰减的波，这就是倏逝波。 其等位相面为一系列平行于z轴的等距直线。 其等振幅面与等位相面垂直。
光课程设计——光波在介质中界面上的反射及透射特性的仿真教学提纲
光课程设计——光波在介质中界面上的反射及透射特性的仿真西安邮电大学光学报告学院：电子工程学生姓名：专业名称：光信息科学与技术班级：光信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。
西安邮电大学光学实验matlab仿真结果分析与程序
光学实验实验报告课程名称：光学实验*名：***学院：电子工程学院系部：光电子技术系专业：电子科学与技术年级：科技1201学号：********指导教师：**2014年12 月24 日光波在介质中界面上的反射及透射特性一．实验目的：1．掌握反射系数及透射系数的概念；2．掌握反射光与透射光振幅和相位的变化规律； 3．掌握布儒斯特角和全反射临界角的概念。
二．实验原理：1 反射定律和折射定律光由一种介质入射到另一种介质时，在界面上将产生反射和折射。
现假设二介质为均匀、透明、各向同性介质，分界面为无穷大的平面，入射、反射和折射光均为平面光波，其电场表示式为)(0r k t i l l l l e E E ⋅--=ω l =i, r, t式中，脚标i 、r 、t 分别代表入射光、反射光和折射光；r 是界面上任意点的矢径，在图2-1所示的坐标情况下，有r=ix+jy图2-1 平面光波在界面上的反射和折射 图2-2 k i 、k r 、k t 三波矢关系根据电磁场的边界条件，可以得到如下关系)(0)(t i r i tr i =⋅-=⋅-==r k k r k k ωωω 这些关系表明：①入射光、反射光和折射光具有相同的频率；②入射光、反射光和折射光均在入射面内，k i 、k r 和k t 波矢关系如图2-2所示。
进一步可得tt i i r r i i sin sin sin sin θθθθk k k k == 或tt i i r r i i sin sin sin sin θθθθn n n n ==即介质界面上的反射定律和折射定律，它们给出了反射光、折射光的方向。
折射定律又称为斯涅耳(Snell)定律。
2 菲涅耳公式 s 分量和p 分量通常把垂直于入射面振动的分量称做s 分量，把平行于入射面振动的分量称做p 分量。
为讨论方便起见，规定s 分量和p 分量的正方向如图2-3所示。
图2-3 s 分量和p 分量的正方向反射系数和透射系数 假设介质中的电场矢量为)(i 0e r k t l l l E E ⋅--=ω l =i, r, t其s 分量和p 分量表示式为)(i 0e r k t lm lm l E E ⋅--=ω m =s,p则定义s 分量、p 分量的反射系数、透射系数分别为tmtm m im rmm E E t E E r 0000==菲涅耳公式假设界面上的入射光、反射光和折射光同相位，根据电磁场的边界条件及s 分量、p 分量的正方向规定，可得ts rs s E E E i =+和2tp 1rp 1ip cos cos cos θθθH H H =-利用E H εμ=，上式变为22ts 11rs is cos cos )(θθn E n E E =-再利用折射定律，消去E ts ，经整理可得)sin()sin(1212is rs θθθθ+-=E E 根据反射系数定义，得到)sin()sin(2121θθθθ+--=s r221111cos cos cos 2θθθn n n t s +=将所得到的表示式写成一个方程组，就是著名的菲涅耳公式：212122112*********tan tan tan tan cos cos cos cos )sin()sin(θθθθθθθθθθθθ+--=+-=+--==n n n n E E r is rs s 2121211221122121002sin 2sin 2sin 2sin cos cos cos cos )tan()tan(θθθθθθθθθθθθ+-=+-=+-==n n n n E E r iprp p 21121121112100221111212100cos cos cos 2)cos()sin(sin cos 2cos cos cos 2)sin(sin cos 2θθθθθθθθθθθθθθθθn n n E E t n n n E E t iptp p is ts s +=-+==+=+==这些系数首先是由菲涅耳用弹性波理论得到的，所以又叫做菲涅耳系数。
光波在介质中界面的反射及透射的特性仿真实验题
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),由此我们可以得出全反射角公式的倒数大，也就是说，在相同变量的情况下它的数值大，从而我们也就说明了临界角大于布儒斯特角。
光波在介质界面上的反射和折射 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
§1-7光在两个介质分界面上的反射和折射
一、 内容回顾 二、折、反射波性质的进一步讨论：
§1-8全反射
一、内容回顾： 1.电磁场的边值关系 是研究光在两个介质分界面上的反射和折 射规律的基础。 电磁场的边值关系总结为：尽管两种介质 的分界面上，电磁场量整个的是不连续的， 但在界面上没有自由电荷和面电流时，磁 感应强度矢量和电位移矢量法向分量与电 场强度和磁场强度的切向分量是连续的。
§1-7光在两个介质分界面上的 反射和折射
在规定了电场、磁场的正方向后可以得 到一组关于入射波、反射波、折射波电 场的振幅之间的关系——菲涅尔公式。
rs ts sin( i t ) sin( i t )
rp tp tg ( i t ) tg ( i t ) 2 sin t cos i sin( i t ) cos( i t )
n1 cos i in2 n1 cos i in2
r
s
exp( i
2
rs
)
n r n
p
2
2
cos
cos
i i
i n1 i n1
r
p
exp( i
rp
)
可知：
tg
rs
2
sin i n cos i
2
tg
rp
2
sin i n
rs n1 cos i n2 cos t n1 cos i n2 cos t
n2 cos
2
r0
t
n1 n2 n1 n2
光波在介质界面上的反射和折射 菲涅耳公式
( ki k r ) r 0 ( ki k t ) r 0
(123) (124)
n1 n2 O
kr ki kt
r
B
分界面
(121) (122)
i
t
A C
2.1 反射定律和折射定律 又因为 k n / c ，可将上二式改写为
H ip cos1 H rp cos1 H tp cos 2 (132)
利用
H E ，上式变为
(Eis Ers )n1cos1 Ets n2 cos 2 (133)
3. 菲涅耳公式 再利用折射定律，并由（131）式和（133）式消去 Ets，经整理可得
Ers sin ( 2 1 ) Eis sin ( 2 1 )
sin (1 2 ) rs =sin (1 2 )
(134)
(Eis Ers )n1cos1 Ets n2 cos 2 (133)
3. 菲涅耳公式 利用类似方法，可以推出 p 分量的反射系数和透射系 数表示式， 这就是著名的菲涅耳公式：
O
kr
2
Ers k t
1.s 分量和 p 分量
E p1
H s1
z
1 Hs
E p2 H s2
y
o
E p1
x
2. 反射系数和透射系数 假设介质中的电场矢量为
El E0l e－i(l t－kl r ) l i, r, t ( 127)
其 s 分量和 p 分量表示式为
Elm E0lm e－i(l t－kl r ) m s, p （ 128）
2. 反射系数和透射系数 则定义 s 分量、p 分量的反射系数、透射系数分别为
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
西安邮电大学光学课程设计课程设计名称：光波在介质中界面上的反射及透射特性的仿真一、课程设计目的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 分量的反射系数、透射系数分别为m E E t E E r imtmm 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。