MATLAB光学模拟仿真实训报告模板.pdfVIP

MATLAB光学模拟仿真实训报告模

板(1)(1)(总24页)

--本页仅作为文档封面，使用时请直接删除即可--

--内页可以根据需求调整合适字体及大小--

MATLAB光学模拟仿真

一、实训目的

1、熟悉matlab绘图和仿真功能

2、复习物理光学等相关知识

3、掌握运用matlab软件的仿真等功能与专业知识相结合

二、实训内容

（一）MATLAB基础训练及光波在介质分界面的反射和折射

1、相关原理

当一个单色平面波射到两种不同介质的分界面时，将分为两个波，一个折

射波一个反射波。从电磁场的边值关系可以求出它们的传播方向和入射波的振

幅关系相位关系。

2、实训任务

已知界面两侧的折射率n2、n1和入射角，绘出在n1n2（光由光疏介质

射向光密介质）和n1n2（光由光密介质射向光疏介质）两种情况下，反射系

数、透射系数随入射角的变化曲线。

程序如下：

clear;%清空

disp(请输入介质折射率n1和n2);%在显示括号内语句

n1=input(n1=);%接受键盘任意输入合适的折射率n1

n2=input(n2=);%接受键盘任意输入合适的折射率n2

theta=0::90;%入射角范围范围0~90，步距

a=theta*pi/180;%角度化为弧度

rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-

(n1/n2*sin(a)).^2));%p分量振幅反射率

rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-

(n1/n2*sin(a)).^2));%s分量振幅反射率

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);%创建一个窗口

2

subplot(1,2,1);%作图rp,rs,|rp|,|rs|随入射角的变化曲线

plot(theta,rp,-,theta,rs,--,theta,abs(rp),:,theta,abs(rs),-.,LineWidth,2);%

用‘-’、‘--’，‘：’、‘-.’符号标注对应曲线

legend(rp,rs,|rp|,|rs|);%标注曲线图例

xlabel(入射角/theta_i);%命名x轴

ylabel(振幅);%命名y轴

title([n_1=,num2str(n1),,n2=,num2str(n2),时反射系数随入射角的变化曲

线]);%命名图像

axis([090-11]);%设定作图区间

gridon;%添加网格

subplot(1,2,2);%tp,ts,|tp|,|ts|随入射角的变化曲线

plot(theta,tp,-,theta,ts,--,theta,abs(tp),:,theta,abs(ts),-,LineWidth,2);%用

-‘’、‘--’，‘：’、‘-.’符号标注对应曲线

legend(tp,ts,|tp|,|ts|);%标注曲线图例

xlabel(入射角/theta_i);%命名x轴

ylabel(振幅);%命名y轴

title([n_1=,num2str(n1),,n2=,num2str(n2),时透射系数随入射角的变化曲

线]);%命名图像

ifn1n2%如果此时从光疏到光密

axis([09001]);%设定作图区间

else%否则