一种强噪声背景下的拉曼信号重构方法，从信号重构的角度，采用基于高斯分析字典的信号稀疏表达从强噪声背景下的原始测量信号中重构出有用的拉曼信号，重构的信号能同时实现噪声去除和基线校正；同时，通过计算字典与残差的内积，能构实现加快信号处理速度和减少内存占用；采用本算法可实现强背景噪声下的拉曼信号有效重构，拉曼峰位移定位准确，峰面积衰减小，信噪比改善明显，且算法处理效率高，可以用于高时空分辨率的拉曼测量系统信号预处理。
【技术实现步骤摘要】
本专利技术涉及信号处理领域，特别是。
技术介绍
拉曼光谱技术是一种强有力的分子结构研究工具，但是，拉曼散射效应非常微弱， 其散射光强度约为入射光强度的IO 8-IO 6。为了获得高时空分辨率，常常采用很小的积分 时间（如0.1 s)情况下，此时，微弱的有效拉曼信号会淹没在强噪声背景下。 常用的信号去噪方法有移动窗口平均平滑法、EMD (Empirical Mode Decomposition经验模态分解），小波变换法等。 移动窗口平均平滑法，基于信号和噪声统计特性之间的差异，其基本假设是噪声 为零均值噪声，通过对原始信号取平均值达到提高信噪比的目的，是消除噪声最常用方法。 通过选择一个宽度为奇数2 ω +1的平滑窗口，以中心波长点k为参照点从左到右移动窗口， 将窗体覆盖区域内所有测量的平均值代替中心波长点对应的测量值，直至完成所有点的平 滑，该方法平滑窗口宽度ω影响平滑结果，宽度太小，平滑效果不佳，宽度太大，则平滑掉 特征峰信息，造成光谱失真；且存在边界问题。对于宽度为2 ω+1的平滑窗口，光谱左右两 端各有ω个点不能被处理； EMD将信号自适应的分解成有限个由高到低、时间尺度由小到大的頂F(instrinsi cmode function本征模态函数）分量及单一趋势项，但是缺少严格的数学根据、停止准则 难以确定，且存在端点效应、模态混叠效应，算法稳定性不好。 小波变换法经由小波基的伸缩和平移，达到信号时频分析局域化，能同时保留信 号的时域特征和频域特征，在适当的尺度内，非平稳信号的有效成分能暴露出与噪声不同 特点，运用信号和噪声在多尺度范围内的各种传递特性，可在噪声干扰背景提取出有效信 号，但是小波基的选取对去噪结果干扰较大，同时在低信噪比的情况下，效果较差。 移动窗口平均平滑法、EMD、小波变换法等传统去噪方法处理低信噪比拉曼信号时 效果并不理想，不能实现强背景噪声下拉曼信号的有效还原。
技术实现思路
本专利技术的主要目的在于克服现有技术中的上述缺陷，提出一种强噪声背景下的拉 曼信号重构方法，其算法处理效率高、信噪比改善明显、同时实现噪声去除和基线校正。 本专利技术采用如下技术方案： ，其特征在于，预先定义初始辅助内积为 a °= D Ty，y为原始拉曼测量信号，D为高斯字典；辅助矩阵G = DtD ;初始平方误差ε ° = yTy ;目标平方误差为ε ;残差r = (y-D γ)，γ为y的稀疏表达；包括如下步骤： 1)初始化：令迭代过程所选原子序列I:=()，对辅助矩阵Glil的Cholesky分解 L赋初值L: = ,辅助内积α : = α°,迭代过程中间变量δ°: =〇,初始迭代序号η:= 1 ; 2)计算高斯字典原子和残差信号的辅助内积α，选出与当前残差信号最匹配的 原子索引判断当前迭代序号是否大于丨，若是，则进入步骤 3);否则，跳转 到步骤4); 3)通过引入中间变量更新L ; 4)将选出的最匹配的原子索引 5)根据式LlZc = W，《:；'为D/y即为当前更新的原子序列组成的对应的高斯字典 D的子矩阵D1与被测原始拉曼信号y的内积，求解c，c为更新后的当前高斯字典对原始拉 曼测量信号的稀疏表达，即γΙ; 6)根据更新后的原子序列计算第η次迭代后的迭代中间变量β =G1Y1，其中G1 为对应更新后的原子序列辅助矩阵G = DtD的子矩阵，并令当前辅助内积α等于初始辅助 内积，减去得到的迭代中间变量β ; 7)根据更新后的原子序列求得第η次迭代后的另一中间变量，^为 更新的原子序列所对应的迭代中间变量β的子向量，更新第η次迭代平方误差εη = 8)判断是否满足终止条件，若是，则输出原始拉曼测量信号y稀疏表达γ，重构的 有效拉曼信号X = D γ，否则η = η+1，返回步骤2)。 优选的，所述的步骤3)具体为，将求得的序号为k原子添加到高斯字典D最后一 列，更新辅助矩阵G 1J: = !>,/?,,,，再根据Liy = Grf求取中间变量ω，根据求取的ω及式，.更新L。 优选的，在步骤8)中，所述的终止条件为第η次迭代的平方误差ε η〈 ε。 由上述对本专利技术的描述可知，与现有技术相比，本专利技术具有如下有益效果： 本专利技术基于信号重构的角度从原始测量信号提取有效拉曼信号，重构的光谱曲线 光滑，能同时实现噪声去除和基线校正，拉曼峰位移定位准确、峰面积衰减小，信噪比改善 明显，且算法处理效率高。【附图说明】 图1为本专利技术的流程图； 图2为理想仿真信号； 图3为对图2的理想仿真信号加噪后的信号； 图4为Batch-OMP处理后的仿真加噪信号与理想仿真信号的对比； 图5为分别采用移动窗口平均、小波变换与Batch-OMP处理原始加噪信号的仿真 的效果对比图； 图6为别采用移动窗口平均、小波变换与Batch-OMP处理原始信号的仿真的效果 对比图。【具体实施方式】 以下通过【具体实施方式】对本专利技术作进一步的描述。 本专利技术的，其采用的原理是通过每步挑选 出内积最大的候选原子，再将所选原子进行施密特正交化处理，将信号投影在这些正交原 子所构成的空间上，得到信号在候选原子上的分量及残余，其后采用相同的方法分解残差 分量，进过有限次分解，得到原信号由有限个原子组成的线性。施密特正交化处理的步骤保 证了单次迭代的最优性，但是最小二乘法的迭代过程较为复杂，因此OMP算法存在求解速 度较慢和存储量较大的问题。本专利技术借助Cholesky分解和计算内积辅助内积D Tr代替直 接计算残差r，提高了算法处理效率，本专利技术Batch-OMP具体实施流程框图如图1所示。稀 疏表达字典可分为基于特定函数的分析字典和基于样本学习的学习字典。由于有用拉曼信 号淹没在强背景噪声下，若采用样本学习字典，噪声会被误作为有用信号进行表达，因此我 们选用分析字典。由于拉曼峰可近似为高斯峰，且高斯函数各参数（峰强度、峰位置和标准 差）物理意义明确，因此我们选用高斯函数作为分析字典对强噪声背景下的拉曼信号进行 稀疏表达。 定义第η次迭代结束时残差信号和信号y的稀疏表达分别为rn，γ n，高斯字典为 D，由于每次迭代残差都与当前所选原子正交，所以对于所有的迭代序号n，满足（rn)TDy n =0 定义平方误差S = IriI，并引入辅助矩阵G = DTD，由上式可得 引入中间变量δη= (γ n)TGyn，因此第η次迭代后的平方误差更新步骤可表 示为ε η= ε η Ιδη+δ11 \定义中间变量β = Gy ;由于信号y的稀疏表达γ = D+y = (DtD) Vy，引入Cholesky分解求解第η次迭代后原始拉曼测量信号稀疏表达γι，进而提高 算法计算效率。 本专利技术的当前第1页1 2 本文档来自技高网...
【技术保护点】
一种强噪声背景下的拉曼信号重构方法，其特征在于，预先定义初始辅助内积为α0＝DTy，y为原始拉曼测量信号，D为高斯字典；辅助矩阵G＝DTD；初始平方误差ε0＝yTy；目标平方误差为ε；残差r＝(y‑Dγ)，γ为y的稀疏表达；包括如下步骤：1)初始化：令迭代过程所选原子序列I:＝(),对辅助矩阵GI,I的Cholesky分解L赋初值L:＝[1],辅助内积α:＝α0,迭代过程中间变量δ0:＝0，初始迭代序号n:＝1；2)计算高斯字典原子和残差信号的辅助内积α，选出与当前残差信号最匹配的原子索引判断当前迭代序号是否大于1，若是，则进入步骤3)；否则，跳转到步骤4)；3)通过引入中间变量更新L；4)将选出的最匹配的原子索引加入集合5)根据式为DITy即为当前更新的原子序列组成的对应的高斯字典D的子矩阵DI与被测原始拉曼信号y的内积，求解c，c为更新后的当前高斯字典对原始拉曼测量信号的稀疏表达，即γI；6)根据更新后的原子序列计算第n次迭代后的迭代中间变量β＝GIγI，其中GI为对应更新后的原子序列辅助矩阵G＝DTD的子矩阵，并令当前辅助内积α等于初始辅助内积α0减去得到的迭代中间变量β；7)根据更新后的原子序列求得第n次迭代后的另一中间变量βI为更新的原子序列所对应的迭代中间变量β的子向量，更新第n次迭代平方误差εn＝εn‑1‑δn+δn‑1；8)判断是否满足终止条件，若是，则输出原始拉曼测量信号y稀疏表达γ，重构的有效拉曼信号x＝Dγ，否则n＝n+1，返回步骤2)。...
【技术特征摘要】
【专利技术属性】
技术研发人员：王昕，王秀芬，范贤光，胡振邦，何浩，王小东，阙靖，汤明，李韦，
申请(专利权)人：厦门大学，
类型：发明
国别省市：福建;35
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。