南邮光电综合设计报告2

《南邮光电综合设计报告2》由会员分享，可在线阅读，更多相关《南邮光电综合设计报告2（17页珍藏版）》请在装配图网上搜索。

1、课程设计报告B01060702邓心惟A类：课题2.设计任务：有一平面镜和曲率半径为 R的凹面镜，画出光束发散角与腔长 一.课题要求：1.有输入输出界面；2.可输入不同凹面镜曲率半径值，查看结果。参考：激光原理第二章L的关系曲线。二.课题分析及设计思路：1 .问题分析：根据激光原理，一般稳定球面腔基模远场发散角为： 2(2L - 尺 - &)24胡-3U2L(Ri-L)(R-L)(R+R2-L) =2gi g2gig22l gig2(igg2)14其中L为腔长，R为半径又由课题条件，为平-凹腔，故r1 = ：可简化公式为:其中变量为L, R2 ,，常量为二2设计思路：要求发散角与腔长 L的的关系

2、，即需要给定 & , 的值，R2要求根据输入确定，而 可以为程序内含或者外部输入。考虑到一般激光器的波长不是任意值，而为了使设计单元体现不同波长对发散角的影 响，这里利用分离选择项作为波长输入，一般我们使用的激光器为CO激光器，波长10.6um；氦氖激光器，波长 632.8nm。因此选择支为两个。采用MATLAB用户界面(GUI)工具设计，输入变量为两个，R为编辑输入，单位cm,缺省值1m；波长为选择输入，10.6um或者632.8nm，缺省值为10.6um。三.模型创建与编程：本题较为直观，除了在 MATLAB计中注意矩阵元与数值的差别外，没有难点。 直接给出相应处理的程序部分：global

