matlab中透射率&反射率计算单层膜（采用特征矩阵方法进行计算）

1. 这里提供的代码是采用膜层的干涉特征矩阵来进行反射率和透过率计算，其中介质1和介质3的折射率和膜层2的厚度需要进行输入；该代码计算的是无吸收的膜层反射率和透过率。

 代码如下：

clear all;
close all;
clc
h=input('输入膜层厚度：');
n1=input('输入膜层上方介质折射率：');
n3=input('输入膜层下方介质折射率：');
data1=importdata('D:\数据\nk值\PMMA.txt');
data2=data1.data;%将data1中的数值赋值给data2；
x=1000.*data2(:,1);%折射率对应波长；
y1=data2(:,2);%n值；
p=polyfit(x,y1,3);%3次多项式拟合，获取对应拟合系数；
theta=0:1:90;%入射角定义，0-90°，step 1°；
a=theta*pi/180;%角度转换为弧度；
w=380:1:780;%波长定义，380-780nm，step 1nm;
for i=1:401
    n2(i)=polyval(p,w(i));%材料不同波长下折射率实部拟合值；
    for m=1:91
        c2(i,m)=sqrt(1-(n1/n2(i)*sin(a(m)))^2);%介质中折射角度的余弦值cosθ1计算；
        p2(i,m)=2*pi*h/w(i)*n2(i)*c2