2009～2010 学年第一学期 MATLAB 课程 研究生考题： 1、用 MATLAB 语言，绘制 SPR 理论中 rp、rs、R 对于入射角和折率比值 n=n1/n2的二维和三维图形,其中 n 的范围为（0.6，1.5），并结合每一个图形分析其特点及物理意义。 2211211222211211cos1sincoscoscoscoscos1sinpnnnnrnnnn 2211112222112211cos1sincoscoscoscoscos1sinsnnnnrnnnn (其中 rp 为反射光中平行分量的反射系数；rs 为反射光中垂直分量的反射系数； R-反射率；n1、n2 分别表示两种不同介质的折射率) 答案： 程序： （1）%反射系数r，和反射率R与折射率之比n的关系； clear; clc; clf; n=0.6:0.05:1.5; zeta1=pi/10; %入射角 zeta2=real(asin(n.*sin(zeta1))); %折射角 rpz=-n.*cos(zeta2)+cos(zeta1); %平行分量反射部分分子 rpm=n.*cos(zeta2)+cos(zeta1); %平行分量反射部分分母 rp=rpz./rpm; %平行分量反射系数 rsz=n.*cos(zeta1)-cos(zeta2); %垂直分量反射部分分子 rsm=n.*cos(zeta1)+cos(zeta2); %垂直分量反射部分分母 rs=rsz./rsm; %垂直分量反射系数 Rp=rp.^2; %平行分量反射率 Rs=rs.^2; %垂直分量反射率 t=0.6:0.005:1.5; y=0; y1=-0.4:0.005:0; t1=1; plot(n,rp,'c',n,rs,'*:',n,Rp,'k-',n,Rs,'g',t,y,'-',t1,y1,'-') text(1.3,-0.1,'rp'); text(1.2,0.09,'rs'); text(0.65,0.023,'Rp'); text(0.65,0.065,'Rs'); text(0.9,-0.2,'(n<1)'); text(1.2,-0.2,'(n>1)'); title('relationships between rp,rs,Rp,Rs and n=n1/n2'); xlabel('n=n1/n2'); ylabel('r or R'); （1）分析其物理意义：（取入射角pi/10，小于Brew ster 角，小于临界角） 当n<1（光疏到光密介质）时，反射光中垂直分量（s）的反射系数rs<0，说明反射光中s 分量与入射光中的s 分量相位相反；而p 分量，即平行分量，其反射系数rp>0，说明反射光中p 分量与入射光中的p 分量相位相同。 当n=1 时，相当于光束在同一种介质中传输，因此无反射，反射系数为零； 当n>1（光密到光疏介质）时，反射光中rs>0，说明反射光中s 分量与入射光中的s 分量相位相同；而p 分量rp<0，说明p 分量在两种介质中相位相反。 （2） （2-1）%反射系数r，和反射率R与入射角x的关系（n<1）； clear; clc; clf; n=0.6:0.05:1.5; zeta1=pi/10; %入射角 zeta2=real(a...
