matlab模擬菲涅爾波帶片

1、4.菲涅爾波帶片（1）偶數(shù)波帶片：將波帶片所在屏幕分為1001*1001個(gè)點(diǎn)，依次求出各點(diǎn)處的半波帶數(shù)，為奇數(shù)則涂黑，為偶數(shù)則不涂黑表示透光。設(shè)波帶片中心為原點(diǎn)。clear; lam=600e-6; %設(shè)置波長為600nmR=3; %設(shè)置菲涅爾波帶片半徑為3mmf=1000; %設(shè)置焦距1mym=R;xm=R;y=linspace(-ym,ym,1001); x=linspace(-xm,xm,1010); %設(shè)置程序中x，y向量，將屏幕分為1001*1001個(gè)點(diǎn)for m=1:1001 for n=1:1001 %設(shè)置二重循環(huán)依次求菲涅爾波帶片屏幕上各點(diǎn) p=sqrt(x(m).2+y(n)

2、.2); %求各點(diǎn)所在圓半徑 k=fix(p.2./(lam.*f); %求各點(diǎn)半波點(diǎn)數(shù) if pR%如果屏幕上點(diǎn)大于波帶片半徑 I(m,n)=0;% 則用灰色表示背景，不是波帶片部分 else if mod(k,2)=1%判斷半波帶數(shù)奇偶，為奇數(shù)則涂黑 I(m,n)=0; else I(m,n)=1; %為偶數(shù)則不涂黑表示透光 end end endendcolormap(gray(255)N=255; %設(shè)置灰度等級Ir=I*N; image(Ir); %繪制波帶片圖像運(yùn)行程序得到如下圖像，為菲涅爾偶數(shù)波帶片：（2）奇數(shù)波帶片：將波帶片所在屏幕分為1001*1001個(gè)點(diǎn)，依次求出各點(diǎn)處的半

3、波帶數(shù)，為偶數(shù)則涂黑，為奇數(shù)則不涂黑表示透光。設(shè)波帶片中心為原點(diǎn)。clear; %設(shè)置波長為600nmlam=600e-6; %設(shè)置波長為600nmR=3; %設(shè)置菲涅爾波帶片半徑為3mmf=1000; %設(shè)置焦距1mym=R;xm=R;y=linspace(-ym,ym,1001); x=linspace(-xm,xm,1010); %設(shè)置程序中x，y向量，將屏幕分為1001*1001個(gè)點(diǎn)for m=1:1001 for n=1:1001 %設(shè)置二重循環(huán)依次求菲涅爾波帶片屏幕上各點(diǎn) p=sqrt(x(m).2+y(n).2); %求各點(diǎn)所在圓半徑 k=fix(p.2./(lam.*f); %求各點(diǎn)半波點(diǎn)數(shù) if pR%如果屏幕上點(diǎn)大于波帶片半徑 I(m,n)=0.5;% 則用灰色表示背景，不是波帶片部分 else if mod(k,2)=1%判斷半波帶數(shù)奇偶，為奇數(shù)則不圖黑表示透光 I(m,n)=1; else I(m,n)=0; %為偶數(shù)則涂黑 end end endendcolormap(gray(255)N=255;