匿名

用戶注銷 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1282.6
散金: 853
紅花: 1
帖子: 172
在線: 85小時(shí)
蟲號(hào): 0
注冊: 2012-03-25
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

本帖僅樓主可見

» 本帖已獲得的紅花（最新10朵）

banuo2005

» 猜你喜歡

291求調(diào)劑已經(jīng)有25人回復(fù)
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京科技，085601總分305求調(diào)劑已經(jīng)有10人回復(fù)
化學(xué)工程專碩324分，一志愿中國礦業(yè)大學(xué)求調(diào)劑已經(jīng)有6人回復(fù)
279求調(diào)劑已經(jīng)有3人回復(fù)
26年申博已經(jīng)有4人回復(fù)
材料工程322分已經(jīng)有5人回復(fù)
318求調(diào)劑，計(jì)算材料方向已經(jīng)有7人回復(fù)
292分，材料與化工，申請調(diào)劑已經(jīng)有20人回復(fù)
085600 295分求調(diào)劑已經(jīng)有18人回復(fù)

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

析出相衍射斑點(diǎn)如何標(biāo)定已經(jīng)有6人回復(fù)
MATLAB入門經(jīng)典--word 版已經(jīng)有116人回復(fù)
誰能講講x射線三軸衍射及reciprocal space mapping技術(shù)��？已經(jīng)有14人回復(fù)
請教向量Q滿足高斯分布，用matlab如何表示？已經(jīng)有11人回復(fù)
Fluent中壓力入口邊界條件怎么設(shè)置已經(jīng)有21人回復(fù)
matlab處于busy 已經(jīng)有7人回復(fù)
TEM（透射電鏡）多晶體衍射環(huán)分析已經(jīng)有11人回復(fù)
matlab反求未知輸入?yún)?shù) 已經(jīng)有4人回復(fù)
新人求指教啊請問怎么從matlab圖中提取600組數(shù)據(jù) 已經(jīng)有11人回復(fù)
matlab解微分方程已經(jīng)有10人回復(fù)
請教關(guān)于matlab矩陣轉(zhuǎn)換問題已經(jīng)有6人回復(fù)
fluent中的體熱源加載問題（實(shí)際上是熱流密度，不是熱生成率）已經(jīng)有19人回復(fù)
ansys 13.0 網(wǎng)格劃分和邊界定義怎么設(shè)置，拜求高手已經(jīng)有13人回復(fù)
求matlab 解方程的問題已經(jīng)有5人回復(fù)
求助大俠幫忙，有關(guān)matlab模擬數(shù)據(jù) 已經(jīng)有29人回復(fù)
【求助】判斷菲涅爾區(qū)及夫瑯和費(fèi)衍射區(qū)的距離，菲涅爾數(shù) 已經(jīng)有5人回復(fù)

1樓 2012-04-04 00:34:27

已閱同方向廣播申請物理EPI 回復(fù)此樓編輯查看我的主頁

匿名

用戶注銷 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1282.6
散金: 853
紅花: 1
帖子: 172
在線: 85小時(shí)
蟲號(hào): 0
注冊: 2012-03-25
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

本帖僅樓主可見

贊一下

2樓2012-04-04 00:35:55

已閱申請物理EPI 回復(fù)此樓編輯查看我的主頁

匿名

用戶注銷 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1282.6
散金: 853
紅花: 1
帖子: 172
在線: 85小時(shí)
蟲號(hào): 0
注冊: 2012-03-25
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

本帖僅樓主可見

贊一下

3樓2012-04-04 00:37:31

已閱申請物理EPI 回復(fù)此樓編輯查看我的主頁

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
glddoyss: 金幣+5, ★有幫助, 雖然沒指出我的錯(cuò)誤，不過提供的程序也有用啦 2012-04-04 15:24:08

利用光波衍射的菲涅爾近似－菲涅爾衍射公式來計(jì)算雙縫的衍射：