3、 r;r=str2double(get(hObject,Stri ng);%r为半径变量，从界面处获得输入半径值global bocha ng;contents = get(hObject,Stri ng);bocha ng=str2double(c onten tsget(hObject,Value);%bochang为波长变量，从界面处获得输入波长值，因为设计实现时直接选择的就%是波长值，因此直接将其转化成双精度值即可。也可以选择激光器类型，用字符%串变量接收然后再在程序中转换。l=0.01*r:0.99*r;猥使平凹腔稳定，腔长必须为比半径小的值，设定曲线横坐标范围an g=2*sqrt

4、(sqrt(bocha ngA2./(piA2.*(r-l).*l);9发散角函数，注意矩阵量的运算axes(ha ndles.axesl);cla;plot(l,a ng);%乍图四仿真调试与结果分析：1 仿真调试：模拟确省值 R=1m九=10.6um,改变半径值 R=2m入不变 输入后输出图形分别为：腔长保持100cm,改变输入波长，输出图形:2 结果分析：首先可以看到，不同输入半径，不同波长，输出 曲线形状大体一致，都是凹形曲线，存在一个极小值 的点，在腔长为半径一半处，此时球面发射镜焦点落 在平面镜与轴交点，即稳定平凹腔的等价共焦腔的一 球面与R重合，获得最小发散角。不同R的值影响发散

5、角，R越大，发散角越小同样，不同波长值也影响发散角，波长越小，发散角越小 这一点可以通过函数理论得出：1/42L(R2 -L)发散角二是R,，的单调函数；当L=R-L时，即L=R/2时，分母获得最大，相应的发散角最小。课题4.设计任务：实验测得半导体发光二极管结偏压与注入电流的一组数据为:V(v)1.002.002.012.022.032.04Ij(A)0.0000.0490.0730.1100.1690.265V(v)2.052.062.072.082.092.10Ij(A)0.4220.6841.1261.8753.1575.360试根据该组数据，建立适当的数学模型， 可根据结电压计算相应

6、注入电流，并画出发光二极管的电流一电压特性曲线。一. 课题要求：（1）具有输入输出界面；调节输入值。得到不同结果。（2 ）画出曲线，并分析电流 一电压特性。 参考：光电子器件与 OEIC模拟第二章、第三章二. 课题分析及设计思路：本题中给出的电压值均为正，因此我们要讨论的发光二极管的电压电流特性也是正向偏压时的情况。根据所给的12个点数据，需要做出曲线，并根据曲线读值。这种情况可以通过插值拟和作出相应曲线，MATLAB提供了丰富的插值函数，实现中采用一维插值函数in terl进行拟和，再根据拟和后的曲线读出不同偏压下的电流值。实现的关键即为拟和。我们注意到所给的12个点中，除第一个点外均匀分布

7、与x轴，而第一第二个点跨距大，通过实践看到，如果图形中包含第一个点，会在1,2区间产生极大的抖动，不能反映整体的趋势故使横坐标从靠近 2开始以避免出现这样的干扰，图中不显示（1, 0）点，但是可以通过曲线趋势看出。考虑到输出时的问题，由于在1 , 2区段电流值非常小， 故可考虑用线形简单拟和原曲线，而获得在1 , 2区段较为可信的注入电流值。三. 模型创建与编程：问题实现分为两步：1. 曲线拟和输入12个值的数据放入矩阵变量中，利用一维插值函数拟和后显示。实现函数部分为：global v;global Vj;global Ij;v=1;%输入确省的电压值Vj=1.00 2.00 2.01 2.

8、02 2.03 2.04 2.05 2.06 2.07 2.08 2.09 2.10;Ij=0.000 0.049 0.073 0.110 0.169 0.265 0.422 0.684 1.126 1.875 3.157 5.306; %输入要插值的数据x=1:0.01:2.12;% 插值横坐标范围 y=interp1(Vj,Ij,x,spline);%三次样条插值axes(handles.axes2);cla;plot(Vj,Ij,*,x,y,m),axis(1.99 2.12 0 8); %显示图形，包括原来输入的点和拟和后的曲线，定横纵轴范围，避免出现抖动 grid on;% 网格显示

9、2 输入偏压值输出注入电流一个输入量U, 个输出量I实现函数部分为：global v; v=str2double(get(hObject,String);%v为接收输入偏压值的变量if (v=1) result=0.049*(v-1);elseresult=interp1(Vj,Ij,v,spline);end;%result为根据拟和确定的相应偏压v的注入电流，当v在区间1 , 2时，线形拟和输出值set(handles.edit7,String,num2str(result);%将电流值显示在界面输出框中axes(handles.axes2);cla; plot(Vj，Ij，x，y，m，v

10、，result，ob)，axis(1.99 2.12 0 8);grid on;if (v=1) legend(The first point (1.00， 0.000) is omitted.)else if (v1.0) legend(电流随电压变化较小 )else if (v2.05) legend(电流随电压变化很大 )else legend(invalid，0)end; end; end;%根据不同输入情况显示不同的特性四仿真调试与结果分析：1 输出的二极管电流 -电压曲线 2 输入不同的偏压，计算注入电流：在界面中输入框中输入U=1，输出框显示为0；输入U=1.4，显示0.016;

11、输入 U=2.09,显示 3.157 ;3.结果分析：该拟和图验证了二极管典型的伏安特性，正向偏压时只有超过某一数值时才有明显的 电流，该电压称为导通电压，反映在具体数据可以看做为2V,超过2V后，电流随电压迅速增长。在图形显示区域内， 呈现出类似指数变化规律；在图形末段已显示出近似线形的变化规律。B类：课题7.光电子器件参量计算：(1)设计任务：计算光电探测器的光电流。一. 设计要求：1.具有输入输出界面；2. 参量可任意输入；3. 模拟输入一组数据(数据值应与实际相当)，给出结果。参考：课本P356二. 课题分析：该设计单元进行数值计算，不涉及图象，故可编辑相应数量的输入框以及相应数量的输

12、 出框即可。光电探测器光电流的计算涉及11个变量，分别为：二极管横截面积 A, P区掺杂浓度Na, n区掺杂浓度Nd电子扩散系数Dn,空穴扩散系数Dp,少数电子载流子寿命.n反向偏压V,温度Tp-n结基质，影响本征载流子浓度计算中得出四个中间结果，为： 电子扩散长度Ln，空穴扩散长度Lp，内建电压Vbi,耗尽层宽度W最终结果：光电流I其中有判断过程，当计算出的Ln和Lp与W相比很小时可将光电流作为瞬时电流，在此不作判断，得出精确结果。相应的公式为: = , Dn n ；Vb也ln（罂）e m ）（VbiNaNdVr）1/2I 二 eAGL（W Ln Lp）三设计思路：T不同，将影响本征载流子浓

13、度以及内建电压的值，但为了简化问题，这里只计算 温度为300K时的情况，即常温下的情况。考虑到不同基质将有不同的本征载流子浓度，而GaAs也是重要的探测器物质，因此设计单元选择不同物质来获得相应的载流子浓度，给出Si,Ge.GaAs三种常见基质。并给每个输入变量设定相应的缺省值。四模型创建与编程：该问题只涉及计算且较为简单，故直接给出部分函数实现: a=1.0e4;na=2.0e16;nd=1.0e16;dn=20;dp=12;tn=1.0e-8;tp=1.0e-8;gl=1.0e22;v=2;ni=0;液定各参数缺省值global t;contents = get(hObject,Stri

14、ng);t=co nte ntsget(hObject,Value);if (t=Ge) ni=2.4e13else if (t=Si) ni=1.5e10else if (t=As) ni=1.1e7else ni=0en d;e nd;e nd;%从下拉框中获得被选择的基质，转换成响应本征浓度ln=sqrt(d n*t n)*1.0e4;lp=sqrt(dp*tp)*1.Oe4;vbi=0.026*log( na*n d/niA2.);w=sqrt(1.32e7*( na+nd)/na/nd)*(vbi+v)*1.0e4;i=(1.6e-19)*a*(1.0e-12)*gl*(w+l n+

15、lp)*1.0e3;%各步计算set(ha ndles.edit10,Stri ng,n um2str(l n);set(ha ndles.edit11,Stri ng,n um2str(lp);set(ha ndles.edit12,Stri ng, nu m2str(w);set(ha ndles.edit13,Stri ng, nu m2str(vbi);set(ha ndles.edit14,Stri ng, nu m2str(i);en d;%在输出框中输出相应结果五仿真调试与结果分析：(1)选择Si及所有缺省值，得到Ln=4.4721umLp=3.4641umW=0.73324um

16、Vbi=0.71534Vl=0.13871mA选择GeLn=4.4721umLp=3.4641umW=0.67947umVbi=0.3317VI=0.13785mA选择GeAsLn=4.4721umLp=3.4641umW=0.78227umVbi=1.0907VI=0.1395mA可以看到，在这组条件的情况下，虽然本征载流子浓度差别引起内建电压的较大不同，由于耗尽层相比扩散区宽度很小，因此差别不大，且光电流很小(2)缩短载流子寿命，令其为 0.1 ns;得到三个长度(Si):Ln=0.44721um Lp=0.34641umW=0.73324um 已经在同一数量级上但由于未增大 W光电流减小

17、18-318-3 增大掺杂浓度：令 Na=2*10 cm ,Nd=10 cm结果(Si) Ln=4.4721umLp=3.4641umW=0.076489um Vbi=0.95481Vl=0.1282mA可以看到，增大浓度虽然增加了内建电压，但会很大幅度的减小耗尽层宽度，相应的会减小光电流增大反向偏压，Si的情况下，令 V=4V其他条件不变，W=0.96625um,l=0.14244 因此增大偏压可以增加一定的光电流，但并不十分显著。(5)很明显，加大截面积 A和提高产生率 G将直接增大光电流，每提高一个数量级都将直接提高光电流一个数量级是最有效的提高光电流的方法。扩散系数一般为较固定的值，因

18、 此模拟中就尽量保持其不变，同时 通过不同情况下的模拟发现，不同 基质虽然具有本征浓度几个数量级 上的差别，但在各种情况中相差很 小，只有在掺杂很低的时候反映出 一些差别，但在实际应用的条件下， 光生电流的差别非常小。主要原因 在于本征载流子浓度数量级上的差 别反映到内建电场时已转化为系数 差别，且通常小于外加电压。diode area A光电流的ll翼ctiicuicrtectiicuicrten-side doping HdE 1 Del6cm-3Bledron diHusionOn |20hole difkisiDncoeHidgnt Dp1 12cmZ/Bslectran minofi

19、t/C4rn*rhM1ii11* th1.0e-Bthole minorrtyearner Melime Ip- |L0e-6pide doping | 2 Del6M-3 |G&RESULT:electron ditfu&io niI鑑吨巾 Ln | -4721 umhale diHu&tDn length Lpa | 3 肠即umdeplelionwiditiW- | 0797umhuih-ufi voltage Vbi- | 0.3J17 VFlertran_hol0 pair G1 | 1 OH 3 Cffl郭*pholocwrent I- I 1mAreverse hiai V- |

20、(2)设计任务：计算LED的调制带宽。.课题要求：1.具有输入输出界面；2参量可任意输入；3.参见课本P470,输入并根据p452的图9-6中数据，由n/ r求出复合时间 r，再求调制带宽4. 模拟输入两组数据(数据值应与实际相当)，给出结果。课题分析：r与注入浓度n (=p)该问题涉及读图取值，由于未找到相应的直接函数反映复合时间 的关系，因此无法精确作出该函数图，只能采用先取值后拟和的方法作出近似图。产生拟和曲线后，由于计算出的已知量为横纵坐标的商，因此并不是直接用函数 in terpl读值就可以，涉及满足条件的点的寻找定位问题。找出该点后，已知r，按照公式计算，即可以获得调制带宽。由于9

21、-6的图为300K时的GaAs半导体，因此本单元的应用也限于300K GaAs,共有两个变量，注入电流密度 J，有源区宽度d;输出变量为3db截止频率。这里假设时间响应主要受电子空穴复合时间限制且效率为常数。公式为：end J =red r三设计思路：该问题有两个思考点。1. 对于图的取值拟和看到图9-6为双对数坐标，简单拟和无法产生原来图形，首先考虑的是先将其横纵坐标映射到等距坐标，即取对数，再简单模拟即可。当然，MATLAB提供了 loglog函数绘制双对数坐标的二维图，但原理是相同的。但这里涉及对于读值的两种理解。其一为简单读值，也是我们习惯的读法。在一格中，例如1016, 1017，认

22、为其间一半处为 5*1016,但左右半边不对称，左边为 1-5 （单 位1016）,间距为4个单位，右边为5-10，间距为5个单位，那么必然左边与右边映射 误差就不相同，进一步造成计算误差。其二为对数读值，即在上面的例子中，认为其间一半处为1016 1这样虽然不如其一直观，但是满足坐标的对数关系，且使映射达到统 一。这里考虑采用第二种，也正因为这个差别，可能造成结果上在同一数量级中的差别。2. 满足条件r的寻找定位。由于已知值将是横纵坐标的商，虽然不能直观的用interp1函数读值，但是由于该问题中商是横坐标的单调递增函数，因此问题得到适当简化。最直接的是采用遍历搜索，但由于步长的不同，相应的

23、精度设定也要合理。为了提高效率，考虑运用二分搜索，逐步逼近，也要注意对于不同要求精度的设定。四.模型创建与编程：问题实现分为两步：1. 取值拟和，不需界面输入，部分函数实现为：global nj;global tj;nj=14 14.5 15 15.5 16 16.5 17 17.5 18 18.5 19;tj=-5.5 -6.1 -6.69 -7.19 -7.62 -8.00 -8.31 -8.55 -8.74 -8.89 -9;%按 9-6图取值按取对数后输入，横坐标间距为 0.5ni=14:0.1:19;ti=i nterp1(nj,tj, ni,spli ne);%以和axes(ha

24、ndles.axes1);cla;plot( ni ,ti,m,nj,tj,p),xlabel( n=p(cm-3),ylabel(tr(s),axis(14 -5),lege nd(semico nductor GaAs at 300K);grid on;%显示2. 定位计算，部分函数实现为：global temp;global temp2;temp=j/(1.6e-19)/d*1.0e4;%temp变量保存根据输入计算的一=low=14;high=19;flag=1;while (high-low)0.1)&( flag=1)x=(low+high)/2;y=i nterp1(nj,tj,

25、x,spli ne);temp2=10A(x-y);if temp=temp2 flag=0;elseif (temptemp2)high=x;else low=x;en d;en d;n=10Ax;tr=10Ay;炉分搜索，判断条件为值相等的点找到，或者寻找范围缩小在 set(ha ndles.edit3,Stri ng, nu m2str( n);set(ha ndles.edit5,Stri ng,n um2str(tr);fc=1/2/pi/tr/1000000;set(ha ndles.edit4,Stri ng,n um2str(fc);殆输出框输出ni=14:0.1:19;ti=

26、i nterp1(nj,tj, ni,spli ne);axes(ha ndles.axes2);cla;plot(x,y,ob,ni,ti,m),axis(14 19 -10 -5),legend(interpl plot); grid on;19-100.1以内%显示定位的点五仿真调试与结果分析:1.拟和图形：因而可以作出对数坐标下两者的线形关系，但这里就省2.模拟输入(1)输入 J=100A/cm2d=1.0um结果 n=5.424691e+016r =7.0032e-009fc=22.7259MHz输入 J=500A/cm2结果 n=1.596339e+017.r =4.4498e-0

27、09fc=35.7664MHz输入 J=1000A/cm2,fc=43.1512MHz 输入 J=100A/cm：d=0.1um 结果 fc=43.1512MHz2 、J=500A/cm ,d 不变结果 fc=67.7026MHz2 亠J=1000A/cm ,d 不变结果 fc=87.7234MHz(见图)可以看到，fc与J/d成正比, 略不作了。C类：课题3设计任务： 对图所示的平行光入射的光学系统追迹计算：已知： UO=0 h=10mm 对输入的 R1， R2， R3，R4， d1, d 2, d 3, n 1, n 2值， 计算： (1) 光线 1 在各光学面的入射高度 hi ，像方孔径

28、角 ui ；(2) 像面距最后一面的距离 l 4一课题要求： 1. 有输入输出界面；2. 可根据输入的不同参数值，查看结果 参考：应用光学二课题分析： 对于多透镜组合共轴球面系统直接的方法就是利用光路折射计算公式以及转面公式， 顺 序计算，求得每个参量。 对于该近轴近似问题，可以使用近轴光路计算公式。但是为了精确 计算，仍然利用实际光路计算公式：sinI=(L-r)/r*sinUsinI =n/n sinI U=U+I-I L=r+r*sinI /sinU 转面公式：U2=U1L2=Li -di对于题给具体问题：由于 U0=0，h=10mm；初始时 sinI=h/r然后按上组公式计算。为了使上

29、述公式可以运用于不同光线和球面几何位置，需要注意符号的设定。 涉及该题的符号规则：1 线段以左至右为正，下至上为正。 截距为球面顶点计算到光线与光轴交点，半径为球面顶点计算到球心。为了符合用户习惯，输入时均认为其为正，在程序中进行转换，R2, R3为负，R , R4为正。2 角度顺时针为正，逆时针为负U,U为由光轴转到光线I,I 为由光线转到法线三设计思路： 1由于不同的输入初始值将获得不同坐标轴上的图形，虽然它不影响追迹计算，但是 为了最大程度的模拟实际问题，根据输入画出模拟图，标出相应标签。 将第一个透镜中心定位在坐标( 1， 0)，其余根据输入变化。输入变量为 R1， R2， R3， R

30、4， d1, d 2, d 3, n 1, n 2，画图用到前七个变量，透镜上的 点以圆方程定位。2追击计算，计算中已经考虑到输入半径的正负值，结果也按照符号规定来理解。 需要注意的时对于入射高度的计算，采用了近似公式h=l*tan(u); 而最精确的计算应为圆方程与直线方程的交点，但为了简化程序，就使用近似值。然后根据计算值画出追迹光路图。四模型创建与编程： 设计模型分为三个模块： 1输入模块：输入九个可变参数，并根据输入画出初始模拟图 2输出模块：输出八个结果3追迹模块：画出追击光路2 中特别考虑。实现中考虑到有折射率为 1 时的特殊情况，为了也显示该特殊情况，在模块相应实现部分程序：1

31、模块 1x1=(4-d1/2):0.01:4;y1=sqrt(r1A2-(x1-r1-4+d1/2).A2.);y2=-(sqrt(r1A2-(x1-r1-4+d1/2).A2.);x2=4:0.01:(4+d1/2);y3=sqrt(r2A2-(x2-4-d1/2+r2).A2.);y4=-(sqrt(r2A2-(x2-4-d1/2+r2).A2.);%第一透镜设置x3=(4+d1/2+d2-0.5):0.01:(4+d1/2+d2); y5=sqrt(r3A2-(x3-4-d1/2-d2+r3).A2.);y6=-(sqrt(r3A2-(x3-4-d1/2-d2+r3).A2.); x4=

32、(4+d1/2+d2+d3):0.01:(4+d1/2+d2+d3+0.5); y7=sqrt(r4A2-(x4-4-d1/2-d2-d3-r4).A2.);y8=-(sqrt(r4A2-(x4-4-d1/2-d2-d3-r4).A2.);%第二透镜设置axes(handles.axes1);cla;plot(x1,y1,x1,y2,b,x2,y3,b,x2,y4,b,x3,y5,b,x3,y6,b,x4,y7,b,x4,y8,b ),axis(0 14 -2 2),legend(光学面 1,2,3,4);grid on;%作出两透镜text(3.8,0,n1 ,fontsize,12);te

33、xt(4+d1/2+d2,0,n2 ,fontsize,12);text(3.8,sqrt(r1A2-(3.8-r1-4+d1/2)A2),R1rightarrow,HorizontalAlignment,right,fontsize,12);text(4.2,sqrt(r2A2-(4.2-4-d1/2+r2)A2),leftarrow -R2 ,fontsize,12); text(4+d1/2+d2-0.2,sqrt(r3A2-(4+d1/2+d2-0.2-4-d1/2-d2+r3)A2),-R3 rightarrow,HorizontalAlignment,right,fontsize,

34、12);text(4+d1/2+d2+d3+0.2,sqrt(r4A2-(4+d1/2+d2+d3+0.2-4-d1/2-d2-d3-r4)A2),leftarr ow R4 ,fontsize,12);text(4-d1/2,-0.5,leftarrowd1rightarrow); text(4+d1/2+d2,-0.5,d3);text(4+d1/2,-0.5,leftarrowd2rightarrow);%根据输入值做出标签2 模块 2 i10=asin(1/r1); i1=asin(1/r1/n1); u1=i10-i1;l1=r1+r1*sin(i1)/sin(u1); l20=l1

35、-d1;h1=l20*tan(u1); i20=asin(-(l20+r2)/r2*sin(u1); i2=asin(n1*sin(i20);u2=u1+i20-i2; l2=-r2-r2*sin(i2)/sin(u2);l30=l2-d2; h2=l30*tan(u2);i30=asin(-(l30+r3)/r3*sin(u2); i3=asin(sin(i30)/n2);u3=u2+i30-i3; l3=-r3-r3*sin(i3)/sin(u3);l40=l3-d3; h3=l40*tan(u3); i40=asin(l40-r4)/r4*sin(u3); i4=asin(n2*sin(

36、i40);u4=u3+i40-i4; l4=r4+r4*sin(i4)/sin(u4);%根据实际光路公式计算各变量 if u1=0h1=1;h2=1;h3=1;u2=0;u3=0;u4=0;end;%当u1=0时，会出现分母为0的情况，特别设定其值set(handles.edit10,String,num2str(h1);set(handles.edit11,String,num2str(h2);set(handles.edit12,String,num2str(h3);set(handles.edit13,String,num2str(u1);set(handles.edit14,Stri

37、ng,num2str(u2);set(handles.edit15,String,num2str(u3);set(handles.edit16,String,num2str(u4);set(handles.edit17,String,num2str(l4);%输出到输出框3 模块 3temp1=r1+4-d1/2-sqrt(r1A2-1);temp2=temp1+(1-h1)/tan(u1);temp3=temp2+(h1-h2)/tan(u2);temp4=temp3+(h2-h3)/tan(u3);temp5=h3-(11.5-temp4)*ta n(u4);% 十算各段光路横坐标axes

38、(ha ndles.axesl);cla;if u1=O plot(0.5,11.5,1,1);grid on;hold on;%寸于介质折射率与空气相同的情况，特别作图elseplot(0.5,temp1,1,1,temp1,temp2,1,h1,temp2,temp3,h1,h2,temp3,temp4 ,h2,h3,temp4,11.5,h3,temp5),legend(空气,介质 1,空气,介质 2,空气);grid on;hold on;en d;9追迹光路生成，加上标签五仿真调试与结果分析：-.1 HLaniW3K!WOHmm mm j近SB北將追遵计尊EM： X h*1m模块1的

39、实现结果（根据不同输入的输出图形）输入参数 R=3.5 R 2=3.5 R 3=3.5 R 4=3.5 d 1=1 d 2=3.2 d 3=0.6 n 1=1.55（玻璃）n 2=1.55 模块二及模块三的实现结果：（使用以上参数组）WI0MT1WIU貧背辭特殊情况下n=1时EFUUH近ifi弋囲追或计亶严| DiW I一dWJ HIDUTUZO五帝记帘追越计算M： UaO IrinpRSoil R jmjl pr 尹E3ID12jiHSi沁1ril7丁 丁I 1/Slai_I o追迹光路有仿真中发现，当输入偏离常用近似值的时候，由于所采用的高度近似公式,一点的误差，但是在一般的实际运用中，误

40、差很小，可以接受。改变折射率或者改变半径都将使光线发生较大偏折，改变距离未看到明显改变。课程设计小结：由于对于MATLAB各种函数功能不是非常熟悉，在编程中不断学习发现有直接的函数可 以实现操作，例如对数坐标的显示， 考虑用了先取对数再处理的方式，但是实际上有Loglog函数可以直接实现，较为直观。设计中尽量采用精确算法及定位， 以期获得最为精确的结果。但是在最后一个课题实现 时，采用了一个近似公式，理论上用解方程组可以获得精确值， 但需要在外面编辑一个函数。 仿真结果可见这样的近似在合理参数值情况下是可以接受的。课题实现中涉及插值拟合的过程，因而有精度上的差别，在使用二分搜索的过程中也涉及精度的控制，可以通过改变步长达到更高的精度。且对于拟和数值的不均匀所造成的拟和抖动，考虑外加线形逼近模拟抖动区函数值。通过这次课程设计，很大程度上增加了对与MATLAB用户界面的编辑的熟悉程度，且在对不同特殊条件的考虑分析中，进一步加深对于原来课程的理解，收获良多。