| 8 | 1/1 | 返回列表 |
| 查看: 1632 | 回復(fù): 7 | ||
wupeng2320金蟲 (小有名氣)
|
[求助]
菲涅爾衍射數(shù)值模擬 已有2人參與
|
|
最近在做菲涅爾衍射數(shù)值模擬,但在模擬上遇到了問題,希望光學(xué)加matlab的大神能不嗇賜教! 模擬中,物體是圓孔,然后菲涅爾衍射公式采用的是卷積形式:u(x,y)=exp(jkz)/(j*lamda*z)*【u0(x,y)和exp(jk/(2z)*(x^2+y^2)】。u0(x,y)是物函數(shù),模擬中用上述圓孔做物。二維離散卷積的實(shí)現(xiàn)過程是通過conv2進(jìn)行。但是結(jié)果出現(xiàn)之后不是菲涅爾,所以猜想是菲涅爾公式的實(shí)現(xiàn)出現(xiàn)了問題,還是說取樣的問題?期待各位給解答下,萬分感謝。 程序如下(附件是.m文件): clear all; %做小孔。R=50,圓心(128,128) for x=1:256 for y=1:256 if ((x-128).^2+(y-128).^2<=2500) t(x,y)=1; else t(x,y)=0; end end end %顯示圓孔圖像 subplot(121) imshow(t,[]) %菲涅爾公式參量設(shè)置 z=200; lamda=500*10^(-3); k=2*pi/lamda; %采用的菲涅爾公式形式是:u(x,y)=exp(jkz)/(j*lamda*z)*【u0(x,y)和exp(jk/(2z)*(x^2+y^2)】的卷積 %u0(x,y)是物函數(shù),本程序利用上述圓孔做物函數(shù)。 for x=1:256 for y=1:256 f(x,y)=exp(i*k*(x.^2+y.^2)/2/z); end end m=conv2(t,f) u=exp(i*k*z)/(i*lamda*z).*m %展示衍射圖 subplot(122) imshow(u.*conj(u),[]) |
至尊木蟲 (正式寫手)
火星領(lǐng)導(dǎo)

至尊木蟲 (正式寫手)
火星領(lǐng)導(dǎo)

金蟲 (小有名氣)
金蟲 (小有名氣)
新蟲 (初入文壇)
新蟲 (初入文壇)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|