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

jxxl_1314

已于 2023-06-24 22:36:42 修改

阅读量1.7k

点赞数 1

分类专栏： matlab学习文章标签： matlab 矩阵数学建模

于 2023-06-24 07:29:56 首次发布

本文链接：https://blog.youkuaiyun.com/jxxl_1314/article/details/131357222

版权

这段代码使用特征矩阵方法来计算无吸收单层膜的反射率和透射率，需要输入介质1、介质3的折射率以及膜层2的厚度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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(