fan_she_lv_yu_tou_she_lv.zip_MATLAB光反射_光反射_反射率_透射_透射率

clear %清空内存空间 close all%关闭所有作图界面 n1=1; n2=1.45;%介质折射率 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分量振幅反射率 tp=(2*n1*cos(a))./(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分量振幅反射率 ts=(2*n1*cos(a))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));%s分量振幅透射率 figure(1); subplot(1,2,1);%作图区分割成1行2列，第一列 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.', 'LineWidth',2)%作图'r_p','r_s','|r_p|','|r_s|' legend('r_p','r_s','|r_p|','|r_s|') %添加图例 xlabel('\theta_i') %x坐标 ylabel('Amplitude') %y坐标 title(['n_1=',num2str(n1),',n_2=',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)%作图't_p','t_s','|t_p|','|t_s|' legend('t_p','t_s','|t_p|','|t_s|')%添加图例 xlabel('\theta_i')%x坐标 ylabel('Amplitude') %y坐标 title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) %标题 axis([0 90 0 1])%设定作图区域 grid on %为图添加栅格