版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請盡快對帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

草_小草

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 148.8
散金: 5
帖子: 94
在線: 75小時(shí)
蟲號(hào): 2669777
注冊: 2013-09-21
專業(yè): 計(jì)算機(jī)軟件

[求助] 大神們，求幫我看下我這個(gè)數(shù)字菲涅爾全息圖的代碼有什么事問題已有1人參與

我這個(gè)主要思路是，用點(diǎn)云算法，myfresnel函數(shù)做菲涅爾變化。參考光為平面波，源圖片的各個(gè)點(diǎn)模擬成點(diǎn)光源。算出來的結(jié)果不太對。原圖片和結(jié)果圖片在附件里
大家?guī)臀铱纯�，如果覺得太費(fèi)眼，給我指點(diǎn)一些用matlab編寫全息的書或者論文資料，也感激不盡。信息光學(xué)的書我看過了，不過那上面大多是大段原理公式，想找些有實(shí)例的書和資料。謝謝
% -------------------------------------------------------------------------
% 功能：讀入物光光波振幅分布
[fname,pname] = uigetfile('*.*','讀入物光光波振幅分布圖象');

if pname
name = strcat(pname,fname);
Picture =im2bw( imread(name));
figure,imshow(Picture);title('物光光波振幅分布圖像');
picture = Picture(:,

;
O = picture;
end
clear fname pname name picture Picture;
% 物光光波傳播至全息面上的復(fù)振幅分布 (菲涅耳衍射）
wl = 632.8e-9;       % 單色光波長(m)
z = 0.5;                   % 全息記錄距離(m)
k = 2*pi/wl;             %%波數(shù)
detaO = 2.65e-4;          % 物光波場像素間隔(m)
detaH = 1.325e-4;  % 衍射場光波像素間隔
y=repmat((512:-1:1)',1,512)-256.5;
R=10*exp(j*k*(z*cos(pi*10/180)+y*detaH*sin(pi*10/180)));%%平面參考光復(fù)振幅
tic;
uo=myfresnel(O,detaO,detaH,z,wl);% 物光光波傳播至全息面上的復(fù)振幅分布(菲涅耳衍射）
toc
u=uo+R;       %%在全息圖上，物光波和參考光的復(fù)振幅相加
I=abs(u).^2;
figure,imshow(I,[]);
% -------------------------------------------------------------------------
function [ U ] = myfresnel(Object,DETA_O,DETA_H,Z,WL )
%%菲涅爾全息
%坐標(biāo)軸和矩陣索引的轉(zhuǎn)換：U(256.5-y,x+256.5)等于求xy坐標(biāo)系中（n-256.5,256.5-m）

k = 2*pi/WL;
U=zeros(512,512);
Uo=0;
for x=-255.5:255.5
for y=-255.5:255.5
      for x0=-127.5:127.5
         for y0=-127.5:127.5
      r=((DETA_H*x-DETA_O*x0)^2+(DETA_H*y-DETA_O*y0)^2+Z^2)^(1/2);
      Uo=Uo+(Object(128.5-y0,128.5+x0)/r)*exp(j*k*r);
         end
      end
U(256.5-y,x+256.5)=Uo;
end
end

ABCDEF.jpg

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 結(jié)果.fig

2014-05-22 10:15:50, 1.61 M

» 猜你喜歡

求材料調(diào)劑已經(jīng)有7人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有3人回復(fù)
085600材料與化工已經(jīng)有5人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有6人回復(fù)
286求調(diào)劑已經(jīng)有9人回復(fù)
焦慮已經(jīng)有12人回復(fù)
344求調(diào)劑已經(jīng)有6人回復(fù)
266求調(diào)劑已經(jīng)有9人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有18人回復(fù)
314求調(diào)劑已經(jīng)有8人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

最近卡在一步酰胺化合成很久，求大神們幫幫忙已經(jīng)有27人回復(fù)
求大神們，幫我看看出啥問題了已經(jīng)有1人回復(fù)
求親們幫我看看這個(gè)到底怎么反應(yīng)，愁死我了。已經(jīng)有13人回復(fù)
急求答案，可是好窮沒金幣分給大家，大神們幫我看下這道題該怎么解？小女在此謝過啦！已經(jīng)有3人回復(fù)
求大神幫檢查下一個(gè)簡單的UDF 已經(jīng)有13人回復(fù)
想組裝一臺(tái)電腦，求懂電腦的高手-大神幫幫忙已經(jīng)有11人回復(fù)
關(guān)于2H-MoS2模型建立，文獻(xiàn)中不是很清楚，求大神進(jìn)來幫幫我已經(jīng)有0人回復(fù)
酚醛樹脂發(fā)泡的問題，求大神指導(dǎo)！已經(jīng)有1人回復(fù)
15考研，精細(xì)或藥學(xué)未定，浙江上海未定，本科化學(xué)，求幫忙分析我該考哪所學(xué)校��！已經(jīng)有8人回復(fù)
我想問下，，我是一名的動(dòng)物藥學(xué)的學(xué)生，就要考研了。求有經(jīng)驗(yàn)的大神指導(dǎo)。。。已經(jīng)有6人回復(fù)
過來看過來觀，片片大餐接著上~請求各位植物大神大仙們幫我看看這些是什么植物呢~~ 已經(jīng)有20人回復(fù)
求各位大神幫我看下一個(gè)小小的程序，只需要你一點(diǎn)兒時(shí)間。就會(huì)幫我很大的忙已經(jīng)有15人回復(fù)
求大神幫我看看，我的計(jì)算思路有沒有問題？已經(jīng)有5人回復(fù)
請大神們看下我這個(gè)反應(yīng)，有可能得到不？已經(jīng)有8人回復(fù)

1樓 2014-05-22 10:16:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

smy1982

木蟲 (正式寫手)

物理EPI: 1
應(yīng)助: 21 (小學(xué)生)
金幣: 3128.7
散金: 6748
紅花: 12
帖子: 853
在線: 450.4小時(shí)
蟲號(hào): 305139
注冊: 2006-12-09
性別: GG
專業(yè): 光學(xué)信息獲取與處理

myfresnel
是個(gè)啥函數(shù)？

贊一下

回復(fù)此樓

如果你曾歌頌黎明，那么也請你擁抱黑夜。

2樓2014-05-23 10:16:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

smy1982

木蟲 (正式寫手)

物理EPI: 1
應(yīng)助: 21 (小學(xué)生)
金幣: 3128.7
散金: 6748
紅花: 12
帖子: 853
在線: 450.4小時(shí)
蟲號(hào): 305139
注冊: 2006-12-09
性別: GG
專業(yè): 光學(xué)信息獲取與處理

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

clear
clc
Obj=imread('Obj.bmp');
Obj=double(Obj);
[M,N]=size(Obj); % M,N are the pixel number of the object.
Lamda=0.6328; % um, Wavelength of the input light
DDist=150.0e+3; % um, Diffraction distance to the output plane
% M=512; N=512; % Number of the sample array
dx=10.0; dy=10.0; % um, Sampling space
figure(1);
colormap(gray);
imshow(Obj);%物體

% Set coordinate of the sample array
[x,y]=meshgrid(-M*dx/2:dx

M-1)*dx/2,-N*dy/2:dy

N-1)*dy/2);

SphFunct=exp(i*pi*(x.^2+y.^2)/(Lamda*DDist)); % The spherical wave function
% Simulate the Fresnel diffraction
Obj_F=fftshift(fft2(fftshift(Obj))); %Fourier transform of the slit function.
SphFunct_F=fftshift(fft2(fftshift(SphFunct))); %Fourier transform of the spherical function.
FrDiffract=fftshift(ifft2(fftshift(Obj_F.*SphFunct_F))); % Inverse Fourier transform
figure(2);
FrDiffract_I=abs(FrDiffract).^2;
FrDiffract_I=FrDiffract_I/max(max(FrDiffract_I));
imshow(FrDiffract_I,[0 0.5]);%物光波
% Display the diffraction pattern
a=40*pi/180;b=90*pi/180;
R=1000*exp(i*2*pi/Lamda*DDist*sqrt(1-cos(a).^2-cos(b).^2)).*exp(i*2*pi/Lamda*(x*cos(a)+y*cos(b)));%%平面參考光復(fù)振幅
Hol=abs(R+FrDiffract).^2;
% FrDiffract_I=FrDiffract.*conj(FrDiffract); % Intensity distribution
Hol=Hol./max(max(Hol));
figure(3);
imshow(Hol,[0 0.5]); %全息圖
%%

大神們，求幫我看下我這個(gè)數(shù)字菲涅爾全息圖的代碼有什么事問題-1

全息圖.jpg

大神們，求幫我看下我這個(gè)數(shù)字菲涅爾全息圖的代碼有什么事問題-2

物體.jpg

大神們，求幫我看下我這個(gè)數(shù)字菲涅爾全息圖的代碼有什么事問題-3

物體150毫米的衍射圖像.jpg

贊一下

回復(fù)此樓

如果你曾歌頌黎明，那么也請你擁抱黑夜。

3樓2014-05-23 11:14:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

smy1982

木蟲 (正式寫手)

物理EPI: 1
應(yīng)助: 21 (小學(xué)生)
金幣: 3128.7
散金: 6748
紅花: 12
帖子: 853
在線: 450.4小時(shí)
蟲號(hào): 305139
注冊: 2006-12-09
性別: GG
專業(yè): 光學(xué)信息獲取與處理

引用回帖:

3樓: Originally posted by smy1982 at 2014-05-23 11:14:13
clear
clc
Obj=imread('Obj.bmp');
Obj=double(Obj);
=size(Obj); % M,N are the pixel number of the object.
Lamda=0.6328; % um, Wavelength of the input light
DDist=150.0e+3; % um, Diffraction dis ...

clear
clc
Obj=imread('Obj.jpg');
Obj=double(Obj(:,:,1));
[M,N]=size(Obj); % M,N are the pixel number of the object.
Lamda=0.6328; % um, Wavelength of the input light
DDist=100.0e+3; % um, Diffraction distance to the output plane
% M=512; N=512; % Number of the sample array
dx=10.0; dy=10.0; % um, Sampling space
figure(1);
colormap(gray);
imshow(Obj);%物體

% Set coordinate of the sample array
[x,y]=meshgrid(-M*dx/2:dx

M-1)*dx/2,-N*dy/2:dy

N-1)*dy/2);

