工程光學(xué)等厚干涉Matlab仿真

1、第五題，并設(shè)置相關(guān)地參數(shù)即四個(gè)圖均為等厚干涉因此相位差，光強(qiáng)地式子均相同，只要更改光程差可.丄I =4cos2()21.平面楔板(1):=2nh 2(2)參數(shù)設(shè)置波長(zhǎng) =500nm 楔角 =0.1rad折射率n =1(3) Matlab仿真程序Sgioit12131415clc:clear all: close all： Lanibda=500e-9: thet a-1 e-1 :n=l;vect orl=linspace (-0, 01f 0. Olj 400); vect or2=linspace (-Q. 0】，0* Olj 4OQ): x, y=meshEridtvect or L v

2、ector2): A=x+t an(t het a):Dm Lt a= 2 *n* A-+L amb da. 2 ;Phi=2#pi*E?lta. /Lambda:I=4*os(Phi/2). "2 mesh(x,yj I);阪設(shè)置垃長(zhǎng)騁設(shè)置楔角暫設(shè)蚤折射率場(chǎng)設(shè)置衷方向向里空間殳畫y方向向重空間喘枸成備y二維向里空間%計(jì)算光程差 鷗計(jì)負(fù)相位差 驚計(jì)尊光強(qiáng)甌根據(jù)光強(qiáng)繪制干涉條紋%在網(wǎng)格片內(nèi)采用顏色插值處理%設(shè)蠱顏色為灰度pcolor'Xj yj IJ ; shading interp colormap gray;(4)仿真結(jié)果2.柱面楔板(1) 光程差分析72 2 2 2 2

3、y x R 二 y- R_xA = y -h 二.R2 -x2 士Z:二 2nhA 2(2)參數(shù)設(shè)置波長(zhǎng) =500nm平行平板距柱面圓心透鏡距離h=0.5mm柱面透鏡半徑R=4.5mm折射率n =1(3)Matlab仿真程序2 -clc;clear all; closs all:3 -Lanbda=SOOe-9;舄設(shè)畫波長(zhǎng)4 -R=4 5：%設(shè)蠱楔角5 -h=0.5：顒殳豊平行平板到柱面團(tuán)心的距廃6 -rl;%設(shè)養(yǎng)折射率7 一vectorl=linspace ' 0. Olj 0. 01, 400 J ;箱設(shè)SxhlnJId里空間8 二vect or2=liriEpace (-0, 0

4、1# 0. 01, 400);9(設(shè)置y方向向里空間9 -x,y =meshgrid(v亡ctorlj vector' ;%構(gòu)咸My】二維向里空間10 -"乃 *0.5-h;11 一D e 11 a= 2 * A n+L amb da/ 2 ;爲(wèi)計(jì)算光程差12 -Phi=2*pi*Delt a. /Lambda:鼎計(jì)算相位差13 -I=4*cos(Phi/2). 2軸計(jì)尊光強(qiáng)14 -mesh(x, y, I);附艮據(jù)光強(qiáng)絵制干涉粲紋15 -pcolor(x,y,I) ; shading irrt亡rp%在網(wǎng)格片內(nèi)采用顏色插值處理16 一colormap gray;%設(shè)置顏色為

5、頁度(4)仿真結(jié)果3.球面楔板(1) 光程差分析72222rR =(R-h°) r 二 h。：RA 二 R - h0 - hAk:=2nhA - _2(2)參數(shù)設(shè)置波長(zhǎng)=500nm-2lu121315160H,04.不規(guī)則柱面楔板(1) 光程差分析以柱面楔板干涉為基礎(chǔ)，將以常數(shù)R為半徑地柱面更改為R = (1O0x)2 0.1x地二次平行平板距球面球心透鏡距離h=0.5mm球面透鏡半徑R=4.5mm折射率n二1(3)Matlab仿真程序(4)仿真結(jié)果殳畫楔角咕殳養(yǎng)平行平板到球面球己的距離 舄設(shè)査折射軍場(chǎng)i殳蠱艮方向向里空間舄i殳蠱y方向向里空間場(chǎng)構(gòu)成xj y_維向量空間 舄將平面直角

6、坐標(biāo)轉(zhuǎn)換為扱坐標(biāo)賠計(jì)算光程差 騎計(jì)算相位差 覊計(jì)算光強(qiáng) 黑根據(jù)光強(qiáng)繪制干涉條紋 黑在網(wǎng)格片內(nèi)采用顏色插值處理 舄將顏色設(shè)畫為換度clc;clear all; close all:Lajnbda=500e9;n=l. 5vect orllinspace (-0, 01, 0. 01, 400);v«ctor2=linspace(-0. Olj 0(H, 400); be, yleshEridCvectcrlj vect or2);r=abs(x+i*y): A=R-r. rt2/R/2-h'*A+n+Lajib da/ 2 ;Phi=2*pi*D?lta,/Lambda: ll

7、4*cos(Phi/2). "2neshI);pcolor(x, y, I):shading interpcolormap gray;0.010.0080.0060.0D40.002-u.002-0 004*0.006-0.0080.002 0.004 D.OOB 0.0080.01°%01 -0WB -0 006 -0.004 -0.002(2)曲線，即可滿足題目要求 參數(shù)設(shè)置波長(zhǎng)=500nm平行平板距柱面圓心透鏡距離h=0.5mm折射率n =1(3)Matlab仿真程序2 - tic:clear all: close all:3 -LauLbda=500?-9 ：4 -

8、h=0- 5;5 -n=l;vectorl=linspace(-0. 01, 0. 01400)vect or2=lmspace (-0. 01t Q. 01400 J ; x, y=jneshgrid( vectorl, vect or2?;R=<100*k), *2+0. l*x：A=(R 2-xr n2),0r5-h;Dlt a-2jKA4n+Lambda/2 : Phi=2*pi*D亡It乩/Lambda:I=4*cob(Phi/2)- "2nesh(k, y, I);p co Lor (x, y, I); shading interp colormap gray;%設(shè)置玻