




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、精選優(yōu)質(zhì)文檔-傾情為你奉上 1、 課程設(shè)計題目: 用matlab仿真光束的傳輸特性。2、 任務(wù)和要求 用matlab仿真光束通過光學(xué)元件的變換。 設(shè)透鏡材料為k9玻璃,對1064nm波長的折射率為1.5062,鏡片中心厚度為3mm,凸面曲率半徑,設(shè)為100mm,初始光線距離透鏡平面20mm。用matlab仿真近軸光線(至少10條)經(jīng)過平凸透鏡的焦距,與理論焦距值進(jìn)行對比,得出誤差大小。 已知透鏡的結(jié)構(gòu)參數(shù)為,(K9玻璃),物點A距第一面頂點的距離為100,由A點計算三條沿光軸夾角分別為10、20、30的光線的成像。試用Matlab對以上三條光線光路和近軸光線光路進(jìn)行仿真,并得出實際光線的球差大
2、小。 設(shè)半徑為1mm的平面波經(jīng)凸面曲率半徑為25mm,中心厚度3mm的平凸透鏡。用matlab仿真平面波在透鏡幾何焦平面上的聚焦光斑強度分布,計算光斑半徑。并與理論光斑半徑值進(jìn)行對比,得出誤差大小。(方法:采用波動理論,利用基爾霍夫菲涅爾衍射積分公式。)2、用MATLAB仿真平行光束的衍射強度分布圖樣。(夫朗和費矩形孔衍射、夫朗和費圓孔衍射、夫朗和費單縫和多縫衍射。)3、用MATLAB仿真厄米高斯光束在真空中的傳輸過程。(包括三維強度分布和平面的灰度圖。)4、(補充題)查找文獻(xiàn),掌握各類空心光束的表達(dá)式,采用費更斯-菲涅爾原理推導(dǎo)各類空心光束在真空中傳輸?shù)墓鈴姳磉_(dá)式。用matlab對不同傳輸距
3、離處的光強進(jìn)行仿真。三、理論推導(dǎo)部分將坐標(biāo)原點選在透鏡中心處,1=arcsin(y1/r),由n1*sin1=n2*sin2可得出2=arcsin(n1/n2)*(y1/r),由幾何關(guān)系可得到=21,則出射光線的斜率k=tan(21),當(dāng)入射直線y=y1時,x1=d(r),并設(shè)出射直線為y=k*x+b;由直線經(jīng)過(x1,y1)即可求出b值,從而就可以求出射直線。由單透鏡焦點計算公式1/f=(n-1)*(1/r11/r2)可求得f=193.6858。利用近軸光學(xué)公式i1=(l1-r1)*u1/r1,i11=n1*i/n11,u11=u1+i1-i11l11=r1+r1*i11/u11和轉(zhuǎn)面公式u
4、2=u11,l2=l11-d1可以求得u11、u22、l22、h2等。入射光線的夾角為u1,設(shè)入射光線為y1=k1*x1+b1其中的斜率k1=u1又由于入射光線經(jīng)過經(jīng)過(-100,0)就可以求出b1。由h1=l1*u1即為y1,當(dāng)y1為定值時就可以得到第一個橫坐標(biāo)x0,再利用最后的出射光線公式y(tǒng)3=k3*x3+b3,k3=u22,又因為最終出射經(jīng)過(d+l22,0)可求出b3,利用轉(zhuǎn)面公式h2=h1d*u11,即為y3可求出第二個橫坐標(biāo)x00。再求在透鏡中的直線斜率k2=(h2h1)/(x00x0),y2=k2*x2+b2經(jīng)過(x0,h1)即可求得b2值,從而即可求得三條直線。實際光束求法同理
5、。利用菲涅耳近似公式求衍射面上的光強要對孔徑上的點求積分可以轉(zhuǎn)換成對其x1,y1的微分求和,其中公式中的z1=f。2. (1)夫朗和費矩形孔衍射若衍射孔為矩形則在透鏡焦平面上得到的衍射圖樣如圖,衍射圖樣的主要特征為衍射亮斑集中分布在兩個相互垂直的方向上,并且x軸上的亮斑寬度與y軸亮斑寬度之比,恰與矩形孔在兩個軸上的寬度相反。 其中的為x,同樣的中的為y,利用x=x/f,y=y/f進(jìn)行求解。(2)夫朗和費圓形孔衍射夫朗和費圓孔衍射的討論方法和矩形孔衍射的討論方法相同,只是由于圓孔的幾何對稱性,采用極坐標(biāo)更為方便。=ka(3)夫朗和費單縫衍射對于前面討論的夫朗和費矩形孔衍射,如果矩形的一個方向的尺
6、寸比另一個方向大得多,則該矩形孔衍射就變成單縫衍射(如圖),這時沿y方向的衍射效應(yīng)不明顯,只在x方向有亮暗變化的衍射圖樣。實驗中通過利用=x/f進(jìn)行求解(4)夫朗和費多縫衍射夫朗和費多縫衍射裝置如圖,其每條狹縫均平行于y1方向,沿x1方向的縫寬為a,相鄰狹縫的間距為d,在研究多縫衍射時,由于后透鏡的存在使衍射屏上每個單縫的衍射條紋位置與位置無關(guān)。因此,用平行光照射多縫時,其每一個單縫都要產(chǎn)生自己的衍射,形成各自一套衍射條紋。當(dāng)每個單縫等寬時,各套衍射條紋在透鏡焦平面上完全重疊,其總光強分布為它們的干涉疊加。四、Matlab仿真部分clear allr=100;n1=1.5163;n2=1;%透
7、鏡的曲率半徑為100mm,透鏡的折射率n1=1.5,空氣的折射率n2=1d=3;%x=77:0.1:320;figure(1)for n=-5:5 y1=0.1*n; %hold on; %plot(x1,y1); a1=asin(y1/r);%入射角 a2=asin(n1/n2*(y1/r);%折射角 a=a2-a1; k=tan(a);%出射光線的斜率 x1=sqrt(r2-y12); x2=x1-r+d; b=y1+k*x2; %出射光線經(jīng)過(x2,y1) x=-20:0.01:x2;%零坐標(biāo)選在透鏡中心,入射光線距透鏡20mm,故x=-20 hold on plot(x,y1);%平行
8、光束 x3=x2:0.01:300; y=-k*x3+b;%出射光線 hold on plot(x3,y);Endclear all%透鏡的結(jié)構(gòu)參數(shù)r1=10;r2=-50;l1=-100;L1=-100;n1=1.0;d1=5;n11=1.563n2=1.563;n22=1.0;figure(1)for n=-3:-1%沿光軸分別為1、2、3度的光線進(jìn)行入射 %近軸光學(xué)成像公式第一個面 u1=n;i1=(l1-r1)*u1/r1i11=n1*i/n11;u11=u1+i1-i11;l11=r1+r1*i11/u11;%轉(zhuǎn)面公式u2=u11;l2=l11-d1;%近軸光學(xué)成像公式第二個面i2=
9、(l2-r2)*u2/r2;i22= n2*i2/n22;u22=u2+i2-i22;l22=r2+r2*i22/u22;%入射光線與第一個透鏡交點的縱坐標(biāo),坐標(biāo)原點選在第一個透鏡的頂點處h1=l1*(u1*pi/180);k1=-u1*pi/180;%入射光線的斜率b1=100*k1;%因為入射光線經(jīng)過(-100,0)點x0=(h1-b1)/k1;%入射光線與第一個透鏡交點的橫坐標(biāo)x1=-100:0.01:x0;y1=k1*x1+b1;hold onplot(x1,y1);%輸出入射光線k3=-u22*pi/180;%第二次折射后出射光線的斜率b3=-k3*(d1+l22);%因為第二次折射
10、后出射光線經(jīng)過(d1+l22,0)點h2=h1-d1*(u11*pi/180);%第一次折射后入射到第二個透鏡的縱坐標(biāo)x00=(h2-b3)/k3;%第一次折射后入射到第二個透鏡的橫坐標(biāo)k2=(h2-h1)/(x00-x0);%第一次折射后光線的斜率b2=h1-k2*x0;%因為第一次折射后光線經(jīng)過(x0,h1)點x2=x0:0.01:x00;y2=k2*x2+b2;hold onplot(x2,y2);%輸出第一次折射在兩個透鏡中的光線x3=x00:0.01:30;%選在30是為了將輸出圖形看得更清晰些y3=k3*x3+b3;hold onplot(x3,y3);%輸出經(jīng)過第二個透鏡后的輸出
11、光線%實際光路U1=n*pi/180;I1=asin(L1-r1)*sin(U1)/r1);I11=asin(n1*sin(I1)/n11);U11=U1+I1-I11;L11=r1+r1*sin(I11)/sin(U11);%轉(zhuǎn)面公式U2=U11;L2=L11-d1;%實際光學(xué)成像公式第二個面I2=asin(L2-r2)*sin(U2)/r2);I22=asin(n2*sin(I2)/n22);U22=U2+I2-I22;L22=r2+r2*sin(I22)/sin(U22);%入射光線與第一個透鏡交點的縱坐標(biāo),坐標(biāo)原點選在第一個透鏡的頂點處h3=L1*tan(U1);k4=-tan(U1)
12、;%入射光線的斜率b4=100*k4;%因為入射光線經(jīng)過(-100,0)點x01=(h3-b4)/k4;%入射光線與第一個透鏡交點的橫坐標(biāo)x4=-100:0.01:x01;y4=k4*x4+b4;hold onplot(x4,y4,'r');%輸出入射光線k6=-tan(U22);b6=-k6*(d1+L22);%因為第二次折射后出射光線經(jīng)過(d1+L22,0)點h4=h3-d1*tan(U11);%第一次折射后入射到第二個透鏡的縱坐標(biāo)x02=(h4-b6)/k6;%第一次折射后入射到第二個透鏡的橫坐標(biāo)k5=(h4-h3)/(x02-x01);%第一次折射后光線的斜率b5=h4
13、-k5*x02;%因為第一次折射后光線經(jīng)過(x02,h4)點x5=x01:0.01:x02y5=k5*x5+b5;hold onplot(x5,y5,'r');%輸出第一次折射在兩個透鏡中的光線x6=x02:0.01:30;%選在30是為了將輸出圖形看得更清晰些x6=x02:0.01:30;y6=k6*x6+b6;hold onplot(x6,y6,'r');%輸出經(jīng)過第二個透鏡后的輸出光線%球差m=(L22+d1)-(l22+d1);endclear alln=1.5062;%K9玻璃的折射率d=3;%透鏡的中心厚度R=25;%透鏡凸面曲率半徑f=R/(n-1
14、);%透鏡焦距R0=1;%入射光束半徑lambda=1.064e-3;%波長k=2*pi/lambda;phy=lambda*0.61/R0;%角半徑w0=sqrt(f*lambda/pi);%實際光斑半徑data=w0-f*phy; %誤差z=f;rmax=3*f*phy;%艾利斑半徑r=linspace(0,rmax,100);%產(chǎn)生從0到rmax之間的100點行矢量 將衍射半徑100等分eta=linspace(0,2*pi,100);%將0到2*pi100等分rho,theta=meshgrid(r,eta);%生成繪制3D圖形所需的網(wǎng)格數(shù)據(jù)x,y=pol2cart(theta,rho
15、);%衍射斑某點的坐標(biāo)轉(zhuǎn)換極坐標(biāo)到直角坐標(biāo)r0=linspace(0,R0,100);%將入射光束半徑100等分eta0=linspace(0,2*pi,100);rho0,theta0=meshgrid(r0,eta0);x0,y0=pol2cart(theta0,rho0);for dx=1:100%都是為了建立網(wǎng)格 for dy=1:100 Ep=-i/(lambda*z)*exp(i*k*z)*exp(i*k*(x-x0(dx,dy).2+(y-y0(dx,dy).2)/(2*z); E2(dx,dy)=sum(Ep(:);%積分公式的求和表達(dá) endendIe=conj(E2).*E
16、2;%光強表達(dá)式figure(1);surf(x,y,Ie);figure(2)plot(x(50,:),Ie(50,:);2. (1)夫朗和費矩形孔衍射clear all;lamda=500e-9;a=1e-3;b=1e-3;f=1;m=500;ym=8000*lamda*f;ys=linspace(-ym,ym,m)xs=ys;n=255;for i=1:m sinth2=ys./sqrt(ys.2+f2);%相當(dāng)于x/f sinth1=xs(i)/sqrt(xs(i).2+f2);%xs(i)作用每給一個ys值,要遍歷到所有的x值 angleA=pi*a*sinth1/lamda;%相當(dāng)
17、于書上的alfa=kax/2f k=2*pi/lamda angleB=pi*b*sinth2./lamda;B(:,i)=(sin(angleA).2.*sin(angleB).2.*5000./(angleA.2.*angleB.2);%光強度公式endsubplot(1,2,1)image(xs,ys,B)colormap(gray(n)subplot(1,2,2)plot(B(m/2,:),ys)(2)夫朗和費圓孔衍射clearlam=500e-9a=1e-3f=1m=300;ym=5*0.61*lam*f/a;%取愛里光斑半徑的5倍ys=linspace(-ym,ym,m);xs=y
18、s;n=200;for i=1:m r=xs(i)2+ys.2;%相當(dāng)于r的平方 sinth=sqrt(r./(r+f2);%角度 fai=2*pi*a*sinth./lam;%fai=k*a*sinth hh=(2*BESSELJ(1,fai).2./fai.2;%貝塞爾函數(shù) b(:,i)=hh.*5000;endsubplot(1,2,1)image(xs,ys,b)colormap(gray(n)subplot(1,2,2)b(:,m/2)plot(ys,b(:,m/2)(3)夫朗和費單縫衍射clear alllam=500e-9;a=3;f=1;xm=3*lam*f/a;nx=50;x
19、=linspace(-xm,xm,nx);ny=50;y=linspace(0,a,ny);for i=1:nysinphi=x/f;%角af=(pi*a*sin(sinphi)/lam; I(i,:)=5*(sin(af)./af).2;endN=255;%確定灰度等級Br=(I/max(I(1,:)*N;subplot(1,2,1)image(x,y,Br);colormap(gray(N);%顏色subplot(1,2,2)plot(x,I(1,:);(4)夫朗和費多縫衍射clear all;lamda=500e-9; %波長N=2; %縫數(shù),可以隨意更改變換a=2e-4;f=5;d=5
20、*a;ym=2*lamda*f/a;%選擇坐標(biāo)范圍xs=ym;n=1001;ys=linspace(-ym,ym,n);for i=1:nsinphi=ys(i)/f;alpha=pi*a*sinphi/lamda;fai=2*pi*d*sinphi/lamda;I1=(sin(alpha)./alpha).2;%單縫衍射因子B(i,:)=I1*(sin(N*fai/2)./sin(fai/2).2;%多縫衍射光強的計算公式B1=B/max(B);%歸一化光強endNC=256; %確定灰度的等級Br=(B/max(B)*NC;subplot(1,2,1)image(xs,ys,Br);col
21、ormap(gray(NC); %色調(diào)處理subplot(1,2,2)plot(B1,ys,'k');五、畫出仿真圖形2. (1)夫朗和費矩形孔(2) 夫朗和費圓孔衍射(3)夫朗和費單縫衍射(4)夫朗和費多縫衍射六、仿真過程中發(fā)生的問題及解決的方法1. 在用matlab寫程序時一定要保證是在英文狀態(tài)下輸入標(biāo)點符號的不然程序運行時會報錯。2. 在求第一個圖形y=-k*x+b時由于忘記加負(fù)號,使出射光線是發(fā)散的。由于在算b值,忘記了加負(fù)號將b=y-k*x使入射光線和出射光線接不上。3. 在仿真的工程中由于正、余弦函數(shù)內(nèi)只對弧度值識別,在計算的忘記了導(dǎo)致出現(xiàn)的圖形出現(xiàn)的是錯誤的,后來
22、經(jīng)過老師指導(dǎo)將其改正輸出圖形正確。4.在編寫程序的過程中由于不懂原理而不會編寫,經(jīng)老師講解重新編寫。七、分析和總結(jié)在實驗時有好多不明白的地方,通過查看資料在用Matlab進(jìn)行仿真,能夠讓我們更清楚的理解每個課題的實驗原理,也讓我們體會到了Matlab的重要性,同時通過對程序進(jìn)行改動和觀察現(xiàn)象能夠讓我們更好的掌握光學(xué)原理及其現(xiàn)象。這次實驗也讓我深刻的體會到理論與實踐結(jié)合的重要性,我們不僅要學(xué),而且要學(xué)著把知識用到實踐中去,這才使學(xué)得有意義、有價值。而且我覺得不懂就要問,這次不懂,下次就懂了,一點一點的弄懂了也就會得多了才會有信心去做更難的題目。八、參考文獻(xiàn)Matlab講義及參考書物理光學(xué)與應(yīng)用光學(xué)(第二版)石順祥、王學(xué)恩、劉勁松編應(yīng)用光學(xué)(第四版)
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年旅游管理專業(yè)實務(wù)考卷及答案
- 學(xué)區(qū)房學(xué)位鎖定與購房保障合同
- 青春偶像劇改編影視作品授權(quán)合同
- 機場航站樓精裝修房驗收標(biāo)準(zhǔn)與保障協(xié)議
- 藝術(shù)考級項目品牌推廣與承包合作協(xié)議
- 留學(xué)歸國人員國際科技合作項目執(zhí)行聘用合同
- 零售企業(yè)智能倉儲物流系統(tǒng)開發(fā)與維護(hù)服務(wù)合同
- 網(wǎng)絡(luò)信息安全防護(hù)方案設(shè)計與實施協(xié)議
- 東南亞熱帶雨林度假村租賃及生態(tài)保護(hù)協(xié)議
- 跨國貿(mào)易進(jìn)口貨物違約金計算及結(jié)算協(xié)議
- 基于PLC交流變頻調(diào)速系統(tǒng)的設(shè)計 畢業(yè)設(shè)計(論文)
- DB11-490-2007地鐵工程監(jiān)控量測技術(shù)規(guī)程高清最新版
- PROE中英文對照表
- 漢字六書課件
- 勞動爭議仲裁調(diào)解法課件
- 職場利器-陽光心態(tài)課件
- 發(fā)展少數(shù)民族醫(yī)藥實施方案
- GB∕T 20394-2019 體育用人造草
- 中考物理總復(fù)習(xí)幾點建議和思考(ppt課件)
- 乳腺的影像診斷課件
- 三峽庫區(qū)榨菜廢水污染治理技術(shù)導(dǎo)則
評論
0/150
提交評論