天津大学 小波分析及应用第三次作业 小波变换在thz时域光谱去噪方面的 .doc

天津大学 小波分析及应用 第三次作业 小波变换在THz时域光谱去噪方面的应用 学生: 苏晓强 学号: 1012202009 学院:精密仪器与光电子工程学院 专业: 光学工程（博士） 日期: 2013年6月3日  小波变换在THz时域光谱去噪方面的应用 1引言 小波分析具有时频分析，多分辨率分析的特点，传统分析信号的理想工具是傅里叶变换，但是对于非平稳信号，傅里叶变换变得有些无能为力，因为傅里叶变换是对整个时域到频域的变换，只能分析频域信息，而对于时域则无法描述。而信号的时频域性质正是非平稳信号的关键所在。小波变换能够进行时频分析，通过小波母函数的伸缩和平移运算能够对信号进行多尺度分解，然后从各尺度上处理和分析信息，被誉为“数学显微镜”，在调和分析的发展史上具有里程碑式的重要作用。太赫兹频段是指处于毫米波与红外光之间的电磁辐射区域，频率从 0.1THz到 10THz。太赫兹波在电磁波谱中的位置是属于微波和红外线的过渡，对科学和研究都有很重要的意义，THz时域光谱技术（THz-TDS）[1-2]是研究THz科学的主要手段之一，因其大带宽、高信噪比、高时间分辨率和相位相干探测的优势，该技术得到了越来越广泛的应用。THz-TDS本质上是相位相干电磁辐射探测技术,它能直接探测THz波与物质相互作用的电场信息,即同时测得电磁波的振幅和相位。因此,可精确测量物质的色散和吸收特性。这种新兴的光谱技术,在气体、液体、电介质、半导体、超导介质,及生物样品等光谱研究中显示出独到的技术特点。然而，太赫兹波在检测中对环境的要求比较苛刻，因为水蒸气对太赫兹波有着较为强烈的吸收，很容易对实际检测的样品信号引入噪声，使得样品对太赫兹波真实的吸收和响应无法得到准确的分析，目前使用的THz时域光谱技术主要是基于Fourier变换，Fourier变换是对时域上整体的分析，无法对THz时域或频域上的任意细节进行分析，如果利用小波分析对时域上某个细节的信号进行分析和变换，并对它进行去噪处理，从而大大优化和改善对实际样品信号的分析结果。 2小波变换对THz时域光谱进行去噪 利用小波分析的方法对信号进行去噪，主要是利用了信号和噪声在小波分解 后有不同的特征表现。分解后的噪声小波系数会随着尺度的增加衰减的很快，噪声都集中在粗尺度上，并且在各尺度间也没有相关性可言。但是真实信号的幅值随着尺度的增加几乎是没有什么变化。小波去噪，也就是先把信号分解为小波系数，然后对分解出来的小波系数做对应处理，实质是希望减小或者完全将噪声产生的系数去除掉，同时尽量做到不损失真实的信号系数，然后用处理过的系数重构回复原信号。 小波去噪的特点为： (1) 首先，它不是平滑。平滑是去除高频信息而保留低频部分，而小波去噪是试图去除所有的噪声，保留所有信号，并不考虑它们的频率范围。 (2) 其次，它是小波变换域对小波系数进行非线性处理。 (3) 滤波过程一般有三个步骤实现：①小波变换；②对小波系数进行非线性处理，以滤除噪声；③小波逆变换。 (4) 小波去噪是一种非参数方法。在参数方法中，必须要对一个特定的先验模型进行参数估计。 在THz时域光谱测量中，信噪比的重要性是显而易见的，它决定了该测量系统的动态范围。吸收峰可以分辨的必要条件是该吸收峰对太赫兹光谱的调制度必须大于系统的测量噪声；只有当系统的信噪比与太赫兹波经过样品损耗的乘积大于吸收峰对太赫兹波调制度的倒数时，该吸收峰才能被太赫兹时域光谱辨认出来。 这里采用了空域相关去噪的方法[3-5]，对采集的THz时域光谱信号进行了去噪的处理和分析。空域相关方法的原理：小波分解后，信号在各个尺度上的突变点都会有一个比较大的峰值，而且在每个尺度上出现的位置几乎都是相同的。但是分解后的噪声却随着尺度的增加，能量有明显的衰减现象。所以，对两个相邻的尺度上分解的小波系数相乘来做相关运算，这样能够使得信号的边缘更加的清晰，提高了定位的精度，重要的特征点都能提取出来，这样的同时也就抑制了噪声。 信号经小波变换后，其小波系数在各尺度上有较强的相关性，因此，如果分解的每个尺度上对应点处都有边缘点，则可以判断这个点是信号系数，反之则是噪声系数，予以剔除。这样处理后的小波系数重建，就将噪声去除了。 设分解的最大尺度为 J。Wf(j,n)表示尺度 j 上位置 n 处含噪信号 f 的离散小波变换值，取相邻尺度的变换值进行相关计算，定义： 式中 Corr2(j,n)是尺度 j 上点 n