下面的這個(gè)程序你明白了的話，你這個(gè)問題就解決了，其他物體的衍射模擬你也就會(huì)了。
（1）實(shí)驗(yàn)參數(shù)設(shè)定：
Lamda=0.6328; % um, Wavelength of the input light
DDist=200.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
SWidth=4; % pixels, Width of the Slits
SDist=30; % pixels, Distance between two slits
（2）建立衍射屏函數(shù)（Slit）的抽樣點(diǎn)陣
DScreen=zeros(M,N); % define a M*N matrix with zero elements;
DScreen(:,(N-SDist-SWidth)/2

N-SDist+SWidth)/2)=1.0; % Slit one
DScreen(:,(N+SDist-SWidth)/2

N+SDist+SWidth)/2)=1.0; % Slit two
%也可用下列語句調(diào)入實(shí)現(xiàn)設(shè)計(jì)好的衍射屏圖樣（bmp格式512×512）。
DScreen=double(imread('DScreen1.bmp'));
（3）衍射屏函數(shù)的圖形顯示：
figure(1); % 激活1號(hào)顯示窗口為當(dāng)前窗口
colormap(gray); ％設(shè)置灰度圖形顯示格式
imshow(Slit); ％在當(dāng)前圖形窗口顯示衍射屏函數(shù)
（4）建立高斯球面波函數(shù) 的抽樣點(diǎn)陣：
% 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
  （5）通過計(jì)算菲涅爾衍射積分求衍射光場：
% Simulate the Fresnel diffraction
DScreen_F=fftshift(fft2(fftshift(DScreen))); %Fourier transform of the slit function.
SphFunct_F=fftshift(fft2(fftshift(SphFunct))); %Fourier transform of the spherical function.
FrDiffract=fftshift(ifft2(fftshift(DScreen_F.*SphFunct_F))); % Inverse Fourier transform
  （6）顯示輸出平面上的光強(qiáng)度分布：
% Display the diffraction pattern
FrDiffract_I=FrDiffract.*conj(FrDiffract); % Intensity distribution
FrDiffract_I=FrDiffract_I/max(max(FrDiffract_I));
imshow(FrDiffract_I/max(max(FrDiffract_I)));
  （7）保存衍射圖樣

贊一下(1人)

回復(fù)此樓

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

4樓2012-04-04 09:53:16

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

匿名

用戶注銷 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1282.6
散金: 853
紅花: 1
帖子: 172
在線: 85小時(shí)
蟲號(hào): 0
注冊: 2012-03-25
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

本帖僅樓主可見

贊一下

5樓2012-04-04 12:52:11

已閱申請物理EPI 回復(fù)此樓編輯查看我的主頁

smy1982

木蟲 (正式寫手)

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

引用回帖:

5樓: Originally posted by glddoyss at 2012-04-04 12:52:11:
這個(gè)不很懂啊。學(xué)matlab不久。。能幫忙看看我函數(shù)定義哪里出了問題嗎？
M文件這個(gè)
function CS=fresnelC(a)
CS=quad(@fun1,0,a);
function fun=fun1(u)
fun=cos(0.5.*pi.*u.^2)+i.*sin(0.5.*pi.*u.^2);
...

你這個(gè)程序不是用run來運(yùn)行的，你的是個(gè)函數(shù)。得用函數(shù)命令。
在matlab命令行里直接輸入比如：CS=fresnelC(200e-3)；
但是你這個(gè)程序里面，“l(fā)(n)=abs(F(x1(n))+sqrt(0.5)*exp(i.*pi./4));%表示的是菲涅爾積分，f(n)=F(x1(n));%得到菲涅爾積分的列陣
” 這個(gè)F（）是什么函數(shù)��？沒有給出來��？

贊一下

回復(fù)此樓

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

6樓2012-04-04 20:07:42

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

zhfeng

木蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 4592.3
散金: 6
帖子: 202
在線: 237.8小時(shí)
蟲號(hào): 623364
注冊: 2008-10-11
性別: GG
專業(yè): 光學(xué)

glddoyss: 回帖置頂 2012-04-06 19:28:32

對于quad函數(shù)，積分限應(yīng)該是常數(shù)，a必須賦值。可以看一下幫助文件

贊一下(1人)

回復(fù)此樓

7樓2012-04-04 23:20:30

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

banuo2005

送鮮花一朵

