原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
光纤光学大作业
T903-220090030231
程序如下：
%有限差分光束传输法(FDBPM)
%以微米为单位
%常数的定义
clc,clear
lambda=1.55;%波长
k0=2*pi/lambda;
neff=(1.4457+1.4378)/2;
a1=k0*neff;
delta_x=0.1;delta_y=0.1;delta_z=2*i;
N=300;M=300;%lengthisN,widthisM;
x0=N/2;y0=M/2;
%定义网格各点的折射率
n=1.4378*ones(N,M);
fori1=1:N
forj=1:M
if(i1-x0)^2+(j-y0)^2=2500
n(i1,j)=1.4457;
end
end
end
%初始光场
A=1.2;
[x,y]=meshgrid(1:N,1:M);
fai0=exp(-(x-x0).^2/100/A).*exp(-(y-y0).^2/100/A);
mesh(0.1*x,0.1*y,fai0)
title(初始光场);xlabel(x);ylabel(y);
Z=400;
fai=zeros(N,M,Z)+eps;
fai(:,:,1)=fai0;
P=2-k0^2*delta_x^2/2*(n.^2-neff^2)+4*i*a1*delta_x^2/delta_z;
Q=-2+k0^2*delta_x^2/2*(n.^2-neff^2)+4*i*a1*delta_x^2/delta_z;
%%%%%%%%%
kl=i/delta_x*log(fai(:,1,1)./fai(:,2,1));
k_left=abs(real(kl))+i*imag(kl);
kr=i/delta_x*log(fai(:,N,1)./fai(:,N-1,1));
k_right=abs(real(kr))+i*imag(kr);
%%%%%%%%%%%
a=[zeros(N,1),-1*ones(N,N-3)];
光纤光学大作业T903‑2宋志
鹏20090030231
程序如下：
％有限差分光束传输法（fdbpm）％以微米为单位％
常数的定义clc，clearlambda=1.55;％波长k0=2
*pi/lambda;neff=（1。4457+1.4378）/2;a1=
___
k0*neff;deltax=0.1;deltay=0.1;deltaz=2
*i;n=300;m=300;百分比is，宽度为m；x0=n/2;y
0=w/2;％定义网格各点的折射率n=1.4378*一个
（n，m）;forI1=1：nforj=1：mif（i1‑x0）^2+
（j‑y0）^2=2500n（i1，j）=1.4457;端端％初始光
场a=1.2;[x，y]=meshgrid（1：n，1：m）;fai0=
exp（‑（x‑x0）。^2/100/a）。*exp（‑（y‑y0）。
^2/100/a）;网格（0.1*x，0.1*y，fai0）标题（初始
光场）;Xlabel（x）;Ylabel（y）;z=400;fai=
zeros（n，m，z）+eps;fai（：，：，1）=fai0;
_
p=2‑k0^2*deltax^2/2*（n。^2‑nff^2）+4*i*a1
文档评论（0）