SphFunct=exp(i*pi*(x.^2+y.^2)/(Lamda*DDist)); % The spherical wave function
% Simulate the Fresnel diffraction
Obj_F=fftshift(fft2(fftshift(Obj))); %Fourier transform of the slit function.
SphFunct_F=fftshift(fft2(fftshift(SphFunct))); %Fourier transform of the spherical function.
FrDiffract=fftshift(ifft2(fftshift(Obj_F.*SphFunct_F))); % Inverse Fourier transform
figure(2);
FrDiffract_I=abs(FrDiffract).^2;
FrDiffract_I=FrDiffract_I/max(max(FrDiffract_I));
imshow(FrDiffract_I,[0 0.5]);%物光波
% Display the diffraction pattern
a=40*pi/180;b=90*pi/180;
R=100000*exp(i*2*pi/Lamda*DDist*sqrt(1-cos(a).^2-cos(b).^2)).*exp(i*2*pi/Lamda*(x*cos(a)+y*cos(b)));%%平面參考光復(fù)振幅
Hol=abs(R+FrDiffract).^2;
% FrDiffract_I=FrDiffract.*conj(FrDiffract); % Intensity distribution
Hol=Hol./max(max(Hol));
figure(3);
imshow(Hol,[0 0.5]); %全息圖
%%%

你的物體最好保存為灰度的圖像》

大神們，求幫我看下我這個(gè)數(shù)字菲涅爾全息圖的代碼有什么事問題-4

obj.jpg

大神們，求幫我看下我這個(gè)數(shù)字菲涅爾全息圖的代碼有什么事問題-5

全息圖1.jpg

贊一下

回復(fù)此樓

如果你曾歌頌黎明，那么也請你擁抱黑夜。

4樓2014-05-23 11:25:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主草_小草的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 265求調(diào)劑 +3	梁梁校校 2026-03-17	3/150	2026-03-18 09:12 by zhukairuo
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 本人考085602 化學(xué)工程專碩 +16	不知道叫什么！ 2026-03-15	18/900	2026-03-17 17:05 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[基金申請] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

草_小草

[求助] 大神們，求幫我看下我這個(gè)數(shù)字菲涅爾全息圖的代碼有什么事問題 已有1人參與

» 本帖附件資源列表

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

smy1982

smy1982

【答案】應(yīng)助回帖

smy1982

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 大神們，求幫我看下我這個(gè)數(shù)字菲涅爾全息圖的代碼有什么事問題已有1人參與