用matlab绘制不同角度下的透射反射强度图

代码分别讨论了 玻璃到空气、空气到玻璃的两种情况。为了满足精度需要,还在关键地方使用了变步

clc;clear;close all;
%入射 th1 折射 th2
syms rs ts rp tp
th1 = 0:0.001:pi/2;
n1 = 1;n2 = 1.5; n = n1/n2;
% 当从空气中入射玻璃时
sinth2 = sin(th1).*n;
costh2 = sqrt(1-sinth2.^2);
ts = 2*n1.*cos(th1)./(n1.*cos(th1)+n2.*costh2);
rs = ts-1;
figure;
plot(th1,ts, 'LineWidth',2);hold on; 
plot(th1,rs, 'LineWidth',2);hold on; 

tp = 2*n1.*cos(th1)./(n2.*cos(th1)+n1.*costh2);
rp = (n2/n1).*tp-1;
plot(th1,tp, 'LineWidth',2);hold on; 
plot(th1,rp, 'LineWidth',2);hold on; 
% 辅助线绘制
legend('ts','rs','tp','rp');legend('boxoff');
title('从空气中入射玻璃');
x = [0,1.6];
y = [rp(end),rp(end)];
plot(x,y,'HandleVisibility','off','Color','k','LineStyle','--');hold on;
x = [0,1.6];
y = [tp(end),tp(end)];
plot(x,y,'HandleVisibility','o
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值