案例1：光疏到光密介质

平面光波从空气（折射率为）入射到石英玻璃中（折射率为），用 MATLAB作出p、s分量的振幅反射率和振幅透射率以及它们的绝对值随入射角度的变化曲线。

MATLAB代码：

在 MATLAB中分别调用平面光波的电矢量平行入射面的情况下的菲涅尔公式，该情况也常常被称为TE波或者p-偏振入射：

再调用平面光波的电矢量垂直入射面的情况下的菲涅尔公式，该情况也常常被称为TM波或者s-偏振入射：

计算出它们的绝对值，然后作图即可。

MATLAB代码如下：

clear; %清空内存空间
close all; %关闭所有的作图页面
n1 = 1; n2 = 1.45; %介质折射率
theta=0:0.1:90; %入射角范围0—90度，步距0.1度
a=theta*pi/180; %角度化为弧度
rp = (n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./... %p分量振幅反射率(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
rs = (n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./... %s分量振幅反射率(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)); %p分量振幅透射率
ts = 2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); %s分量振幅透射率
figure(1);
subplot(1,2,1); %作图rp、rs、｜rpl、｜rs｜随入射角的变化曲线
plot(theta,rp,' -',theta,rs,'--',theta,abs(rp),':',....theta,abs (rs) , ' -.' ,'LineWidth',2)
legend('r_p','r_s','|r_p|','|r_s|')
xlabel('\theta_i');
ylabel('Amplitude');
title( ['n_1=',num2str(n1),', n_2=',num2str(n2)]);
axis([0 90 -1 1]) %设定作图区间
grid on; %作图加栅格
subplot(1,2,2); %作图rp、rs、｜rp｜、｜rs｜随入射角的变化曲线
plot (theta,tp,'-' ,theta,ts,'--',theta,abs(tp),':',...theta,abs (ts) , '-.' ,'LineWidth',2)
legend('t_p','t_s','|t_pl','|t_s|');
xlabel('\theta_i');
ylabel('Amplitude');
title (['n_1=',num2str(n1),', n_2=', num2str(n2)]);
axis([0 90 0 1]);
grid on;

结果：

源代码下载地址：菲涅尔公式MATLAB作图：振幅反射率和振幅透射率随入射角度的变化曲线资源-CSDN文库

案例2：光密到光疏介质

平面光波从石英玻璃（折射率为）入射到空气中（折射率为），用 MATLAB作出p、s分量的振幅反射率和振幅透射率以及它们的绝对值随入射角度的变化曲线。

MATLAB代码：

将上述代码中的和数值改动一下，调整一下图形的y轴范围，重新作图，代码如下：

clear       %清空内存空间
close all   %关闭所有的作图页面n1 = 1.45; n2 = 1;  %介质折射率
theta = 0:0.1:90;   %入射角范围0-90度，步距0.1度
a = theta*pi/180;   %角度化为弧度rp = (n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./...      %p分量振幅反射率(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
rs = (n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./...      %s分量振幅反射率(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));%p分量振幅透射率
ts = 2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));%s分量振幅透射率figure(1);
subplot(1,2,1);      %作图rp、rs、|rp|、|rs|随入射角的变化曲线
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',... theta,abs(rs),'-.','LineWidth',2)
legend('r_p','r_s','|r_p|','|r_s|')
xlabel('\theta_i')
ylabel('Amplitude')
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
axis([0 90 -1 1.2])   %设定作图区间
grid on             %作图加栅格subplot(1,2,2);     %作图rp、rs、|rp|、|rs|随入射角的变化曲线
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',...theta,abs(ts),'-.','LineWidth',2)
legend('t_p','t_s','|t_p|','|t_s|')
xlabel('\theta_{i_1}')
ylabel('Amplitude')
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
axis([0 90 0 3])
grid on

结果：

本文未经允许，不得转载。

关注CSDN“光学码农”！

MATLAB：菲涅尔公式（反射/透射公式）相关推荐

Matlab 菲涅尔系数计算
1 简介 Matlab 菲涅尔系数计算 2 部分代码 function varargout = Fresnel(varargin) % FRESNEL MATLAB code for Fresnel ...
matlab菲涅尔衍射_有问必答——SYNOPSYS安装体验课堂——可以设计菲涅尔透镜吗?...
问:SYNOPSYS可以设计菲涅尔透镜吗? 答:在USS中有多种菲涅尔面型,用户输入参数即可. 问:SYNOPSYS中具有的输入方式? 答:大家总是有个误区,以为SYNOPSYS需要输入命令运行,其实 ...
matlab 菲涅尔环,菲涅耳区讨论 - 通信原理与基础 - 通信人家园 - Powered by C114
1．菲涅耳区当需要计算传播主区的几何尺寸时,要应用惠更斯-菲涅尔原理.惠更斯-菲涅尔原理认为,波在传播过程中,波面上的每一点都是一个进行二次辐射球面波(子波)的波源,而下一个波面,就是前一个波面所辐 ...
庄懂的TA笔记（九）＜菲涅尔 + MatCap + CubeMap＞
庄懂的TA笔记(九)<菲涅尔 + MatCap + CubeMap> 课程内容: 1.菲涅尔 0.效果展示.公式: 菲涅尔现象: Fresnel(菲涅尔) 用到的两个向量 NdirWS ...
学习笔记（二）：WiFi-无接触感知与菲涅尔区
学习笔记(二):WiFi-无接触感知与菲涅尔区 <毫米级的 Wi-Fi 无接触感知:从模式到模型>笔记整体总结: 本文在理论上揭示和评估什么样的人体行为可被检测以及何时被检测到.具体 ...
菲涅尔区，前后比，VSWR
1.菲涅尔区对于垂直化天线系统,波从奇数菲涅尔区反射时产生的信号更强,偶数则弱. 对于水平化天线系统,波从奇数菲涅尔区反射时产生的信号更弱,偶数则强. 2.后瓣辐射方向与主瓣相反 3.前后比主瓣 ...
菲涅尔公式实现边缘光效果
Shader"review/EdgeEmission" {///标准光照模型光照数据=自发光(E)+漫反射光(D)+环境光(A)+高光(S)Properties{_MainTex ...
基于Matlab——夫琅禾夫衍射以及菲涅尔衍射
我以往在学习F分析的时候,编写了一个作业代码,在此附上供大家学习交流. 引言: 在傅里叶光学信息基础中,主要研究的是光在传播过程携带的信息如何去检测得到.如果光在自由空间(均匀 ...
菲涅尔定理及MATLAB实现
** 菲涅尔定理及MATLAB实现本文是学习胡章芳<MATLAB仿真及其在光学课程中的应用>和赵凯华<光学>这两本书时所做笔记. 1. 菲尼尔反射与折射公式一般来说 ...

MATLAB：菲涅尔公式（反射/透射公式）

案例1：光疏到光密介质

案例2：光密到光疏介质

MATLAB：菲涅尔公式（反射/透射公式）相关推荐

最新文章

热门文章