MATLAB仿真光的衍射（五）

最新推荐文章于 2023-06-04 12:06:13 发布

matlab练习生

最新推荐文章于 2023-06-04 12:06:13 发布

阅读量1.1k

点赞数

分类专栏： # MATLAB光学建模仿真文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/fendoudecheng/article/details/130404582

版权

MATLAB光学建模仿真专栏收录该内容

25 篇文章 97 订阅

订阅专栏

直边的菲涅尔衍射。

之前考虑的衍射均为光屏离物体（如小孔）较远的情况，也就是夫琅禾费衍射情况，

对于距离较近的菲涅尔衍射，我们很难得到确定的解析解。

不过，特殊的菲涅尔衍射情况比如：圆孔、圆盘、直边是可以解出的。

这里以直边的菲涅尔衍射为例，计算结果。

光强的表达式如下：

$I=\frac{1}{2}I_0\left\{\left[\frac{1}{2}-C(v_0)\right]^2+\left[\frac{1}{2}-S(v_0)\right]^2\right\}$

其中：

$C\left (v_0 \right) =\int_ {0} ^{v_0} {cos(\frac{\pi t^2}{2}) dt}$

$S\left (v_0 \right) =\int_ {0} ^{v_0} {sin(\frac{\pi t^2}{2}) dt}$

$v_0=\sqrt{\frac{2}{\lambda d}}x$

菲涅尔正弦和余弦积分可以直接调用matlab中的函数进行计算。

clear
close all
ScreenX = 1048;
ScreenY = 350;
lambda = 632.8e-9;
d = 10;
x = 0.02;

upsilon0 = linspace(-sqrt(2/(lambda*d))*x,sqrt(2/(lambda*d))*x,ScreenX);
I = 1/2*((1/2-fresnelc(-upsilon0)).^2+(1/2-fresnels(-upsilon0)).^2);
II = repmat(I,[ScreenY 1]);

subplot(2,1,1)
imshow(II);
subplot(2,1,2)
plot(linspace(-x/2,x/2,ScreenX),I)
grid on
ylabel('相对光强');
xlabel('x/m');

结果：