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