8樓2012-04-05 20:47:18

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

匿名

用戶注銷 (小有名氣)

本帖僅樓主可見

9樓2012-04-06 19:28:49

已閱申請物理EPI 回復(fù)此樓編輯查看我的主頁

卡洛比

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 88
帖子: 8
在線: 12.1小時(shí)
蟲號(hào): 1614147
注冊: 2012-02-13
專業(yè): 光學(xué)

??? Error using ==> times
Matrix dimensions must agree.

Error in ==> Fresnelyanshe at 25
FrDiffract=fftshift(ifft2(fftshift(DScreen_F.*SphFunct_F))); % Inverse Fourier transform傅里葉逆變換（球面函數(shù)及狹縫函數(shù)的乘積）

??? Error using ==> times
Matrix dimensions must agree.

Error in ==> Fresnelyanshe at 25
FrDiffract=fftshift(ifft2(fftshift(DScreen_F.*SphFunct_F))); % Inverse Fourier transform傅里葉逆變換（球面函數(shù)及狹縫函數(shù)的乘積）
同學(xué) 你這個(gè)程序運(yùn)行下去出現(xiàn)了如上錯(cuò)誤，是不是你球面函數(shù)和狹縫函數(shù)的取樣尺寸不同導(dǎo)致此上錯(cuò)誤呢？

贊一下

回復(fù)此樓

10樓2012-04-25 10:36:58

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

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

[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +12	睿08 2026-03-27	12/600	2026-04-01 22:36 by Dyhoer
[考研] 337求調(diào)劑 +9	《樹》 2026-03-29	9/450	2026-04-01 18:05 by 記事本2026
[考研] 生物與醫(yī)藥考研調(diào)劑 +5	鐵憨憨123425 2026-03-31	5/250	2026-04-01 18:01 by syh9288
[考研] 324分 085600材料與化工 +17	呆鵝oor 2026-03-27	17/850	2026-04-01 17:51 by JYD2011
[考研] 086502化學(xué)工程342求調(diào)劑 +7	阿姨復(fù)古不過 2026-03-27	7/350	2026-04-01 16:14 by yanflower7133
[考研] 調(diào)劑 +5	好好讀書。 2026-03-28	7/350	2026-04-01 15:32 by 王亮_大連醫(yī)科大
[考研] 0856初試324分求調(diào)劑 +5	想上學(xué)求調(diào) 2026-04-01	5/250	2026-04-01 15:28 by 考研上岸呀一定?/a>
[考研] 材料0856 英一數(shù)二 323 求調(diào)劑 +9	袁sy 2026-04-01	9/450	2026-04-01 14:30 by wangjy2002
[考研] 0855機(jī)械初試280求調(diào)劑 +3	kazenotori 2026-03-31	3/150	2026-04-01 10:08 by JourneyLucky
[考研] 349求調(diào)劑 +6	zwjjjjjj 2026-03-31	6/300	2026-04-01 09:16 by JourneyLucky
[考研] 初試301，代碼085701環(huán)境工程，本碩一致，四六級已過，有二區(qū)一作，共發(fā)表5篇論文 +3	axibli 2026-04-01	3/150	2026-04-01 08:43 by i_cooler
[考研] 085601 329分調(diào)劑 +6	yzsa12 2026-03-31	6/300	2026-03-31 15:23 by yanflower7133
[考研] 274求調(diào)劑 +6	xiao愛同學(xué) 2026-03-30	6/300	2026-03-31 10:04 by cal0306
[考研] 085600 286分材料求調(diào)劑 +11	麻辣魷魚 2026-03-27	12/600	2026-03-30 19:33 by Wang200018
[考研] 348求調(diào)劑 +6	小懶蟲不懶了 2026-03-28	6/300	2026-03-30 10:29 by Evan_Liu
[考研] 求調(diào)劑 +10	張zz111 2026-03-27	11/550	2026-03-30 09:17 by 無際的草原
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7	前行必有光 2026-03-28	8/400	2026-03-28 23:22 by 小木蟲tim
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 08開頭275求調(diào)劑 +4	拉誰不重要 2026-03-26	4/200	2026-03-27 14:12 by Delta2012