matlab大作业题题单,matlab大作业第一题

matlab大作业第一题

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)运行结果：,