简介：这是一份面向光学初学者、物理实验教学及MATLAB入门者的透镜成像原理可视化学习工具，通过GUI交互式仿真帮助理解物距、像距、焦距等核心概念与薄透镜成像规律。资源包含17个文件，以13幅BMP格式的成像过程对比图（如Image_before_Phase200.bmp、Image_after_abs600.bmp等）直观呈现不同参数下的光场相位/振幅演化，辅以3个关键M文件（fun_my_lens.m、my_lens.m、fun_Fresnel.m）实现计算逻辑，以及1个FIG界面文件（my_lens.fig）构成完整可运行GUI系统，压缩包仅1.25MB，轻量易部署。已有1222人学习下载，适合课堂教学演示、课设实践或自学探究。用户可直接运行GUI，拖动滑块实时调节透镜焦距与物距，同步观察光线追迹效果、成像位置变化及前后场图像差异，配套代码结构清晰、注释充分，便于二次开发与原理验证。
1. 这不是“画个透镜再标几条线”的演示——它是一套可调参、可验证、可复现的几何光学仿真闭环
你见过太多“透镜成像 demo”：静态图配公式，箭头从物点出发，经透镜折射后交于像点，最后贴一句“符合高斯公式”。但真实教学与工程预研需要的是——当物距从 300mm 拉到 120mm，焦距从 200mm 降到 80mm，像高怎么非线性跳变？像距是否真在焦点外连续收敛？相位延迟引入后，衍射效应如何撕裂理想高斯像？这个透镜成像模拟.zip不是幻灯片，而是一个完整闭环：GUI 输入 → 几何光路计算 → 相位/振幅双域建模 → Fresnel 衍射传播 → 多层图像输出（Image_before_Phase200.bmp到Image_after_abs600.bmp），所有中间态都落盘可查。它用 MATLAB 原生 GUIDE 构建，不依赖任何第三方工具箱，核心逻辑封装在fun_my_lens.m和fun_Fresnel.m中，连苏庆帅.bmp这张测试图都带进项目——说明作者做过实测校准。适合高校光学实验课教师快速部署课堂交互环境，也适合光电工程师在无 Zemax 许可时做初步系统参数扫掠。如果你正被“为什么理论像距和仿真结果差 5mm”这类问题卡住，这个包里my_lens.fig的控件绑定逻辑和output文件夹的逐帧图像，就是你的第一手排错依据。
2. 从 GUIDE 界面到物理模型：MATLAB GUI 光学仿真的四层数据流设计
2.1 GUIDE 界面结构解析：控件-回调-参数映射关系
该 GUI 使用 MATLAB 经典 GUIDE（而非较新的 App Designer）构建，主界面文件为my_lens.fig，对应代码文件my_lens.m。打开.fig文件后可见 7 类控件：
- 3 个滑块（uicontrol, 'Style', 'slider'）：分别控制
Object Distance (mm)、Focal Length (mm)、Lens Diameter (mm)； - 2 个文本框（uicontrol, 'Style', 'edit'）：实时显示当前物距/像距数值，支持手动输入覆盖滑块值；
- 1 个下拉菜单（uicontrol, 'Style', 'popupmenu'）：切换成像模式：
Geometric Optics（纯几何追迹）、Fresnel Diffraction（含衍射）、Phase Only（仅相位调制）； - 1 个按钮（uicontrol, 'Style', 'pushbutton'）：
Run Simulation触发核心计算流程。
提示：所有滑块的
Min/Max/Value属性在 GUIDE 编辑器中硬编码，例如物距滑块范围设为[100, 1000]，步长10，这直接约束了可探索的物理区间——若需研究微距成像（物距 < 100mm），必须手动修改my_lens.m中set(hObject, 'Min', 50)并重置初始值。
2.2 回调函数链：从用户操作到物理计算的触发路径
GUI 的响应逻辑由my_lens.m中的回调函数驱动。关键链路如下：
- 滑块移动触发
slider1_Callback→ 调用update_display()更新文本框； Run Simulation按钮触发pushbutton1_Callback→ 读取当前所有控件值 → 调用fun_my_lens()主计算函数；fun_my_lens.m接收参数后，根据下拉菜单选择调用不同子模块：if strcmp(mode, 'Geometric Optics') [image_out, obj_pos, img_pos] = geometric_ray_trace(obj_dist, f_len, lens_diam); elseif strcmp(mode, 'Fresnel Diffraction') [image_out, phase_map] = fun_Fresnel(obj_dist, f_len, lens_diam, 'wavelength', 532e-9); end
2.2.1 几何光学模块：严格遵循符号规则的高斯公式实现
geometric_ray_trace()并非简单套用1/f = 1/u + 1/v，而是实现带符号约定的完整追迹：
- 物距
u以透镜光心为原点，实物为负值（MATLAB 坐标系惯例）； - 像距
v计算后自动判断实像（v > 0，像在透镜右侧）或虚像（v < 0，像在左侧）； - 同时计算横向放大率
m = -v/u，并据此缩放苏庆帅.bmp的像素尺寸生成Image_on_Phase.bmp（即理想像平面图）。
function [img_out, obj_pos, img_pos] = geometric_ray_trace(u, f, diam) % u: object distance (negative for real object) % f: focal length (positive for converging lens) v = 1/(1/f - 1/u); % Gaussian formula with sign convention m = -v/u; % magnification % Load and resize source image src_img = imread('苏庆帅.bmp'); [h, w, ~] = size(src_img); new_h = round(h * abs(m)); new_w = round(w * abs(m)); img_out = imresize(src_img, [new_h, new_w]); obj_pos = [0, -abs(u)]; % object at (0, -|u|) img_pos = [0, v]; % image at (0, v) end注意：此处
imresize使用双线性插值，若需保留锐利边缘（如测试图中的文字），应在调用前添加'Method','nearest'参数，否则小字体可能模糊。
2.3 Fresnel 衍射模块：从相位屏到强度图的物理引擎
当选择Fresnel Diffraction模式时，fun_Fresnel.m成为计算核心。它不模拟单条光线，而是将透镜建模为相位调制器：
- 输入平面：
苏庆帅.bmp作为物平面振幅分布； - 透镜作用：在空间频域施加二次相位因子
exp(j*k*(x^2+y^2)/(2*f))； - 传播：采用角谱法（Angular Spectrum Method）或直接 Fresnel 积分，本项目使用后者（更直观）：
function [intensity, phase_map] = fun_Fresnel(u, f, diam, varargin) % Parse wavelength wl = 532e-9; % default green laser p = inputParser; addParameter(p,'wavelength',wl); parse(p,varargin{:}); % Define sampling grid N = 512; dx = 10e-6; % 10um pixel pitch [X,Y] = meshgrid((-N/2:N/2-1)*dx, (-N/2:N/2-1)*dx); % Object plane: amplitude from bitmap, phase=0 obj_amp = double(imread('苏庆帅.bmp'))/255; % Lens phase mask: quadratic phase delay lens_phase = exp(1j*2*pi/wl * (X.^2+Y.^2)/(2*f)); % Propagate to image plane (Fresnel integral) z = abs(u) + f; % object-to-lens + lens-to-image distance H = exp(1j*pi/(wl*z)*(X.^2+Y.^2)) / (1j*wl*z); field_out = ifft2(fft2(obj_amp.*lens_phase) .* fft2(H)); intensity = abs(field_out).^2; phase_map = angle(field_out); end
2.3.1 输出图像命名逻辑：理解Image_after_Phase200.bmp的物理含义
output/目录下 14 张 BMP 文件并非随意生成，其命名严格对应计算流程节点：
| 文件名 | 物理含义 | 生成时机 |
|---|---|---|
Image_before_Phase200.bmp | 物平面振幅图（200mm 物距） | fun_Fresnel输入前 |
Image_after_Phase200.bmp | 经透镜相位调制后的场分布 | obj_amp.*lens_phase后 |
Image_after_abs200.bmp | 200mm 物距下的最终强度像 | abs(field_out).^2输出 |
Image_on_abs.bmp | 几何光学模式下的理想像 | geometric_ray_trace输出 |
提示：
Phase200/Phase400中的数字指物距（单位 mm），而非相位值。若需分析其他物距，只需修改滑块值并点击Run Simulation，新图像会自动覆盖同名文件——因此务必在运行前备份output/目录。
3. 实战调试：三类典型异常的定位与修复方法
3.1 “图像全黑”问题：振幅归一化与数值溢出排查
现象：选择Fresnel Diffraction后，Image_after_abs200.bmp显示纯黑，但Image_before_Phase200.bmp正常。
原因：fun_Fresnel.m中intensity = abs(field_out).^2计算后未归一化，原始强度值可能远小于 1（如1e-12），保存为 uint8 BMP 时被截断为 0。
修复步骤：
- 在
fun_Fresnel.m末尾添加归一化：intensity = abs(field_out).^2; intensity = uint8(255 * (intensity - min(intensity(:))) / (max(intensity(:)) - min(intensity(:)) + eps)); - 验证：在命令行运行
fun_Fresnel(200,200,10)，检查intensity矩阵最大值是否接近 255。 - 若仍偏暗，检查
dx（采样间隔）是否过大——dx=10e-6对应 10μm 像素，若物图细节小于此尺度，需减小dx并增大N（如N=1024），但计算量呈平方增长。
3.2 “像位置偏移”问题：坐标系原点与透镜中心对齐校验
现象：几何光学模式下，Image_on_Phase.bmp中像的位置与理论v=200mm不符，偏差达 ±15mm。
原因：苏庆帅.bmp的像素坐标原点（左上角）未映射到物理坐标系原点（透镜光心）。imresize改变尺寸后，图像中心未对齐光轴。
修复步骤：
- 修改
geometric_ray_trace()中图像放置逻辑：% After resizing, center the image on optical axis pad_h = ceil((N - new_h)/2); pad_w = ceil((N - new_w)/2); img_padded = zeros(N,N,3); img_padded(pad_h:pad_h+new_h-1, pad_w:pad_w+new_w-1, :) = img_out; img_out = img_padded; - 在
my_lens.m的pushbutton1_Callback中，确保N（输出图像尺寸）与fun_my_lens内部一致（默认 512）。 - 验证：用
plot(obj_pos(1),obj_pos(2),'ro'); hold on; plot(img_pos(1),img_pos(2),'bx')可视化物/像坐标，确认二者纵坐标差等于v。
3.3 “GUI 响应迟滞”问题：计算耗时与界面冻结优化
现象：拖动滑块时界面卡死 2~3 秒，Run Simulation按钮无反馈。
原因：fun_Fresnel默认使用N=512的二维 FFT，每次调用耗时约 1.2s（i5-8250U 测试），而 GUIDE 默认在主线程执行回调，阻塞 UI 更新。
优化方案（二选一）：
- 轻量级加速（推荐）：在
pushbutton1_Callback开头添加waitbar防止误操作：hwb = waitbar(0, 'Calculating...','Name','Lens Simulation'); % ... call fun_my_lens ... delete(hwb); - 异步计算（高级）：改用
parfeval将fun_Fresnel提交至后台池：future = parfeval(@fun_Fresnel, 1, u, f, diam, 'wavelength', wl); waitbar(0.5, hwb, 'Propagating...'); intensity = fetchOutputs(future);
注意：异步方案需确保 MATLAB Parallel Computing Toolbox 已安装，且
fun_Fresnel不依赖 GUI 句柄（本项目满足）。
4. 进阶应用：用输出图像反推系统参数与验证衍射极限
4.1 从Image_after_abs200.bmp提取实际分辨率——瑞利判据实证
理论瑞利判据给出最小可分辨间距：δx = 0.61 * λ * f / D，其中D为透镜直径。本项目lens_diam默认 10mm，f=200mm，λ=532nm，理论δx ≈ 6.5μm。验证步骤：
- 用
imshow打开Image_after_abs200.bmp，观察苏庆帅.bmp中最细笔画（如“帅”字撇捺交点）是否分离； - 导入图像至 Image Processing Toolbox，用
regionprops提取连通区域：bw = imbinarize(imread('Image_after_abs200.bmp'), 'adaptive'); stats = regionprops(bw, 'Area','Centroid','BoundingBox'); % 计算相邻小区域中心距 dists = pdist([stats.Centroid]); min_dist = min(dists(dists>0)); - 将
min_dist（像素）乘以dx=10e-6得物理距离，对比理论δx。若实测值 >δx，说明系统受像差或离焦主导；若接近，则衍射极限生效。
4.2 相位图分析：识别球差与彗差的指纹特征
Image_after_Phase200.bmp存储的是复数场的相位angle(field_out)，其等高线形态揭示像差类型：
- 理想透镜：相位等高线为同心圆（二次相位项）；
- 球差：边缘相位超前，等高线在径向呈“S”形弯曲；
- 彗差：相位沿某方向不对称，等高线呈泪滴状。
操作指令：
phase_img = imread('Image_after_Phase200.bmp'); % 注意：此文件存储的是 phase map % 若为灰度图，需还原为 [-pi, pi] 范围 phase_unwrapped = (double(phase_img)/255)*2*pi - pi; figure; imagesc(phase_unwrapped); colorbar; title('Unwrapped Phase Map - Look for S-shaped distortion at edges');4.2.1 快速球差量化：径向相位残差计算
提取相位图中心 100×100 区域，拟合二次曲面a*r^2，计算残差 RMS：
[r,~] = meshgrid((-50:49)*dx, (-50:49)*dx); r2 = r.^2; p_fit = polyfit(r2(:), phase_unwrapped(257:356,257:356)(:), 1); residual = phase_unwrapped(257:356,257:356) - p_fit(1)*r2; rms_aberration = rms(residual(:)); fprintf('RMS spherical aberration: %.3f rad\n', rms_aberration);提示：
rms_aberration > 0.2 rad表明球差显著，此时应减小lens_diam或增加f以抑制边缘光线贡献。
4.3 多物距扫掠自动化：批量生成output/图像集
手动拖动滑块效率低下。以下脚本自动生成物距u=150:50:600的全套图像：
u_list = 150:50:600; f_val = 200; diam_val = 10; for i = 1:length(u_list) u = u_list(i); % Simulate and save [inten, ~] = fun_Fresnel(u, f_val, diam_val, 'wavelength', 532e-9); imwrite(uint8(255*mat2gray(inten)), sprintf('output/Image_after_abs%d.bmp',u)); fprintf('Saved Image_after_abs%d.bmp\n', u); end运行后，用montage命令拼接对比图：
files = dir('output/Image_after_abs*.bmp'); montage({files.name}, 'Size',[3,5], 'BackgroundColor','w');观察像尺寸收缩趋势与背景噪声变化，即可直观掌握景深特性——这正是光学系统设计的核心洞察。
本文还有配套的精品资源，点击获取
