Demo entry 6745148

平面光波从空气(折射率n1=1)入射到石英玻璃中(折射率n2=1.45),做出p,s分量的振幅反射率和振幅透射率及它们的绝对值随入射角度的变化曲线。

   

Submitted by 贾天旭 on May 26, 2018 at 04:48
Language: Matlab. Code size: 956 Bytes.

clear
close all
clc
n1 = 1; n2 = 1.45; 
theta = 0:0.1: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));           
rs = (n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./... 
    (n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
tp = 2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
ts = 2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));

figure(1);     
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',...
    theta,abs(rs),'-.','LineWidth',2)
legend('r_p','r_s','|r_p|','|r_s|',2)
xlabel('\theta_i')
ylabel(' Reflectivity ')
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
axis([0 90 -1 1])
grid on

figure(2)
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',...
    theta,abs(ts),'-.','LineWidth',2)
legend('t_p','t_s','|t_p|','|t_s|')
xlabel('\theta_i')
ylabel(' Transmission ')
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
axis([0 90 0 1])
grid on

This snippet took 0.01 seconds to highlight.

Back to the Entry List or Home.

Delete this entry (admin only).