clear;             %清空内存空间
disp('请输入介质折射率n1和n2')
n1=input('n1='); %接受键盘任意输入合适的折射率n1
n2=input('n2='); %接受键盘任意输入合适的折射率n2
theta = 0:0.1:90;   %入射角范围0-90°,步距0.1°
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);
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');
ylabel('振幅');
title(['n1=',num2str(n1),',n2=',num2str(n2),'时反射系数随入射角的变化曲线']);
axis([0 90 -1 1]);    %设定作图区间
grid on;              %作图加栅格
subplot(1,2,2);
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',...theta,abs(ts),'-.','LineWidth',2);
legend('tp','ts','|tp|','|ts|');
xlabel('入射角\theta_i');
ylabel('振幅');
title(['n1=',num2str(n1),',n2=',num2str(n2),'时透射系数随入射角的变化曲线']);
if n1<n2
axis([0 90 0 1]);
else
axis([0 90 0 3.5]);
end
grid on;``````handlebars
>> li_3_1_1
请输入介质折射率n1和n2
n1=1
n2=2
li_3_1_1
请输入介质折射率n1和n2
n1=3
n2=4

MATLAB----光学之分量振幅反射率和分量振幅透射率相关推荐

  1. 读取4:2:0格式YUV序列的Y分量、U分量以及V分量,并分别保存为.yuv格式(matlab实现)

    读取4:2:0格式YUV序列的Y分量.U分量以及V分量,并分别保存为.yuv格式,不多说,具体程序如下: read_Y_U_V.m %% 读取4:2:0格式YUV的Y.U和V分量,然后将三个分量分别保 ...

  2. MATLAB相干成像系统,Matlab光学仿真课程设计-基于Matlb相干与非相干照明成像系统的仿真.docx...

    Matlab光学仿真课程设计-基于Matlb相干与非相干照明成像系统的仿真 东 北 石 油 大 学课 程 设 计课课 程 Matlab光学仿真课程设计 题 目 基于Matlab相干与非相干照明 成像系 ...

  3. 反人类的定义《电力系统分析》——电压的纵分量与横分量

    反人类的定义<电力系统分析>--电压的纵分量与横分量 纵分量与横分量是以"导线的视角"命名的.而不是笛卡尔坐标系.沿着导线方向为纵,垂直于导线方向为横.

  4. matlab将yuv文件的luma分量矩阵转换为图像显示

    问题描述: 我将视频文件中的luma分量提取到txt文本中后,在matlab中存到了数值矩阵中,但当用imshow函数进行显示时居然是白色.其原因如下:(转载于别人的博客:http://blog.si ...

  5. 颜色空间直方图matlab,使用Matlab绘制图像的rgb颜色空间和Lab颜色空间分量图和分量直方图 | 学步园...

    %author sagittariusboy All right reserved 2013年7月28日 14:58:11 org = imread('laaa.jpg'); %org = imrea ...

  6. matlab中lab颜色空间,使用Matlab绘制图像的rgb颜色空间和Lab颜色空间分量图和分量直方图...

    %author sagittariusboy All right reserved 2013年7月28日 14:58:11 org = imread('laaa.jpg'); %org = imrea ...

  7. 小波功率谱 matlab,小波变换后的各频率分量的功率谱,

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 下面程序是用来提取脑电波信号的,利用小波提取四种频率分量的波α.β.δ.θ,现在不明白如何求频率分量的功率谱? 其中一段代码是: %dalt fs1=ff ...

  8. matlab光学毕业论文,光学信息处理实验的Matlab仿真.doc

    摘要:光学实验系统往往比较精密和复杂,理论课堂教学中加入相关实验的演示往往比较费时,而且也不一定能得到正确的实验结果.随着计算机技术的发展,强大的运算性能和图形处理能力可以很好的解决这个问题,满足课堂 ...

  9. 从光波叠加到条纹分布的matlab仿真,matlab光学仿真实验 13171019

    光学仿真实验 一.前言 此次光学仿真实验,是基于matlab来进行的.在这仿真 的一系列过程中,对于光学现象出现的条件,以及干涉. 衍射是光波叠加的本质都有了更深的认识.还从中学习 了matlab这一 ...

最新文章

  1. python调用webservice接口实例_python调用webservice接口的实现
  2. Eclipse中的集成Git插件删除线上远程分支
  3. [设计模式] ------ 原型模式(浅拷贝和深拷贝)
  4. SpringBoot2 集成 xxl-job任务调度中心_路由策略
  5. Python判断不可变对象(字符串,整数,浮点数,数组)相等的办法以及其底层实现原理
  6. QT每日一练day5:QLabel和按钮窗口打印功能
  7. JAXB实现java对象与xml之间转换
  8. 应用程序工程文件组织重要性
  9. LeetCode-326. Power of Three
  10. python模块文件的扩展名不一定是py_Python导入:导入没有.py扩展名的模块?
  11. @EnableWebMVC注解理解
  12. 2022最新版VMware虚拟机及CentOS-7安装教程
  13. CleanMyPC比360管家好用N倍的电脑清理软件
  14. 算法的两大性能指标:时间复杂度和空间复杂度
  15. 当你的浏览器页面被篡改时,我们该怎么办
  16. android捕获全局异常lin,全局获取 (Activity)Context,实现全局弹出 Dialog
  17. 差分 --算法竞赛专题解析(32)
  18. 人工智能时代,如何让机器狗听懂你说的话?
  19. JAVA第二次作业《胖瘦程度计算》
  20. 3椭圆曲线密码学:ECDH和ECDSA

热门文章

  1. 史上最简单的openshift免费空间上传代码教程!没有之一!
  2. 转 为什么数码相机可以拍出彩色照片?
  3. c语言如何反复执行一段程序,C语言中重复执行程序的问题
  4. 矩阵在游戏开发中的应用
  5. comparison lemma
  6. 数字助手:Siri、Google和Alexa之间的PK
  7. Unity获取组件的几种方式(拖拽法、标签法、名字法)
  8. 盛世昊通解析为什么你买不到想要的新能源汽车
  9. ORACLE US7ASCII编码 读取时乱码问题
  10. 四个角不是直角的四边形_四边形的特点是有四条直的边和四个直角对吗