CN105027519B - A kind of signal processing method and device - Google Patents
A kind of signal processing method and device Download PDFInfo
- Publication number
- CN105027519B CN105027519B CN201380033634.0A CN201380033634A CN105027519B CN 105027519 B CN105027519 B CN 105027519B CN 201380033634 A CN201380033634 A CN 201380033634A CN 105027519 B CN105027519 B CN 105027519B
- Authority
- CN
- China
- Prior art keywords
- signal
- iteration
- sparsity
- value
- iterative
- 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
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/0222—Estimation of channel variability, e.g. coherence bandwidth, coherence time, fading frequency
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Complex Calculations (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明的实施例提供一种信号处理方法和装置,涉及通信技术领域,能够解决稀疏度先验信息受限导致信号重建性能降低的问题,并降低或基本消除重建信号产生的偏差,提高单比特压缩感知技术的实用性与准确度。其方法为:通过获取输入参数与编码端数据,在对过程参数与相关参数值初始化后,进行自适应稀疏度估计操作,并在获取了目标稀疏度后,根据目标稀疏度进行信号重建得到原始稀疏信号的重建信号。本发明的实施例用于单比特压缩感知技术场景下的信号处理。
Embodiments of the present invention provide a signal processing method and device, which relate to the field of communication technology, and can solve the problem of reduced signal reconstruction performance due to limited sparsity prior information, reduce or basically eliminate the deviation generated by the reconstructed signal, and improve the single bit rate. Practicality and accuracy of compressive sensing techniques. The method is as follows: by obtaining input parameters and encoding end data, after initializing the process parameters and related parameter values, performing an adaptive sparsity estimation operation, and after obtaining the target sparsity, reconstructing the signal according to the target sparsity to obtain the original Reconstruction of sparse signals. Embodiments of the present invention are used for signal processing in the scenario of single-bit compressed sensing technology.
Description
技术领域technical field
本发明涉及通信技术领域,尤其涉及一种信号处理方法和装置。The present invention relates to the field of communication technology, in particular to a signal processing method and device.
背景技术Background technique
在传统信号处理理论中,依据香农采样定理,采集速率应至少等于两倍信号带宽才可以无失真地恢复原信号,但现有硬件设备的模数转换和信号处理能力已无法满足不断增长的信号高速采样需求,且高能耗的海量数据采集也并不是必不可少的,故最新的信号处理方法采用压缩感知(Compressive Sensing,CS)技术,以较低的采样速率来采集数据仍可以准确地恢复原信号。In traditional signal processing theory, according to Shannon sampling theorem, the acquisition rate should be at least twice the signal bandwidth to restore the original signal without distortion, but the analog-to-digital conversion and signal processing capabilities of existing hardware equipment can no longer meet the growing signal High-speed sampling is required, and massive data acquisition with high energy consumption is not essential. Therefore, the latest signal processing method uses Compressive Sensing (CS) technology, and the data collected at a lower sampling rate can still be accurately restored. original signal.
CS技术基于信号具有稀疏性这一前提,对输入信号进行低速采样后获得降维的采样信号,译码过程根据降维的采样信号重建原输入信号,实现了通过较少的低维采样数据来重建原始的高维信号,显著降低信号获取开销。应用中可扩展采用单比特CS技术,在对输入信号进行低速采样后,对采样信号进行符号量化,以便模拟域的输入信号转换到数字域进行后续处理、传输、存储等操作,可降低硬件设备模数转换实现的复杂度、信息接收和获取的复杂度、系统传输和存储的数据量,且数据鲁棒性好,使CS技术更具实用性。CS technology is based on the premise that the signal is sparse. After low-speed sampling of the input signal, the dimension-reduced sampling signal is obtained. The decoding process reconstructs the original input signal according to the dimension-reducing sampling signal. Reconstruct the original high-dimensional signal, significantly reducing signal acquisition overhead. In the application, the single-bit CS technology can be extended. After the input signal is sampled at a low speed, the sampling signal is quantized so that the input signal in the analog domain can be converted to the digital domain for subsequent processing, transmission, storage, etc., which can reduce hardware equipment. The complexity of analog-to-digital conversion, the complexity of information reception and acquisition, the amount of data transmitted and stored by the system, and the data robustness make CS technology more practical.
但是,现有技术在上述重建信号的过程中,需要预知原始信号的稀疏度,并以此作为输入参数,然而在实际应用中,由于原始输入稀疏信号的动态变化等因素,要确定实时且准确的稀疏度这一先验信息通常是受限的,即使通过对信号的长时间统计观测设置最大稀疏度,也会与实际稀疏度存在偏差,导致信号的重建性能受到严重的影响,获取的重建信号产生较大偏差。However, in the prior art, in the above-mentioned process of reconstructing the signal, it is necessary to predict the sparsity of the original signal and use it as an input parameter. The prior information of sparsity is usually limited. Even if the maximum sparsity is set through long-term statistical observation of the signal, there will be deviations from the actual sparsity, which will seriously affect the reconstruction performance of the signal. The obtained reconstruction The signal has a large deviation.
发明内容Contents of the invention
本发明的实施例提供一种信号处理方法和装置,用于单比特压缩感知技术的场景下,解决了稀疏度先验信息受限导致信号重建性能降低的问题,降低或基本消除了重建信号产生的偏差,提高了单比特压缩感知技术的实用性与准确度。Embodiments of the present invention provide a signal processing method and device, which are used in the scenario of single-bit compressed sensing technology, which solves the problem of reduced signal reconstruction performance due to limited sparsity prior information, and reduces or basically eliminates the generation of reconstructed signals. The deviation improves the practicability and accuracy of single-bit compressed sensing technology.
为达到上述目的,本发明的实施例采用如下技术方案:In order to achieve the above object, embodiments of the present invention adopt the following technical solutions:
第一方面,提供一种信号处理方法,所述方法包括:In a first aspect, a signal processing method is provided, the method comprising:
获取输入参数与编码端数据,并对过程参数与第一迭代信号进行初始化,所述编码端数据包括编码端对原始稀疏信号进行低速采样处理及符号量化处理获取的测量信号;Acquiring input parameters and encoding end data, and initializing the process parameters and the first iteration signal, the encoding end data including the measurement signal obtained by the encoding end performing low-speed sampling processing and symbol quantization processing on the original sparse signal;
根据所述输入参数、所述编码端数据、初始化后的所述过程参数和所述第一迭代信号进行自适应稀疏度估计操作,获得目标稀疏度;performing an adaptive sparsity estimation operation according to the input parameters, the encoding end data, the initialized process parameters, and the first iteration signal to obtain a target sparsity;
根据所述目标稀疏度,对所述测量信号进行信号重建得到所述原始稀疏信号的重建信号。Perform signal reconstruction on the measurement signal according to the target sparsity to obtain a reconstructed signal of the original sparse signal.
结合第一方面,在第一种可能的实现方式中,所述获取输入参数与编码端数据包括:With reference to the first aspect, in a first possible implementation manner, the acquiring input parameters and encoding end data includes:
获取所述输入参数,所述输入参数包括最大稀疏度、判决门限;Acquiring the input parameters, the input parameters include a maximum sparsity and a decision threshold;
接收所述编码端传递的所述编码端数据,所述编码端数据包括所述测量信号、采样矩阵;所述测量信号为所述编码端结合所述采样矩阵对所述原始稀疏信号进行低速采样处理获取采样信号后,再对所述采样信号进行符号量化处理获得的。Receive the encoding end data transmitted by the encoding end, the encoding end data includes the measurement signal and the sampling matrix; the measurement signal is the low speed sampling of the original sparse signal by the encoding end in combination with the sampling matrix After the sampling signal is acquired, the sampling signal is obtained by performing sign quantization processing on the sampling signal.
结合第一方面或第一方面的第一种可能的实现方式,在第二种可能的实现方式中,所述过程参数包括所述目标稀疏度的第一估计值、第一使用值,所述对过程参数与第一迭代信号进行初始化包括:With reference to the first aspect or the first possible implementation manner of the first aspect, in a second possible implementation manner, the process parameters include a first estimated value and a first used value of the target sparsity, and the Initializing process parameters and first iteration signals involves:
将所述第一估计值赋值为所述最大稀疏度,将所述第一使用值赋值为所述最大稀疏度;assigning the first estimated value as the maximum sparsity, and assigning the first use value as the maximum sparsity;
将所述第一迭代信号的初始值赋值为零向量。Assigning the initial value of the first iteration signal as a zero vector.
结合第一方面至第一方面的第二种可能的实现方式中的任一种,在第三种可能的实现方式中,所述根据所述输入参数、所述编码端数据、初始化后的所述过程参数和所述第一迭代信号进行自适应稀疏度估计操作,获得目标稀疏度包括:With reference to any one of the first aspect to the second possible implementation manner of the first aspect, in a third possible implementation manner, the Performing an adaptive sparsity estimation operation on the process parameters and the first iteration signal, and obtaining the target sparsity includes:
根据所述测量信号与所述采样矩阵,获取所述第一迭代信号经低速采样及符号量化后与所述测量信号的差值,并对所述差值左乘所述采样矩阵的转置矩阵,得到第一迭代梯度;According to the measurement signal and the sampling matrix, obtain the difference between the first iteration signal and the measurement signal after low-speed sampling and symbol quantization, and multiply the difference by the transpose matrix of the sampling matrix on the left , get the gradient of the first iteration;
通过梯度下降法获取过程信号,所述过程信号为所述第一迭代信号与所述第一迭代梯度的差值;Obtaining a process signal by a gradient descent method, where the process signal is a difference between the first iteration signal and the first iteration gradient;
对所述第一迭代梯度是否趋近于零向量进行判断,若所述第一迭代梯度趋近于零向量,则将所述目标稀疏度的第二使用值赋值为所述第一估计值;或,若所述第一迭代梯度不趋近于零向量,则将所述第二使用值赋值为所述第一使用值;Judging whether the first iterative gradient approaches a zero vector, and if the first iterative gradient approaches a zero vector, assigning the second used value of the target sparsity as the first estimated value; Or, if the first iterative gradient does not approach the zero vector, assigning the second use value as the first use value;
根据所述过程信号与所述第二使用值,通过门限函数获取第二迭代信号,所述门限函数用于将所述过程信号中元素幅值最大的所述第二使用值个元素值保留,同时将所述过程信号中除元素幅值最大的所述第二使用值个元素之外的其他所有元素值置为零,并将所述门限函数的处理结果赋值给所述第二迭代信号;According to the process signal and the second use value, a second iterative signal is obtained through a threshold function, and the threshold function is used to reserve the second use value and element values with the largest element amplitude in the process signal, At the same time, all other element values in the process signal except the second use value elements with the largest element amplitude are set to zero, and the processing result of the threshold function is assigned to the second iteration signal;
对所述第二迭代信号进行单位归一化,得到归一化迭代信号,统计所述归一化迭代信号中元素幅值绝对值超过所述判决门限的非零元素个数,并将所述元素个数赋值给所述目标稀疏度的第二估计值;Carrying out unit normalization on the second iterative signal to obtain a normalized iterative signal, counting the number of non-zero elements whose absolute value of the element amplitude in the normalized iterative signal exceeds the decision threshold, and calculating the Assigning the number of elements to the second estimated value of the target sparsity;
若所述第一迭代梯度不趋近于零向量或所述第二估计值与所述第一使用值不相等,则将所述第一估计值赋值为所述第二估计值、所述第一使用值赋值为所述第二使用值、所述第一迭代信号赋值为所述第二迭代信号,并重新执行所述自适应稀疏度估计操作;或,若所述第一迭代梯度趋近于零向量且所述第二估计值与所述第一使用值相等,则将所述目标稀疏度赋值为所述第二估计值,得到所述目标稀疏度。If the first iterative gradient does not approach the zero vector or the second estimated value is not equal to the first used value, assign the first estimated value to the second estimated value, the first estimated value A use value is assigned as the second use value, the first iteration signal is assigned as the second iteration signal, and the adaptive sparsity estimation operation is re-executed; or, if the gradient of the first iteration approaches If it is a zero vector and the second estimated value is equal to the first used value, assign the target sparsity as the second estimated value to obtain the target sparsity.
结合第一方面至第一方面的第三种可能的实现方式中的任一种,在第四种可能的实现方式中,所述根据所述目标稀疏度,对所述测量信号进行信号重建得到所述原始稀疏信号的重建信号包括:With reference to any one of the first aspect to the third possible implementation manner of the first aspect, in a fourth possible implementation manner, performing signal reconstruction on the measurement signal according to the target sparsity to obtain The reconstructed signal of the original sparse signal includes:
根据所述编码端数据、所述第一迭代信号与所述目标稀疏度,执行信号重建操作,所述信号重建操作具体包括:Perform a signal reconstruction operation according to the encoding end data, the first iteration signal, and the target sparsity, and the signal reconstruction operation specifically includes:
根据所述测量信号与所述采样矩阵,获取所述第一迭代信号经低速采样及符号量化后与所述测量信号的差值,并对所述差值左乘所述采样矩阵的转置矩阵,得到第一迭代梯度;According to the measurement signal and the sampling matrix, obtain the difference between the first iteration signal and the measurement signal after low-speed sampling and symbol quantization, and multiply the difference by the transpose matrix of the sampling matrix on the left , get the gradient of the first iteration;
通过梯度下降法获取过程信号,所述过程信号为所述第一迭代信号与所述第一迭代梯度的差值;Obtaining a process signal by a gradient descent method, where the process signal is a difference between the first iteration signal and the first iteration gradient;
结合所述过程信号与所述目标稀疏度,通过门限函数获取第二迭代信号,所述门限函数用于将所述过程信号中元素幅值最大的所述目标稀疏度个元素值保留,同时将所述过程信号中除元素幅值最大的所述目标稀疏度个元素之外的其他所有元素值置为零,并将所述门限函数的处理结果赋值给所述第二迭代信号;Combining the process signal and the target sparsity, a second iterative signal is obtained through a threshold function, the threshold function is used to retain the target sparsity element values with the largest element amplitude in the process signal, and at the same time In the process signal, all element values except the target sparsity elements with the largest element amplitude are set to zero, and the processing result of the threshold function is assigned to the second iteration signal;
根据所述测量信号、所述采样矩阵与所述第二迭代信号,通过非零项统计函数获取汉明距离;所述非零项统计函数用于,获取所述第二迭代信号经低速采样处理及符号量化处理后与所述测量信号的差值,并统计所述差值中非零元素的个数,将所述非零元素的个数赋值给所述汉明距离;所述汉明距离为所述第二迭代信号经低速采样处理及符号量化处理后与所述测量信号对应位置元素值不相同的元素个数;According to the measurement signal, the sampling matrix and the second iteration signal, the Hamming distance is obtained through a non-zero statistical function; the non-zero statistical function is used to obtain the second iteration signal through low-speed sampling processing and the difference between the measurement signal after sign quantization processing and counting the number of non-zero elements in the difference, assigning the number of non-zero elements to the Hamming distance; the Hamming distance is the number of elements whose value is different from that of the measurement signal at the corresponding position of the second iterative signal after low-speed sampling processing and sign quantization processing;
若所述汉明距离大于预设门限值,则将所述第一迭代信号赋值为所述第二迭代信号,并重新执行所述信号重建操作;或,若所述汉明距离小于或等于所述预设门限值,则对所述第二迭代信号进行单位归一化处理,得到所述原始稀疏信号的所述重建信号。If the Hamming distance is greater than a preset threshold value, assign the first iterative signal to the second iterative signal, and re-execute the signal reconstruction operation; or, if the Hamming distance is less than or equal to For the preset threshold value, unit normalization processing is performed on the second iteration signal to obtain the reconstructed signal of the original sparse signal.
第二方面,提供一种信号处理装置,其特征在于,所述装置包括:In a second aspect, a signal processing device is provided, wherein the device includes:
参数获取单元,用于获取输入参数与编码端数据,并对过程参数与第一迭代信号进行初始化,所述编码端数据包括编码端对原始稀疏信号进行低速采样处理及符号量化处理获取的测量信号;A parameter acquisition unit, configured to acquire input parameters and encoding end data, and initialize process parameters and the first iteration signal. The encoding end data includes measurement signals acquired by the encoding end through low-speed sampling processing and symbol quantization processing on the original sparse signal ;
稀疏度估计单元,用于根据所述输入参数、所述编码端数据、初始化后的所述过程参数和所述第一迭代信号进行自适应稀疏度估计操作,获得目标稀疏度;A sparsity estimation unit, configured to perform an adaptive sparsity estimation operation according to the input parameters, the encoding end data, the initialized process parameters, and the first iteration signal to obtain a target sparsity;
信号重建单元,用于根据所述目标稀疏度,对所述测量信号进行信号重建得到所述原始稀疏信号的重建信号。The signal reconstruction unit is configured to perform signal reconstruction on the measurement signal according to the target sparsity to obtain a reconstructed signal of the original sparse signal.
结合第二方面,在第一种可能的实现方式中,所述参数获取单元具体用于:With reference to the second aspect, in a first possible implementation manner, the parameter acquisition unit is specifically configured to:
获取所述输入参数,所述输入参数包括最大稀疏度、判决门限;Acquiring the input parameters, the input parameters include a maximum sparsity and a decision threshold;
接收所述编码端传递的所述编码端数据,所述编码端数据包括所述测量信号、采样矩阵;所述测量信号为所述编码端结合所述采样矩阵对所述原始稀疏信号进行低速采样处理获取采样信号后,再对所述采样信号进行符号量化处理获得的。Receive the encoding end data transmitted by the encoding end, the encoding end data includes the measurement signal and the sampling matrix; the measurement signal is the low speed sampling of the original sparse signal by the encoding end in combination with the sampling matrix After the sampling signal is acquired, the sampling signal is obtained by performing sign quantization processing on the sampling signal.
结合第二方面或第二方面的第一种可能的实现方式,在第二种可能的实现方式中,所述过程参数包括所述目标稀疏度的第一估计值、第一使用值,所述参数获取单元还包括初始化单元,具体用于:With reference to the second aspect or the first possible implementation manner of the second aspect, in the second possible implementation manner, the process parameters include a first estimated value and a first used value of the target sparsity, and the The parameter acquisition unit also includes an initialization unit, which is specifically used for:
将所述第一估计值赋值为所述最大稀疏度,将所述第一使用值赋值为所述最大稀疏度;assigning the first estimated value as the maximum sparsity, and assigning the first use value as the maximum sparsity;
将所述第一迭代信号的初始值赋值为零向量。Assigning the initial value of the first iteration signal as a zero vector.
结合第二方面至第二方面的第二种可能的实现方式中的任一种,在第三种可能的实现方式中,所述稀疏度估计单元用于执行所述自适应稀疏度估计操作,具体包括:With reference to any one of the second aspect to the second possible implementation manner of the second aspect, in a third possible implementation manner, the sparsity estimation unit is configured to perform the adaptive sparsity estimation operation, Specifically include:
根据所述测量信号与所述采样矩阵,获取所述第一迭代信号经低速采样及符号量化后与所述测量信号的差值,并对所述差值左乘所述采样矩阵的转置矩阵,得到第一迭代梯度;According to the measurement signal and the sampling matrix, obtain the difference between the first iteration signal and the measurement signal after low-speed sampling and symbol quantization, and multiply the difference by the transpose matrix of the sampling matrix on the left , get the gradient of the first iteration;
通过梯度下降法获取过程信号,所述过程信号为所述第一迭代信号与所述第一迭代梯度的差值;Obtaining a process signal by a gradient descent method, where the process signal is a difference between the first iteration signal and the first iteration gradient;
对所述第一迭代梯度是否趋近于零向量进行判断,若所述第一迭代梯度趋近于零向量,则将所述目标稀疏度的第二使用值赋值为所述第一估计值;或,若所述第一迭代梯度不趋近于零向量,则将所述第二使用值赋值为所述第一使用值;Judging whether the first iterative gradient approaches a zero vector, and if the first iterative gradient approaches a zero vector, assigning the second used value of the target sparsity as the first estimated value; Or, if the first iterative gradient does not approach the zero vector, assigning the second use value as the first use value;
根据所述过程信号与所述第二使用值,通过门限函数获取第二迭代信号,所述门限函数用于将所述过程信号中元素幅值最大的所述第二使用值个元素值保留,同时将所述过程信号中除元素幅值最大的所述第二使用值个元素之外的其他所有元素值置为零,并将所述门限函数的处理结果赋值给所述第二迭代信号;According to the process signal and the second use value, a second iterative signal is obtained through a threshold function, and the threshold function is used to reserve the second use value and element values with the largest element amplitude in the process signal, At the same time, all other element values in the process signal except the second use value elements with the largest element amplitude are set to zero, and the processing result of the threshold function is assigned to the second iteration signal;
对所述第二迭代信号进行单位归一化,得到归一化迭代信号,统计所述归一化迭代信号中元素幅值绝对值超过所述判决门限的非零元素个数,并将所述元素个数赋值给所述目标稀疏度的第二估计值;Carrying out unit normalization on the second iterative signal to obtain a normalized iterative signal, counting the number of non-zero elements whose absolute value of the element amplitude in the normalized iterative signal exceeds the decision threshold, and calculating the Assigning the number of elements to the second estimated value of the target sparsity;
若所述第一迭代梯度不趋近于零向量或所述第二估计值与所述第一使用值不相等,则将所述第一估计值赋值为所述第二估计值、所述第一使用值赋值为所述第二使用值、所述第一迭代信号赋值为所述第二迭代信号,并重新执行所述自适应稀疏度估计操作;或,若所述第一迭代梯度趋近于零向量且所述第二估计值与所述第一使用值相等,则将所述目标稀疏度赋值为所述第二估计值,得到所述目标稀疏度。If the first iterative gradient does not approach the zero vector or the second estimated value is not equal to the first used value, assign the first estimated value to the second estimated value, the first estimated value A use value is assigned as the second use value, the first iteration signal is assigned as the second iteration signal, and the adaptive sparsity estimation operation is re-executed; or, if the gradient of the first iteration approaches If it is a zero vector and the second estimated value is equal to the first used value, assign the target sparsity as the second estimated value to obtain the target sparsity.
结合第二方面至第二方面的第三种可能的实现方式中的任一种,在第四种可能的实现方式中,所述信号重建单元具体用于:With reference to any one of the second aspect to the third possible implementation manner of the second aspect, in a fourth possible implementation manner, the signal reconstruction unit is specifically configured to:
根据所述编码端数据、所述第一迭代信号与所述目标稀疏度,执行信号重建操作,所述信号重建操作具体包括:Perform a signal reconstruction operation according to the encoding end data, the first iteration signal, and the target sparsity, and the signal reconstruction operation specifically includes:
根据所述测量信号与所述采样矩阵,获取所述第一迭代信号经低速采样及符号量化后与所述测量信号的差值,并对所述差值左乘所述采样矩阵的转置矩阵,得到第一迭代梯度;According to the measurement signal and the sampling matrix, obtain the difference between the first iteration signal and the measurement signal after low-speed sampling and symbol quantization, and multiply the difference by the transpose matrix of the sampling matrix on the left , get the gradient of the first iteration;
通过梯度下降法获取过程信号,所述过程信号为所述第一迭代信号与所述第一迭代梯度的差值;Obtaining a process signal by a gradient descent method, where the process signal is a difference between the first iteration signal and the first iteration gradient;
结合所述过程信号与所述目标稀疏度,通过门限函数获取第二迭代信号,所述门限函数用于将所述过程信号中元素幅值最大的所述目标稀疏度个元素值保留,同时将所述过程信号中除元素幅值最大的所述目标稀疏度个元素之外的其他所有元素值置为零,并将所述门限函数的处理结果赋值给所述第二迭代信号;Combining the process signal and the target sparsity, a second iterative signal is obtained through a threshold function, the threshold function is used to retain the target sparsity element values with the largest element amplitude in the process signal, and at the same time In the process signal, all element values except the target sparsity elements with the largest element amplitude are set to zero, and the processing result of the threshold function is assigned to the second iteration signal;
根据所述测量信号、所述采样矩阵与所述第二迭代信号,通过非零项统计函数获取汉明距离;所述非零项统计函数用于,获取所述第二迭代信号经低速采样处理及符号量化处理后与所述测量信号的差值,并统计所述差值中非零元素的个数,将所述非零元素的个数赋值给所述汉明距离;所述汉明距离为所述第二迭代信号经低速采样处理及符号量化处理后与所述测量信号对应位置元素值不相同的元素个数;According to the measurement signal, the sampling matrix and the second iteration signal, the Hamming distance is obtained through a non-zero statistical function; the non-zero statistical function is used to obtain the second iteration signal through low-speed sampling processing and the difference between the measurement signal after sign quantization processing and counting the number of non-zero elements in the difference, assigning the number of non-zero elements to the Hamming distance; the Hamming distance is the number of elements whose value is different from that of the measurement signal at the corresponding position of the second iterative signal after low-speed sampling processing and sign quantization processing;
若所述汉明距离大于预设门限值,则将所述第一迭代信号赋值为所述第二迭代信号,并重新执行所述信号重建操作;或,若所述汉明距离小于或等于所述预设门限值,则对所述第二迭代信号进行单位归一化处理,得到所述原始稀疏信号的所述重建信号。If the Hamming distance is greater than a preset threshold value, assign the first iterative signal to the second iterative signal, and re-execute the signal reconstruction operation; or, if the Hamming distance is less than or equal to For the preset threshold value, unit normalization processing is performed on the second iteration signal to obtain the reconstructed signal of the original sparse signal.
第三方面,提供一种信号处理装置,所述信号处理装置包括:总线,以及连接到所述总线的处理器、存储器和接口,其中所述接口用于与外部设备进行通信;所述存储器用于存储指令,所述处理器执行所述指令用于:In a third aspect, a signal processing device is provided, and the signal processing device includes: a bus, a processor connected to the bus, a memory, and an interface, wherein the interface is used for communicating with external devices; the memory uses In order to store instructions, the processor executes the instructions for:
获取输入参数与编码端数据,并对过程参数与第一迭代信号进行初始化,所述编码端数据包括编码端对原始稀疏信号进行低速采样处理及符号量化处理获取的测量信号;Acquiring input parameters and encoding end data, and initializing the process parameters and the first iteration signal, the encoding end data including the measurement signal obtained by the encoding end performing low-speed sampling processing and symbol quantization processing on the original sparse signal;
根据所述输入参数、所述编码端数据、初始化后的所述过程参数和所述第一迭代信号进行自适应稀疏度估计操作,获得目标稀疏度;performing an adaptive sparsity estimation operation according to the input parameters, the encoding end data, the initialized process parameters, and the first iteration signal to obtain a target sparsity;
根据所述目标稀疏度,对所述测量信号进行信号重建得到所述原始稀疏信号的重建信号。Perform signal reconstruction on the measurement signal according to the target sparsity to obtain a reconstructed signal of the original sparse signal.
结合第三方面,在第一种可能的实现方式中,所述处理器执行所述指令具体用于:With reference to the third aspect, in a first possible implementation manner, the processor executes the instruction specifically for:
获取所述输入参数,所述输入参数包括最大稀疏度、判决门限;Acquiring the input parameters, the input parameters include a maximum sparsity and a decision threshold;
接收所述编码端传递的所述编码端数据,所述编码端数据包括所述测量信号、采样矩阵;所述测量信号为所述编码端结合所述采样矩阵对所述原始稀疏信号进行低速采样处理获取采样信号后,再对所述采样信号进行符号量化处理获得的。Receive the encoding end data transmitted by the encoding end, the encoding end data includes the measurement signal and the sampling matrix; the measurement signal is the low speed sampling of the original sparse signal by the encoding end in combination with the sampling matrix After the sampling signal is acquired, the sampling signal is obtained by performing sign quantization processing on the sampling signal.
结合第三方面或第三方面的第一种可能的实现方式,在第二种可能的实现方式中,所述处理器执行所述指令具体用于:With reference to the third aspect or the first possible implementation manner of the third aspect, in a second possible implementation manner, the processor executes the instruction specifically for:
将所述第一估计值赋值为所述最大稀疏度,将所述第一使用值赋值为所述最大稀疏度;assigning the first estimated value as the maximum sparsity, and assigning the first use value as the maximum sparsity;
将所述第一迭代信号的初始值赋值为零向量。Assigning the initial value of the first iteration signal as a zero vector.
结合第三方面至第三方面的第二种可能的实现方式中的任一种,在第三种可能的实现方式中,所述处理器执行所述指令还具体用于:With reference to any one of the third aspect to the second possible implementation manner of the third aspect, in a third possible implementation manner, the processor executing the instruction is further specifically used for:
根据所述测量信号与所述采样矩阵,获取所述第一迭代信号经低速采样及符号量化后与所述测量信号的差值,并对所述差值左乘所述采样矩阵的转置矩阵,得到第一迭代梯度;According to the measurement signal and the sampling matrix, obtain the difference between the first iteration signal and the measurement signal after low-speed sampling and symbol quantization, and multiply the difference by the transpose matrix of the sampling matrix on the left , get the gradient of the first iteration;
通过梯度下降法获取过程信号,所述过程信号为所述第一迭代信号与所述第一迭代梯度的差值;Obtaining a process signal by a gradient descent method, where the process signal is a difference between the first iteration signal and the first iteration gradient;
对所述第一迭代梯度是否趋近于零向量进行判断,若所述第一迭代梯度趋近于零向量,则将所述目标稀疏度的第二使用值赋值为所述第一估计值;或,若所述第一迭代梯度不趋近于零向量,则将所述第二使用值赋值为所述第一使用值;Judging whether the first iterative gradient approaches a zero vector, and if the first iterative gradient approaches a zero vector, assigning the second used value of the target sparsity as the first estimated value; Or, if the first iterative gradient does not approach the zero vector, assigning the second use value as the first use value;
根据所述过程信号与所述第二使用值,通过门限函数获取第二迭代信号,所述门限函数用于将所述过程信号中元素幅值最大的所述第二使用值个元素值保留,同时将所述过程信号中除元素幅值最大的所述第二使用值个元素之外的其他所有元素值置为零,并将所述门限函数的处理结果赋值给所述第二迭代信号;According to the process signal and the second use value, a second iterative signal is obtained through a threshold function, and the threshold function is used to reserve the second use value and element values with the largest element amplitude in the process signal, At the same time, all other element values in the process signal except the second use value elements with the largest element amplitude are set to zero, and the processing result of the threshold function is assigned to the second iteration signal;
对所述第二迭代信号进行单位归一化,得到归一化迭代信号,统计所述归一化迭代信号中元素幅值绝对值超过所述判决门限的非零元素个数,并将所述元素个数赋值给所述目标稀疏度的第二估计值;Carrying out unit normalization on the second iterative signal to obtain a normalized iterative signal, counting the number of non-zero elements whose absolute value of the element amplitude in the normalized iterative signal exceeds the decision threshold, and calculating the Assigning the number of elements to the second estimated value of the target sparsity;
若所述第一迭代梯度不趋近于零向量或所述第二估计值与所述第一使用值不相等,则将所述第一估计值赋值为所述第二估计值、所述第一使用值赋值为所述第二使用值、所述第一迭代信号赋值为所述第二迭代信号,并重新执行所述自适应稀疏度估计操作;或,若所述第一迭代梯度趋近于零向量且所述第二估计值与所述第一使用值相等,则将所述目标稀疏度赋值为所述第二估计值,得到所述目标稀疏度。If the first iterative gradient does not approach the zero vector or the second estimated value is not equal to the first used value, assign the first estimated value to the second estimated value, the first estimated value A use value is assigned as the second use value, the first iteration signal is assigned as the second iteration signal, and the adaptive sparsity estimation operation is re-executed; or, if the gradient of the first iteration approaches If it is a zero vector and the second estimated value is equal to the first used value, assign the target sparsity as the second estimated value to obtain the target sparsity.
结合第三方面至第三方面的第三种可能的实现方式中的任一种,在第四种可能的实现方式中,所述处理器执行所述指令还具体用于:With reference to any one of the third aspect to the third possible implementation manner of the third aspect, in a fourth possible implementation manner, the processor executing the instruction is further specifically used for:
根据所述编码端数据、所述第一迭代信号与所述目标稀疏度,执行信号重建操作,所述信号重建操作具体包括:Perform a signal reconstruction operation according to the encoding end data, the first iteration signal, and the target sparsity, and the signal reconstruction operation specifically includes:
根据所述测量信号与所述采样矩阵,获取所述第一迭代信号经低速采样及符号量化后与所述测量信号的差值,并对所述差值左乘所述采样矩阵的转置矩阵,得到第一迭代梯度;According to the measurement signal and the sampling matrix, obtain the difference between the first iteration signal and the measurement signal after low-speed sampling and symbol quantization, and multiply the difference by the transpose matrix of the sampling matrix on the left , get the gradient of the first iteration;
通过梯度下降法获取过程信号,所述过程信号为所述第一迭代信号与所述第一迭代梯度的差值;Obtaining a process signal by a gradient descent method, where the process signal is a difference between the first iteration signal and the first iteration gradient;
结合所述过程信号与所述目标稀疏度,通过门限函数获取第二迭代信号,所述门限函数用于将所述过程信号中元素幅值最大的所述目标稀疏度个元素值保留,同时将所述过程信号中除元素幅值最大的所述目标稀疏度个元素之外的其他所有元素值置为零,并将所述门限函数的处理结果赋值给所述第二迭代信号;Combining the process signal and the target sparsity, a second iterative signal is obtained through a threshold function, the threshold function is used to retain the target sparsity element values with the largest element amplitude in the process signal, and at the same time In the process signal, all element values except the target sparsity elements with the largest element amplitude are set to zero, and the processing result of the threshold function is assigned to the second iteration signal;
根据所述测量信号、所述采样矩阵与所述第二迭代信号,通过非零项统计函数获取汉明距离;所述非零项统计函数用于,获取所述第二迭代信号经低速采样处理及符号量化处理后与所述测量信号的差值,并统计所述差值中非零元素的个数,将所述非零元素的个数赋值给所述汉明距离;所述汉明距离为所述第二迭代信号经低速采样处理及符号量化处理后与所述测量信号对应位置元素值不相同的元素个数;According to the measurement signal, the sampling matrix and the second iteration signal, the Hamming distance is obtained through a non-zero statistical function; the non-zero statistical function is used to obtain the second iteration signal through low-speed sampling processing and the difference between the measurement signal after sign quantization processing and counting the number of non-zero elements in the difference, assigning the number of non-zero elements to the Hamming distance; the Hamming distance is the number of elements whose value is different from that of the measurement signal at the corresponding position of the second iterative signal after low-speed sampling processing and sign quantization processing;
若所述汉明距离大于预设门限值,则将所述第一迭代信号赋值为所述第二迭代信号,并重新执行所述信号重建操作;或,若所述汉明距离小于或等于所述预设门限值,则对所述第二迭代信号进行单位归一化处理,得到所述原始稀疏信号的所述重建信号。If the Hamming distance is greater than a preset threshold value, assign the first iterative signal to the second iterative signal, and re-execute the signal reconstruction operation; or, if the Hamming distance is less than or equal to For the preset threshold value, unit normalization processing is performed on the second iteration signal to obtain the reconstructed signal of the original sparse signal.
本发明的实施例提供一种信号处理方法和装置,通过获取输入参数与编码端数据,在对过程参数与相关参数值初始化后,进行自适应稀疏度估计操作,并在获取了目标稀疏度后,根据目标稀疏度进行信号重建得到原始稀疏信号的重建信号。这样,解决了稀疏度先验信息受限导致信号重建性能降低的问题,降低或基本消除了重建信号产生的偏差,提高了单比特压缩感知技术的实用性与准确度。Embodiments of the present invention provide a signal processing method and device. By acquiring input parameters and encoding end data, after initializing process parameters and related parameter values, an adaptive sparsity estimation operation is performed, and after acquiring the target sparsity , reconstruct the signal according to the target sparsity to obtain the reconstructed signal of the original sparse signal. In this way, the problem of signal reconstruction performance degradation due to limited sparsity prior information is solved, the deviation caused by reconstruction signal is reduced or basically eliminated, and the practicability and accuracy of single-bit compressed sensing technology are improved.
附图说明Description of drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following will briefly introduce the drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only These are some embodiments of the present invention. Those skilled in the art can also obtain other drawings based on these drawings without creative work.
图1为本发明实施例提供的一种信号处理方法的流程示意图;FIG. 1 is a schematic flowchart of a signal processing method provided by an embodiment of the present invention;
图2为本发明实施例提供的另一种信号处理方法的流程示意图一;FIG. 2 is a first schematic flowchart of another signal processing method provided by an embodiment of the present invention;
图3为本发明实施例提供的另一种信号处理方法的流程示意图二;FIG. 3 is a second schematic flow diagram of another signal processing method provided by an embodiment of the present invention;
图4为本发明实施例提供的另一种信号处理方法的流程示意图三;FIG. 4 is a third schematic flowchart of another signal processing method provided by an embodiment of the present invention;
图5为本发明实施例提供的一种信号处理装置的结构示意图一;FIG. 5 is a first structural schematic diagram of a signal processing device provided by an embodiment of the present invention;
图6为为本发明实施例提供的一种信号处理装置的结构示意图二;FIG. 6 is a second schematic structural diagram of a signal processing device provided by an embodiment of the present invention;
图7为为本发明实施例提供的一种信号处理装置的结构示意图三。FIG. 7 is a schematic structural diagram III of a signal processing device provided by an embodiment of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some, not all, embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without creative efforts fall within the protection scope of the present invention.
本发明的实施例提供一种信号处理方法,如图1所示,该方法包括:Embodiments of the present invention provide a signal processing method, as shown in FIG. 1, the method includes:
S101、获取输入参数与编码端数据,并对过程参数与第一迭代信号进行初始化;S101. Obtain input parameters and encoding end data, and initialize the process parameters and the first iteration signal;
其中,输入参数包括最大稀疏度、判决门限;编码端数据包括测量信号、采样矩阵;过程参数包括目标稀疏度的第一估计值、第一使用值。Wherein, the input parameters include the maximum sparsity and the decision threshold; the data at the encoding end include the measurement signal and the sampling matrix; the process parameters include the first estimated value and the first use value of the target sparsity.
最大稀疏度、判决门限与采样矩阵为与目标稀疏度估计相关的预给定计算参数;第一迭代信号、第一估计值与第一使用值为参与计算的中间变量;测量信号为编码端对原始稀疏信号进行低速采样处理及符号量化处理后获得的。The maximum sparsity, decision threshold, and sampling matrix are predetermined calculation parameters related to target sparsity estimation; the first iteration signal, the first estimated value, and the first used value are intermediate variables involved in the calculation; the measurement signal is the encoding end pair The original sparse signal is obtained after low-speed sampling processing and symbol quantization processing.
值得一提的,原始稀疏信号可以为视频信号、音频信号、图像处理信号、信道估计信号、无线传感器网络信号、认知无线电频谱检测信号等,且上述列举的原始稀疏信号仅为示例性的,包括但不限于此。It is worth mentioning that the original sparse signal can be video signal, audio signal, image processing signal, channel estimation signal, wireless sensor network signal, cognitive radio spectrum detection signal, etc., and the original sparse signal listed above is only exemplary, Including but not limited to this.
S102、根据输入参数、编码端数据、初始化后的过程参数和第一迭代信号进行自适应稀疏度估计操作,获取目标稀疏度;S102. Perform an adaptive sparsity estimation operation according to the input parameters, encoding end data, initialized process parameters, and the first iteration signal to obtain the target sparsity;
S103、根据目标稀疏度,对测量信号进行信号重建得到原始稀疏信号的重建信号。S103. Perform signal reconstruction on the measurement signal according to the target sparsity to obtain a reconstructed signal of the original sparse signal.
本发明的实施例提供一种信号处理方法,通过获取输入参数与编码端数据,在对过程参数与相关参数值初始化后,进行自适应稀疏度估计操作,并在获取了目标稀疏度后,根据目标稀疏度进行信号重建得到原始稀疏信号的重建信号。这样,解决了稀疏度先验信息受限导致信号重建性能降低的问题,降低或基本消除了重建信号产生的偏差,提高了单比特压缩感知技术的实用性与准确度。Embodiments of the present invention provide a signal processing method. By acquiring input parameters and encoding end data, after initializing process parameters and related parameter values, an adaptive sparsity estimation operation is performed, and after acquiring the target sparsity, according to The target sparsity is used to reconstruct the signal to obtain the reconstructed signal of the original sparse signal. In this way, the problem of signal reconstruction performance degradation due to limited sparsity prior information is solved, the deviation caused by reconstruction signal is reduced or basically eliminated, and the practicability and accuracy of single-bit compressed sensing technology are improved.
为了使本领域技术人员能够更清楚地理解本发明实施例提供的技术方案,下面通过具体的实施例,对本发明实施例提供的另一种信号处理方法进行详细说明,如图2所示,该方法包括:In order to enable those skilled in the art to more clearly understand the technical solutions provided by the embodiments of the present invention, another signal processing method provided by the embodiments of the present invention will be described in detail below through specific embodiments, as shown in FIG. 2 , the Methods include:
S201、获取输入参数与编码端数据。S201. Obtain input parameters and encoding end data.
具体的,获取的输入参数包括最大稀疏度Km、判决门限λ;其中,最大稀疏度Km与判决门限λ可以为系统缺省设置,也可以为外部输入给定,此处不做限定。Specifically, the acquired input parameters include the maximum sparsity Km and the decision threshold λ; wherein, the maximum sparsity Km and the decision threshold λ can be the default settings of the system, or can be given by external input, which are not limited here.
接收编码端传递的编码端数据,编码端数据包括测量信号y、采样矩阵Φ;其中,测量信号y为编码端在对原始稀疏信号θ进行低速采样处理,并对得到的采样信号x进行符号量化处理后获得的,包含了可用单比特表示的符号信息,计算公式为:Receive the encoding end data transmitted by the encoding end. The encoding end data includes the measurement signal y and the sampling matrix Φ; where the measurement signal y is the low-speed sampling process of the original sparse signal θ at the encoding end, and the symbol quantization is performed on the obtained sampling signal x After processing, it contains symbol information that can be represented by a single bit, and the calculation formula is:
y=sign(Φθ),y=sign(Φθ),
其中,θ为原始高维输入信号,表示为一个N×1的列向量,N表示信号维度;值得一提的,作为应用压缩感知技术的前提,认定该输入信号θ具有稀疏性,即输入信号θ中除仅有的K个非零元素外其他元素均为零,且K值远远小于信号维度N。采样矩阵Φ表示为一个M×N的矩阵,且行数M小于列数N。输入信号θ左乘以采样矩阵Φ为,可实现低速采样的作用,将输入信号θ从原有的N维降维映射到M维,得到降维后的采样信号x=Φθ,然后对采样信号x进行单比特量化得到测量信号y,具体的,sign()为单比特符号量化函数,用于对采样信号x进行单比特符号量化处理,将采样信号x中的正值量化为1、负值量化为-1,并分别用比特1、比特0表示,这样测量信号y中的每一个测量值都可通过一个比特来表示。Among them, θ is the original high-dimensional input signal, expressed as an N×1 column vector, and N represents the signal dimension; it is worth mentioning that, as a prerequisite for applying compressed sensing technology, the input signal θ is considered to be sparse, that is, the input signal Except for the only K non-zero elements in θ, the other elements are all zero, and the K value is much smaller than the signal dimension N. The sampling matrix Φ is expressed as an M×N matrix, and the number of rows M is smaller than the number of columns N. The input signal θ is multiplied by the sampling matrix Φ to the left, which can realize the function of low-speed sampling. The input signal θ is mapped from the original N-dimensional dimension reduction to M-dimensional, and the dimension-reduced sampling signal x=Φθ is obtained, and then the sampling signal Perform single-bit quantization on x to obtain measurement signal y. Specifically, sign() is a single-bit sign quantization function, which is used to perform single-bit sign quantization processing on sampled signal x, and quantize positive values in sampled signal x to 1 and negative values. Quantized to -1, and represented by bit 1 and bit 0 respectively, so that each measurement value in the measurement signal y can be represented by one bit.
示例性的,上述接收编码端传递的编码端数据可以为译码端接收编码端传递的编码端数据,其中,译码端与编码端可以为相互独立的两个设备,以有线和/或无线的方式进行传输,也可以为一个设备中的两个功能模块,以有线的方式进行传输,此处不做限定。Exemplarily, the encoding end data transmitted by the receiving encoding end may be the encoding end data transmitted by the decoding end receiving the encoding end, wherein the decoding end and the encoding end may be two independent devices, wired and/or wireless It can also be transmitted in a wired manner for two functional modules in one device, which is not limited here.
值得一提的,原始稀疏信号θ可以为视频信号、音频信号、图像处理信号、信道估计信号、无线传感器网络信号、认知无线电频谱检测信号等,且上述列举的原始稀疏信号仅为示例性的,包括但不限于此。It is worth mentioning that the original sparse signal θ can be video signal, audio signal, image processing signal, channel estimation signal, wireless sensor network signal, cognitive radio spectrum detection signal, etc., and the original sparse signal listed above is only exemplary , including but not limited to.
S202、对过程参数与第一迭代信号的初始值进行初始化。S202. Initialize process parameters and initial values of the first iteration signal.
具体的,过程参数包括目标稀疏度K的第一估计值pre_Est_S、第一使用值pre_Use_S。初始化时,将第一估计值pre_Est_S与第一使用值pre_Use_S都赋值为输入参数中的最大稀疏度Km,将第一迭代信号的初始值赋值为零向量。Specifically, the process parameters include a first estimated value pre_Est_S of the target sparsity K and a first usage value pre_Use_S. During initialization, both the first estimated value pre_Est_S and the first use value pre_Use_S are assigned as the maximum sparsity K m in the input parameters, and the first iteration signal initial value of Assignment is a zero vector.
其中,上述对过程参数与第一迭代信号进行初始化时的初始值仅为示例性的,还可根据实际应用选择其它可用初始值进行初始化,包括但不限于此。Among them, the above pair of process parameters and the first iteration signal The initial value during initialization is only exemplary, and other available initial values can also be selected for initialization according to actual applications, including but not limited to this.
S203、根据输入参数、编码端数据、初始化后的过程参数和第一迭代信号进行自适应稀疏度估计操作,获取目标稀疏度。S203. Perform an adaptive sparsity estimation operation according to the input parameters, the encoding end data, the initialized process parameters, and the first iteration signal to obtain a target sparsity.
具体可以为,如图3所示,执行自适应稀疏度估计操作,包括:Specifically, as shown in FIG. 3, performing an adaptive sparsity estimation operation includes:
S2031、计算获取第一迭代梯度。S2031. Calculate and obtain the first iteration gradient.
具体的,根据S201中接收的测量信号y与采样矩阵Φ,以及第一迭代信号计算当前迭代过程对应的第一迭代梯度其中,第一迭代信号为上一次迭代过程更新得到的信号列向量,若为第一次进行该迭代过程,则第一迭代信号即为S202中进行初始化后的第一迭代信号初始值 Specifically, according to the measurement signal y and the sampling matrix Φ received in S201, and the first iteration signal Calculate the first iteration gradient corresponding to the current iteration process Among them, the first iteration signal It is the signal column vector obtained by updating the last iteration process. If the iteration process is performed for the first time, the first iteration signal That is, the initial value of the first iteration signal after initialization in S202
第一迭代梯度为反映迭代更新过程中迭代信号变化率的参数,计算公式为:first iteration gradient In order to reflect the parameters of the iterative signal change rate in the iterative update process, the calculation formula is:
具体的,结合采样矩阵Φ,在对第一迭代信号进行低速采样处理及符号量化处理后,获取与测量信号y的差值,并对该差值左乘采样矩阵Φ的转置矩阵,得到第一迭代梯度 Specifically, combined with the sampling matrix Φ, in the first iterative signal After low-speed sampling processing and symbol quantization processing, the difference with the measurement signal y is obtained, and the difference is multiplied by the transpose matrix of the sampling matrix Φ to the left to obtain the first iteration gradient
S2032、计算获取过程信号。S2032. Calculate and acquire the process signal.
具体的,通过梯度下降法,根据第一迭代梯度与第一迭代信号沿梯度下降的方向更新迭代中间过程信号at,计算公式为:Specifically, through the gradient descent method, according to the first iteration gradient with the first iteration signal Update the iterative intermediate process signal at along the direction of gradient descent , the calculation formula is:
S2033、判断第一迭代梯度是否趋近于零向量。S2033. Determine whether the gradient of the first iteration approaches the zero vector.
具体的,若第一迭代梯度趋近于零,则执行S2034;或,若第一迭代梯度不趋近于零,则执行S2035。Specifically, if the first iteration gradient tends to zero, execute S2034; or, if the gradient of the first iteration If it is not close to zero, execute S2035.
S2034、将第二使用值更新为第一估计值。S2034. Update the second usage value to the first estimated value.
具体的,目标稀疏度K的第二使用值cur_Use_S用于在后续更新迭代信号过程中,作为门限函数的阈值。Specifically, the second use value cur_Use_S of the target sparsity K is used as a threshold of the threshold function in the subsequent process of updating the iterative signal.
若第一迭代梯度趋近于零,则将第二使用值cur_Use_S更新为上一次迭代过程中稀疏度估计的结果,即为目标稀疏度K的第一估计值pre_Est_S;若为第一次进行该迭代过程,则第一估计值pre_Est_S等于输入参数给定的最大稀疏度Km。If the first iteration gradient approaching zero, update the second use value cur_Use_S to the result of the sparsity estimation in the last iteration process, which is the first estimated value pre_Est_S of the target sparsity K; if this iteration process is performed for the first time, the second An estimated value pre_Est_S is equal to the maximum sparsity K m given by the input parameters.
S2035、将第二使用值更新为第一使用值。S2035. Update the second usage value to the first usage value.
具体的,目标稀疏度K的第二使用值cur_Use_S用于在后续更新迭代信号过程中,作为门限函数的阈值。Specifically, the second use value cur_Use_S of the target sparsity K is used as a threshold of the threshold function in the subsequent process of updating the iterative signal.
若第一迭代梯度不趋近于零,则第二使用值cur_Use_S沿用上一次迭代过程中门限函数使用的阈值,即将第二使用值cur_Use_S赋值为目标稀疏度K的第一使用值pre_Use_S;若为第一次进行该迭代过程,则第一使用值pre_Use_S等于输入参数给定的最大稀疏度Km。If the first iteration gradient If it is not close to zero, the second usage value cur_Use_S follows the threshold value used by the threshold function in the last iteration, that is, the second usage value cur_Use_S is assigned as the first usage value pre_Use_S of the target sparsity K; In an iterative process, the first use value pre_Use_S is equal to the maximum sparsity K m given by the input parameter.
S2036、通过门限函数执行硬阈值操作获取第二迭代信号。S2036. Acquire a second iteration signal by performing a hard threshold operation through a threshold function.
具体的,根据过程信号at与第二使用值cur_Use_S,通过门限函数获取第二迭代信号计算公式为:Specifically, according to the process signal at and the second use value cur_Use_S , the second iteration signal is obtained through the threshold function The calculation formula is:
其中,Γ()表示门限函数,执行硬阈值操作用于将过程信号at中元素幅值最大的第二使用值cur_Use_S个元素值保留,同时将过程信号at中除元素幅值最大的第二使用值cur_Use_S个元素之外的其他所有元素值置为零,并将门限函数Γ()的处理结果赋值给第二迭代信号 Among them, Γ() represents the threshold function, and the hard threshold operation is used to retain the second use value cur_Use_S element values with the largest element amplitude in the process signal at, and at the same time to remove the first element with the largest amplitude in the process signal at 2. Set the value of all other elements except cur_Use_S elements to zero, and assign the processing result of the threshold function Γ() to the second iteration signal
S2037、获取目标稀疏度的第二估计值。S2037. Acquire a second estimated value of the target sparsity.
具体的,对第二迭代信号进行单位归一化,得到归一化迭代信号计算公式为:Specifically, for the second iteration signal Perform unit normalization to obtain normalized iterative signals The calculation formula is:
其中,||||2表示向量的二范数;然后根据输入参数给定的目标稀疏度判决门限λ,比较并统计归一化迭代信号中元素幅值绝对值超过判决门限λ的非零元素个数,并将该非零元素个数记录为目标稀疏度K的第二估计值cur_Est_S,可表示为:Among them, |||| 2 represents the two-norm of the vector; then, according to the target sparsity decision threshold λ given by the input parameters, compare and count the normalized iteration signals The number of non-zero elements in which the absolute value of the element amplitude exceeds the decision threshold λ, and the number of non-zero elements is recorded as the second estimated value cur_Est_S of the target sparsity K, which can be expressed as:
其中,count()表示计数函数,用于记录向量中满足条件的元素的个数。Among them, count() represents a counting function, which is used to record the number of elements in the vector that satisfy the condition.
S2038、判断第一迭代梯度是否趋近于零且第二估计值与第一使用值是否相等。S2038. Determine whether the first iteration gradient approaches zero and whether the second estimated value is equal to the first used value.
具体的,若第一迭代梯度不趋近于零向量或第二估计值cur_Est_S与第一使用值pre_Use_S不相等,则执行S2039;Specifically, if the first iteration gradient If the vector does not approach zero or the second estimated value cur_Est_S is not equal to the first use value pre_Use_S, execute S2039;
或,若第一迭代梯度趋近于零向量且第二估计值cur_Est_S与第一使用值pre_Use_S相等,则执行S2030。Or, if the first iteration gradient The vector approaches to zero and the second estimated value cur_Est_S is equal to the first usage value pre_Use_S, then execute S2030.
S2039、更新过程参数,继续进行后续迭代。S2039 , update the process parameters, and continue to perform subsequent iterations.
具体的,若第一迭代梯度不趋近于零向量或第二估计值cur_Est_S与第一使用值pre_Use_S不相等,则将第一估计值pre_Est_S赋值为第二估计值cur_Est_S、第一使用值pre_Use_S赋值为第二使用值cur_Use_S、将迭代信号的下角标t加1,即将当前迭代过程中的第二迭代信号作为下一次迭代更新过程中的第一迭代信号参与计算,并重新执行S2031至2038的操作。Specifically, if the first iteration gradient If the vector does not approach zero or the second estimated value cur_Est_S is not equal to the first use value pre_Use_S, then the first estimate value pre_Est_S is assigned as the second estimate value cur_Est_S, the first use value pre_Use_S is assigned as the second use value cur_Use_S, and Add 1 to the subscript t of the iterative signal, that is, the second iterative signal in the current iterative process as the first iteration signal during the next iterative update Participate in the calculation, and re-execute the operations from S2031 to S2038.
S2030、获取目标稀疏度。S2030. Obtain the target sparsity.
具体的,若第一迭代梯度趋近于零向量且第二估计值cur_Est_S与第一使用值pre_Use_S相等,则将目标稀疏度K赋值为第二估计值cur_Est_S,即可得到目标稀疏度K。Specifically, if the first iteration gradient If the vector is close to zero and the second estimated value cur_Est_S is equal to the first use value pre_Use_S, then assign the target sparsity K to the second estimated value cur_Est_S to obtain the target sparsity K.
S204、根据目标稀疏度,对测量信号进行信号重建得到原始稀疏信号的重建信号。S204. Perform signal reconstruction on the measurement signal according to the target sparsity to obtain a reconstructed signal of the original sparse signal.
具体的,根据测量信号y、采样矩阵Φ、第一迭代信号目标稀疏度K,执行信号重建操作获取重建信号,如图4所示,包括:Specifically, according to the measurement signal y, the sampling matrix Φ, the first iteration signal The target sparsity is K, and the signal reconstruction operation is performed to obtain the reconstructed signal, as shown in Figure 4, including:
S2041、计算获取第一迭代梯度。S2041. Calculate and obtain the first iteration gradient.
具体的,根据S201中接收的测量信号y与采样矩阵Φ,以及第一迭代信号计算当前迭代过程对应的第一迭代梯度 Specifically, according to the measurement signal y and the sampling matrix Φ received in S201, and the first iteration signal Calculate the first iteration gradient corresponding to the current iteration process
值得一提的,第一迭代信号可以为S203中最后一次迭代过程结束后保留的值,也可以为S202中进行初始化后的第一迭代信号初始值 It is worth mentioning that the first iteration signals It can be the value retained after the last iteration process in S203, or the initial value of the first iteration signal after initialization in S202
第一迭代梯度为反应迭代更新过程中迭代信号变化率的参数,计算公式为:first iteration gradient In order to reflect the parameters of the iterative signal change rate in the iterative update process, the calculation formula is:
具体的,结合采样矩阵Φ,在对第一迭代信号进行低速采样处理及符号量化处理后,获取与测量信号y的差值,并对该差值左乘以采样矩阵Φ的转置矩阵,得到第一迭代梯度 Specifically, combined with the sampling matrix Φ, in the first iterative signal After low-speed sampling processing and symbol quantization processing, the difference with the measurement signal y is obtained, and the difference is multiplied by the transpose matrix of the sampling matrix Φ to the left to obtain the first iteration gradient
S2042、计算获取过程信号。S2042. Calculate and acquire the process signal.
具体的,通过梯度下降法,根据第一迭代梯度与第一迭代信号沿梯度下降的方向更新迭代中间过程信号at,计算公式为:Specifically, through the gradient descent method, according to the first iteration gradient with the first iteration signal Update the iterative intermediate process signal at along the direction of gradient descent , the calculation formula is:
S2043、根据过程信号与目标稀疏度获取第二迭代信号。S2043. Acquire a second iteration signal according to the process signal and the target sparsity.
具体的,根据S203中获取的目标稀疏度K,与S2042中获取的过程信号at,通过门限函数执行硬阈值操作,获取第二迭代信号计算公式为:Specifically, according to the target sparsity K obtained in S203 and the process signal at obtained in S2042 , the hard threshold operation is performed through the threshold function to obtain the second iteration signal The calculation formula is:
其中,Γ()表示门限函数,执行硬阈值操作用于将过程信号at中元素幅值最大的目标稀疏度K个元素值保留,同时将过程信号at中除元素幅值最大的目标稀疏度K个元素之外的其他所有元素值置为零,并将门限函数Γ()的处理结果赋值给第二迭代信号 Among them, Γ() represents the threshold function, and the hard threshold operation is used to retain the K element values of the target sparsity with the largest element amplitude in the process signal at, and at the same time divide the target sparseness with the largest element amplitude in the process signal at Set the value of all other elements except K elements to zero, and assign the processing result of the threshold function Γ() to the second iteration signal
S2044、计算获取汉明距离。S2044. Calculate and obtain the Hamming distance.
具体的,根据测量信号y、采样矩阵Φ,对第二迭代信号进行与编码端相同的信号处理操作,即先进行低速采样处理再进行符号量化处理,进一步计算与测量信号y的差值,并统计差值中非零元素的个数作为汉明距离hd,计算公式为:Specifically, according to the measurement signal y and the sampling matrix Φ, the second iteration signal Carry out the same signal processing operation as that of the encoder, that is, first perform low-speed sampling processing and then perform symbol quantization processing, and further calculate the difference with the measured signal y, and count the number of non-zero elements in the difference as the Hamming distance hd, calculate The formula is:
其中,函数nnz()表示求解向量中非零元素的个数,汉明距离hd为第二迭代信号经低速采样处理及符号量化处理后得到的与测量信号y对应位置元素值不相同的元素个数,为一个非负整数。Among them, the function nnz() represents the number of non-zero elements in the solution vector, and the Hamming distance hd is the second iteration signal After low-speed sampling processing and symbol quantization processing The number of elements whose value is different from that of the corresponding position element of the measurement signal y, which is a non-negative integer.
S2045、判断汉明距离是否大于预设门限值。S2045. Determine whether the Hamming distance is greater than a preset threshold.
具体的,对汉明距离是否大于预设门限值进行判断,若汉明距离hd大于预设门限值,则执行S2046;或,若汉明距离hd小于或等于预设门限值,则执行S2047。Specifically, it is judged whether the Hamming distance is greater than the preset threshold value, if the Hamming distance hd is greater than the preset threshold value, then execute S2046; or, if the Hamming distance hd is less than or equal to the preset threshold value, then Execute S2047.
其中,预设门限值可以为系统缺省设置,也可以为外部输入给定,具体可根据对重建信号准确度的要求来进行设置及后续调节,此处不做限制。Wherein, the preset threshold value can be a default setting of the system, or can be given by an external input, which can be set and subsequently adjusted according to the requirements for the accuracy of the reconstructed signal, and there is no limitation here.
S2046、跳转执行下一次迭代更新过程。S2046. Jump to execute the next iterative update process.
若汉明距离hd大于预设门限值,则将迭代信号的下角标t加1,即将当前迭代过程中的第二迭代信号作为下一次迭代更新过程中的第一迭代信号参与计算,并重新执行S2041至S2045的操作。If the Hamming distance hd is greater than the preset threshold value, add 1 to the subscript t of the iterative signal, that is, the second iterative signal in the current iterative process as the first iteration signal during the next iterative update Participate in the calculation, and re-execute the operations from S2041 to S2045.
S2047、获取重建信号。S2047. Acquire a reconstruction signal.
具体的,若汉明距离小于或等于预设门限值,则对第二迭代信号进行单位归一化处理,得到原始稀疏信号的重建信号计算公式为:Specifically, if the Hamming distance is less than or equal to the preset threshold value, then for the second iteration signal Perform unit normalization to obtain the reconstructed signal of the original sparse signal The calculation formula is:
S205、将获取的重建信号作为最终结果输出。S205. Output the acquired reconstructed signal as a final result.
本发明的实施例提供一种信号处理方法,通过获取输入参数与编码端数据,在对过程参数与相关参数值初始化后,进行自适应稀疏度估计操作,并在获取了目标稀疏度后,根据目标稀疏度进行信号重建得到原始稀疏信号的重建信号。这样,解决了稀疏度先验信息受限导致信号重建性能降低的问题,降低或基本消除了重建信号产生的偏差,提高了单比特压缩感知技术的实用性与准确度。Embodiments of the present invention provide a signal processing method. By acquiring input parameters and encoding end data, after initializing process parameters and related parameter values, an adaptive sparsity estimation operation is performed, and after acquiring the target sparsity, according to The target sparsity is used to reconstruct the signal to obtain the reconstructed signal of the original sparse signal. In this way, the problem of signal reconstruction performance degradation due to limited sparsity prior information is solved, the deviation caused by reconstruction signal is reduced or basically eliminated, and the practicability and accuracy of single-bit compressed sensing technology are improved.
本发明实施例还提供一种信号处理装置00,如图5,该信号处理装置00包括:The embodiment of the present invention also provides a signal processing device 00, as shown in Figure 5, the signal processing device 00 includes:
参数获取单元001,用于获取输入参数与编码端数据,并对过程参数与第一迭代信号进行初始化,编码端数据包括编码端对原始稀疏信号进行低速采样处理及符号量化处理获取的测量信号;The parameter acquisition unit 001 is used to acquire input parameters and encoding end data, and initialize the process parameters and the first iteration signal. The encoding end data includes the measurement signal obtained by the encoding end performing low-speed sampling processing and symbol quantization processing on the original sparse signal;
稀疏度估计单元002,用于根据输入参数、编码端数据、初始化后的过程参数和第一迭代信号进行自适应稀疏度估计操作,获得目标稀疏度;The sparsity estimation unit 002 is used to perform an adaptive sparsity estimation operation according to the input parameters, the encoding end data, the initialized process parameters and the first iteration signal to obtain the target sparsity;
信号重建单元003,用于根据目标稀疏度,对测量信号进行信号重建得到原始稀疏信号的重建信号。The signal reconstruction unit 003 is configured to perform signal reconstruction on the measurement signal according to the target sparsity to obtain a reconstructed signal of the original sparse signal.
可选的,参数获取单元001具体用于:Optionally, the parameter acquisition unit 001 is specifically used for:
获取输入参数,输入参数包括最大稀疏度Km、判决门限λ;Obtain input parameters, including the maximum sparsity K m and the decision threshold λ;
接收编码端传递的编码端数据,编码端数据包括测量信号y、采样矩阵Φ;测量信号y为编码端结合采样矩阵Φ对原始稀疏信号θ进行低速采样处理获取采样信号x后,再对采样信号x进行符号量化处理获得的。Receive the encoding end data transmitted by the encoding end, the encoding end data includes the measurement signal y, the sampling matrix Φ; the measurement signal y is the encoding end combined with the sampling matrix Φ to perform low-speed sampling processing on the original sparse signal θ to obtain the sampling signal x, and then the sampling signal x is obtained by sign quantization.
可选的,过程参数包括目标稀疏度K的第一估计值pre_Est_S、第一使用值pre_Use_S;Optionally, the process parameters include the first estimated value pre_Est_S and the first use value pre_Use_S of the target sparsity K;
参数获取单元001还包括初始化单元004,如图6所示,具体用于:The parameter acquisition unit 001 also includes an initialization unit 004, as shown in Figure 6, specifically for:
将第一估计值pre_Est_s赋值为最大稀疏度Km,将第一使用值pre_Use_S赋值为最大稀疏度Km;Assign the first estimated value pre_Est_s as the maximum sparsity K m , and assign the first use value pre_Use_S as the maximum sparsity K m ;
将第一迭代信号的初始值赋值为零向量。The first iteration signal initial value of Assignment is a zero vector.
可选的,稀疏度估计单元002具体用于执行自适应稀疏度估计操作,其中,自适应稀疏度估计操作包括:Optionally, the sparsity estimation unit 002 is specifically configured to perform an adaptive sparsity estimation operation, where the adaptive sparsity estimation operation includes:
根据测量信号y与采样矩阵Φ,获取第一迭代信号经低速采样处理及符号量化处理后与测量信号y的差值,并对该差值左乘采样矩阵Φ的转置矩阵,得到第一迭代梯度 According to the measurement signal y and the sampling matrix Φ, the first iteration signal is obtained The difference between the measured signal y and the measured signal y after low-speed sampling processing and symbol quantization processing, and the difference is multiplied by the transpose matrix of the sampling matrix Φ to the left to obtain the first iteration gradient
通过梯度下降法获取过程信号at,过程信号at为第一迭代信号与第一迭代梯度的差值; Obtain the process signal at by gradient descent method, and the process signal at is the first iteration signal with the first iteration gradient the difference;
对第一迭代梯度是否趋近于零向量进行判断,若第一迭代梯度趋近于零向量,则将目标稀疏度K的第二使用值cur_Use_S赋值为第一估计值pre_Est_S;或,若第一迭代梯度不趋近于零向量,则将第二使用值cur_Use_S赋值为第一使用值pre_Use_S;For the first iteration gradient It is judged whether it is close to the zero vector, if the gradient of the first iteration approaching the zero vector, assign the second use value cur_Use_S of the target sparsity K to the first estimated value pre_Est_S; or, if the first iteration gradient If it does not approach the zero vector, assign the second use value cur_Use_S to the first use value pre_Use_S;
根据过程信号at与第二使用值cur_Use_S,通过门限函数Γ()获取第二迭代信号门限函数Γ()用于将过程信号at中元素幅值最大的第二使用值cur_Use_S个元素值保留,同时将过程信号at中除元素幅值最大的第二使用值cur_Use_S个元素之外的其他所有元素值置为零,并将门限函数Γ()的处理结果赋值给第二迭代信号 According to the process signal at and the second use value cur_Use_S , the second iteration signal is obtained through the threshold function Γ() The threshold function Γ() is used to keep the second use value cur_Use_S element values with the largest element amplitude in the process signal at, and at the same time save the process signal at except the second use value cur_Use_S elements with the largest element amplitude All other element values of are set to zero, and the processing result of the threshold function Γ() is assigned to the second iteration signal
对第二迭代信号进行单位归一化,得到归一化迭代信号统计归一化迭代信号中元素幅值绝对值超过判决门限λ的非零元素个数,并将该非零元素个数赋值给目标稀疏度K的第二估计值cur_Est_S;For the second iteration signal Perform unit normalization to obtain normalized iterative signals Statistically normalized iterative signal The number of non-zero elements whose absolute value of the element amplitude exceeds the decision threshold λ, and assign the number of non-zero elements to the second estimated value cur_Est_S of the target sparsity K;
若第一迭代梯度不趋近于零向量或第二估计值cur_Est_S与第一使用值pre_Use_S不相等,则将第一估计值pre_Est_S赋值为第二估计值cur_Est_S、第一使用值pre_Use_S赋值为第二使用值cur_Use_S、将迭代信号的下角标t加1,即将当前迭代过程中的第二迭代信号作为下一次迭代更新过程中的第一迭代信号参与计算,并重新执行自适应稀疏度估计操作;If the first iteration gradient If the vector does not approach zero or the second estimated value cur_Est_S is not equal to the first use value pre_Use_S, then the first estimate value pre_Est_S is assigned as the second estimate value cur_Est_S, the first use value pre_Use_S is assigned as the second use value cur_Use_S, and Add 1 to the subscript t of the iterative signal, that is, the second iterative signal in the current iterative process as the first iteration signal during the next iterative update Participate in the calculation and re-execute the adaptive sparsity estimation operation;
或,若第一迭代梯度趋近于零向量且第二估计值cur_Est_S与第一使用值pre_Use_S相等,则将目标稀疏度K赋值为第二估计值cur_Est_S,得到目标稀疏度K。Or, if the first iteration gradient is approaching the zero vector and the second estimated value cur_Est_S is equal to the first use value pre_Use_S, then the target sparsity K is assigned as the second estimated value cur_Est_S to obtain the target sparsity K.
可选的,信号重建单元003具体用于:Optionally, the signal reconstruction unit 003 is specifically used for:
根据编码端数据、第一迭代信号与目标稀疏度K,执行信号重建操作,信号重建操作具体包括:According to the encoder data, the first iteration signal With the target sparsity K, the signal reconstruction operation is performed, and the signal reconstruction operation specifically includes:
根据测量信号y与采样矩阵Φ,获取第一迭代信号经低速采样处理及符号量化处理后与测量信号y的差值,并对该差值左乘采样矩阵Φ的转置矩阵,得到第一迭代梯度 According to the measurement signal y and the sampling matrix Φ, the first iteration signal is obtained The difference between the measured signal y and the measured signal y after low-speed sampling processing and symbol quantization processing, and the difference is multiplied by the transpose matrix of the sampling matrix Φ to the left to obtain the first iteration gradient
通过梯度下降法获取过程信号at,过程信号at为第一迭代信号与第一迭代梯度的差值; Obtain the process signal at by gradient descent method, and the process signal at is the first iteration signal with the first iteration gradient the difference;
结合过程信号at与目标稀疏度K,通过门限函数Γ()获取第二迭代信号门限函数Γ()用于将过程信号at中元素幅值最大的目标稀疏度K个元素值保留,同时将过程信号at中除元素幅值最大的目标稀疏度K个元素之外的其他所有元素值置为零,并将门限函数Γ()的处理结果赋值给第二迭代信号 Combining the process signal at and the target sparsity K, the second iteration signal is obtained through the threshold function Γ() The threshold function Γ() is used to keep the target sparsity K elements with the largest element amplitude in the process signal a t , and at the same time save the other elements in the process signal a t except the target sparsity K elements with the largest element amplitude All element values are set to zero, and the processing result of the threshold function Γ() is assigned to the second iteration signal
根据测量信号y、采样矩阵Φ与第二迭代信号通过非零项统计函数获取汉明距离hd;非零项统计函数用于,获取第二迭代信号经低速采样处理及符号量化处理后与测量信号y的差值,并统计该差值中非零元素的个数,将非零元素的个数赋值给汉明距离hd;汉明距离hd为第二迭代信号经低速采样及符号量化处理后与测量信号y对应位置元素不相同的元素个数;According to the measurement signal y, the sampling matrix Φ and the second iteration signal Obtain the Hamming distance hd through the non-zero statistical function; the non-zero statistical function is used to obtain the second iteration signal The difference between the measured signal y and the measured signal y after low-speed sampling processing and symbol quantization processing, and the number of non-zero elements in the difference are counted, and the number of non-zero elements is assigned to the Hamming distance hd; the Hamming distance hd is the first Second iteration signal After low-speed sampling and symbol quantization processing, the number of elements is different from the corresponding position element of the measurement signal y;
若汉明距离hd大于预设门限值,则将迭代信号的下角标t加1,即将当前迭代过程中的第二迭代信号作为下一次迭代更新过程中的第一迭代信号参与计算,并重新执行信号重建操作;或,若汉明距离hd小于或等于预设门限值,则对第二迭代信号进行单位归一化处理,得到原始稀疏信号的重建信号 If the Hamming distance hd is greater than the preset threshold value, add 1 to the subscript t of the iterative signal, that is, the second iterative signal in the current iterative process as the first iteration signal during the next iterative update Participate in the calculation, and re-execute the signal reconstruction operation; or, if the Hamming distance hd is less than or equal to the preset threshold value, then for the second iteration signal Perform unit normalization to obtain the reconstructed signal of the original sparse signal
本发明的实施例提供一种信号处理装置,通过获取输入参数与编码端数据,在对过程参数与相关参数值初始化后,进行自适应稀疏度估计操作,并在获取了目标稀疏度后,根据目标稀疏度进行信号重建得到原始稀疏信号的重建信号。这样,解决了稀疏度先验信息受限导致信号重建性能降低的问题,降低或基本消除了重建信号产生的偏差,提高了单比特压缩感知技术的实用性与准确度。An embodiment of the present invention provides a signal processing device, which performs an adaptive sparsity estimation operation after initializing the process parameters and related parameter values by acquiring input parameters and encoding end data, and after acquiring the target sparsity, according to The target sparsity is used to reconstruct the signal to obtain the reconstructed signal of the original sparse signal. In this way, the problem of signal reconstruction performance degradation due to limited sparsity prior information is solved, the deviation caused by reconstruction signal is reduced or basically eliminated, and the practicability and accuracy of single-bit compressed sensing technology are improved.
本发明实施例还提供了一种信号处理装置01,如图7所示,该信号处理装置01包括:The embodiment of the present invention also provides a signal processing device 01, as shown in Figure 7, the signal processing device 01 includes:
总线011、以及连接到总线011的处理器012、存储器013和接口014,其中接口014用于与外部设备进行通信;A bus 011, and a processor 012 connected to the bus 011, a memory 013 and an interface 014, wherein the interface 014 is used to communicate with external devices;
该存储器013用于存储指令,该处理器012用于执行该指令用于获取输入参数与编码端数据,并对过程参数与第一迭代信号进行初始化,编码端数据包括编码端对原始稀疏信号进行低速采样处理及符号量化处理获取的测量信号;The memory 013 is used to store instructions, and the processor 012 is used to execute the instructions to obtain input parameters and encoding end data, and initialize process parameters and the first iteration signal, and the encoding end data includes the original sparse signal performed by the encoding end Measurement signals obtained by low-speed sampling processing and symbol quantization processing;
该处理器012执行该指令还用于根据输入参数、编码端数据、初始化后的过程参数和第一迭代信号进行自适应稀疏度估计操作,获得目标稀疏度;The processor 012 executes the instruction and is also used to perform an adaptive sparsity estimation operation according to the input parameters, the encoding end data, the initialized process parameters and the first iteration signal, so as to obtain the target sparsity;
该处理器012执行该指令还用于根据目标稀疏度,对测量信号进行信号重建得到原始稀疏信号的重建信号。The processor 012 executing the instruction is also used to perform signal reconstruction on the measurement signal according to the target sparsity to obtain a reconstructed signal of the original sparse signal.
在本发明实施例中,可选的,该处理器012执行该指令可以具体用于获取输入参数,输入参数包括最大稀疏度Km、判决门限λ;In the embodiment of the present invention, optionally, the execution of the instruction by the processor 012 may be specifically used to obtain input parameters, the input parameters include the maximum sparsity K m and the decision threshold λ;
接收编码端传递的编码端数据,编码端数据包括测量信号y、采样矩阵Φ;测量信号y为编码端结合采样矩阵Φ对原始稀疏信号θ经过低速采样处理获取采样信号x后,再对采样信号x进行符号量化处理后获得的。Receive the encoding end data transmitted by the encoding end, the encoding end data includes the measurement signal y, the sampling matrix Φ; the measurement signal y is the encoding end combined with the sampling matrix Φ to the original sparse signal θ after low-speed sampling processing to obtain the sampling signal x, and then the sampling signal x is obtained after sign quantization processing.
在本发明实施例中,可选的,过程参数包括目标稀疏度K的第一估计值pre_Est_S、第一使用值pre_Use_S,该处理器012执行该指令可以具体用于:In the embodiment of the present invention, optionally, the process parameters include the first estimated value pre_Est_S of the target sparsity K and the first use value pre_Use_S, and the processor 012 executing the instruction can be specifically used for:
将第一估计值pre_Est_S赋值为最大稀疏度Km,将第一使用值pre_Use_S赋值为最大稀疏度Km;Assign the first estimated value pre_Est_S as the maximum sparsity K m , and assign the first use value pre_Use_S as the maximum sparsity K m ;
将第一迭代信号的初始值赋值为零向量。The first iteration signal initial value of Assignment is a zero vector.
在本发明实施例中,可选的,该处理器012执行该指令可以具体用于执行自适应稀疏度估计操作,其中,自适应稀疏度估计操作包括:In this embodiment of the present invention, optionally, the processor 012 executing the instruction may be specifically used to perform an adaptive sparsity estimation operation, where the adaptive sparsity estimation operation includes:
根据测量信号y与采样矩阵Φ,获取第一迭代信号经低速采样处理及符号量化处理后与测量信号y的差值,并对该差值左乘采样矩阵Φ的转置矩阵,得到第一迭代梯度 According to the measurement signal y and the sampling matrix Φ, the first iteration signal is obtained The difference between the measured signal y and the measured signal y after low-speed sampling processing and symbol quantization processing, and the difference is multiplied by the transpose matrix of the sampling matrix Φ to the left to obtain the first iteration gradient
通过梯度下降法获取过程信号at,过程信号at为第一迭代信号与第一迭代梯度的差值; Obtain the process signal at by gradient descent method, and the process signal at is the first iteration signal with the first iteration gradient the difference;
对第一迭代梯度是否趋近于零向量进行判断,若第一迭代梯度趋近于零向量,则将目标稀疏度K的第二使用值cur_Use_S赋值为第一估计值pre_Est_S;或,若第一迭代梯度不趋近于零向量,则将第二使用值cur_Use_S赋值为第一使用值pre_Use_S;For the first iteration gradient It is judged whether it is close to the zero vector, if the gradient of the first iteration approaching the zero vector, assign the second use value cur_Use_S of the target sparsity K to the first estimated value pre_Est_S; or, if the first iteration gradient If it does not approach the zero vector, assign the second use value cur_Use_S to the first use value pre_Use_S;
根据过程信号at与第二使用值cur_Use_S,通过门限函数Γ()获取第二迭代信号门限函数Γ()用于将过程信号at中元素幅值最大的第二使用值cur_Use_S个元素值保留,同时将过程信号at中除元素幅值最大的第二使用值cur_Use_S个元素之外的其他所有元素值置为零,并将门限函数Γ()的处理结果赋值给第二迭代信号 According to the process signal at and the second use value cur_Use_S , the second iteration signal is obtained through the threshold function Γ() The threshold function Γ() is used to keep the second use value cur_Use_S element values with the largest element amplitude in the process signal at, and at the same time save the process signal at except the second use value cur_Use_S elements with the largest element amplitude All other element values of are set to zero, and the processing result of the threshold function Γ() is assigned to the second iteration signal
对第二迭代信号进行单位归一化,得到归一化迭代信号统计归一化迭代信号中元素幅值绝对值超过判决门限λ的非零元素个数,并将该非零元素个数赋值给目标稀疏度K的第二估计值cur_Est_S;For the second iteration signal Perform unit normalization to obtain normalized iterative signals Statistically normalized iterative signal The number of non-zero elements whose absolute value of the element amplitude exceeds the decision threshold λ, and assign the number of non-zero elements to the second estimated value cur_Est_S of the target sparsity K;
若第一迭代梯度不趋近于零向量或第二估计值cur_Est_S与第一使用值pre_Use_S不相等,则将第一估计值pre_Est_S赋值为第二估计值cur_Est_S、第一使用值pre_Use_S赋值为第二使用值cur_Use_S、将迭代信号的下角标t加1,即将当前迭代过程中的第二迭代信号作为下一次迭代更新过程中的第一迭代信号参与计算,并重新执行自适应稀疏度估计操作;If the first iteration gradient If the vector does not approach zero or the second estimated value cur_Est_S is not equal to the first use value pre_Use_S, then the first estimate value pre_Est_S is assigned as the second estimate value cur_Est_S, the first use value pre_Use_S is assigned as the second use value cur_Use_S, and Add 1 to the subscript t of the iterative signal, that is, the second iterative signal in the current iterative process as the first iteration signal during the next iterative update Participate in the calculation and re-execute the adaptive sparsity estimation operation;
或,若第一迭代梯度趋近于零向量且第二估计值cur_Est_S与第一使用值pre_Use_S相等,则将目标稀疏度K赋值为第二估计值cur_Est_S,得到目标稀疏度K。Or, if the first iteration gradient is approaching the zero vector and the second estimated value cur_Est_S is equal to the first use value pre_Use_S, then the target sparsity K is assigned as the second estimated value cur_Est_S to obtain the target sparsity K.
在本发明实施例中,可选的,该处理器012执行该指令可以具体用于根据编码端数据、第一迭代信号与目标稀疏度K,执行信号重建操作,信号重建操作具体包括:In this embodiment of the present invention, optionally, the processor 012 executing the instruction may be specifically used to With the target sparsity K, the signal reconstruction operation is performed, and the signal reconstruction operation specifically includes:
根据测量信号y与采样矩阵Φ,获取第一迭代信号经低速采样处理及符号量化处理后与测量信号y的差值,并对该差值左乘采样矩阵Φ的转置矩阵,得到第一迭代梯度 According to the measurement signal y and the sampling matrix Φ, the first iteration signal is obtained The difference between the measured signal y and the measured signal y after low-speed sampling processing and symbol quantization processing, and the difference is multiplied by the transpose matrix of the sampling matrix Φ to the left to obtain the first iteration gradient
通过梯度下降法获取过程信号at,过程信号at为第一迭代信号与第一迭代梯度的差值; Obtain the process signal at by gradient descent method, and the process signal at is the first iteration signal with the first iteration gradient the difference;
结合过程信号at与目标稀疏度K,通过门限函数Γ()获取第二迭代信号门限函数Γ()用于将过程信号at中元素幅值最大的目标稀疏度K个元素值保留,同时将过程信号at中除元素幅值最大的目标稀疏度K个元素之外的其他所有元素值置为零,并将门限函数Γ()的处理结果赋值给第二迭代信号 Combining the process signal at and the target sparsity K, the second iteration signal is obtained through the threshold function Γ() The threshold function Γ() is used to keep the target sparsity K elements with the largest element amplitude in the process signal a t , and at the same time save the other elements in the process signal a t except the target sparsity K elements with the largest element amplitude All element values are set to zero, and the processing result of the threshold function Γ() is assigned to the second iteration signal
根据测量信号y、采样矩阵Φ与第二迭代信号通过非零项统计函数获取汉明距离hd;非零项统计函数用于,获取第二迭代信号经低速采样处理及符号量化处理后与测量信号y的差值,并统计该差值中非零元素的个数,将非零元素的个数赋值给汉明距离hd;汉明距离hd为第二迭代信号经低速采样及符号量化处理后与测量信号y对应位置元素不相同的元素个数;According to the measurement signal y, the sampling matrix Φ and the second iteration signal Obtain the Hamming distance hd through the non-zero statistical function; the non-zero statistical function is used to obtain the second iteration signal The difference between the measured signal y and the measured signal y after low-speed sampling processing and symbol quantization processing, and the number of non-zero elements in the difference are counted, and the number of non-zero elements is assigned to the Hamming distance hd; the Hamming distance hd is the first Second iteration signal After low-speed sampling and symbol quantization processing, the number of elements is different from the corresponding position element of the measurement signal y;
若汉明距离hd大于预设门限值,则将迭代信号的下角标t加1,即将当前迭代过程中的第二迭代信号作为下一次迭代更新过程中的第一迭代信号参与计算,并重新执行信号重建操作;或,若汉明距离hd小于或等于预设门限值,则对第二迭代信号进行单位归一化处理,得到原始稀疏信号的重建信号 If the Hamming distance hd is greater than the preset threshold value, add 1 to the subscript t of the iterative signal, that is, the second iterative signal in the current iterative process as the first iteration signal during the next iterative update Participate in the calculation, and re-execute the signal reconstruction operation; or, if the Hamming distance hd is less than or equal to the preset threshold value, then for the second iteration signal Perform unit normalization to obtain the reconstructed signal of the original sparse signal
本发明的实施例提供一种信号处理装置,通过获取输入参数与编码端数据,在对过程参数与相关参数值初始化后,进行自适应稀疏度估计操作,并在获取了目标稀疏度后,根据目标稀疏度进行信号重建得到原始稀疏信号的重建信号。这样,解决了稀疏度先验信息受限导致信号重建性能降低的问题,降低或基本消除了重建信号产生的偏差,提高了单比特压缩感知技术的实用性与准确度。An embodiment of the present invention provides a signal processing device, which performs an adaptive sparsity estimation operation after initializing the process parameters and related parameter values by acquiring input parameters and encoding end data, and after acquiring the target sparsity, according to The target sparsity is used to reconstruct the signal to obtain the reconstructed signal of the original sparse signal. In this way, the problem of signal reconstruction performance degradation due to limited sparsity prior information is solved, the deviation caused by reconstruction signal is reduced or basically eliminated, and the practicability and accuracy of single-bit compressed sensing technology are improved.
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。The above is only a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Anyone skilled in the art can easily think of changes or substitutions within the technical scope disclosed in the present invention. Should be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.
Claims (12)
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2013/090961 WO2015100559A1 (en) | 2013-12-30 | 2013-12-30 | Signal processing method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105027519A CN105027519A (en) | 2015-11-04 |
CN105027519B true CN105027519B (en) | 2018-11-13 |
Family
ID=53492918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380033634.0A Active CN105027519B (en) | 2013-12-30 | 2013-12-30 | A kind of signal processing method and device |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN105027519B (en) |
WO (1) | WO2015100559A1 (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106059971A (en) * | 2016-07-07 | 2016-10-26 | 西北工业大学 | Sparse reconstruction based correlation detection method under signal correlation attenuation condition |
CN108614252B (en) * | 2018-07-02 | 2022-04-15 | 南京理工大学 | Single-bit compressed sensing radar target time delay estimation method based on parameterized sparse representation |
CN111523656B (en) * | 2019-02-03 | 2024-03-26 | 上海寒武纪信息科技有限公司 | Processing devices and methods |
CN110007284B (en) * | 2019-04-10 | 2023-01-31 | 南京航空航天大学 | A Dimensionality Reduction Method for Nonlinear Target Reconstruction Problem in Impulse System 1-bit Radar |
CN112906046B (en) * | 2021-01-27 | 2024-04-19 | 清华大学 | A model training method and device using single-bit compressed sensing technology |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5673210A (en) * | 1995-09-29 | 1997-09-30 | Lucent Technologies Inc. | Signal restoration using left-sided and right-sided autoregressive parameters |
CN101908889A (en) * | 2010-07-30 | 2010-12-08 | 哈尔滨工业大学 | A Compressive Sensing Reconstruction Method for Sparse Signals with Unknown Block Sparsity |
CN101908890A (en) * | 2010-07-30 | 2010-12-08 | 哈尔滨工业大学 | Blind reconstructing method of block sparse signal with unknown block size |
CN102045118A (en) * | 2010-10-22 | 2011-05-04 | 清华大学 | Compressed sensing reconstruction algorithm based on pseudo-inverse multiplication |
CN102497337A (en) * | 2011-12-11 | 2012-06-13 | 天津大学 | Compressed sensing wireless communication channel estimation method based on sparsity self-adapting |
CN102624399A (en) * | 2012-03-30 | 2012-08-01 | 北京邮电大学 | A Reconstruction Method for Compressed Sensing Signals |
CN102984711A (en) * | 2012-11-21 | 2013-03-20 | 北京邮电大学 | Multi-user collaborative spectrum sensing method based on single bit compression sensing technology |
CN103107815A (en) * | 2013-01-09 | 2013-05-15 | 华为技术有限公司 | Method and device for signal reconstruction |
CN103250352A (en) * | 2011-01-10 | 2013-08-14 | 上海贝尔股份有限公司 | Method and apparatus for measuring and recovering sparse signals |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8855190B2 (en) * | 2009-10-21 | 2014-10-07 | Futurewei Technologies, Inc. | Communication system with compressive sensing |
US8219360B2 (en) * | 2009-10-30 | 2012-07-10 | Mitsubishi Electric Research Laboratories, Inc. | Method for reconstructing sparse signals from sign measurements using greedy search |
CN103036573B (en) * | 2012-12-03 | 2015-10-28 | 南京工业大学 | compressed sensing signal acquisition method based on filtering |
-
2013
- 2013-12-30 CN CN201380033634.0A patent/CN105027519B/en active Active
- 2013-12-30 WO PCT/CN2013/090961 patent/WO2015100559A1/en active Application Filing
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5673210A (en) * | 1995-09-29 | 1997-09-30 | Lucent Technologies Inc. | Signal restoration using left-sided and right-sided autoregressive parameters |
CN101908889A (en) * | 2010-07-30 | 2010-12-08 | 哈尔滨工业大学 | A Compressive Sensing Reconstruction Method for Sparse Signals with Unknown Block Sparsity |
CN101908890A (en) * | 2010-07-30 | 2010-12-08 | 哈尔滨工业大学 | Blind reconstructing method of block sparse signal with unknown block size |
CN102045118A (en) * | 2010-10-22 | 2011-05-04 | 清华大学 | Compressed sensing reconstruction algorithm based on pseudo-inverse multiplication |
CN103250352A (en) * | 2011-01-10 | 2013-08-14 | 上海贝尔股份有限公司 | Method and apparatus for measuring and recovering sparse signals |
CN102497337A (en) * | 2011-12-11 | 2012-06-13 | 天津大学 | Compressed sensing wireless communication channel estimation method based on sparsity self-adapting |
CN102624399A (en) * | 2012-03-30 | 2012-08-01 | 北京邮电大学 | A Reconstruction Method for Compressed Sensing Signals |
CN102984711A (en) * | 2012-11-21 | 2013-03-20 | 北京邮电大学 | Multi-user collaborative spectrum sensing method based on single bit compression sensing technology |
CN103107815A (en) * | 2013-01-09 | 2013-05-15 | 华为技术有限公司 | Method and device for signal reconstruction |
Also Published As
Publication number | Publication date |
---|---|
CN105027519A (en) | 2015-11-04 |
WO2015100559A1 (en) | 2015-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105027519B (en) | A kind of signal processing method and device | |
CN102938649B (en) | Power quality data self-adapting reconstruction decompression method based on compressive sensing theory | |
CN101908890B (en) | Blind reconstructing method of block sparse signal with unknown block size | |
CN109255429B (en) | A Parameter Decompression Method for Sparse Neural Network Models | |
CN106570559A (en) | Data processing method and device based on neural network | |
CN103368578A (en) | Compressed-sensing-based signal sampling method for distributed wireless sensor network nodes | |
CN110175641A (en) | Image-recognizing method, device, equipment and storage medium | |
CN102253117A (en) | Acoustic signal collection method based on compressed sensing | |
CN110717958A (en) | Image reconstruction method, device, equipment and medium | |
CN104767535A (en) | A Low Complexity Block Sparse Signal Reconstruction Method | |
CN103036574A (en) | Self-check sparseness self-adaption matching pursuit arithmetic based on compressive sensing | |
Mangia et al. | Zeroing for HW-efficient compressed sensing architectures targeting data compression in wireless sensor networks | |
CN106841865B (en) | The single-bit of Short Duration Power Quality Disturbance signal samples and reconstructing method | |
CN110309877B (en) | Feature map data quantization method and device, electronic equipment and storage medium | |
CN110728732A (en) | An image reconstruction method, device, equipment and medium | |
CN106301384A (en) | A kind of signal reconfiguring method based on splits' positions perception | |
CN104901704B (en) | Has the body Sensor Network signal reconfiguring method of time-space correlative characteristics | |
CN105846826B (en) | Compressed sensing signal reconfiguring method based on approximate smooth L0 norm | |
Zamani et al. | On the use of compressive sensing (CS) exploiting block sparsity for neural spike recording | |
CN106374936B (en) | A kind of power grid real-time control method based on compression sensing technology | |
CN103997347B (en) | Fast estimation method for sparseness of signals during signal reconstruction | |
CN110233624A (en) | Data reconstruction method based on distributed compression perception | |
CN114049514A (en) | Image classification network compression method based on parameter reinitialization | |
CN115329971A (en) | Method and apparatus for eliminating amplitude damping noise, electronic device, and medium | |
CN104767533B (en) | A kind of signal reconstruction method and device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20191218 Address after: 314400 No.11, Weisan Road, Nongfa District, Chang'an Town, Haining City, Jiaxing City, Zhejiang Province Patentee after: Haining hi tech Zone Science and Innovation Center Co.,Ltd. Address before: 510000 unit 2414-2416, building, No. five, No. 371, Tianhe District, Guangdong, China Patentee before: GUANGDONG GAOHANG INTELLECTUAL PROPERTY OPERATION Co.,Ltd. Effective date of registration: 20191218 Address after: 510000 unit 2414-2416, building, No. five, No. 371, Tianhe District, Guangdong, China Patentee after: GUANGDONG GAOHANG INTELLECTUAL PROPERTY OPERATION Co.,Ltd. Address before: 518129 Bantian HUAWEI headquarters building, Longgang District, Shenzhen City, Guangdong Province, Guangdong, Shenzhen Patentee before: HUAWEI TECHNOLOGIES Co.,Ltd. |
|
TR01 | Transfer of patent right |