用matlab研究光栅衍射

clear all
lamd=650e-9;
f=5;
Ns=900;
str{1}=‘2缝d/a=5‘;
str{2}=‘5缝d/a=4‘;
str{3}=‘8缝d/a=3‘;
M=[258];
a=2e-4;
d=[5*a4*a3*a];
ym=3*lamd*f/a;
ys=linspace（-ymymNs）;
for j=1:3
  for i=1:Ns
      sinTHE=ys（i）/f;
      phi=2*pi*a*sinTHE/（lamd*Ns）;
      num=sin（Ns*phi/2）;
      deno=sin（phi/2）;
      v=pi*d（j）*sinTHE