信息光学第三章第三讲
信息光学第03章

输入是指施加于系统的作用,称为系统的输入激励(excitation);输出是要求系统完成的功能,称为系统的 输出响应(response)。可见,系统的特性决定对某一输入激励会产生什么样的输出响应。当研究一个系统的 性质时,不必过多地关心系统内部的结构,只需知道其输入端和输出端的性质就行了。 分析一个系统,首先要对系统建立数学模型,然后运用数学方法进行求解,最后又回到系统,对结果 做出物理解释,并赋予其物理意义。所谓系统的模型是指系统物理特性的抽象,以数学表达式或具有物理 特性的符号图形来表征系统特性。系统的分类比较复杂,从数学模型的差异来看,可划分为: (1) 连续系统和离散系统:输入和输出均为连续信号的系统称为连续系统,输入和输出均为离散信号 的系统称为离散系统。 (2) 线性系统和非线性系统:线性系统是指具有线性特性的系统。所谓线性(linearity) 特性是指齐次性 与叠加性。若系统输入增加 k 倍,输出也增加 k 倍,这就是齐次性(homogeneity)。若有几个输入同时作用 于系统,而系统总的输出等于每一个单独作用所引起的输出之和,这就是叠加性(superposition property)。 系统同时具有齐次性和叠加性便呈现线性特性。一般线性系统性必须具有以下特性:① 分解性 (decomposition);② 零输入线性;③ 零状态线性。凡不具备上述特性的系统则为非线性系统。 (3) 空间不变系统和空间变系统:只要初始状态不变,系统的输出仅取决于输入而与输入的起始作用 点无关,这种特性称为空间不变性。具有空间不变特性的系统为空间不变系统(space invariant system),不 具有空不变特性的系统为空间变系统(space varying system)。 (4) 因果系统和非因果系统:因果系统(causal system)是指其响应不会超前激励的系统,非因果系统 (noncausal system) 是指响应能领先于激励的系统,它的输出取决于输入的将来值。 为了用简洁的语言来分析物理系统,最常用的方法是寻找一个数学模型,使其在数学意义上能恰当地 描述该系统的性质和状态。在傅里叶光学中,常常采用一种算符把光学系统的激励与对此产生的响应联系 起来,系统的作用就是完成数学上的某种变换或运算。如图 3.1.1 所示,算符 L{} 表示系统的作用,激励 函数 f ( x, y ) 通过系统后变成相应的响应函数 g ( x, y ) ,两函数之间满足下列关系:
信息光学第三章

2 d o di
P( x, y) exp{ j
2 [(xi ~ xo ) x ( yi ~ yo ) y ]}dxdy d i
h( xi ~ xo , yi ~ yo )
1 2 do di
P( x, y) exp{ j
2 [( xi ~ xo ) x ( yi ~ yo ) y]}dxdy di
(1)当 d 2 f,即后焦面为观察平面
U ( x, y) jf e
j k d (1 1 )( x 2 y 2 ) 2f f
其中:
1 1 1 d1 d 2 f
j
dx0 dy0
U 0 ( x0 , y0 )e
j 2 ( x x y0 y ) f 0
2 d o di
1
P( x, y) exp{ jk[(
P( x, y ) exp{ j
xi xo y y ) x ( i o ) y ]}dxdy di d o di d o
(3.34)
2 [( xi Mxo ) x ( yi Myo ) y ]}dxdy d i
(3-4-11)
三、相干照明下衍射受限系统的成像规律
U i ( xi , yi ) U o (xo , yo )h( xi ~ xo , yi ~ yo )dxo dyo
1 2 M
~ xo U o ( M
~ yo , )h( xi ~ xo , yi ~ yo )d~ xo d~ yo M
( x2 y 2 ) jk 2q
t ( x0 , y0 )e
信息光学_第三章

