| 查看: 1319 | 回复: 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-21 10:32:09, 712.53 K
» 猜你喜欢
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
面上合作单位盖章
已经有8人回复
-
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
找到一些相关的精华帖子，希望有用哦~
关于安徽省2015年度自然科学基金拟立项项目的公示
已经有20人回复
- 2014年度浙江省自然科学基金拟资助项目清单
已经有14人回复
- 2014年天津市第一批科技计划项目立项公告
已经有19人回复
- 场发射性质可以用在哪些方面呢？
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
5