matlab-對菲涅爾公式畫圖

1、2010～2011學(xué)年第一學(xué)期MATLAB課程1、用MATLAB語言，繪制SPR理論中rp、rs、R對于入射角和折率比值n=n1n2的二維和三維圖形其中n的范圍為（0.6，1.5），并結(jié)合每一個圖形分析其特點及物理意義。2211211222211211cos1sincoscoscoscoscos1sinpnnnnrnnnn????????????????2211112222112211cos1sincoscoscoscoscos1sin

2、snnnnrnnnn????????????????(其中rp為反射光中平行分量的反射系數(shù)；rs為反射光中垂直分量的反射系數(shù)；R反射率；n1、n2分別表示兩種不同介質(zhì)的折射率)答案：程序：（1）%反射系數(shù)反射系數(shù)r，和反射率，和反射率R與折射率之比與折射率之比n的關(guān)系；的關(guān)系；clearclcclfn=0.6:0.05:1.5zeta1=pi10%入射角zeta2=real(asin(n.sin(zeta1)))%折射角rpz=n.co

3、s(zeta2)cos(zeta1)%平行分量反射部分分子rpm=n.cos(zeta2)cos(zeta1)%平行分量反射部分分母rp=rpz.rpm%平行分量反射系數(shù)rsz=n.cos(zeta1)cos(zeta2)%垂直分量反射部分分子rsm=n.cos(zeta1)cos(zeta2)%垂直分量反射部分分母rs=rsz.rsm%垂直分量反射系數(shù)Rp=rp.^2%平行分量反射率Rs=rs.^2%垂直分量反射率t=0.6:0.00

4、5:1.5y=0y1=0.4:0.005:0t1=1plot(nrpcnrs:nRpknRsgtyt1y1)text(1.30.1rp)text(1.20.09rs)text(0.650.023Rp)text(0.650.065Rs)text(0.90.2(n1))title(relationshipsbetweenrprsRpRsn=n1n2)Rs=rs.^2%垂直分量反射率Brewster=acot(n)180pi%布魯斯特角y=0

5、y1=1:0.005:0plot(xrpxrsxRpxRsxyBrewstery1c:)text(100.25rp)text(400.2rs)text(770.1Rp)text(720.3Rs)text(28.50.6Brewsangle56.3099)xlabel(xangle)ylabel(rR)title(TherelationshipsofrprsRsRpincidentanglex||(n0，說明此范圍內(nèi)，反射光與入射光中的p