角谱 matlab,关于角谱法实现数字全息 - 程序语言 - MATLAB/Mathematica

通过角谱法实现数字全息，程序语法没有错误，但就是重现出来的图是黑色一片，不知道什么问题，希望能有能人给我解答！最好能直接给出改进后的程序！！！

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)

角谱 matlab,关于角谱法实现数字全息 - 程序语言 - MATLAB/Mathematica - 小木虫论坛-学术科研互动平台...相关推荐

matlab寻峰代码,寻峰的函数！！ - 程序语言 - MATLAB/Mathematica - 小木虫论坛-学术科研互动平台...
我这里的数据是pgm的,我将其处理成多个高斯拟合的形式,现阶段只能将其最大的那个拟合出来,其他的高斯拟合我需要找到其峰值的位置! 我把前边的语句先列举上: function [ OutArr ] = ...
matlab迭代算法实例sor,SOR迭代 - 程序语言 - MATLAB/Mathematica - 小木虫论坛-学术科研互动平台...
方法一:建立了SOR.m的脚本文件,实现的是SOR迭代,程序语言如下: %SOR迭代 clear; clc; format long; i=1; n=6; H=hilb(n); X=ones(n,1) ...
matlab模糊控制m函数,模糊控制m文件运行出错 - 程序语言 - MATLAB/Mathematica - 小木虫论坛-学术科研互动平台...
Error using parsrule (line 182) Output MF index is too high Error in readfis (line 231) out=parsrule ...
用MATLAB做结构灵敏度分析怎么做,如何做灵敏度分析 - 数学 - 小木虫 - 学术科研互动社区...
我有以下数据年份平均气温最高气温最低气温相对湿度日照时数实际风速日均ET 1953 ...
matlab双线性z变换法设计数字低通滤波器
matlab双线性z变换法设计数字低通滤波器双线性z变换法利用了正切函数的非线性特点,将整个jΩ轴压缩到了单位圆的一周上. 低通: clear; close all; clc; fp=100;fs= ...
豪斯荷尔德变换 matlab,隐式QR法求实矩阵的全部特征值matlab实现
隐式QR法求实矩阵的全部特征值matlab实现隐式QR法求实矩阵的全部特征值matlab实现要求:用matlab编写通用子程序,利用隐式QR法求实矩阵的全部特征值和特征向量. 思想:隐式QR法实质 ...
MATLAB计算杨氏模量,四阶弹性模量Cijkl如何在matlab里表示啊？ - 计算模拟 - 小木虫 - 学术科研互动社区...
matlab 四元数运算计算包就可以了吧 Matlab 四元数操作函数 2012-06-03 21:02:55| 分类: MATLAB&Mathemati | 标签:四元数 quater ...
matlab decomposition filters,MATLAB小波去噪求助（附算法和显示图片）！不知自己哪个地方出了问题，求指点！ - 信息科学 - 小木虫 - 学术科研互动社区...
star2883 不是非常明白楼主的意思.我的理解是: 问题1:为什么中间两幅图,显示的都是空白? 答:你是指Figure2和Figure3的显示中,全是白色吗?这个是因为你用了colormap(ma ...
matlab中并行用不了,请教matlab在linux下的并行问题 - 程序语言 - 小木虫 - 学术科研互动社区...
刚在网上搜了一下,matlabpool open local 8 设置节点数,但是必须在matlab下面才行,比如我打开linux, >>matlab 进入matlab系统,然后matla ...
matlab比Excel好在哪,matlab计算结果和excel完全不一样 - 仿真模拟 - 小木虫 - 学术科研互动社区...
请教各位: 我有一段程序通过matlab计算出的结果和用Excel中VB计算的结果完全不一样,跪求指导.....(已经确定是程序中这段的问题,实在找不出2者的区别) 两种语言中中Ndata=n, i ...

角谱 matlab,关于角谱法实现数字全息 - 程序语言 - MATLAB/Mathematica - 小木虫论坛-学术科研互动平台...

角谱 matlab,关于角谱法实现数字全息 - 程序语言 - MATLAB/Mathematica - 小木虫论坛-学术科研互动平台...相关推荐

最新文章

热门文章