CN109582917B - 一种基于ssa的信号去噪方法、装置、终端设备及存储介质 - Google Patents
一种基于ssa的信号去噪方法、装置、终端设备及存储介质 Download PDFInfo
- Publication number
- CN109582917B CN109582917B CN201811242280.1A CN201811242280A CN109582917B CN 109582917 B CN109582917 B CN 109582917B CN 201811242280 A CN201811242280 A CN 201811242280A CN 109582917 B CN109582917 B CN 109582917B
- Authority
- CN
- China
- Prior art keywords
- value
- error
- ssa
- vector
- signal
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于SSA的信号去噪方法、装置、终端设备及存储介质,方法包括,采集待去噪信号,并将待去噪信号转化为汉克尔矩阵;将汉克尔矩阵与其转置矩阵相乘后进行奇异值分解,获得若干个二维SSA矩阵;通过L1范数准则建立与每一反对角线对应的误差矢量优化问题;其中,误差矢量优化问题为每一反对角线矢量与每一待求代表值之间的误差的函数,待求代表值为用于代表每一反对角矢量的未知量;计算每一误差矢量优化问题的最优解,根据每一最优解得到每一待求代表值的值;将所有待求代表值的值相加,获得去噪后的信号。通过实施本发明实施例能够降低去噪后的信号的失真度。
Description
技术领域
本发明涉及信号处理领域,具体涉及一种基于SSA的信号去噪方法、装置、终端设备及存储介质。
背景技术
SSA是一种非线性和自适应时频分析工具,用于将信号表示为一维SSA矢量的总和。它广泛应用于去除已知信号的噪声,做信号的基本趋势提取和对信号进行估计。在SSA中,信号首先表示为汉克尔矩阵。然后通过对汉克尔矩阵及其转置矩阵的乘积执行奇异值分解(SVD),获得特征值集合和相应的特征向量集合。然后,导出该组二维SSA矩阵。这里,汉克尔矩阵可以表示为这些二维SSA矩阵的总和。
为了获得一维SSA向量的集合,需要执行去汉克化过程。经典的去汉克尔过程是基于对角平均方法把汉克尔矩阵中的每个二维SSA矩阵变为一维向量,即一个代表值,将每个二维SSA矩阵的一维向量相加后,得到去噪后的信号。具体的,对角平均法指的是计算每个二维SSA矩阵的每个非对角线的平均值。所述平均值为代表值与非对角矢量中的元素之间的均方误差的最小值。式子可以写成
发明内容
本发明提供一种基于SSA的信号去噪方法、装置、终端设备及存储介质,能降低去噪后的信号的失真度。
本发明第一实施例提供了一种SSA的信号去噪方法,包括:
采集待去噪信号,并将所述待去噪信号转化为汉克尔矩阵;
将所述汉克尔矩阵与其转置矩阵相乘后进行奇异值分解,获得若干个二维SSA矩阵;其中,每一所述二维SSA矩阵包括若干反对角线,每一反对角线对应一反对角线矢量;
通过L1范数准则建立与每一所述反对角线对应的误差矢量优化问题;其中,所述误差矢量优化问题为每一所述反对角线矢量与每一待求代表值之间的误差的函数,所述待求代表值为用于代表每一所述反对角矢量的未知量;
计算每一所述误差矢量优化问题的最优解,根据每一所述最优解得到每一所述待求代表值的值;
将所有所述待求代表值的值相加,获得去噪后的信号。
进一步的,通过以下表达式,建立所述通过L1范数准则建立每一所述反对角线对应的误差矢量优化问题:
进一步的,所述计算每一所述误差矢量优化问题的最优解,根据每一所述最优解得到每一所述待求代表值的值,具体为;
Anzn≤x′n;
进一步的,在计算每一所述误差矢量优化问题的最优解之后,还包括:
计算每一所述最优解与对应的每一所述反对角线矢量的误差值;
若所述误差值与0的差值小于预设值,则将所述误差值设为0,得到第二误差值,并根据所述第二误差值求得每一所述待求代表值的值。
相应的在本发明第一实施例的基础上提供了第二实施例
本发明第二实施例提供了一种基于SSA的信号去噪装置,包括:汉克尔矩阵转换模块、奇异值分解模块、误差矢量优化问题生成模块、待求代表值求解模块和信号重构模块;
其中,所述汉克尔矩阵转换模块,用于采集待去噪信号,并将所述待去噪信号转化为汉克尔矩阵;
所述奇异值分解模块,用于将所述汉克尔矩阵与其转置矩阵相乘后进行奇异值分解,获得若干个二维SSA矩阵;其中,每一所述二维SSA矩阵包括若干反对角线,每一反对角线对应一反对角线矢量;
所述误差矢量优化问题生成模块,用于通过L1范数准则建立每一所述反对角线对应的误差矢量优化问题;其中,所述误差矢量优化问题为每一所述反对角线矢量与每一待求代表值之间的误差的函数,所述待求代表值为用于代表每一所述反对角矢量的未知量;
所述待求代表值求解模块,用于计算每一所述误差矢量优化问题的最优解,根据每一所述最优解得到每一所述待求代表值的值;
所述信号重构模块,用于将所有所述待求代表值的值相加,获得去噪后的信号。
进一步的,还包括最优解优化模块,其中所述最优解优化模块,用于计算每一所述最优解与对应的每一所述反对角线矢量的误差值;
若所述误差值与0的差值小于预设值,则将所述误差值设为0,得到第二误差值,并根据所述第二误差值求得每一所述待求代表值的值。
相应的本发明第三实施例提供了一种基于SSA的信号去噪终端设备,包括处理器、存储器以及存储在所述存储器且被配置为由所述处理器执行的计算机程序,所述处理器执行所述计算机程序时实现本发明第一实施例所述的基于SSA的信号去噪方法。
相应的本发明第四实施例提供了种计算机可读存储介质,所述计算机可读存储介质包括存储的计算机程序,其中,在所述计算机程序运行时控制所述计算机可读存储介质所在的设备执行本发明第一实施例中所述的基于SSA的信号去噪方法。
通过实施本发明的实施例,有如下有益效果:
本发明实施例提供了一种基于SSA的信号去噪方法、装置、终端设备及存储介质,在将待去噪信号通过SSA时频分析工具转换为汉克尔矩阵,并在执行奇异值分解后得到,若干个二维的SSA矩阵,紧接着在去汉克尔化过程中,将传统的利用L2范则建立的误差矢量优化问题,改为以L1范则建立的误差矢量优化问题,可以使得求出来的代表值与反对角矢量之间的误差矢量的稀疏度更高,即所获得的代表值更加的准确,从而降低去噪后的信号的失真度。
附图说明
图1是本发明第一实施例提供的基于SSA的信号去噪方法流程示意图;
图2为在L2范则下所求得的误差示意图;
图3为本发明第一实施例提供的基于SSA的信号去噪方法中求得的误差示意图;
图4是本发明第二实施例提供的基于EMD信号重构装置的结构示意图;
附图标记说明:201、汉克尔矩阵转换模块;202、奇异值分解模块;203、误差矢量优化问题生成模块;204、待求代表值求解模块;205、信号重构模块。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1,本发明第一实施例提供的一种基于SSA的信号去噪方法包括:
S101、采集待去噪信号,并将待去噪信号转化为汉克尔矩阵。
S102、将汉克尔矩阵与其转置矩阵相乘后进行奇异值分解,获得若干个二维SSA矩阵;其中,每一二维SSA矩阵包括若干反对角线,每一反对角线对应一反对角线矢量。
S103、通过L1范数准则建立与每一反对角线对应的误差矢量优化问题;其中,误差矢量优化问题为每一反对角线矢量与每一待求代表值之间的误差的函数,待求代表值为用于代表每一反对角矢量的未知量。
S104、计算每一误差矢量优化问题的最优解,根据每一最优解得到每一待求代表值的值。
S105、将所有待求代表值的值相加,获得去噪后的信号。
对于步骤S101、步骤S102将待去噪信号,转化为汉克尔矩阵,并再与其转置矩阵相乘之后进行奇异值分解,得到若干个二维SSA矩阵,为本领域惯用的技术手段,再此不在赘述。
对于步骤S103、在原有的去汉克尔化的过程是通过使用L2范数规则建立优化问题,但这样求出来的代表值与反对角线向量中的许多元素所表示的值不同,因此提出利用L1范数准则来寻找SSA中去汉克化过程中的代表值,建立L1范数准则下的优化问题;具体的,通过以下表达式,建立通过L1范数准则建立每一反对角线对应的误差矢量优化问题:
Anzn≤x′n;
对于步骤S105、在求出代表值后,把所有待求代表值的值相加,获得去噪后的信号。
优选的,在求出最优解后,因为代表误差已经变得稀疏,即有用信号变得相对分散,此时把代表误差的较小值设为0,从而得到一个新的代表误差。再利用这个新的代表误差去求出一个新的代表值,利用这个新的代表值来重构出信号,即在计算每一误差矢量优化问题的最优解之后,还包括:
计算每一最优解与对应的每一反对角线矢量的误差值;
若误差值与0的差值小于预设值,则将误差值设为0,得到第二误差值,并根据第二误差值求得每一待求代表值的值。
把优化问题改写成L1的形式以后,得到的代表误差是稀疏的。利用传统方法去汉克化,小于最大值的1%-5%之间数的数目少于利用L1的方法。这意味着,比传统去汉克化的方法所获得的误差向量的稀疏度,使用L1范数获得的误差向量实现了更高的稀疏度。换句话说,通过我们提出的方法中使用的L1范数标准获得的大多数代表值比通过SSA中的常规去汉克化过程中使用的L2范数标准获得的值更准确。因此,使用改进方法比使用传统方法所获得的代表值更加准确。如图2和图3所示,可以明显看到,在基于L1的情况下,误差是稀疏的。
以下用1%这种情况来做例子:在使用L1的情况下时,有2406个点小于最大误差值的1%。相反地,在使用传统方法的情况下,只有1531个点小于最大误差的1%。所以使用L1方法可以更准确地找到一个匹配汉克矩阵反对角线的值。并且利用该方法可以更好地让信号与噪声分离开,不用混叠在一起。
相应的在本发明第一实施例的基础上提供了第二实施例。
如图4所示,本发明第二实施例提供了一种基于SSA的信号去噪装置,包括:汉克尔矩阵转换模块、奇异值分解模块、误差矢量优化问题生成模块、待求代表值求解模块和信号重构模块;
其中,汉克尔矩阵转换模块,用于采集待去噪信号,并将待去噪信号转化为汉克尔矩阵;
奇异值分解模块,用于将汉克尔矩阵与其转置矩阵相乘后进行奇异值分解,获得若干个二维SSA矩阵;其中,每一二维SSA矩阵包括若干反对角线,每一反对角线对应一反对角线矢量;
误差矢量优化问题生成模块,用于通过L1范数准则建立每一反对角线对应的误差矢量优化问题;其中,误差矢量优化问题为每一反对角线矢量与每一待求代表值之间的误差的函数,待求代表值为用于代表每一反对角矢量的未知量;
待求代表值求解模块,用于计算每一误差矢量优化问题的最优解,根据每一最优解得到每一待求代表值的值;
信号重构模块,用于将所有待求代表值的值相加,获得去噪后的信号。
优选的,还包括最优解优化模块,其中最优解优化模块,用于计算每一最优解与对应的每一反对角线矢量的误差值;
若误差值与0的差值小于预设值,则将误差值设为0,得到第二误差值,并根据第二误差值求得每一待求代表值的值。
相应的在本发明第一实施例的基础上提供了第三实施例。
本发明第三实施例提供了一种SSA的信号去噪终端设备,包括:处理器、存储器以及存储在存储器且被配置为由处理器执行的计算机程序,处理器执行计算机程序时实现本法明第一实施例中的基于SSA的信号去噪方法。
相应的在本发明第一实施例的基础上提供了第四实施例
相应的本发明第四实施例提供了种计算机可读存储介质,计算机可读存储介质包括存储的计算机程序,其中,在计算机程序运行时控制计算机可读存储介质所在的设备执行本发明第一实施例中的SSA的信号去噪方法。
通过实施本发明的实施例,有如下有益效果:
本发明实施例提供了一种基于SSA的信号去噪方法、装置、终端设备及存储介质,在将待去噪信号通过SSA时频分析工具转换为汉克尔矩阵,并在执行奇异值分解后得到,若干个二维的SSA矩阵,紧接着在去汉克尔化过程中,将传统的利用L2范则建立的误差矢量优化问题,改为以L1范则建立的误差矢量优化问题,可以使得求出来的代表值与反对角矢量之间的误差矢量的稀疏度更高,即所获得的代表值更加的准确,从而降低去噪后的信号的失真度。
需说明的是,以上所描述的装置实施例仅仅是示意性的,其中作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。另外,本发明提供的装置实施例附图中,模块之间的连接关系表示它们之间具有通信连接,具体可以实现为一条或多条通信总线或信号线。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。
Claims (6)
1.一种基于SSA的信号去噪方法,其特征在于,包括:
采集待去噪信号,并将所述待去噪信号转化为汉克尔矩阵;
将所述汉克尔矩阵与其转置矩阵相乘后进行奇异值分解,获得若干个二维SSA矩阵;其中,每一所述二维SSA矩阵包括若干反对角线,每一反对角线对应一反对角线矢量;
通过以下表达式,建立所述通过L1范数准则建立每一所述反对角线对应的误差矢量优化问题:
其中,所述误差矢量优化问题为每一所述反对角线矢量与每一待求代表值之间的误差的函数,所述待求代表值为用于代表每一所述反对角矢量的未知量;为矩阵中的第n个反对角线向量,矩阵为SSA算法中的第m个部件,un,m为所述待求代表值,用以表示反对角线向量ιn,m是一条和等长的全1向量;
Anzn≤x′n;
将所有所述待求代表值的值相加,获得去噪后的信号。
2.如权利要求1所述的基于SSA的信号去噪方法,其特征在于,在计算每一所述误差矢量优化问题的最优解之后,还包括:
计算每一所述最优解与对应的每一所述反对角线矢量的误差值;
若所述误差值与0的差值小于预设值,则将所述误差值设为0,得到第二误差值,并根据所述第二误差值求得每一所述待求代表值的值。
3.一种基于SSA的信号去噪装置,其特征在于,包括:汉克尔矩阵转换模块、奇异值分解模块、误差矢量优化问题生成模块、待求代表值求解模块和信号重构模块;
其中,所述汉克尔矩阵转换模块,用于采集待去噪信号,并将所述待去噪信号转化为汉克尔矩阵;
所述奇异值分解模块,用于将所述汉克尔矩阵与其转置矩阵相乘后进行奇异值分解,获得若干个二维SSA矩阵;其中,每一所述二维SSA矩阵包括若干反对角线,每一反对角线对应一反对角线矢量;
所述误差矢量优化问题生成模块,用于
通过以下表达式,建立所述通过L1范数准则建立每一所述反对角线对应的误差矢量优化问题:
其中,所述误差矢量优化问题为每一所述反对角线矢量与每一待求代表值之间的误差的函数,所述待求代表值为用于代表每一所述反对角矢量的未知量;为矩阵中的第n个反对角线向量,矩阵为SSA算法中的第m个部件,un,m为所述待求代表值,用以表示反对角线向量ιn,m是一条和等长的全1向量;所述待求代表值求解模块,用于将每一所述误差矢量优化问题,利用上镜图的形式改写为一个具有两个不等式约束的线性规划问题即:具体公式如下:
Anzn≤x′n;
所述信号重构模块,用于将所有所述待求代表值的值相加,获得去噪后的信号。
4.如权利要求3所述的基于SSA的信号去噪装置,其特征在于,还包括最优解优化模块,其中所述最优解优化模块,用于计算每一所述最优解与对应的每一所述反对角线矢量的误差值;
若所述误差值与0的差值小于预设值,则将所述误差值设为0,得到第二误差值,并根据所述第二误差值求得每一所述待求代表值的值。
5.一种基于SSA的信号去噪终端设备,包括处理器、存储器以及存储在所述存储器且被配置为由所述处理器执行的计算机程序,所述处理器执行所述计算机程序时实现如权利要求1-2中任意一项所述的基于SSA的信号去噪方法。
6.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质包括存储的计算机程序,其中,在所述计算机程序运行时控制所述计算机可读存储介质所在的设备执行如权利要求1-2任意一项所述的基于SSA的信号去噪方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811242280.1A CN109582917B (zh) | 2018-10-24 | 2018-10-24 | 一种基于ssa的信号去噪方法、装置、终端设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811242280.1A CN109582917B (zh) | 2018-10-24 | 2018-10-24 | 一种基于ssa的信号去噪方法、装置、终端设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109582917A CN109582917A (zh) | 2019-04-05 |
CN109582917B true CN109582917B (zh) | 2023-04-14 |
Family
ID=65920483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811242280.1A Active CN109582917B (zh) | 2018-10-24 | 2018-10-24 | 一种基于ssa的信号去噪方法、装置、终端设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109582917B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110292374B (zh) * | 2019-05-31 | 2022-05-17 | 辽宁师范大学 | 基于奇异谱分析和变分模态分解的心电信号去基线漂移方法 |
CN112561306B (zh) * | 2020-12-11 | 2023-12-08 | 领伟创新智能系统(浙江)有限公司 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010088190A2 (en) * | 2009-01-30 | 2010-08-05 | Geco Technology B.V. | Deriving a function that represents data points |
CN105807241A (zh) * | 2016-03-23 | 2016-07-27 | 厦门大学 | 一种利用先验信息的指数信号去噪方法 |
CN106526665A (zh) * | 2016-09-20 | 2017-03-22 | 中国石油天然气股份有限公司 | 一种随机噪声衰减的方法及装置 |
-
2018
- 2018-10-24 CN CN201811242280.1A patent/CN109582917B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010088190A2 (en) * | 2009-01-30 | 2010-08-05 | Geco Technology B.V. | Deriving a function that represents data points |
CN105807241A (zh) * | 2016-03-23 | 2016-07-27 | 厦门大学 | 一种利用先验信息的指数信号去噪方法 |
CN106526665A (zh) * | 2016-09-20 | 2017-03-22 | 中国石油天然气股份有限公司 | 一种随机噪声衰减的方法及装置 |
Non-Patent Citations (3)
Title |
---|
基于奇异值分解的压缩感知观测矩阵优化算法;李周等;《计算机应用》;20180210(第02期);全文 * |
基于改进奇异谱分析的信号去噪方法;戴豪民等;《北京理工大学学报》;20160715(第07期);全文 * |
音乐信号稀疏分析方法研究;丰上等;《复旦学报(自然科学版)》;20180615(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109582917A (zh) | 2019-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bioucas-Dias | Bayesian wavelet-based image deconvolution: A GEM algorithm exploiting a class of heavy-tailed priors | |
CN109582917B (zh) | 一种基于ssa的信号去噪方法、装置、终端设备及存储介质 | |
JP2009037597A (ja) | 入力画像をフィルタリングして出力画像を作成する方法 | |
CN113945967A (zh) | 绕射波的分离方法及装置 | |
CN112906335B (zh) | 一种用于集成电路系统的无源性校正方法及装置 | |
Gupta et al. | On the approximate discrete KLT of fractional Brownian motion and applications | |
CN109756740B (zh) | 基于最优测量矩阵的半张量图像压缩方法和图像恢复方法 | |
CN110929862A (zh) | 定点化的神经网络模型量化装置和方法 | |
CN116012263A (zh) | 一种图像噪声去除方法、装置、存储介质及电子设备 | |
CN109239006B (zh) | 一种基于湿度补偿模型的物质识别方法、装置及存储介质 | |
CN114047499B (zh) | 传感器阵列信号幅度不一致时的信号归一化方法及设备 | |
Mamouni et al. | A new convolutive source separation approach for independent/dependent source components | |
AU761504B2 (en) | Signal processing | |
CN112241740B (zh) | 一种特征提取方法及装置 | |
CN116805370A (zh) | 图像识别方法、装置及相关设备 | |
Novosadová et al. | Piecewise-polynomial signal segmentation using reweighted convex optimization | |
Fernández et al. | Auto-adaptive Laplacian Pyramids. | |
Zhao et al. | Blind source separation based on dictionary learning: a singularity-aware approach | |
CN108597534B (zh) | 基于卷积框架的语音信号稀疏表示方法 | |
CN113037523B (zh) | 网络流量流向预测方法、装置以及存储介质 | |
US10748551B2 (en) | Noise suppression system, noise suppression method, and recording medium storing program | |
CN114386463A (zh) | 基于字典学习的核磁共振回波数据处理方法及装置 | |
Im et al. | A least squares-type density estimator using a polynomial function | |
CN116644265B (zh) | 一种非线性信号的补偿方法、补偿装置及终端设备 | |
Lilienthal | On the method of probability weighted moments in regional frequency analysis |
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 |