基于频谱分析的倾斜干涉条纹的解调
一、模拟生成干涉条纹图
本文仿真模拟一个典型的干涉条纹,并使用频谱分析技术对该条纹进行解调处理。
假设有这样一幅干涉图,它由一个倾斜角度为 的平行干涉条纹组成,其强度可表示为:
I ( x , y ) = a + b sin [ k 0 ( x cos θ + y sin θ ) ] I\left( x,y \right)=a+b\sin \left[ {{k}_{0}}\left( x\cos \theta +y\sin \theta \right) \right] I(x,y)=a+bsin[k0(xcosθ+ysinθ)] (1)
其中的a与b为常数,这些条纹代表频域中的载波频率或者一个脉冲。在光学干涉仪中,这种倾斜的条纹表示参考光束与物光束之间存在倾斜。仿真得到的倾斜干涉条纹图,如图1所示。
图 1 仿真的倾斜干涉条纹图
二、载频干涉条纹的解调
上述干涉强度分布 I ( x , y ) I\left( x,y \right) I(x,y)可以写成:
I ( x , y ) = a ( x , y ) + b ( x , y ) cos [ ( k x 0 x + k y 0 y ) + ϕ ( x , y ) ] I\left( x,y \right)=a(x,y)+b(x,y)\cos \left[ ({{k}_{x0}}x+{{k}_{y0}}y)+\phi (x,y) \right] I(x,y)=a(x,y)+b(x,y)cos[(kx0x+ky0y)+ϕ(x,y)] (2)
其中, k x 0 , y 0 {{k}_{x0,y0}} kx0,y0是空间载波频率, a ( x , y ) a(x,y) a(x,y)是背景噪声, b ( x , y ) b(x,y) b(x,y)与条纹对比度有关。 ϕ ( x , y ) \phi (x,y) ϕ(x,y) 是要从 I ( x , y ) I\left( x,y \right) I(x,y)计算的干涉相位信息,式(2)可写为:
其中, b c ( x , y ) = ( b ( x , y ) / 2 ) exp [ j ϕ ( x , y ) ] {{b}_{c}}(x,y)=\left( {b(x,y)}/{2}\; \right)\exp \left[ j\phi (x,y) \right] bc(x,y)=(b(x,y)/2)exp[jϕ(x,y)],对式(3)进行傅里叶变换,则有:
其中,~表示频域,假设背景强度 a ( x , y ) a(x,y) a(x,y)与条纹间距相比变化缓慢。 I ~ ( k x , k y ) \tilde{I}\left( {{k}_{x}},{{k}_{y}} \right) I~(kx,ky)的振幅谱将是一个三峰函数,其中 a ~ ( k x , k y ) \tilde{a}({{k}_{x}},{{k}_{y}}) a~(kx,ky)是零峰值, b ~ c {{\tilde{b}}_{c}} b~c与 b ~ ∗ c {{\tilde{b}}^{*}}_{c} b~∗c相对于原点对称放置,如图2(a)所示。
图 2 频谱分析与频谱滤波
接着,消除零频谱和 b ~ c {{\tilde{b}}_{c}} b~c与 b ~ ∗ c {{\tilde{b}}^{*}}_{c} b~∗c的任意一个,得到的新谱不再对称,如图2(b)所示,空域函数不再是实函数而是复函数,因而,式(4)变为
这是一个滤波后的频谱,其以零频为中心。然后对式(5)进行逆傅里叶变换,可得:
I ′ ( x , y ) = b c ( x , y ) = 1 2 b ( x , y ) exp [ j ϕ ( x , y ) ] {I}'\left( x,y \right)={{b}_{c}}(x,y)=\frac{1}{2}b(x,y)\exp \left[ j\phi (x,y) \right] I′(x,y)=bc(x,y)=21b(x,y)exp[jϕ(x,y)] (6)
通过下式即可得到包裹相位:
ϕ ( x , y ) = arctan { Im [ b c ( x , y ) ] Re [ b c ( x , y ) ] } \phi (x,y)=\arctan \left\{ \frac{\operatorname{Im}\left[ {{b}_{c}}(x,y) \right]}{\operatorname{Re}\left[ {{b}_{c}}(x,y) \right]} \right\} ϕ(x,y)=arctan{Re[bc(x,y)]Im[bc(x,y)]} (7)
得到的包裹相位如图3所示:
图 3 包裹相位图
三、相位解包裹
采用最小二乘法对其进行相位解包裹,解包裹相位如下图所示:
图 4 解包裹相位
上述解包裹相位以弧度为单位,可根据波长将其进行转化为微米为单位,转化后的相位如下图所示。
图 5 三维深度图
四、资源获取
本文提供上述技术源程序,感兴趣的可从以下链接处获取:
https://download.csdn.net/download/qq_36584460/84866400
基于频谱分析的倾斜干涉条纹的解调相关推荐
- 基于EEMD方法的降噪与解调在轴承故障诊断中的应用实证研究
基于EEMD方法的降噪与解调在轴承故障诊断中的应用实证研究 摘要 1 引言 2 EEMD简介 2.1 EMD简介 2.2 EEMD简介 3 基于阈值处理的EEMD降噪方法 3.1 相关系数 3.2 基 ...
- matlab显示2dpsk误码率,基于MATLAB的2DPSK调制与解调系统的分析.doc
您所在位置:网站首页 > 海量文档  > 计算机 > matlab 基于MATLAB的2DPSK调制与解调系统的分析. ...
- fsk调制解调实验报告 matlab,基于MATLAB的-FSK调制与解调-通信原理实验
<基于MATLAB的-FSK调制与解调-通信原理实验>由会员分享,可在线阅读,更多相关<基于MATLAB的-FSK调制与解调-通信原理实验(2页珍藏版)>请在人人文库网上搜索. ...
- matlab进行动力吸振器设计,干摩擦动力吸振器的MATLAB仿真及其基于频谱分析的优化设计...
干摩擦动力吸振器的MATLAB仿真及其基于频谱分析的优化设计 第 2 l 卷第 2期 振动与冲击 J 0URNA L OF V I BR ATI ON AN D S HO CK 干摩擦 动力 吸振器 ...
- 【基于matlab的mqam调制与解调系统】
提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档 文章目录 前言 一.正交幅度调制原理 二.QAM的解调和判决 三.16QAM调制解调系统的实现与仿真 总结 前言 提示:这里可以添加本 ...
- matlab调制解调开题报告,基于MATLB的数字调制与解调设计与仿真开题报告.doc
基于MATLB的数字调制与解调设计与仿真开题报告 毕业(设计)论文 开 题 报 告 系 别 电子信息工程系 专 业 电子信息工程专业 班 级 122102 学生姓名 朱椰子 学 号 111861 指导 ...
- 信号的解调与调制matlab,基于MATLAB常用数字信号调制与解调
内容简介: 本科毕业设计 基于MATLAB常用数字信号调制与解调 数字通信|MATLAB|ASK|PSK|FSK|调制与解调|毕业设计 文件格式:word+PPT 本科毕业设计 基于MATLAB常用数 ...
- 信号的调制与解调matlab仿真,基于MATLAB对信号调制与解调的仿真
基于MATLAB对信号调制与解调的仿真 2.3 2PSK的基本原理和调制解调实现 数字相移键控,记作PSK(Phase shift-keying ),二进制相位键控记作2PSK. 2PSK就是根据基带 ...
- matlab 2dpsk调制频谱,基于MATLAB的2DPSK调制与解调系统的分析
摘要 MATLAB集成环境下的Simulink仿真平台,设计一个2DPSK调制与解调系统.用示波器观察调制前后的信号波形;用频谱分析模块观察调制前后信号频谱的变化;加上噪声源,用误码测试模块测量误码率 ...
最新文章
- 使用条件注释完成浏览器兼容
- 机器学习——图解SVM中gamma和c参数的作用
- ContextMap详解
- Java获得不同格式的系统时间
- java 身份证地址提取籍贯_函数 | 身份证提取籍贯信息
- 1.1 《数据库系统概论》之数据库系统概述(数据、数据库、数据库管理系统、数据库系统、数据库的产生和发展)
- 网站都变成灰色,一行代码就搞定了!
- URL Schemer
- 菜鸟教程C#高级教程
- NLP机器翻译任务中,如何用Bleu score评价翻译质量(学习心得)
- python 拼音库_python有没有拼音库python进阶之socket详解
- tf2 + keras学习
- Linux下修改虚拟机的root密码
- 计算机丢失quartz.dll什么意思,电脑系统提示没有找到quartz.dll该怎么办
- treeShaking和sideEffects详解
- [转]秒杀系统架构分析与实战
- 黑马程序员——Java基础(十一)之GUI
- iOS常用数学函数(公式)
- 示波器的探头补偿方法
- Android中启动其它应用的几种方式