P1
P2
s o1 o2
s’
p
q
U1(x,y) t(x,y) U1‘(x,y)
U
'
1
(
x,
y)
U1
(
x,
y)t
(
x,
y)
透镜的复振幅透过率为:
t(x, y) U1' (x, y) U1(x, y)
在傍轴近似下,单色点光源S发出的发散
球面光波在P1平面上造成的光场分布为
U1(x,
y)
Aexp(
jkp)
➢ 孔径光栏、入瞳、出瞳由系统元件参数及相对位置决定。 ➢ 对整个光学系统而言,入瞳和出瞳保持物像共轭关系。由入瞳限 制的物方光束必定能全部通过系统,成为被出瞳所限制的像方光束
2)“黑箱”模型
系统成像分三部分:
1、物平面到光入瞳平面2、入 瞳平面到出瞳平面3、出瞳平
面到像平面
➢ 在第1、3部分:光波传播可按菲涅耳衍射处理 ➢ 在第2部分,在等晕条件下,可把它当作一个黑箱来处理,黑 箱两端分别是入瞳和出瞳,只要能够确定黑箱两个边端的性质,
y0 ;
x,
y)
P(x,
y) exp
jk
x2 2f
y2
dUl
( x0 ,
y0 ;
x,
y)
透镜后表面xi,yi平面: 再次运用菲涅耳衍射公式:
h(x0
,
y0 ;
xi
,
yi
)
exp( jkdi
jdi
)
dUl(x0
,
y0
;
x,
y)
exp
jk
(xi
x)2 ( 2di
信息光学导论第三章

线性系统概论◆引言在信息光学系统中光学装置被看成收集、传递或变换信息的系统。
一个光学系统的理想成像,就是将无空间的物体信息传递、变换物空间,在像面上形成不变的物体的像。
这样的理想光学系统显然是一线性系统。
虽然实际光学成像系统由于不可避免的存在相差,总会产生失真,是非线性的,但在把研究的问题看成线性的而不会引起明显误差,或只在某个小范围满足现行性质时,就可以将其当作现行未提来处理。
所以线性系统理论与傅立叶分析方法一样,是研究信息光学中成像系统和信息处理系统的重要理论基础。
本章主要介绍线性系统特别是空间不变线性系统的定义、特点和分析方法。
3.1线性系统的基本概念◆系统及其分类所谓系统,是指一组相互关联的事物构成的总体。
这样的系统可分为物理系统和非物理系统。
这里仅讨论物理系统。
如图所示一个物理系统,它是这样的装置,当对其作用一个激励时,他就产生一个响应。
从数学上着眼,很多现象都可抽象为使函数)(x f 通过一定的变换,形成)(x g 函数的运算过程.这种实现函数变换的运算过程称为系统.这种意义下的系统,既可以是特定功能的 元器件组合,例如电子线路、光学透镜组等,也可以是与实际元件无关的物理现象,如光学系统,通讯系统,管理系统和指挥系统等。
系统论的引入,使得我们在研究一个光学系统时,所关心的是系统对于给定的激励产生什么样的响应,而不去考虑系统内部的具体结构和具体工作原理。
线性系统理论是从总体上研究系统输入输出之间的对应关系和他们的共同特性。
◆线性系统的定义及其算符表示假设一个激励)(1x f 作用于某系统产生的响应为)(1x g ,而激励)(2x f 作用于某系统产生的响应为)(2x g ,用符号表示为)()(),()(2211x g x f x g x f →→如果系统满足可加性)()()()(2121x g x g x f x f +=+和奇次性(均匀性))()(),()(22221111x g c x f c x g c x f c →→则这样的系统为线性系统。
信息光学课件第三章

相干系统的点扩散函数 可看成是复振幅透过率 的光瞳被 半径为di的球面波照明后所得的分布。 称广义光瞳。 就是广义光瞳 的傅里叶变换。
相干传递函数定义为相干点扩散函数的傅里叶变换
由
得
(无像差)
有像差系统的通频带没有变化,截止频率也没有变化,但在通频 带内引入了与频率有关的位相畸变,使像质变坏。 非相干光照明下强度点扩散函数仍然是相干点扩散函数模的平方 但峰值减小。 Strehl Ratio
的频谱函数(相干传递函数)H(ξ ,η )
描述系统的变换特性更为方便。
3.3.1相干传递函数
相干成像系统的物像卷积关系
是几何关系理想像的复振幅分布。ĥ是系统的脉冲响应。 从频域上看,对上式进行傅里叶变换,可得到系统对各种频率成 分的传递特性。
系统的输入频谱 输出频谱 相干传递函数 CTF
已知
说明相干传递函数等于光瞳函数,只是将空域坐标变换为频域坐标 (-λ diξ ,-λ diη ),通常光瞳都具有中心对称性,正负号无关紧要, 忽略负号后取
因hI是实函数,H是厄密型的,即
因此模是偶函数
辐角是奇函数
3.6相干与非相干成像系统的比较
各有优缺点。 3.6.1截止频率 OTF 的截止频率是CTF的2倍。但是 OTF是随空间频率增大而降低的。而CTF 是在空间频率小于某值前均为1,大于某 值时突变为0。
相干传递函数
3.6.2像强度的频谱
利用卷积定理和自相关定理得到像强度频谱
D为出瞳直径。
相干照明时,两点源产生的艾利斑按复振幅叠加。因而各点的 相位关系对强度分辨影响很大。
Φ =0,两点源位相相同,I(x)没有凹陷两点完全不能分辨。 Φ =π /2 与非相干光完全相同。 Φ =π 时,两点源位相相反。 两点源能否分辨与点源位相有关。
信息光学原理第3章

1
焦面场是透镜前端场的傅里叶变换(空间频谱)。 如上图所示,距离透镜前端有一物体,其透过率为t(x0,y0)。若用振幅为A 的平面波垂直照明物体,则物体的透射光场为:
U0 x0 , y0 A t x0 , y0
根据角谱理论,透镜前端场的角谱为:
F U1 x, y F U 0 x0 , y0 H f x , f y
U l x, y U x, y
在傍轴近似下,忽略透镜对光波振幅的影响,紧靠透镜前后的平面上产生的 复振幅分布为
k U l x, y A exp jkd 0 exp j x 2 y 2 2d 0
k U l x, y A exp jkd i exp j x 2 y 2 2d i
x 1
2
y2 R12
x 1
x2 y 2 1 2 R2 2 2 x y 2 x, y 02 R2 R2 2 x 2 y 2 02 R2 1 1 2 仅考虑傍轴光 R2
f
f
f
j f
2f
2
2
f
f
k 2 U 2 x, y exp j x 2 y 2 exp j xx f yy f dxdy 2f f
?
3.2 透镜的傅里叶变换性质
后焦面上的场分布为
3.1 透镜的位相调制作用
则透镜复振幅透过率表示为:
k A exp jkdi exp j x 2 y 2 U x, y 2d i tl x, y l U l x, y k 2 2 A exp jkd 0 exp j x y 2d 0
信息光学教程李俊昌第三章

然而 , 上面对连续函数被无穷δ 序列取样离散后的频谱研究只是一个理论结果 , 因为实际
频域均是周期离散函数的傅里叶变换问题 , 还要将 离 散 函 数 截 断 及 延 拓 才 能 满 足 要 求 . 因 此,
r x)= L x(
( ) ] 得到具有 Nx 个点的离散分布 [ 图3 1 2 e 1
{
是, 对于实际给定的衍射问题 , 能够直接从傅里叶变 换 求 出 解 析 表 达 式 的 函 数 非 常 有 限 , 在研 究实际问题时 , 不得不将函数按一定规律在二维空间进行取样及延拓 , 变成该函数的周期离散 分布作离散傅里叶变换 . 然而 , 离散傅里叶变换的计算仍然十分繁杂 , 如果没有计算机 , 事实上
由图可见 , 由于矩形窗函数的 频 谱 RL (f 卷积运算的结 x ) 具 有 较 大 的 起 伏 变 化 的 旁 瓣, x
( ) ] 函数 [ 图3 1 2 1 . g 空域及频域离散函数均以 Nx 为周期 , 我们只要分别知道一个周期内的离散值或样本点便
菲涅耳衍射积分可以表示为傅里叶变换及卷积两种形 式 , 存在一次快速傅里叶变换及快速卷
㊃4 8㊃
信息光学教程
上不可能作取样点为无限多的数值计算 . 并且 , 由于离散傅里叶变换事实上讨论的是在空域及 ( ) ] ( ) ] 将空域非周期的离散函数 [ 图3 先通过下述矩形窗函数 [ 图3 截断 : 1 2 c 1 1 2 d 1 / 1, -Tx/ 2< x < L 2 x -T x 0
n=- ∞
∞
( 由于梳状函数δ 可以表示为傅里叶级数 x)为周期 Tx 的δ 函数 , T x ( δ x)= T x 式中 , j=
( ) 3 1 1
Ak = -1 ,
信息光学(傅里叶光学)Chap3-1

x
, f y ) exp[ j 2p ( f x x f y y )]df x df y
即: 把U(x,y)看作频率不同的复指数分量的线性组合, 各分 量的权重因子是A(fx, fy).
A( f x , f y ) U ( x, y ) exp[ j 2p ( f x x f y y )]dxdy
fx
X
l
;
fy
Y
s单色平面波 在xy 平面的复振幅分布可以表示为
U ( x, y ) A exp[ j 2p ( f x x f y y )]
#
光波的数学描述
平面波的空间频率-信息光学中最基本的概念
练习 1
单位振幅的单色平面波, 波矢量k与x轴夹 角为30, 与y轴夹角为60. (1)画出z = z1平面上间隔为2p的等相线族, 并求出Tx、 Ty、T 和fx 、fy和 f。 (2)画出y = y1平面上间隔为2p的等相线族, 并求出Tx、 Tz 和fx 、fz.
§3-1 光波的数学描述
单色光波场的复振幅表示
将光场用复数表示,有利于简化运算
u(P,t) = a(P)cos[2pnt - j(P)]} = e{a(P)e-j[2pnt -j(P)] }
复数表示有利于 = e{a(P) e jj(P). e -j2pnt } 将时空变量分开
光场随时间的变化e -j2pnt不重要: n ~1014Hz, 无法探测 n为常数,线性运算后亦不变 对于携带信息的光波, 感兴趣的是其空间变化部分. 故引入复振幅U(P):
l
l
l
l
cosa cos b 称为xy平面上复振幅分布的 A( , )
l
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
•入射平面 波
•近场衍射 •Fresnel
•远场衍射 •Fraunhofer
•小 结
• 菲涅耳衍射公式(空域和频域) • 夫琅和费衍射公式(空域和频域) • 爱里斑
•C
•A’
•阿贝成像原理将成像过程分为两步:
•第一步“分频”; •第二步“合成”.
•阿贝成像真正意义:提供了新的频谱语言描述信息,启发人们用改变
•
频谱的手段来改造信息,即信息光学处理基础
•光信息处理的典型实验:
•物平面
•频谱面
•像平面
•五、菲涅耳衍射与夫琅禾费衍射的定量关系
•很容易实现
•菲涅耳衍射区:
• 矩形孔复振幅透过率函数对应的频谱:
• 采用单位振幅的单色平面波垂直照明孔径,即 : • 孔径透射的复振幅分布为:
• 观察面上得到的夫琅禾费衍射场的复振幅分布为:
• 观察面的强度分布为:
• 矩形孔衍射的MATLAB仿真:
•(1)实验参数设定:
•L=0.2;
%side length (m)
•M=250;
•别 来无
恙
•Contemporaries, but not collaborators (nor competitors).
•一、 菲涅耳衍射与夫琅和费衍射的定性区别
Fraunhofer diffraction:
both incident and diffracted waves may be considered to be planar (i.e. both S and P are far from the aperture)
•中央亮斑称为爱里斑
,其半径为:
•衍射场三维图
•衍射场强度三维图
•爱里斑的应用
•不可分辨
•Rayleigh •分辨极限
•小孔
•可分辨
•不同距离处圆孔的衍射图像
•不同距离处圆孔的衍射图像
•不同距离处圆孔的衍射图像
•From Fresnel to Fraunhofer diffraction
•w2=0.2e-3; %inner radius
•lambda=0.633e-6;%wavelength
•z=50;
%prop distance
•k=2*pi/lambda; %wavenumber
•lz=lambda*z;
•(2)观察平面强度:
•%irradiance •I2=(1/(lambda*z))^2.*((w1^2*jinc(w1/(lamb da*z).*R))(w2^2*jinc(w2/(lambda*z).*R))).^2; • •figure(1) %irradiance image •imagesc(x,y,nthroot(I2,3)); •xlabel('x (m)'); ylabel('y (m)'); •colormap('gray'); •axis square •axis xy; • •figure(2) %x-axis profile •plot(x,I2(M/2+1,:)); •xlabel('x(m)'); ylabel('Irradiance');
•(1)实验参数设定:
•L=0.2;
%side length (m)
•M=250;
%# samples
•dx=L/M; %sample interval
•x=-L/2:dx:L/2-dx; y=x; %coords
•[X,Y]=meshgrid(x,y);
•R=sqrt(X.^2+Y.^2);
•
•w1=1e-3; %outer radius
▪ 这些衍射斑和图象的空间频率一一对应,后焦面就是图象的傅立 叶频谱面,称为傅立叶面,夫琅和费衍射斑称为谱斑。这就是现
代光学对夫琅和费衍射的新认识。
•y
•y’
•F
•x
•x’
•衍射屏 •光学图象
•透镜 •频谱分析器
•夫琅禾费衍射场 •傅立叶频谱面
•出的球面波,经透镜 折射后,各自会聚到它们的像点A,B,C.
•S
•P
Fresnel diffraction: occurs when either S or P are close enough to the aperture that wavefront curvature is not negligible
•六、夫琅和费衍射计算实例及MATLAB仿真
•1、矩形孔衍 射• 矩形孔复振幅透过率:
•S
•P
Fresnel diffraction: occurs when either S or P are close enough to the aperture that wavefront curvature is not negligible
•From Fresnel to Fraunhofer diffraction
•入射平面 波
•近场衍射 •Fresnel
•远场衍射 •Fraunhofer
•二、衍射公式的初步近似
•旁轴近似 :
•则 :
•,所以分母中的
•三、菲涅耳衍射的空域分析
•1、菲涅耳衍射的卷积形式
•指数中的r不能用z近似,必需采用更高一级的近似:
•令
•脉冲响应
:
•此即菲涅耳 衍射的卷积形
式
•2、菲涅耳衍射的傅里叶变换形式
•衍射场三维 图
•衍射强度三维 图
•2、圆孔衍射
• 圆形孔复振幅透过率:
• 圆形孔复振幅透过率函数对应的频谱:
• 采用单位振幅的单色平面波垂直照明孔径,即 : • 孔径透射的复振幅分布为:
• 观察面上得到的夫琅和费衍射场的复振幅分布为:
• 观察面的强度分布为:
• 圆形孔衍射的MATLAB仿真:
•The propagated wave follows the periphery of the wavelets.
•Huygens, just add the wavelets considering interference!
•衍射双雄: Fresnel and Fraunhofer
•别来 无恙
•对于菲涅耳衍射公式:
•指数展开并移项 :
•若观察平面离开孔径平面的距离z进一步增大,使其不仅满足 菲涅耳近似条件,而且还满足
•即 : •则,平方相位因子(上述红色)在整个孔径上近似为1。
•于是,菲涅耳衍射公式可以进一步简化为夫琅和费衍射公式:
•上式表明:•观察平面上的光场分布正比于孔径平面透射光场 分布的傅里叶变换。
•观察平面上衍射图样的强度分布是:
•其中
信息光学(或傅里叶光学)的基本思想
▪ 理想夫琅和费衍射系统起到空间频率分析器的作用.当单色光波 入射到待分析的图象上时,通过夫琅禾费衍射,一定空间频率的信 息就被一特定方向的平面衍射波输送出来。
▪ 这些衍射波在近场彼此交织在一起,到了远场它们彼此分开,从而 达到分频的目的.常用远场分频装置是透镜,将不同方向的平面 波汇聚到后焦面上不同的点上,形成一个个衍射斑。
•像面 •透镜
•图中光线不同 的颜色表示发自 不同的物点.
•
• ••
•
•
•现代光学对透镜成像的新观点
•阿贝成像原理: 物是一系列不同空间频率的集合。入射光经物平面发生 夫琅和费衍射,在透镜焦面(频谱面)上形成一系列衍射光斑,各衍射光 斑发出的球面次波在相面上相干叠加,形成像。
•F
•C’ •A
•B
•B’
信息光学第三章第三讲
•主要内容
• 菲涅耳衍射的空域形式 • 菲涅耳衍射的频域形式 • 夫琅和费衍射的空域形式 • 夫琅和费衍射的MATLAB仿真
•Huygens-Fresnel 原理
•波前上的每一个面元都可以看做是一个次级扰动中心,它们能产生球面子波, 后一时刻的波前位置是所有这些子波波前的包络面。----《论光》C.Huygens
•z=7;
%prop distance
•k=2*pi/lambda; %wavenumber
•lz=lambda*z;
•(2)观察平面强度:
•%irradiance •I2=(4*wx*wy/(lambda*z))^2.*(sinc(2*wx/(la mbda*z).*X).*sinc(2*wy/(lambda*z).*Y)).^2; • •figure(1) %irradiance image •imagesc(x,y,nthroot(I2,3)); •xlabel('x (m)'); ylabel('y (m)'); •colormap('gray'); •axis square •axis xy; • •figure(2) %x-axis profile •plot(x,I2(M/2+1,:)); •xlabel('x(m)'); ylabel('Irradiance');
%# samples
•dx=L/M; %sample interval
•x=-L/2:dx:L/2-dx; y=x; %coords
•[X,Y]=meshgrid(x,y);
•
•wx=0.1e-3; %x half-width
•wy=0.05e-3; %y half-width
•lambda=0.633e-6; %wavelength
•夫琅禾费衍射区 :
•From Fresnel to Fraunhofer diffraction
•入射平面 波
•近场衍射 •Fresnel
•远场衍射 •Fraunhofer