matlab大作业第一题.pptVIP

Matlab大作业第一题 光栅衍射的强度分布函数为 式中，相邻两缝光线的位相差 ； 是光栅常数； 为衍射角。每个单缝衍射的光强可写为 式中， ；a 是缝宽； 是每条缝的光强。试编写一计算机程序，画出8缝光栅（N 8）衍射的相对光强分布图（即 曲线），并观察单缝衍射对多缝干涉的调制作用（设缝宽 ， 光栅常数 ，波长 ）。 （1）分析过程： 由题意求得：相对光强为 式中， 是单缝衍射因子， 是多缝干涉因子。 从此式可知，光栅衍射光强是单缝衍射和多缝干涉因子的乘积，所以 I/I0 为多缝干涉相对光强。 因此，只需画出： （a）画出光栅衍射的相对光强分布图。 （b）画出多缝干涉的相对光强分布图。 这样便能观察单缝衍射对多缝干涉的调制作用。 2 在matlab中用到的函数： a、 linspace x1,x2,N linspace是Matlab中的一个指令，用于产生x1,x2之间的N点行矢量。其中x1、x2、N分别为起始值、终止值、元素个数。 b、 subplot（m,n,p） c、plot（x，y）绘制二维图形。 d、title（图形名称） e、xlabel x轴说明 ，ylable（y轴说明）。 3 程序： clear a 4e-6; d 4e-5; lambda 5e-7;theta linspace -0.01*pi,0.01*pi,1000000 ; u pi*a*sin theta /lambda; v pi*d*sin theta /lambda; I sin u ./u .^2 .* sin 8*v ./sin v .^2 ; I1 sin 8*v ./sin v .^2; subplot 2,1,1 ; plot sin theta ,I ; xlabel sin theta ;ylabel I/I0 ; title N 8时的光栅衍射相对光强分布图 ; subplot 2,1,2 ; plot sin theta ,I1 ; xlabel sin theta ;ylabel I/I0 ; title N 8时的多缝干涉相对光强分布图 ; 4 运行结果：