matlab反射率&透过率仿真（采用膜层干涉特征矩阵）——可计算无吸收和吸收材质

本片代码作为上一篇的补充，可以进行吸收材质的反射率和透过率计算，即在吸收材质中将实数折射率变更为复数折射率（虚部为材料的消光系数）。代码如下：

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