簡介本資源是《現(xiàn)代信息光學(xué)與MATLAB》中英雙語合訂版教材PDF系統(tǒng)講授傅里葉光學(xué)、光學(xué)圖像處理、數(shù)字全息及計(jì)算機(jī)生成全息等核心內(nèi)容特別強(qiáng)化非相干光學(xué)系統(tǒng)建模與前沿數(shù)字全息技術(shù)實(shí)現(xiàn)面向工程與應(yīng)用物理專業(yè)高年級本科生、研究生及光學(xué)信息處理領(lǐng)域研究人員助力理論理解與MATLAB仿真能力同步提升。資源為1個24.87MB的PDF文件完整涵蓋全書章節(jié)、公式推導(dǎo)、MATLAB代碼示例、課后習(xí)題及雙語作者簡介結(jié)構(gòu)清晰、步驟詳盡無須額外安裝或解壓其他依賴。目前已有96人學(xué)習(xí)下載讀者可直接獲取權(quán)威教材原文、配套算法實(shí)現(xiàn)邏輯、光學(xué)現(xiàn)象可視化方法及中英文術(shù)語對照體系尤其適合開展數(shù)字全息實(shí)驗(yàn)設(shè)計(jì)、光學(xué)圖像處理課程設(shè)計(jì)或科研建模參考。1. 為什么光學(xué)實(shí)驗(yàn)不再依賴暗室和激光臺用 MATLAB 把傅里葉變換、衍射仿真、全息圖生成全搬到筆記本上跑通你有沒有試過在實(shí)驗(yàn)室調(diào)一上午光路就為了驗(yàn)證一個簡單的夫瑯禾費(fèi)衍射公式激光器溫漂、透鏡軸向偏移、CCD噪聲干擾——這些物理世界的“玄學(xué)變量”常讓理論推導(dǎo)和實(shí)測結(jié)果差出兩個數(shù)量級。而《現(xiàn)代信息光學(xué)與MATLAB》這個標(biāo)題背后不是教你怎么寫 MATLAB 語法而是提供一套可復(fù)現(xiàn)、可調(diào)試、可量化驗(yàn)證的光學(xué)建模工作流把光波傳播抽象為復(fù)數(shù)矩陣運(yùn)算把空間濾波變成二維卷積核設(shè)計(jì)把全息圖編碼轉(zhuǎn)化為相位迭代優(yōu)化。它適合三類人高校光學(xué)/光電專業(yè)學(xué)生替代部分昂貴實(shí)驗(yàn)課、工業(yè)界成像算法工程師快速驗(yàn)證光學(xué)預(yù)處理方案、跨領(lǐng)域研究者如用數(shù)字全息做微流控粒子追蹤。核心價值在于——所有結(jié)果都帶明確物理量綱單位是 mm、μm、nm、rad所有參數(shù)可調(diào)、所有中間態(tài)可可視化沒有黑匣子。這不是“MATLAB 光學(xué)插件包”而是用基礎(chǔ)函數(shù)fft2、ifft2、fspecial、imfilter搭出來的、經(jīng)得起同行復(fù)現(xiàn)的光學(xué)計(jì)算棧。2. 從標(biāo)量衍射理論出發(fā)為什么必須用復(fù)振幅建模而不是直接處理灰度圖光學(xué)信息處理的本質(zhì)是光場在空間中的線性傳播與非線性調(diào)制。若只把圖像當(dāng)灰度矩陣處理會徹底丟失相位信息——而這恰恰是干涉、全息、相襯成像的物理基礎(chǔ)。MATLAB 的強(qiáng)項(xiàng)正在于對復(fù)數(shù)矩陣的原生支持double類型天然承載real 1i*imagfft2默認(rèn)輸出復(fù)頻譜imshow(abs(U))和imshow(angle(U))可分通道查看振幅與相位。下面以最典型的菲涅爾衍射為例說明如何構(gòu)建符合物理意義的傳播模型。2.1 菲涅爾衍射的離散化實(shí)現(xiàn)繞不開的波長、采樣間隔與傳播距離三重約束菲涅爾近似下的衍射積分可離散化為$$ U(x,y,z) \frac{e^{ikz}}{i\lambda z} e^{i\frac{k}{2z}(x^2y^2)} \iint U_0(\xi,\eta) e^{i\frac{k}{2z}[(x-\xi)^2(y-\eta)^2]} d\xi d\eta $$在 MATLAB 中我們不直接數(shù)值積分而是利用卷積定理將其轉(zhuǎn)為頻域乘法。關(guān)鍵步驟如下function Uz fresnel_propagate(U0, lambda, dx, dy, z, N) % U0: 輸入復(fù)振幅矩陣 (N x N) % lambda: 波長 (單位同dx,dy如 um) % dx, dy: 物平面采樣間隔 (um) % z: 傳播距離 (um) % N: 矩陣尺寸 k 2*pi/lambda; % 計(jì)算頻域傳遞函數(shù) H(u,v) du 1/(N*dx); dv 1/(N*dy); % 頻域采樣間隔 [u, v] meshgrid((-N/2:N/2-1)*du, (-N/2:N/2-1)*dv); H exp(1i*k*z) .* exp(-1i*pi*lambda*z*(u.^2 v.^2)); % 注意此處已含相位因子 % 物平面補(bǔ)零至2N保證無混疊 U0_padded padarray(U0, [N N], post); U0_fft fftshift(fft2(fftshift(U0_padded))); % 頻域相乘 逆變換 Uz_padded ifftshift(ifft2(ifftshift(U0_fft .* H))); Uz Uz_padded(N1:2*N, N1:2*N); % 截取中心N×N區(qū)域 end邏輯說明該函數(shù)嚴(yán)格遵循標(biāo)量衍射理論中菲涅爾近似的數(shù)學(xué)形式。H是頻域傳遞函數(shù)其指數(shù)項(xiàng)中的-pi*lambda*z*(u2v2)來源于對原始積分式中二次相位因子的傅里葉變換結(jié)果。注意fftshift的兩次使用——這是 MATLABfft2默認(rèn)將零頻置于左上角而物理頻譜需零頻居中必須手動搬移。參數(shù)說明dx,dy決定空間帶寬若物平面最小特征尺寸為 5 μm則dx必須 ≤ 2.5 μm滿足奈奎斯特采樣z不能過小當(dāng)z (N*dx)2/(2*lambda)時菲涅爾近似失效應(yīng)改用更精確的角譜法N建議為 2 的整數(shù)次冪如 512、1024避免 FFT 效率損失。2.2 夫瑯禾費(fèi)衍射菲涅爾的極限特例也是最常用的頻譜分析工具當(dāng)傳播距離足夠大z (Δx)2/λ其中Δx為孔徑尺寸菲涅爾相位因子中的一次項(xiàng)主導(dǎo)二次項(xiàng)可忽略此時衍射圖即為輸入光場的傅里葉變換模平方。這正是光學(xué)信息處理中“4f 系統(tǒng)”的理論基礎(chǔ)。function Ufar farfield_diffraction(U0, lambda, dx, dy, f) % U0: 輸入復(fù)振幅 (N x N) % f: 等效傅里葉變換透鏡焦距 (um)用于縮放頻譜尺寸 % 輸出 Ufar 為夫瑯禾費(fèi)衍射復(fù)振幅未取模 N size(U0,1); k 2*pi/lambda; % 空間頻率坐標(biāo) fx (-N/2:N/2-1)/(N*dx); fy (-N/2:N/2-1)/(N*dy); [FX, FY] meshgrid(fx, fy); % 相位因子exp(i*k*f) * exp(-i*pi*lambda*f*(FX.^2FY.^2)) 在遠(yuǎn)場可忽略二次項(xiàng) % 實(shí)際只需做 FFT 并加全局相位 Ufar fftshift(fft2(fftshift(U0))); % 縮放使頻譜中心對應(yīng)零頻單位為 mm?1 → 乘 f 得到像面坐標(biāo) (mm) scale_factor f * lambda; % 單位mm因 fx 單位為 mm?1 Ufar Ufar * scale_factor; % 保持能量守恒的歸一化可選 end關(guān)鍵點(diǎn)此函數(shù)輸出的是復(fù)振幅頻譜而非強(qiáng)度。abs(Ufar).^2才是探測器記錄的夫瑯禾費(fèi)衍射圖樣。很多初學(xué)者誤將fft2結(jié)果直接imshow看到一片黑——因?yàn)槟J(rèn)imshow對復(fù)數(shù)取abs但未歸一化動態(tài)范圍極大。正確做法是imshow(mat2gray(abs(Ufar).^2))或手動設(shè)定clim。3. 光學(xué)系統(tǒng)建模實(shí)戰(zhàn)4f 濾波、阿貝成像原理與空間頻率響應(yīng)可視化光學(xué)信息處理的核心思想是把圖像看作不同空間頻率成分的疊加再通過“濾波器”如光闌、掩模選擇性地通過或阻擋某些頻率。MATLAB 提供了完整的信號處理鏈路從空域輸入 → 頻域分解 → 濾波操作 → 空域重建。本節(jié)以經(jīng)典 4f 系統(tǒng)為例完整走通這一閉環(huán)。3.1 構(gòu)建 4f 系統(tǒng)用兩個傅里葉變換和一個頻域掩模串聯(lián)一個標(biāo)準(zhǔn) 4f 系統(tǒng)由四部分組成輸入面 → 第一透鏡前焦面頻譜面→ 頻譜面放置濾波器 → 第二透鏡后焦面輸出面。在 MATLAB 中這等價于對輸入U_in做 FFT → 得到頻譜U_f將U_f與濾波器H_filter逐點(diǎn)相乘 →U_f_filtered對U_f_filtered做 IFFT → 得到輸出U_out% 示例低通濾波圓形光闌 N 512; U_in imread(letter_A.png); % 讀入二值圖像 U_in im2double(U_in); % 歸一化為 0~1 U_in U_in - 0.5; % 轉(zhuǎn)為 -0.5 ~ 0.5模擬透射率調(diào)制0全吸收1全透過 U_in complex(U_in, zeros(size(U_in))); % 初始化為復(fù)振幅相位為0 % 步驟1前焦面頻譜夫瑯禾費(fèi)衍射 U_f farfield_diffraction(U_in, 0.6328, 10, 10, 1000); % lambda632.8nm, dxdy10um, f1000um % 步驟2設(shè)計(jì)低通濾波器半徑 R 對應(yīng)截止頻率 fc R/(N*dx) R 30; % 頻譜面像素半徑 [H_filter, fx, fy] create_circular_filter(N, R); U_f_filtered U_f .* H_filter; % 步驟3逆傅里葉變換重建輸出 U_out ifftshift(ifft2(ifftshift(U_f_filtered))); U_out real(U_out); % 輸出為實(shí)數(shù)強(qiáng)度分布忽略相位 % 可視化 figure; subplot(1,3,1); imshow(U_in,[]); title(Input Object); subplot(1,3,2); imshow(mat2gray(abs(U_f).^2),[]); title(Fourier Spectrum); subplot(1,3,3); imshow(mat2gray(U_out),[]); title(Filtered Output);function [H, fx, fy] create_circular_filter(N, R) % 生成 N×N 圓形低通濾波器半徑 R像素 fx (-N/2:N/2-1)/N; fy (-N/2:N/2-1)/N; [FX, FY] meshgrid(fx, fy); H sqrt(FX.^2 FY.^2) R/N; % 歸一化半徑保證圓心在中心 end參數(shù)說明R是濾波器在頻譜面上的像素半徑其物理截止頻率為fc R / (N * dx)單位μm?1。例如N512,dx10 μm,R30→fc ≈ 0.00586 μm?1 5.86 mm?1即能分辨最小線寬約1/(2*fc) ≈ 0.085 mm。這個換算關(guān)系是連接代碼參數(shù)與光學(xué)性能指標(biāo)的關(guān)鍵橋梁。3.2 阿貝成像原理驗(yàn)證為什么顯微鏡分辨率受限于物鏡 NA阿貝指出物體只有被物鏡收集的衍射級才能參與成像。因此提高分辨率 收集更高階衍射斑。我們可用上述 4f 框架定量驗(yàn)證不同“數(shù)值孔徑”NA對應(yīng)的截止頻率。濾波器半徑 R對應(yīng) NA近似理論最小可分辨線寬 δ (μm)觀察現(xiàn)象10~0.1~3.2字母邊緣嚴(yán)重模糊細(xì)節(jié)丟失30~0.3~1.1A 的橫杠可見但內(nèi)部紋理仍糊80~0.8~0.4橫杠邊緣銳利內(nèi)部柵格隱約可辨提示NA 與頻譜面半徑的關(guān)系為R ≈ NA * N * dx / lambda。代入lambda0.6328 μm,dx10 μm,N512可得R≈80對應(yīng)NA≈0.8。這正是高端顯微物鏡的典型值。4. 數(shù)字全息圖生成與再現(xiàn)從干涉記錄到三維重構(gòu)的端到端流程全息術(shù)是信息光學(xué)皇冠上的明珠——它不僅記錄光強(qiáng)更記錄光波的全部復(fù)振幅信息。傳統(tǒng)光學(xué)全息需穩(wěn)定干涉環(huán)境而數(shù)字全息Digital Holography將參考光與物光的干涉圖樣用 CCD 記錄再用 MATLAB 數(shù)值再現(xiàn)。本節(jié)聚焦離軸數(shù)字全息最易實(shí)現(xiàn)、抗噪性強(qiáng)完成從合成干涉圖到三維層析的全流程。4.1 合成干涉圖物光 參考光 → 強(qiáng)度疊加 ≠ 簡單相加干涉圖強(qiáng)度為$$ I(x,y) |U_o(x,y)|^2 |U_r(x,y)|^2 2\Re{U_o(x,y)U_r^*(x,y)} $$其中U_o是物光復(fù)振幅含三維形貌信息U_r是傾斜平面波參考光提供載頻分離孿生像。關(guān)鍵在第三項(xiàng)——它攜帶物光全部信息且被參考光的載頻調(diào)制到高頻區(qū)。function I_holo generate_offaxis_hologram(U_obj, lambda, dx, dy, alpha, beta, N) % U_obj: 物光復(fù)振幅 (N x N)可由高度圖生成 % alpha, beta: 參考光波矢方向余弦弧度決定載頻位置 % 返回干涉圖 I_holo (N x N) k 2*pi/lambda; % 生成傾斜參考光U_r exp(i*k*(alpha*x beta*y)) [x, y] meshgrid((-N/2:N/2-1)*dx, (-N/2:N/2-1)*dy); U_ref exp(1i*k*(alpha*x beta*y)); % 干涉圖 |Uo|2 |Ur|2 2*Re(Uo*Ur*) I_holo abs(U_obj).^2 abs(U_ref).^2 2*real(U_obj .* conj(U_ref)); end邏輯說明U_ref的構(gòu)造是核心。alpha和beta控制參考光入射角從而在頻譜中產(chǎn)生一個遠(yuǎn)離零頻的載頻點(diǎn)。典型取值alpha 0.01,beta 0沿 x 方向傾斜 0.57°此時載頻位于(u,v) ≈ (alpha/lambda, 0) ≈ (0.0158, 0) mm?1在 512×512 頻譜中約在第 120 行附近與零頻充分分離。4.2 全息圖再現(xiàn)頻譜濾波 逆傅里葉變換 三維形貌重建再現(xiàn)過程本質(zhì)是對干涉圖I_holo做 FFT → 在頻譜中選取包含物光信息的載頻分量濾除零頻和共軛分量→ IFFT → 得到復(fù)振幅再現(xiàn)波前 → 取angle()得相位 → 轉(zhuǎn)換為高度。function [height_map, U_recon] reconstruct_hologram(I_holo, lambda, dx, dy, z, filter_radius) % I_holo: 干涉圖 % z: 再現(xiàn)距離對應(yīng)虛擬觀察距離 % filter_radius: 頻譜濾波器半徑像素 N size(I_holo,1); % 步驟1干涉圖頻譜 I_fft fftshift(fft2(fftshift(I_holo))); % 步驟2設(shè)計(jì)帶通濾波器中心在載頻位置半徑 filter_radius % 假設(shè)載頻在 (u0,v0) (120, 256) —— 需根據(jù)實(shí)際 alpha/beta 計(jì)算 u0 120; v0 256; [U, V] meshgrid(1:N, 1:N); H_filter ((U-v0).^2 (V-u0).^2) filter_radius^2; % 步驟3濾波 逆變換 I_filtered I_fft .* H_filter; U_recon ifftshift(ifft2(ifftshift(I_filtered))); % 步驟4相位解包裹 高度轉(zhuǎn)換 phase angle(U_recon); phase_unwrapped unwrap(phase, [], 1); % 沿行解包裹 height_map phase_unwrapped * lambda / (2*pi); % 單位μm end參數(shù)說明filter_radius通常取 10~20 像素。太小則物光信息截斷太大則混入噪聲和共軛像。z不影響相位解算因離軸全息再現(xiàn)為菲涅爾傳播此處簡化為直接取相位但若要做深度掃描需對不同z值重復(fù)fresnel_propagate。5. 避坑指南那些讓光學(xué)仿真實(shí)驗(yàn)集體翻車的 4 個硬核陷阱光學(xué)仿真是個“差之毫厘失之千里”的領(lǐng)域。以下是我?guī)н^多個模擬項(xiàng)目 X 時學(xué)生踩得最多、debug 時間最長的 4 類問題。每一條都附帶真實(shí)報錯現(xiàn)象、物理根源和可立即執(zhí)行的檢查清單。5.1 現(xiàn)象衍射圖樣完全對稱但與理論預(yù)測的方位相反如豎條紋變橫條紋原因fft2和ifft2的默認(rèn)零頻位置在左上角而物理頻譜要求零頻在中心。未用fftshift搬移導(dǎo)致相位因子符號錯誤進(jìn)而使空間坐標(biāo)系鏡像翻轉(zhuǎn)。解決強(qiáng)制在所有 FFT/IFFT 前后加fftshift/ifftshift。檢查代碼中fft2出現(xiàn)次數(shù)確保每次都有對應(yīng)fftshift用meshgrid生成的x,y坐標(biāo)也需確認(rèn)是否以(0,0)為中心。自查命令max(abs(fftshift(fft2(fftshift(ones(8,8)))) - ones(8,8)))應(yīng)返回接近 0若為64說明fftshift缺失。5.2 現(xiàn)象全息再現(xiàn)高度圖出現(xiàn)密集條紋狀偽影且隨z變化劇烈抖動原因相位解包裹unwrap失敗。當(dāng)物面高度變化超過λ/2如 316 nm相鄰像素相位跳變 πunwrap無法識別真實(shí)相位連續(xù)性強(qiáng)行“接續(xù)”導(dǎo)致條紋。解決改用二維相位解包裹算法如 Goldstein 算法或先對相位圖做高斯平滑imgaussfilt再unwrap。更可靠的是——用多波長或多步相移法但仿真中可先降低物面斜率減小alpha/beta測試。自查命令max(abs(diff(angle(U_recon),1,1)))若 3.0說明存在大梯度unwrap必然失敗。5.3 現(xiàn)象4f 系統(tǒng)輸出一片全黑或僅中心一個亮點(diǎn)原因輸入U_in未轉(zhuǎn)為復(fù)數(shù)類型或U_in值域錯誤。imread讀入uint8圖像im2double后為[0,1]但光學(xué)透射率模型要求0全吸收對應(yīng)復(fù)振幅01全透過對應(yīng)10i。若直接complex(U_in)則U_in0時振幅為 0但U_in1時振幅為 1看似合理——然而干涉項(xiàng)2*Re(Uo*Ur*)要求Uo有合理相位參考。常見錯誤是U_in全為實(shí)數(shù)且無背景光導(dǎo)致干涉項(xiàng)被零頻淹沒。解決統(tǒng)一初始化U_in complex(im2double(I_obj), zeros(size(I_obj)))并確保I_obj是二值圖0/1或用U_in 0.5 0.5*complex(im2double(I_obj), zeros(...))加直流偏置。自查命令min(abs(U_in(:)))和max(abs(U_in(:)))應(yīng)分別接近 0 和 1若全為 0 或全為 1輸入有誤。5.4 現(xiàn)象菲涅爾傳播后圖像嚴(yán)重畸變出現(xiàn)環(huán)狀偽影且隨z增大而加劇原因采樣定理 violated。當(dāng)傳播距離z過大頻譜收縮高頻成分被fft2的有限帶寬截斷混疊。物理上這相當(dāng)于物平面采樣太粗無法分辨?zhèn)鞑ズ笳箤挼牟ㄇ啊＝鉀Q增大N補(bǔ)零或減小dx提高物面采樣率。經(jīng)驗(yàn)公式z_max ≈ (N*dx)^2 / (2*lambda)。例如N512,dx10 μm,lambda0.6328 μm→z_max ≈ 20.5 mm。若需仿真z50 mm則必須N≥1200或dx≤6.5 μm。自查命令計(jì)算z_actual / z_max若 0.8必須升級采樣。6. 進(jìn)階技巧用 GPU 加速百倍衍射仿真以及如何把 MATLAB 模型部署為 Python 接口當(dāng)你的項(xiàng)目從單張圖擴(kuò)展到視頻序列如動態(tài)全息顯示、或需要遍歷數(shù)百組參數(shù)如優(yōu)化濾波器形狀純 CPU 的fft2會成為瓶頸。MATLAB 自 R2019a 起原生支持 GPU 數(shù)組而無需改寫核心算法——只需將輸入數(shù)據(jù)gpuArray化后續(xù)函數(shù)自動調(diào)用 CUDA。6.1 GPU 加速三行代碼提速 80~120 倍% 原始 CPU 版本耗時 ~1.2s on i7-10875H tic; Uz_cpu fresnel_propagate(U0, 0.6328, 10, 10, 10000, 1024); toc; % Elapsed time is 1.234567 seconds. % GPU 版本耗時 ~0.015s on RTX 3060 U0_gpu gpuArray(U0); % 數(shù)據(jù)搬入 GPU 顯存 Uz_gpu fresnel_propagate(U0_gpu, 0.6328, 10, 10, 10000, 1024); % 函數(shù)內(nèi)所有計(jì)算在 GPU Uz_cpu_from_gpu gather(Uz_gpu); % 搬回 CPU 內(nèi)存 toc; % Elapsed time is 0.014923 seconds.關(guān)鍵點(diǎn)fresnel_propagate函數(shù)本身無需修改MATLAB 會自動將fft2、.*、exp等內(nèi)置函數(shù)映射到 GPU。唯一要求是輸入為gpuArray且所有中間變量如H,U0_fft自動繼承 GPU 屬性。內(nèi)存瓶頸在于顯存容量——1024×1024 復(fù)數(shù)矩陣占約 16 MBRTX 30606GB可同時處理 300 個。6.2 MATLAB 模型轉(zhuǎn) Python 接口用 MATLAB Compiler SDK 生成.so庫工業(yè)界常用 Python 做前端交互與數(shù)據(jù)處理但核心光學(xué)算法仍用 MATLAB 開發(fā)因其信號處理生態(tài)成熟。MATLAB Compiler SDK 可將函數(shù)編譯為 C/C 共享庫供 Python 調(diào)用。步驟概要將fresnel_propagate.m封裝為類OpticalPropagator添加compile方法運(yùn)行compiler.build.cppSharedLibrary(OpticalPropagator, OpticalPropagator.m)生成libOpticalPropagator.so及頭文件Python 中用ctypes加載并傳參import ctypes import numpy as np lib ctypes.CDLL(./libOpticalPropagator.so) lib.fresnel_propagate.argtypes [ np.ctypeslib.ndpointer(dtypenp.complex128, flagsC_CONTIGUOUS), ctypes.c_double, ctypes.c_double, ctypes.c_double, ctypes.c_double, ctypes.c_int ] lib.fresnel_propagate.restype None U0 np.random.rand(1024,1024).astype(np.complex128) Uz np.zeros_like(U0) lib.fresnel_propagate(U0, 0.6328, 10.0, 10.0, 10000.0, 1024)血淚經(jīng)驗(yàn)MATLAB 編譯時務(wù)必用coder.typeof顯式聲明輸入類型如coder.typeof(complex(0), [1024 1024])否則 Python 傳入的ndarray可能被解釋為double而非complex128導(dǎo)致結(jié)果全零。這是我在某公司部署全息渲染引擎時花兩天才定位的坑。最后說句實(shí)在話我堅(jiān)持用 MATLAB 做光學(xué)仿真并非守舊而是它把fft2、ifft2、fspecial這些底層算子封裝得足夠干凈讓你專注物理建模而不是跟 FFTW 的 plan 創(chuàng)建、內(nèi)存對齊、線程綁定死磕。當(dāng)你在深夜調(diào)通一個衍射仿真看到屏幕上清晰浮現(xiàn)的艾里斑那種確定性帶來的踏實(shí)感是任何框架文檔都給不了的。希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取