光纖matlab編程

1、光電綜合課程設(shè)計(jì)報(bào)告姓 名:李方圓 學(xué) 號(hào): 1150730006 專 業(yè): 應(yīng)用物理學(xué) 班 級(jí)：2011級(jí) 指導(dǎo)教師：李丹 實(shí)驗(yàn)二 光的電磁理論及仿真【實(shí)驗(yàn)?zāi)康摹?掌握Bessel函數(shù)及其特性。2掌握MATLAB的一維、二維曲線作圖，并對(duì)單模光纖的特性進(jìn)行分析。【實(shí)驗(yàn)內(nèi)容】一、用MATLAB做出前五階第一類Bessel函數(shù)曲線。分析：利用MATLAB內(nèi)建貝塞爾函數(shù)分別作出0階，1階，2階，3階，4階圖像【實(shí)驗(yàn)程序及相應(yīng)結(jié)果】程序：clear;clc;x=0:0.1:25;y0=besselj(0,x);y1=besselj(1,x);y2=besselj(2,x);y3=besselj(3,

2、x);y4=besselj(4,x);plot(x,y0,x,y1,x,y2,x,y3,x,y4);grid on;title(前五階第一類Bessel函數(shù)曲線)xlabel(x);ylabel(y);gtext(J0);gtext(J1);gtext(J2);gtext(J3);gtext(J4);二、在MATLAB中作圖，給出單模光纖中LP01模在V 分別為0.8、1.6和2.4時(shí)，電場(chǎng)分量E 相對(duì)歸一化直徑Ra的歸一化曲線。分析：給定V的條件下，知道了對(duì)應(yīng)U,W的值，故可以在MATLAB中用hold on分別作出R對(duì)應(yīng)E的曲線【實(shí)驗(yàn)程序及相應(yīng)結(jié)果】程序：clear all;V=0.8 1

3、.6 2.4;U=0.7974 1.3670 1.6453;W=0.0640 0.3815 1.7473;Ra1=(-1:0.1:1);Ra2=(1:0.1:10),(-10:0.1:-1);E1_1=besselj(0,0.7974*Ra1);E1_2=besselj(0,0.7974)*besselk(0,0.0640*Ra2)/besselk(0,0.0640);E2_1=besselj(0,1.3670*Ra1);E2_2=besselj(0,1.3670)*besselk(0,0.3815 *Ra2)/besselk(0,0.3815 );E3_1=besselj(0,1.6453*

4、Ra1);E3_2=besselj(0,1.6453)*besselk(0,1.7473*Ra2)/besselk(0,1.7473);plot(Ra1,E1_1,Ra2,E1_2,Ra1,E2_1,Ra2,E2_2,Ra1,E3_1,Ra2,E3_2);三、在MATLAB中進(jìn)行二維作圖，給出單模光纖中LP01模在V 為2.4時(shí)，歸一化光強(qiáng)I相對(duì)歸一化直徑Ra的二維分布圖。分析：根據(jù)光強(qiáng)公式作出其分布圖【實(shí)驗(yàn)程序及相應(yīng)結(jié)果】程序： clearclose all V=2.4000;U=1.6453;W=1.7473; Npoint=501; x=linspace(-5,5,Npoint);y=l

5、inspace(-5,5,Npoint);X=meshgrid(x,y);Y=meshgrid(y,x);Y=Y;R=sqrt(X.2+Y.2); I1=besselj(0,U*R).2;I2=(besselj(0,U).*besselk(0,W.*R)./besselk(0,W).2;I=I1;pos=find(R=1);I(pos)=I2(pos); imagesc(x,y,I,0 1);colormap(gray);colorbarxlabel(x)ylabel(y)zlabel(I)四、在MATLAB中進(jìn)行二維作圖，給出單模光纖中LP01模在V 為2.4時(shí)，歸一化光強(qiáng)I相對(duì)歸一化直徑Ra的三維分布圖。分析：根據(jù)光強(qiáng)公式作出其分布圖【實(shí)驗(yàn)程序及相應(yīng)結(jié)果】程序： clear all;r1=0:0.01:1;r2=1:0.01:5;theta=0:0.01:6.28;R1,th1=meshgrid(r1,theta);R2,th2=meshgrid(r2,theta);x1=R1.*cos(th1),y1=R1.*sin(th1);x2=R2.*cos(th2),y2=R2.*sin(th2);u=1.6453,w=1.7473;I1=besselj(0,u.*R1).*besselj(0,u.*R1);I2=besselj(0,u)*besselj(0