幾種光柵。正弦光柵，矩形光柵，正弦光柵的強度衍射分布。

NO1

%正弦振幅型光柵

clear all;clc

close all

I=zeros(768,1024);

for i=1:768

for j=1:1024

I(i,j)=128+127*cos(i/32*2*pi+pi/4);%頻率為1/32

%或者I(i,j)=0.5+0.5*cos(2*pi/32+pi/4)

%if I(i,j)<128

% ? I(i,j)=0; %二值化處理

%else

% ? I(i,j)=1;

%end

end

end

I1=mat2gray(I); %實現(xiàn)圖像矩陣的歸一化操作

% '歸一化'是使矩陣的每個元素的值都在0和1之間

imshow(I1,[])

N02.矩形光柵

%矩形光柵

clear all;clc

%close all

I1=zeros(96,1024);

I2=ones(96,1024);%768/4/2=96

I=zeros(192,1024);

I(1:96,:)=I1;

I(97:192,:)=I2;

Iout=repmat(I,4,1);%橫向條紋

I01=zeros(768,128);

I02=ones(768,128);%1024/4/2=128

I0=zeros(768,256);

I0(:,1:128)=I01;

I0(:,129:256)=I02;

I0t=repmat(I0,1,4);%縱向

Itot=Iout+I0t;

%Itot=mod(Itot,2);

Itot=im2bw(Itot,1);% 二值化，1為閾值，小于等于1化為0，大于1化為1

%此句為橫向，縱向條紋‘黑’‘黑’相加

figure(5)

imshow(Itot)

NO.3

%一維正弦相位光柵的強度衍射分布

%傅里葉光學(xué)P110，(3-138)

clear all;clc

close all

lambda=650e-9;%波長m；

f0=8e6;%光柵頻率；

z=0.01;%入射距離m；

L=1e-3;%光柵長度m

m=8;%相位變化幅度

N=10;%貝塞爾函數(shù)的階數(shù)，取有限的階數(shù)

d=0.4;%觀察屏的觀察范圍m

x=-d:0.002:d;

I0=(L^2/lambda/z).^2;

I1=zeros(2*N+1,length(x));

for q=-N:N

for i=1:length(x)

I1(q+N+1,i)=besselj(q,m/2).*sinc(L/lambda/z*(x(i)-q*f0*lambda*z));

end

end

I2=I0.*sum(I1.^2);

%plot(x,I2)

I3=(I2-min(I2))./(max(I2)-min(I2));

figure;

plot(x,I3);%歸一化

I4=zeros(length(x));

for j=1:length(x)

I4(:,j)=I3(j);

end

figure

imshow(I4)

imwrite(I4,'intensity.jpg')

轉(zhuǎn)載本文請聯(lián)系原作者獲取授權(quán)，同時請注明本文來自宗兆玉科學(xué)網(wǎng)博客。

鏈接地址：http://blog.sciencenet.cn/blog-1120694-854529.html

上一篇：望遠鏡的性能簡介

下一篇：?基于zemax的反射式系統(tǒng)的結(jié)構(gòu)設(shè)計

總結(jié)

以上是生活随笔為你收集整理的matlab正弦光栅,科学网—光栅-正弦，矩形 - 宗兆玉的博文的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。