matlab入门 — 反射率、折射率的变化
matlab的入门,程序为书本上摘抄的,注释为自己添加的。软件模拟过程为自己在matlab软件上实际模拟的,目的记录学习过程,与大家分享互相学习。
用简单的代码实现发射率,折射率的变化曲线。平面光波的入射波、反射波、折射波的电矢量可以分成两个矢量:一个平行于入射面,一个垂直于入射面,分别用p和s表示。
平面波从空气(折射率n1 =1)入射到石英玻璃中(折射率n2 = 1.45)。振幅反射率为rp、rs,振幅折射率为tp、ts。代码如下:
clc %清空内存空间
clear all %清理所有数据
close all %关闭所有的作图界面
n1 = 1,n2 = 1.45; %介质折射率
theta = 0: 0.1:90; %入射角范围0-90度,步距0.1度
a = theta *pi /180; %角度化为弧度
rp = (n2cos(a)-n1sqrt(1-(n1/n2sin(a)).^2))./…
(n2cos(a)+n1sqrt(1-(n1/n2sin(a)).^2)); %p分量振幅反射率
rs = (n1cos(a)-n2sqrt(1-(n1/n2sin(a)).^2))./…
(n1cos(a)+n2sqrt(1-(n1/n2sin(a)).^2)); %s分量振幅反射率
tp = 2n1cos(a)./(n2cos(a)+n1sqrt(1-(n1/n2sin(a)).^2)); %p分量振