| 查看: 1029 | 回复: 0 | ||
[求助]
关于角谱法实现数字全息
|
|
通过角谱法实现数字全息，程序语法没有错误，但就是重现出来的图是黑色一片，不知道什么问题，希望能有能人给我解答！最好能直接给出改进后的程序！！！ clc clear I11=(imread('1.bmp25.bmp')); I1=double(I11(:,:,1)); %figure;imshow(I1,[]);axis on I21=(imread('1.bmp65.bmp')); I2=double(I21(:,:,1)); I31=(imread('1.bmp105.bmp')); I3=double(I31(:,:,1)); I41=(imread('1.bmp145.bmp')); I4=double(I41(:,:,1)); Ips=(I1+i*I2-I3-i*I4)/4; %figure;imshow(abs(Ips),[]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Ips=(I1+i*I2-I3-i*I4)/4; Ips=[zeros(79,768);Ips;zeros(113,768)]; L=7.68; %figure;imshow(abs(Ips),[]);axis on If=fftshift(fft2(Ips)); kk=768*3; If1=zeros(kk,kk); If1(kk/2-768/2+1:kk/2+768/2,kk/2-768/2+1:kk/2+768/2)=If; II1=ifft2(ifftshift(If1)); %figure;imshow(abs(II1),[]);axis on lamda=0.6328*10^-3; k=2*pi/lamda; a=1.007; b=1; zr=60; zo=zr-1.38; deltax = L/kk;deltay = L/kk; aa1=linspace(-L/2,L/2,kk); aa2=linspace(-L/2,L/2,kk); [x y]=meshgrid(aa1,aa2); fx = x/(deltax*deltax*kk);fy = y/(deltay*deltay*kk); H = exp(1i*k*zo*(1-lamda^2*(fx.^2+fy.^2)).^0.5); R=exp(i*k*(zr.^2+x.^2/a+y.^2/b).^0.5)./(zr.^2+x.^2/a+y.^2/b).^0.5; O=ifft2((fft2(II1.*R).*H)); figure;imshow((abs(fftshift(O))).^2,[0 4*10^6]);axis on %figure;plot(1:2048,(abs(fftshift(O(:,1000)))).^2) |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 新建文件夹.rar
2015-04-19 16:15:13, 712.53 K
» 猜你喜欢
265求调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有21人回复
-
295材料求调剂，一志愿武汉理工085601专硕
已经有4人回复
-
0856调剂，是学校就去
已经有9人回复
-
298-一志愿中国农业大学-求调剂
已经有9人回复
-
求调剂
已经有4人回复
-
学校已经提交到NSFC，还能修改吗？
已经有7人回复
-
材料学硕297已过四六级求调剂推荐
已经有6人回复
-
295复试调剂
已经有7人回复
-
本人考085602 化学工程 专硕
已经有20人回复
找到一些相关的精华帖子，希望有用哦~
关于安徽省2015年度自然科学基金拟立项项目的公示
已经有20人回复
- 2014年度浙江省自然科学基金拟资助项目清单
已经有14人回复
- 2014年天津市第一批科技计划项目立项公告
已经有19人回复
- 场发射性质可以用在哪些方面呢？
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
