CN110619296A - 一种基于奇异分解的信号降噪方法 - Google Patents
一种基于奇异分解的信号降噪方法 Download PDFInfo
- Publication number
- CN110619296A CN110619296A CN201910851020.2A CN201910851020A CN110619296A CN 110619296 A CN110619296 A CN 110619296A CN 201910851020 A CN201910851020 A CN 201910851020A CN 110619296 A CN110619296 A CN 110619296A
- Authority
- CN
- China
- Prior art keywords
- matrix
- signal
- singular value
- noise reduction
- singular
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明旨在提供一种基于奇异值分解的信号降噪方法,其特征在于,包括以下步骤:将Hankel矩阵作为采样信号的吸引子轨迹矩阵,根据奇异值能量最大准则获取矩阵的最优列,确定最佳矩阵结构;奇异值分解最佳构造矩阵得到非零奇异值,通过改进后的赤池信息准则获取奇异值的有效阶次;利用有效分量开展奇异值分解的逆运算,得到近似信号的构造矩阵;采用平均法对近似信号的构造矩阵进行还原,得到最终降噪信号。本发明采用改进后的赤池信息准则和奇异值分解对信号进行降噪处理,能避免奇异值分解过程因阶次选择不当而导致的过降噪或欠降噪发生,准确完成信号中的噪声分离。方法具有较强的自适应性,对于非线性和非平稳信号具有很好的降噪效果。
Description
技术领域
本发明涉及一种基于奇异分解的信号降噪方法,属于信号处理技术领域。
背景技术
受环境与结构因素影响,被测对象在信号采集时通常会混入一定的干扰噪声。为了更好地提取***特征参数,在分析前需要对信号进行降噪预处理。由于实际采集过程中的信号大多为非线性或非平稳信号,采用传统处理方法存在着明显缺陷。由时频分析发展出的各类降噪技术(如小波阈值去噪WTD、EMD-SG等)虽然得到广泛应用,但仍存在着诸多局限,如阈值和基函数的选择、处理后的相位移动与波形畸变等。基于奇异值分解(SVD)的信号处理方法是一种分析非线性、非平稳信号的有效工具,近年来被用在振动信号和图像的降噪中。由于矩阵奇异值的稳定特性,信号在奇异值分解时能确保奇异值所代表的信号特征具有很强的鲁棒性。
为了实现信号的有效去噪,采用一种基于奇异值分解的信号降噪方法。该方法根据奇异值能量最大准则选取信号在Hankel矩阵下的最优行列数,对最优构造矩阵进行奇异值分解;然后采用改进后的赤池信息准则进行精准定阶,并利用有效奇异分量求得近似信号的构造矩阵。最后通过平均法将近似信号的构造矩阵进行时序还原,从而完成信号的稳定降噪。
发明内容
本发明要解决的技术问题是:信号在预处理过程中存在降噪效果不稳定、自适应差。
为了解决上述技术问题,本发明的技术方案是提供了一种基于奇异分解的信号降噪方法,其特征在于,包括以下步骤:
步骤1)以m×n维Hankel矩阵为采样信号s=[s(1),s(2),…,s(N)]建立重构吸引子轨迹矩阵,对重构吸引子轨迹矩阵进行奇异值分解;根据奇异值能量最大准则获取矩阵最优列确定最佳矩阵结构;
步骤2)对列Hankel矩阵进行奇异值分解,得到最佳构造矩阵分解后的非零奇异值序列σ=[σ1,σ2,…,σq];采用改进的赤池信息准则计算σi相对应的AIC函数值,i=1,2,...,q,确定最小索引k为奇异值的有效阶次;
步骤3)利用前k阶有效分量进行奇异值分解的逆运算,得到近似信号的构造矩阵式中,ui表示Hankel矩阵的左奇异向量,vi表示Hankel矩阵的右奇异向量;
步骤4)采用平均法对近似信号的构造矩阵进行还原,得到最终的时间序列即为降噪信号
优选地,所述步骤1)中,确定最佳矩阵结构包括以下步骤:
步骤101、将所采集的信号序列s=[s(1),s(2),…,s(N)]构建为m×n维Hankel矩阵:
其中,N为信号的采样个数,信号重构吸引子轨迹矩阵的行m=N-n+1、列n=1,2,...,N;
步骤102、对不同n值的Hankel矩阵进行奇异值分解:A=UΣVT,得到构造矩阵奇异值的递减序列σ=[σ1,σ2,…,σq],式中,U表示Hankel矩阵的左奇异矩阵,V表示Hankel矩阵的右奇异矩阵,Σ表示Hankel矩阵的对角矩阵;
步骤103、根据以下公式求解不同列Hankel矩阵的奇异值能量,按奇异值能量最大准则选取构造矩阵的最优列从而确定信号最佳的矩阵结构:
式中,E(n)表示Hankel矩阵的奇异值能量。
优选地,所述步骤2)中确定奇异值的有效阶次包括以下步骤:
对列Hankel矩阵进行奇异值分解,得到最佳构造矩阵分解后的非零奇异值序列σ=[σ1,σ2,…,σq],为适用于有色噪声信号的去噪,采用改进后的赤池信息准则对奇异值进行定阶:首先判断信号中是否存在有色噪声,若不存在,则特征值μi=σi 2;若存在,则奇异值校正后的特征值然后根据改进后的赤池信息准则计算相应的AIC值AIC(d),
选取AIC中最小值所对应的索引k,即为信号奇异值分解的有效阶次。
优选地,所述步骤4)中对近似信号的构造矩阵进行信号还原包括以下步骤:
由前k阶有效分量得到近似信号的构造矩阵为式中,xij表示近似信号的构造矩阵元素,表示近似信号的构造矩阵行数,表示近似信号的构造矩阵列数,采用平均法对近似构造矩阵进行时序还原,计算公式如下:
式中,表示降噪后的时序信号元素,i=1,2,...,N,
通过上式求得的即为初始信号s的降噪信号。
本发明提供的方法吸收了Hankel矩阵在奇异值分解中的优点,通过奇异值能量最大准则选取Hankel矩阵的最优行列数,使奇异值之间具有良好的区分度,避免信号分量间的特征耦合。然后采用改进后的赤池信息准则确定矩阵在最优结构下的有效奇异值,从而完成对信号有效分量和噪声分量的识别与分离,达到稳定、有力的降噪目的。
本发明采用奇异值能量最大准则和改进后的赤池信息准则依次对信号时间序列的构造矩阵结构及奇异值的有效阶次进行确定,避免信号在处理过程中的不定参数选择,能实现精确的信源数识别与噪声分离,适用于复杂的非稳定信号,具有较强的自适应性,降噪效果相对于其他方法更加稳定有效。
附图说明
图1为本发明提供的一种基于奇异值分解的信号降噪方法的流程图;
图2为本发明提供的信号在不同降噪方法下的效果对比图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
本发明提供了一种基于奇异值分解的信号降噪方法,如图1所示,包括以下步骤:
步骤1)以m×n维Hankel矩阵为采样信号s=[s(1),s(2),…,s(N)]建立重构吸引子轨迹矩阵,重构吸引子轨迹矩阵的结构样式为:
其中,N为信号的采样个数,信号重构吸引子轨迹矩阵的行m=N-n+1、列n=1,2,...,N;
对不同n值的Hankel矩阵进行奇异值分解:A=UΣVT,计算奇异值能量大小由奇异值能量最大准则选取最优列从而确定信号重构吸引子轨迹矩阵的最佳结构。
步骤2)对列Hankel矩阵奇异值分解,得到非零奇异值σ=[σ1,σ2,…,σq]。判断信号中是否存在有色噪声,若无有色噪声,则特征值μi=σi 2;若信号中存在有色噪声,则奇异值校正后的特征值为根据改进后的赤池信息准则计算相应的AIC值AIC(d),
选取AIC中最小值所对应的索引k,即为信号奇异值分解的有效阶次。
步骤3)利用前k阶有效分量进行奇异值分解的逆运算,得到近似信号的构造矩阵式中,ui表示Hankel矩阵的左奇异向量,vi表示Hankel矩阵的右奇异向量。
步骤4)采用平均法对近似信号的构造矩阵进行时序还原,计算公式如下:
式中,表示降噪后的时序信号元素,i=1,2,...,N,
通过上式求得的即为初始信号s的降噪信号。
仿真实验验证
1、确定信号构造矩阵
以信号x1为例进行基于奇异值分解的信号降噪仿真实验,x1的表达式为:
x1=e-2t[sin(40πt)+0.5sin(200πt)]
取采样个数N=1000,向纯净信号x1中添加SNR=5dB高斯白噪声得到待处理信号s1,以m×n维Hankel矩阵为s1的吸引子轨迹矩阵,对其进行奇异值分解,通过奇异值总能量最大准则确定构造矩阵的最优行列
2、确定有效奇异值阶次
对信号s1的列Hankel矩阵进行奇异值分解,得到非零奇异值序列σ;根据公式求解相应的AIC函数值,得到最小AIC值索引k=4;
3、计算近似信号构造矩阵
利用前4阶有效分量进行奇异值分解的逆运算,得到近似信号的构造矩阵
4、降噪信号的时序还原
采用平均法公式对近似信号的构造矩阵进行时序还原,得到最终降噪信号图2显示了信号s1在基于奇异值分解的信号降噪方法(AIC-SVD)、小波阈值去噪(WTD)及(EMD-SG)降噪处理下的最终信号,结合定量指标计算得到AIC-SVD、WTD和EMD-SG降噪后的信噪比(SNR)分别为38.3dB、7.2dB、9.9dB,由此看出信号采用AIC-SVD的降噪效果要明显优于后两者。
Claims (4)
1.一种基于奇异分解的信号降噪方法,其特征在于,包括以下步骤:
步骤1)以m×n维Hankel矩阵为采样信号s=[s(1),s(2),…,s(N)]建立重构吸引子轨迹矩阵,对重构吸引子轨迹矩阵进行奇异值分解;根据奇异值能量最大准则获取矩阵最优列确定最佳矩阵结构;
步骤2)对列Hankel矩阵进行奇异值分解,得到最佳构造矩阵分解后的非零奇异值序列σ=[σ1,σ2,…,σq];采用改进的赤池信息准则计算σi相对应的AIC函数值,i=1,2,...,q,确定最小索引k为奇异值的有效阶次;
步骤3)利用前k阶有效分量进行奇异值分解的逆运算,得到近似信号的构造矩阵式中,ui表示Hankel矩阵的左奇异向量,vi表示Hankel矩阵的右奇异向量;
步骤4)采用平均法对近似信号的构造矩阵进行还原,得到最终的时间序列即为降噪信号
2.如权利要求1所述的一种基于奇异值分解的信号降噪方法,其特征在于,所述步骤1)中,确定最佳矩阵结构包括以下步骤:
步骤101、将所采集的信号序列s=[s(1),s(2),…,s(N)]构建为m×n维Hankel矩阵:
其中,N为信号的采样个数,信号重构吸引子轨迹矩阵的行m=N-n+1、列n=1,2,...,N;
步骤102、对不同n值的Hankel矩阵进行奇异值分解:A=UΣVT,得到构造矩阵奇异值的递减序列σ=[σ1,σ2,…,σq],式中,U表示Hankel矩阵的左奇异矩阵,V表示Hankel矩阵的右奇异矩阵,Σ表示m×n维Hankel矩阵的对角矩阵;
步骤103、根据以下公式求解不同列Hankel矩阵的奇异值能量,按奇异值能量最大准则选取构造矩阵的最优列从而确定信号最佳的矩阵结构:
式中,E(n)表示Hankel矩阵的奇异值能量。
3.如权利要求1所述的一种基于奇异值分解的信号降噪方法,其特征在于,所述步骤2)中确定奇异值的有效阶次包括以下步骤:
对列Hankel矩阵进行奇异值分解,得到最佳构造矩阵分解后的非零奇异值序列σ=[σ1,σ2,…,σq],为适用于有色噪声信号的去噪,采用改进后的赤池信息准则对奇异值进行定阶:首先判断信号中是否存在有色噪声,若不存在,则特征值μi=σi 2;若存在,则奇异值校正后的特征值然后根据改进后的赤池信息准则计算相应的AIC值AIC(d),
选取AIC中最小值所对应的索引k,即为信号奇异值分解的有效阶次。
4.如权利要求1所述的一种基于奇异值分解的信号降噪方法,其特征在于,所述步骤4)中对近似信号的构造矩阵进行信号还原包括以下步骤:
由前k阶有效分量得到近似信号的构造矩阵为式中,xij表示近似信号的构造矩阵元素,表示近似信号的构造矩阵行数,表示近似信号的构造矩阵列数,采用平均法对近似构造矩阵进行时序还原,计算公式如下:
式中,表示降噪后的时序信号元素,i=1,2,...,N,
通过上式求得的即为初始信号s的降噪信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910851020.2A CN110619296B (zh) | 2019-09-10 | 2019-09-10 | 一种基于奇异分解的信号降噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910851020.2A CN110619296B (zh) | 2019-09-10 | 2019-09-10 | 一种基于奇异分解的信号降噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110619296A true CN110619296A (zh) | 2019-12-27 |
CN110619296B CN110619296B (zh) | 2023-09-26 |
Family
ID=68923076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910851020.2A Active CN110619296B (zh) | 2019-09-10 | 2019-09-10 | 一种基于奇异分解的信号降噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110619296B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111582205A (zh) * | 2020-05-13 | 2020-08-25 | 陕西燃气集团新能源发展股份有限公司 | 一种基于多分辨率奇异值分解模型的降噪方法 |
CN112305380A (zh) * | 2020-09-01 | 2021-02-02 | 华南理工大学 | 一种基于s变换与谱峭度的局部放电白噪声抑制方法 |
CN112924175A (zh) * | 2021-01-22 | 2021-06-08 | 南京信息工程大学 | 一种基于奇异值分解的轴承故障信号降噪方法 |
CN112989966A (zh) * | 2021-02-25 | 2021-06-18 | 中国人民解放军海军航空大学 | 一种改进的模拟电路信号降噪方法 |
CN113138377A (zh) * | 2020-01-17 | 2021-07-20 | 中国科学院声学研究所 | 一种基于多分辨二分奇异值分解的自适应底混响抑制方法 |
CN113253047A (zh) * | 2021-04-19 | 2021-08-13 | 南京南瑞继保电气有限公司 | 单相接地选线方法、装置、电子设备及存储介质 |
CN113642185A (zh) * | 2021-08-19 | 2021-11-12 | 东华大学 | 一种卫星高频颤振传递路径辨识方法 |
CN113742985A (zh) * | 2021-08-26 | 2021-12-03 | 湖州师范学院 | 一种生态***碳通量测量信号的处理方法、装置及电子设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106845010A (zh) * | 2017-02-16 | 2017-06-13 | 西南交通大学 | 基于改进SVD降噪和Prony的低频振荡主导模式辨识方法 |
CN107045149A (zh) * | 2017-04-17 | 2017-08-15 | 吉林大学 | 一种基于双奇异值分解的全波核磁共振信号噪声滤除方法 |
-
2019
- 2019-09-10 CN CN201910851020.2A patent/CN110619296B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106845010A (zh) * | 2017-02-16 | 2017-06-13 | 西南交通大学 | 基于改进SVD降噪和Prony的低频振荡主导模式辨识方法 |
CN107045149A (zh) * | 2017-04-17 | 2017-08-15 | 吉林大学 | 一种基于双奇异值分解的全波核磁共振信号噪声滤除方法 |
Non-Patent Citations (2)
Title |
---|
刘承兰 等: "基于数据相关的多雷达融合成像相干配准研究", 《***工程与电子技术》 * |
叶锡均: "基于环境激励的大型土木工程结构模态参数识别研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113138377A (zh) * | 2020-01-17 | 2021-07-20 | 中国科学院声学研究所 | 一种基于多分辨二分奇异值分解的自适应底混响抑制方法 |
CN113138377B (zh) * | 2020-01-17 | 2023-05-16 | 中国科学院声学研究所 | 一种基于多分辨二分奇异值分解的自适应底混响抑制方法 |
CN111582205A (zh) * | 2020-05-13 | 2020-08-25 | 陕西燃气集团新能源发展股份有限公司 | 一种基于多分辨率奇异值分解模型的降噪方法 |
CN111582205B (zh) * | 2020-05-13 | 2023-05-12 | 陕西燃气集团新能源发展股份有限公司 | 一种基于多分辨率奇异值分解模型的降噪方法 |
CN112305380A (zh) * | 2020-09-01 | 2021-02-02 | 华南理工大学 | 一种基于s变换与谱峭度的局部放电白噪声抑制方法 |
CN112924175A (zh) * | 2021-01-22 | 2021-06-08 | 南京信息工程大学 | 一种基于奇异值分解的轴承故障信号降噪方法 |
CN112989966A (zh) * | 2021-02-25 | 2021-06-18 | 中国人民解放军海军航空大学 | 一种改进的模拟电路信号降噪方法 |
CN112989966B (zh) * | 2021-02-25 | 2022-06-24 | 中国人民解放军海军航空大学 | 一种改进的模拟电路信号降噪方法 |
CN113253047A (zh) * | 2021-04-19 | 2021-08-13 | 南京南瑞继保电气有限公司 | 单相接地选线方法、装置、电子设备及存储介质 |
CN113642185A (zh) * | 2021-08-19 | 2021-11-12 | 东华大学 | 一种卫星高频颤振传递路径辨识方法 |
CN113742985A (zh) * | 2021-08-26 | 2021-12-03 | 湖州师范学院 | 一种生态***碳通量测量信号的处理方法、装置及电子设备 |
CN113742985B (zh) * | 2021-08-26 | 2024-06-28 | 湖州师范学院 | 一种生态***碳通量测量信号的处理方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN110619296B (zh) | 2023-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110619296A (zh) | 一种基于奇异分解的信号降噪方法 | |
CN107220606B (zh) | 基于一维卷积神经网络的雷达辐射源信号的识别方法 | |
CN110865357B (zh) | 一种基于参数优化vmd的激光雷达回波信号降噪方法 | |
CN109890043B (zh) | 一种基于生成式对抗网络的无线信号降噪方法 | |
CN109872305B (zh) | 一种基于质量图生成网络的无参考立体图像质量评价方法 | |
CN110751612A (zh) | 多通道多尺度卷积神经网络的单幅图像去雨方法 | |
CN107392123B (zh) | 一种基于相参积累消噪的射频指纹特征提取和识别方法 | |
CN111624570B (zh) | 基于二维卷积神经网络的雷达目标识别方法 | |
CN105429719B (zh) | 基于功率谱和多尺度小波变换分析强干扰信号检测方法 | |
CN112426160A (zh) | 一种心电信号类型识别方法及装置 | |
CN105572473B (zh) | 高分辨率线性时频分析方法 | |
CN103093432B (zh) | 基于极化分解和图像块相似性的极化sar图像降斑方法 | |
CN102184451B (zh) | 一种基于多小波融合特征的神经元动作特征提取方法 | |
CN104392086A (zh) | 一种基于皮尔逊秩次变量相关系数的信号检测电路及方法 | |
CN114264953A (zh) | 一种永磁同步电机退磁故障诊断方法和***以及诊断装置 | |
CN103871422A (zh) | 基于三因子迭代联合块对角化的时域混叠盲信号分离方法 | |
CN113723244A (zh) | 一种基于改进变分模态分解的雷达辐射源信号分离方法 | |
CN106448694A (zh) | 一种基于复角检测的欠定盲源分离中的时频单源点提取方法 | |
CN101833752B (zh) | 基于奇异值分解重构的红外小目标预处理方法 | |
CN112689288A (zh) | 一种基于wann的射频指纹提取和识别方法 | |
CN110083865B (zh) | 一种抗脉冲噪声的已知信号检测方法及检测电路 | |
CN108122233B (zh) | 基于局部像素综合特征的彩色图像分割方法 | |
CN107341519B (zh) | 一种基于多分辨率分析的支持向量机识别优化方法 | |
CN110133628A (zh) | 一种短时强干扰下lfm信号的恢复方法 | |
CN115809426A (zh) | 一种辐射源个体识别方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |