3.function [y,t]=mcircle(r) t=0:2*pi/64:2*pi; y=r*ones(size(t));

subplot(121),polar(t,y,'*r') [X,Y]=po12cart(t,y); X=r*cos(t); Y=r*sin(t); subplot(122),plot(X,Y,'*r') axis equal; axis square;

4.t=0:pi/180:2*pi*5; r1=5;

x1=r1*cos(t); y1=r1*sin(t); z=t/(2*pi);

Subplot(121),plot3(x1,y1,z) grid on

r2=linspace(5,0,length(t)); x2=r2.*cos(t);y2=r2.*sin(t); subplot(122),plot3(x2,y2,z) grid on

9

实验小结：本次实验我们掌握二维图形的绘制，在老师的指导下掌握图形的

标注，同时了解三维曲线和曲面图形的绘制，收获颇丰！

10

实验四 离轴傅里叶全息术的光学仿真

实验

（一）实验目的

掌握离轴全息术的基本原理，掌握傅里叶全息图的仿真过程。 （二）实验内容

选择合适的原始图片，通过软件编写代码，仿真离轴傅里叶全息图的制作、再现过程。

程序和图像结果如下： clear;

cm=0.01;um=1e-6; N=1024;%全息图分辨率

% hx=1*cm;hy=1*cm;%原始图大小 % dhx=hx/N;dhy=hy/N;%原始图各方向增量 dhx=17*um;dhy=17*um; x=dhx.*(ones(N,1)*[-N/2:N/2-1]); y=(dhy.*(ones(N,1)*[-N/2:N/2-1]))'; p1=imread('kzz.bmp'); p1=rgb2gray(p1); p1=imresize(p1,[512,512]); p1=double(p1); p0=zeros(1024); p0(1:512,1:512)=p1; % D=zeros(1024);

% L=imread('E:\\动物\\dongwu2_033.jpg'); % figure(1);imshow(L); % p=rgb2gray(L); % %p=double(p); % niao=p(1:768,1:768); % niao=imresize(niao,[512,512]); % D(1:512,1:512)=niao; randnum=rand(N,N);

randphase=exp(i*2*pi*randnum); p=double(p0).*randphase; % p1=double(D);

11

%p=mat2gray(p);

figure(4);imshow(abs(p).^2);%原始图 %z=2;%菲涅尔传输距离 w=0.6328*um;%波长 k=2*pi/w;%波数

%[f,dx1,dy1,x1,y1]=fourier(p1,N,N,dhx,dhy,z,w); f=fftshift(fft2(fftshift(p))); sit=pi/30;%参考光角度

ur=exp(-j*k*y*sin(sit));%参考光的复波函数 f1=f+ur;

% f=f.*exp(i*2*pi*(x+y)); f1=abs(f1).^2;

f1=f1-abs(ur).^2-abs(f).^2; f3=mat2gray(f1); f3=histeq(f3);

% imwrite(f3,'C:\\Documents and Settings\\Administrator\\桌面\\carh.bmp'); %L=(ones(1024,1)*[1:1024])';

% comphase=2*17*sin(pi/18)/0.6328*2*pi*ones(1024,1)*[1:1024]; %comphase=0.66*L; ó=f3+comphase;

%imwrite(f3,'C:\\Documents and Settings\\ahu\\桌面\\fmaSHangXianc.bmp') figure(5);imshow(f3);%全息图 rp1=f1.*conj(ur);%再现光照射 % rp1=f1.*conj(ur); % rp1=f3.*abs(ur).^2;

%[f2,dx2,dy2,x2,y2]=fourier(rp1,N,N,dx1,dy1,z,w); f2=fftshift(ifft2(fftshift(rp1))); f2=abs(f2).^2; f2=mat2gray(f2);

imwrite(f2,'C:\\Documents and Settings\\Administrator\\桌面\\kzz.bmp'); figure(6);imshow(f2);

12