马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
close all;clc;clear;<BR>A=zeros(64);<BR>A(15:20,10:50)=1;<BR>A(20:50,30:35)=1;<BR>figure;imshow(A);<BR>B=fft2(fftshift(A));<BR>C=fftshift(B);<BR>Am=abs(C);<BR>Ph=angle(C);<BR>th=max(max(Am));<BR>Am=Am/th;<BR>Ph=Ph/(pi*2);<BR>cgh=zeros(64*13);<BR>hol=zeros(13);<BR>for m=1:64;<BR> for n=1:64;<BR>if abs(Ph(m,n))<=1/4;<BR>h=round(Am(m,n)*12);<BR>p=round(abs(Ph(m,n))*12);<BR>k=ones(h+1,7);<BR>hol(round(7-h/2):round(7+h/2),4+p:10+p)=k;<BR>else <BR> hol(round(7-h/2):round(7+h/2),4+p:13)=ones(h+1,10-p);<BR> hol(round(7-h/2):round(7+h/2),1:p-2)=ones(h+1,p-2);<BR> end;<BR>cgh((m-1)*13+1:m*13,(n-1)*13+1:n*13)=hol;<BR>end;<BR>end;<BR>figure;imshow(cgh,[]);<BR>Re=ifft2(cgh);<BR>Re=fftshift(Re);<BR>figure;imshow(abs(Re),[]); |