全息图matlab程序.doc_第1页
全息图matlab程序.doc_第2页
全息图matlab程序.doc_第3页
全息图matlab程序.doc_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

%Modified off-axis reference beam CGH of Burch% %读入原始图象 PI=3.14159; Image=imread(lena256.jpg); N=512; A=zeros(N,N); B=zeros(N,N); for I=1:1:256 for J=1:1:256 ImageNum=double(Image(I,J,1); A(I,J)=ImageNum/255; B(I,J)=0; end;end; figure; imshow(A); %为降低全息图的动态范围，乘一随机位相因子% for I=1:1:N for J=1:1:N R=rand(1,1); B(I,J)=A(I,J)*sin(R*2*PI); A(I,J)=A(I,J)*cos(R*2*PI); F(I,J)=A(I,J)+j*B(I,J); end;end; %为降低全息图的动态范围，乘一随机位相因子% %对物函数做二维FFT变换% F=fft2(F); Max=max(max(abs(F); F=F/Max; A=real(F); B=imag(F); %对物函数做二维FFT变换% %定义载波参数 alpha=0.5; for I=1:1:N for J=1:1:N Xcos=(J-1)/127; A1(I,J)=cos(2*PI*alpha*Xcos); B1(I,J)=sin(2*PI*alpha*Xcos); end;end; %定义载波参数 %全息图数据区 for I=1:1:N for J=1:1:N Holodata(I,J)=0.5+0.5*(A(I,J)*A1(I,J)+B(I,J)*B1(I,J); end;end; %全息图数据区 %绘制全息图，制作灰阶全息图像，灰度等级256 M=512;%定义全息图的大小， Hologram=zeros(M,M); S=M/N;%定义每个抽样单元大小为S,S for I=1:1:N for J=1:1:N Xa=(J-1)*S+1; Xb=J*S; Ya=(I-1)*S+1; Yb=I*S; for Ix=Xa:1:Xb for Iy=Ya:1:Yb Hologram(Iy,Ix)=Holodata(I,J); end;end; end;end; Max=max(max(Hologram); Hologram=Hologram/Max; figure; imshow(Hologram); %绘制全息图，制作灰阶全息图像，灰度等级256 %再现 Object=fft2(Hologram); Object=fftshift(Object); Object=abs(Object);
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论