CN105471529A - Spectrum signal sensing method and device - Google Patents

Spectrum signal sensing method and device Download PDF

Info

Publication number
CN105471529A
CN105471529A CN201511021220.3A CN201511021220A CN105471529A CN 105471529 A CN105471529 A CN 105471529A CN 201511021220 A CN201511021220 A CN 201511021220A CN 105471529 A CN105471529 A CN 105471529A
Authority
CN
China
Prior art keywords
spectrum signal
iteration
kth
reconstructed
subsignal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201511021220.3A
Other languages
Chinese (zh)
Other versions
CN105471529B (en
Inventor
徐文波
董超
林家儒
张聪
王一凡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts and Telecommunications
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN201511021220.3A priority Critical patent/CN105471529B/en
Publication of CN105471529A publication Critical patent/CN105471529A/en
Application granted granted Critical
Publication of CN105471529B publication Critical patent/CN105471529B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention provides a spectrum signal sensing method and device. The method comprises the following steps: determining a weighted term v(pi) corresponding to the ith element of a spectrum signal vector according to the ith element in a reconstruction spectrum signal (the formula is described in the specification) obtained by (k-1)th iteration and the minimum distance of the elements in a support set of (the formula is described in the specification); ai expresses the ith column of a known measurement matrix, rk-1 expresses a residual error, adding the imkth column of the known measurement matrix in a measurement matrix reconstruction atom set obtained in the (k-1)th iteration to form a measurement matrix reconstruction atom set updated by k th iteration; determining the reconstruction spectrum signal (the formula is described in the specification) of the k th iteration according to the minimum value of (the formula is described in the specification); and wherein, Y expresses a measured value of the spectrum signal to be reconstructed, and (the formula is described in the specification) expresses the measurement matrix reconstruction atom set updated by the k th iteration. The spectrum signal sensing method and device provided by the invention are convenient to apply and can be used for improving the signal reconstruction performance.

Description

Spectrum signal cognitive method and device
Technical field
The present invention relates to wireless communication field, particularly relate to a kind of spectrum signal cognitive method and device.
Background technology
Along with the growth of radio communication service, available frequency band day is becoming tight, and the problem of frequency spectrum resource scarcity is day by day serious.The existing frequency usage policy in countries in the world adopts license granting mechanisms mostly, and the user secured permission the not round-the-clock licensed band that takies, in some band portion times, user does not use, and separately has some just occupied once in a while.How to improve the availability of frequency spectrum, in zones of different and different time sections, effectively utilize different idle channels, become the technical problem that people pay special attention to.
Cognition wireless network is a kind of new intellectual technology for improving wireless communication spectrum utilance.The Wireless Telecom Equipment with cognitive function can environment around perception, recycling distributes to authorized user, but under a certain specific moment and environment, do not have occupied frequency band, namely dynamically " frequency spectrum cavity-pocket " is recycled, to realize no matter when and where can ensureing the high efficiency that wireless frequency spectrum utilizes.Therefore, perceived spectral environment is the primary task of cognitive radio, and it embodies the most significant feature of cognitive radio: can perception analyze the frequency range of specific region, finds out the spectrum interposition being applicable to communication.
At present, the spectrum signal cognitive method of compressed sensing technology is adopted to obtain more concern.Compressive sensing theory is pointed out: if signal is sparse (or compressible), so just can to sample to signal far below the sample rate of Nyquist (Nyquist) speed, and with high probability Accurate Reconstruction original signal.In compressive sensing theory, can carry out under low sampling rate the sampling of signal, therefore the requirement of sample devices be reduced greatly.
Distributed compression perception utilizes correlation and cross correlation in signal to carry out combined reconstruction to multiple signal, takes full advantage of the temporal correlation of distributed sensor data, decreases transmitted data amount, further reduce communication overhead.Joint sparse model (jointsparsemodel, JSM) is the basis that distributed compression perception theory is set up, and whether contain common sparse part according in all signals, JSM can be divided into three kinds of models, JSM-1, JSM-2, JSM-3.
In JSM-1 model, related sparse spectrum signal is added by common ground and unique portion and forms, namely
s j=z c+z j,j∈{1,2,…,J}
Signal z crepresent all s jcommon ground, signal z jrepresent each s junique portion, and two parts are all sparse.JSM-1 model is highly suitable for the modeling to frequency spectrum.Common ground corresponds to all sensing nodes frequency spectrum that can receive, unique portion corresponding to the frequency spectrum only having certain sensing node to receive, the frequency spectrum namely only used near this sensing node.
The combined signal reconfiguration scheme overwhelming majority under existing JSM-1 model adopts l 1norm algorithm realization, this algorithm quality reconstruction is better, but its speed is slow, is difficult to practical application for solution large scale problem; Under JSM-1 model based on l 0the restructing algorithm speed of norm, but this kind of algorithm needs the prior informations such as known degree of rarefication, and these prior informations are difficult to obtain, and limit its practical application.
Summary of the invention
The embodiment of the present invention provides a kind of spectrum signal cognitive method and device, in order to solve the problem in frequency spectrum sensing method of the prior art, the requirement of signal element nonzero probability prior information being limited to its practical application.
The spectrum signal cognitive method that the embodiment of the present invention provides, comprising:
According to the reconstructed spectrum signal that kth-1 iteration obtains in i-th element and described support set in the minimum range of each element, determine the weighted term v (p that spectral signal vector i-th element is corresponding i); The support set of reconstructed spectrum signal represents the set of nonzero element position in described reconstructed spectrum signal; I and k be more than or equal to 1 integer;
Determine to make reach the value of the i of maximum, in kth time iteration, this value is designated as i mk; Wherein, a irepresent the i-th row of known calculation matrix, r k-1the residual error of the measured value of the described spectrum signal rebuild after representing the measured value of spectrum signal and kth-1 iteration;
By the i determined in described kth time iteration mkvalue joins in the support set of the reconstructed spectrum signal that kth-1 iteration obtains, and forms the support set after upgrading in kth time iteration, and by i-th of described known calculation matrix mkrow add the calculation matrix obtained in kth-1 iteration and rebuild in atom set, and the calculation matrix formed after kth time iteration renewal rebuilds atom set; The support set of described reconstructed spectrum signal represents the set of nonzero element position in described reconstructed spectrum signal;
According to the reconstructed spectrum signal of minimum value determination kth time iteration wherein, Y represents the measured value of spectrum signal to be reconstructed, represent that the calculation matrix after time iteration renewal of described kth rebuilds atom set, the measured value of the described spectrum signal rebuild after representing described kth time iteration.
In another embodiment, described basis the reconstructed spectrum signal of minimum value determination kth time iteration afterwards, also comprise:
According to the i that described kth time iteration obtains mk, described in determining i-th mksubsignal z belonging to individual element j; Wherein z j∈ { z c, z 1..., z j, according to upgrade described subsignal z jsub-support set wherein, n represents described subsignal z jlength, J represents the number of described frequency spectrum subsignal to be reconstructed, J and N be more than or equal to 1 integer; z crepresent the common ground of described J frequency spectrum subsignal to be reconstructed, z 1..., z jrepresent the unique portion of described J frequency spectrum subsignal to be reconstructed;
Determine described subsignal z jsub-support set element number whether reach described signal z jdegree of rarefication K j∈ { K c, K 1..., K j, wherein, degree of rarefication K c, K 1..., K jrepresent subsignal z c, z 1..., z jmiddle nonzero element number; If reached, then kth+1 iteration finds i m (k+1)time do not consider subsignal z jany element.
In another embodiment, the described i obtained according to described kth time iteration mk, described in determining i-th mksubsignal z belonging to individual element jalso comprise: described in determining i-th mkindividual element is at described subsignal z jin position i', for finding i in kth+1 iteration m (k+1)time do not consider i-th in each subsignal ' individual element.
In another embodiment, the described reconstructed spectrum signal obtained according to kth-1 iteration in i-th element and described support set in weighted term corresponding to minimum range determination spectral signal vector i-th element of each element comprise:
According to the reconstructed spectrum signal that described kth-1 iteration obtains in i-th element and the support set of reconstructed spectrum signal that obtains of described kth-1 iteration in the minimum range of each element, determine the nonzero probability of described reconstructed spectrum signal i-th element;
Weighted term v (the p of described reconstructed spectrum signal i-th element is determined according to the nonzero probability of i-th element in described reconstructed spectrum signal i).
In another embodiment, the described reconstructed spectrum signal obtained according to described kth-1 iteration in i-th element and the support set of reconstructed spectrum signal that obtains of described kth-1 iteration in the minimum range of each element, determine the nonzero probability of described reconstructed spectrum signal i-th element, comprising:
According to
p i=C·exp(-n i/α)
Determine the nonzero probability of described reconstructed spectrum signal i-th element;
Wherein, p irepresent the nonzero probability of i-th element, C represents a constant, is positive count, n irepresent the reconstructed spectrum signal that kth-1 iteration obtains the minimum value of the distance of each element in the support set of the reconstructed spectrum signal that i-th element and described kth-1 iteration obtain, α represents decay factor, control p ithe speed of decay.
In another embodiment, the described nonzero probability according to i-th element in described reconstructed spectrum signal determines the weighted term of described reconstructed spectrum signal i-th element, comprising:
According to
v ( p i ) = g 2 ( 2 ( K + 1 - k ) - 1 + 2 ( σ g ) 2 ) l o g p i 1 - p i
Determine the weighted term of described reconstructed spectrum signal i-th element;
Wherein, v (p i) representing the weighted term of described reconstructed spectrum signal i-th element, g is the average of nonzero element in the measured value of described spectrum signal, and σ is the standard deviation of noise, and k represents iterations, and K represents the joint sparse degree of spectrum signal to be reconstructed, wherein, K=K c+ K 1+ ... + K j.
In another embodiment, describedly determine described subsignal z jsub-support set element number whether reach described signal z jdegree of rarefication K j∈ { K c, K 1..., K j, if reached, then next iteration finds i m (k+1)time do not consider subsignal z jany element after, also comprise:
The residual error of the measured value rebuild after upgrading the measured value of described spectrum signal and kth time iteration and iterations k;
Judge whether described iterations k meets k>K, if met, then stop iteration.
In another embodiment, describedly judge whether described iterations k meets k>K, if met, then, after stopping iteration, also comprise:
By described reconstructed spectrum signal matrix be split as reconstruction subsignal obtain the reconstructed results of described frequency spectrum subsignal
The embodiment of the present invention also provides a kind of spectrum signal sensing device, comprising:
Determination module, for:
According to the reconstructed spectrum signal that kth-1 iteration obtains in i-th element and described support set in the minimum range of each element, determine the weighted term that spectral signal vector i-th element is corresponding; I and k be more than or equal to 1 integer;
Determine to make or reach the value of the i of maximum, in kth time iteration, this value is designated as i mk; Wherein, a irepresent the i-th row of known calculation matrix, r k-1the residual error of the measured value rebuild after representing the measured value of spectrum signal and kth-1 iteration, v (p i) represent the weighted term that spectral signal vector i-th element is corresponding; I and k be more than or equal to 1 integer;
Update module, for: by the i determined in described kth time iteration mkvalue joins in the support set of the reconstructed spectrum signal that kth-1 iteration obtains, form the support set after upgrading in kth time iteration, and the i-th row of calculation matrix described in described kth time iteration are added in the calculation matrix reconstruction atom set obtained in kth-1 iteration, the calculation matrix after formation kth time iteration upgrades rebuilds atom set; The support set of described reconstructed spectrum signal represents the set of nonzero element position in described reconstructed spectrum signal;
Described determination module, also for:
According to the reconstructed spectrum signal of minimum value determination kth time iteration wherein, Y represents the measured value of spectrum signal to be reconstructed, represent that the calculation matrix of kth time iteration rebuilds atom set; the measured value of the described spectrum signal rebuild after representing described kth time iteration.
The embodiment of the present invention also provides another spectrum signal sensing device, comprising: processor and memory, and the software program instructions that described processor is used for storing in execute store realizes the method described in said method embodiment.
The frequency spectrum sensing method that the embodiment of the present invention provides and device, the reconstructed spectrum signal obtained by kth-1 iteration in i-th element and described support set in the minimum range of each element, determine the position i of the nonzero element of kth time iteration intermediate frequency spectrum signal vector mk, and by i mkjoin in the support set of the reconstructed spectrum signal that kth-1 iteration obtains, form the support set after upgrading in kth time iteration, and the i-th row of calculation matrix described in described kth time iteration are added in the calculation matrix reconstruction atom set obtained in kth-1 iteration, the calculation matrix after formation kth time iteration upgrades rebuilds atom set; Basis again the reconstructed spectrum signal of minimum value determination kth time iteration wherein, Y represents the measured value of spectrum signal to be reconstructed, represent that the calculation matrix after time iteration renewal of described kth rebuilds atom set, the measured value of the described spectrum signal rebuild after representing described kth time iteration; The frequency spectrum sensing method that the embodiment of the present invention provides, utilizes the minimum range of spectrum signal element and support set element, determines the nonzero probability of spectrum signal element, thus obtains the position of the nonzero element in spectral signal vector, avoids based on l 0to the requirement of signal element nonzero probability prior information in the restructing algorithm of norm, can be used widely, and the performance of signal reconstruction can be improved.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, introduce doing one to the accompanying drawing used required in embodiment or description of the prior art simply below, apparently, accompanying drawing in the following describes is some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the schematic flow sheet of embodiment of the present invention spectrum signal cognitive method;
Fig. 2 is the structural representation of embodiment of the present invention spectrum signal sensing device;
Fig. 3 is the structural representation of another embodiment of the present invention spectrum signal sensing device;
Fig. 4 is the structural representation of yet another embodiment of the invention spectrum signal sensing device.
Embodiment
For making the object of the embodiment of the present invention, technical scheme and advantage clearly, below in conjunction with the accompanying drawing in the embodiment of the present invention, technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.
The spectrum signal cognitive method that the embodiment of the present invention provides and device are for perception in cognition wireless network and analyze the frequency range of specific region, find out the spectrum interposition being applicable to communication, so that Wireless Telecom Equipment can utilize distribute to authorized user, but under a certain particular moment and environment, there is no occupied frequency band, improve the utilance of wireless frequency spectrum.
The spectrum signal cognitive method that the embodiment of the present invention provides and device are in cognition wireless network, utilize J sensing node and 1 central fusion center to carry out collaborative spectrum sensing.For jth (1≤j≤J) individual sensing node, to be the measured value of the spectrum signal of M be its length received known calculation matrix is then y jjs j, wherein the spectrum signal of to be length be N.According to JSM-1 model, any 1 spectrum signal all can be expressed as s j=z c+ z j, use K cand K jrepresent z respectively cand z jdegree of rarefication, then joint sparse degree is K=K c+ ∑ K j.Because joint sparse degree K is less than the degree of rarefication sum ∑ (K of each node institute perceptual signal c+ K j)=NK c+ ∑ K j, so can obviously improve signaling protein14-3-3 performance by integrated restoration.The embodiment of the present invention accurately judges the support set of spectrum signal by frequency spectrum sensing method spectrum signal being carried out to integrated restoration, i.e. nonzero element location sets.
Particularly, the measured value that J sensing node receives by fusion center merges, and obtains and calculation matrix corresponding for J spectrum signal is merged by mode below:
Fusion center is when known Y and Φ, and the spectrum signal reconstructed is to be J spectrum signal to be reconstructed to recombinate the result obtained by common ground and unique portion described Z, meets Y=Φ Z; The support set of described spectrum signal Z by the subsignal z of Z c, z 1..., z jsupport set get union and obtain.Wherein, represent described subsignal z respectively c, z 1..., z jsupport set.
Fig. 1 is the schematic flow sheet of embodiment of the present invention spectrum signal cognitive method.Refer to Fig. 1, the spectrum signal cognitive method that the embodiment of the present invention provides comprises:
S101: the reconstructed spectrum signal obtained according to kth-1 iteration in i-th element and described support set in weighted term v (p corresponding to minimum range determination spectral signal vector i-th element of each element i); The support set of reconstructed spectrum signal represents the set of nonzero element position in described reconstructed spectrum signal; I and k be more than or equal to 1 integer;
Particularly, the described reconstructed spectrum signal obtained according to kth-1 iteration in i-th element and described support set in weighted term corresponding to minimum range determination spectral signal vector i-th element of each element comprise:
According to the reconstructed spectrum signal that described kth-1 iteration obtains in i-th element and the support set of reconstructed spectrum signal that obtains of described kth-1 iteration in the minimum range of each element, determine the nonzero probability of described reconstructed spectrum signal i-th element;
Weighted term v (the p of described reconstructed spectrum signal i-th element is determined according to the nonzero probability of i-th element in described reconstructed spectrum signal i).
The described reconstructed spectrum signal obtained according to described kth-1 iteration in i-th element and the support set of reconstructed spectrum signal that obtains of described kth-1 iteration in the minimum range of each element, determine the nonzero probability of described reconstructed spectrum signal i-th element, comprising:
According to
p i=C·exp(-n i/α)
Determine the nonzero probability of described reconstructed spectrum signal i-th element;
Wherein, p irepresent the nonzero probability of i-th element, C represents a constant, is positive count, n irepresent the reconstructed spectrum signal that kth-1 iteration obtains the minimum value of the distance of each element in the support set of the reconstructed spectrum signal that i-th element and described kth-1 iteration obtain, α represents decay factor, control p ithe speed of decay.
The described nonzero probability according to i-th element in described reconstructed spectrum signal determines the weighted term of described reconstructed spectrum signal i-th element, comprising:
According to
v ( p i ) = g 2 ( 2 ( K + 1 - k ) - 1 + 2 ( σ g ) 2 ) l o g p i 1 - p i
Determine the weighted term of described reconstructed spectrum signal i-th element;
Wherein, v (p i) representing the weighted term of described reconstructed spectrum signal i-th element, g is the average of nonzero element in the measured value of described spectrum signal, and σ is the standard deviation of noise, and k represents iterations, and K represents the joint sparse degree of spectrum signal to be reconstructed, wherein, K=K c+ K 1+ ... + K j.
S102: determine to make reach the value of the i of maximum, in kth time iteration, this value is designated as i mk; Wherein, a irepresent the i-th row of known calculation matrix, r k-1the residual error of the measured value of the described spectrum signal rebuild after representing the measured value of spectrum signal and kth-1 iteration;
S103: by the i determined in described kth time iteration mkvalue joins in the support set of the reconstructed spectrum signal that kth-1 iteration obtains, and forms the support set after upgrading in kth time iteration, and by i-th of described known calculation matrix mkrow add the calculation matrix obtained in kth-1 iteration and rebuild in atom set, and the calculation matrix formed after kth time iteration renewal rebuilds atom set; The support set of described reconstructed spectrum signal represents the set of nonzero element position in described reconstructed spectrum signal;
S104: according to the reconstructed spectrum signal of minimum value determination kth time iteration wherein, Y represents the measured value of spectrum signal to be reconstructed, represent that the calculation matrix after time iteration renewal of described kth rebuilds atom set, the measured value of the described spectrum signal rebuild after representing described kth time iteration;
Particularly, before execution step S101, the sensing node in cognition wireless network first carries out this locality to spectrum signal and detects, and obtains the measured value y of spectrum signal j.A described J sensing node will detect the described measured value y obtained jsend to described fusion center, the measured value y that a described J sensing node receives by described fusion center jmerge, obtain by calculation matrix Φ corresponding for J spectrum signal jmerge into Φ, can above-mentioned steps S101-S104 be performed.
The frequency spectrum sensing method that the embodiment of the present invention provides, the reconstructed spectrum signal obtained by kth-1 iteration in i-th element and described support set in the minimum range of each element, determine the position i of the nonzero element of kth time iteration intermediate frequency spectrum signal vector mk, and by i mkjoin in the support set of the reconstructed spectrum signal that kth-1 iteration obtains, form the support set after upgrading in kth time iteration, and the i-th row of calculation matrix described in described kth time iteration are added in the calculation matrix reconstruction atom set obtained in kth-1 iteration, the calculation matrix after formation kth time iteration upgrades rebuilds atom set; Basis again the reconstructed spectrum signal of minimum value determination kth time iteration wherein, Y represents the measured value of spectrum signal to be reconstructed, represent that the calculation matrix after time iteration renewal of described kth rebuilds atom set, the measured value of the described spectrum signal rebuild after representing described kth time iteration; The frequency spectrum sensing method that the embodiment of the present invention provides, utilizes the minimum range of signal element and support set element, determines the nonzero probability of spectrum signal element, thus obtains the position of the nonzero element in spectral signal vector, avoids based on l 0to the requirement of signal element nonzero probability prior information in the restructing algorithm of norm, and the performance of signal reconstruction can be improved.
Further, in order to ensure that the number of the nonzero element in described reconstructed spectrum signal is no more than the sparse angle value of setting, described basis the reconstructed spectrum signal of minimum value determination kth time iteration afterwards, also comprise:
According to the i that described kth time iteration obtains mk, described in determining i-th mksubsignal z belonging to individual element j; Wherein z j∈ { z c, z 1..., z j, according to upgrade described subsignal z jsub-support set wherein, n represents described subsignal z jlength, J represents the number of described frequency spectrum subsignal to be reconstructed, J and N be more than or equal to 1 integer; z crepresent the common ground of described J frequency spectrum subsignal to be reconstructed, z 1..., z jrepresent the unique portion of described J frequency spectrum subsignal to be reconstructed;
Determine described frequency spectrum subsignal z jsub-support set element number whether reach described signal z jdegree of rarefication K j∈ { K c, K 1..., K j, wherein, degree of rarefication K c, K 1..., K jrepresent subsignal z c, z 1..., z jmiddle nonzero element number; If reached, then kth+1 iteration finds i m (k+1)time do not consider subsignal z jany element position.
Further, in order to improve the efficiency that described frequency spectrum sensing method is rebuild spectrum signal, the described i obtained according to described kth time iteration mk, described in determining i-th mksubsignal z belonging to individual element jalso comprise: described in determining i-th kindividual element is at described subsignal z jin position i', for finding i in kth+1 iteration m (k+1)time do not consider i-th in each subsignal ' individual element.
Further, in order to ensure to continue before described fusion center reaches stopping criterion for iteration the reconstruction that iteration realizes described spectrum signal, and stopping iteration completing after spectrum signal is rebuild, describedly determining described subsignal z jsub-support set element number whether reach described signal z jdegree of rarefication K j∈ { K c, K 1..., K j, if reached, then next iteration finds i m (k+1)time do not consider subsignal z jany element position after, also comprise:
The residual error of the measured value rebuild after upgrading the measured value of described spectrum signal and kth time iteration and iterations k;
Judge whether described iterations k meets k>K, if met, then stop iteration.
Further, describedly judge whether described iterations k meets k>K, if met, then, after stopping iteration, also comprise:
By described reconstructed spectrum signal matrix be split as reconstruction subsignal obtain the reconstructed results of described frequency spectrum subsignal
The embodiment of the present invention also provides a kind of spectrum signal sensing device, and Fig. 2 is the structural representation of embodiment of the present invention spectrum signal sensing device.Refer to Fig. 2, described spectrum signal sensing device at least comprises: determination module 210 and update module 220.
Described determination module 210, for:
According to the reconstructed spectrum signal that kth-1 iteration obtains in i-th element and described support set in the minimum range of each element, determine the weighted term that spectral signal vector i-th element is corresponding; I and k be more than or equal to 1 integer;
Determine to make or reach the value of the i of maximum, in kth time iteration, this value is designated as i mk; Wherein, a irepresent the i-th row of known calculation matrix, r k-1the residual error of the measured value rebuild after representing the measured value of spectrum signal and kth-1 iteration, v (p i) represent the weighted term that spectral signal vector i-th element is corresponding; I and k be more than or equal to 1 integer;
Described update module 220, for: by the i determined in described kth time iteration mkvalue joins in the support set of the reconstructed spectrum signal that kth-1 iteration obtains, form the support set after upgrading in kth time iteration, and the i-th row of calculation matrix described in described kth time iteration are added in the calculation matrix reconstruction atom set obtained in kth-1 iteration, the calculation matrix after formation kth time iteration upgrades rebuilds atom set; The support set of described reconstructed spectrum signal represents the set of nonzero element position in described reconstructed spectrum signal;
Described determination module 210, also for:
According to the reconstructed spectrum signal of minimum value determination kth time iteration wherein, Y represents the measured value of spectrum signal to be reconstructed, represent that the calculation matrix of kth time iteration rebuilds atom set; the measured value of the described spectrum signal rebuild after representing described kth time iteration.
The frequency spectrum sensing device that the embodiment of the present invention provides, the reconstructed spectrum signal obtained according to kth-1 iteration by determination module in i-th element and described support set in the minimum range of each element, determine the position i of the nonzero element of kth time iteration intermediate frequency spectrum signal vector mk, by described update module by i mkjoin in the support set of the reconstructed spectrum signal that kth-1 iteration obtains, form the support set after upgrading in kth time iteration, and the i-th row of calculation matrix described in described kth time iteration are added in the calculation matrix reconstruction atom set obtained in kth-1 iteration, the calculation matrix after formation kth time iteration upgrades rebuilds atom set; Again by described determination module according to minimum value can determine the reconstructed spectrum signal of kth time iteration wherein, Y represents the measured value of spectrum signal to be reconstructed, represent that the calculation matrix after time iteration renewal of described kth rebuilds atom set, the measured value of the described spectrum signal rebuild after representing described kth time iteration; The frequency spectrum sensing device that the embodiment of the present invention provides, utilizes the minimum range of signal element and support set element, determines the nonzero probability of spectrum signal element, thus obtains the position of the nonzero element in spectral signal vector, avoids based on l 0to the requirement of signal element nonzero probability prior information in the restructing algorithm of norm, and the performance of signal reconstruction can be improved.
Fig. 3 is the structural representation of another embodiment of the present invention spectrum signal sensing device.Refer to Fig. 3, the spectrum signal sensing device that another embodiment of the present invention provides comprises: processor 310 and memory 320, and described processor 310 realizes the method described in preceding method embodiment for the software program instructions stored in execute store 320.Its specific implementation principle and technique effect identical with embodiment of the method, do not repeat them here.
Fig. 4 is the structural representation of yet another embodiment of the invention spectrum signal sensing device.Refer to Fig. 4, the spectrum signal sensing device that yet another embodiment of the invention provides comprises: sensing node 410 and fusion center 420.
Described sensing node 410 is for the measured value of received spectrum signal
The measured value y of described fusion center 420 for sensing node is received jmerge, obtain spectrum signal is pressed common ground and unique portion restructuring, be expressed as by calculation matrix Φ corresponding for J spectrum signal jmerge into Φ, and according to the method in preceding method embodiment, described spectrum signal is rebuild.Its specific implementation principle and technique effect identical with embodiment of the method, do not repeat them here.
One of ordinary skill in the art will appreciate that: all or part of step realizing above-mentioned each embodiment of the method can have been come by the hardware that program command is relevant.Aforesaid program can be stored in a computer read/write memory medium.This program, when performing, performs the step comprising above-mentioned each embodiment of the method; And aforesaid storage medium comprises: ROM, RAM, magnetic disc or CD etc. various can be program code stored medium.
Last it is noted that above each embodiment is only in order to illustrate technical scheme of the present invention, be not intended to limit; Although with reference to foregoing embodiments to invention has been detailed description, those of ordinary skill in the art is to be understood that: it still can be modified to the technical scheme described in foregoing embodiments, or carries out equivalent replacement to wherein some or all of technical characteristic; And these amendments or replacement, do not make the essence of appropriate technical solution depart from the scope of various embodiments of the present invention technical scheme.

Claims (9)

1. a spectrum signal cognitive method, is characterized in that, comprising:
According to the reconstructed spectrum signal that kth-1 iteration obtains in i-th element and described support set in the minimum range of each element, determine the weighted term v (p that spectral signal vector i-th element is corresponding i); The support set of reconstructed spectrum signal represents the set of nonzero element position in described reconstructed spectrum signal; I and k be more than or equal to 1 integer;
Determine to make reach the value of the i of maximum, in kth time iteration, this value is designated as i mk; Wherein, a irepresent the i-th row of known calculation matrix, r k-1the residual error of the measured value of the described spectrum signal rebuild after representing the measured value of spectrum signal and kth-1 iteration;
By the i determined in described kth time iteration mkvalue joins in the support set of the reconstructed spectrum signal that kth-1 iteration obtains, and forms the support set after upgrading in kth time iteration, and by i-th of described known calculation matrix mkrow add the calculation matrix obtained in kth-1 iteration and rebuild in atom set, and the calculation matrix formed after kth time iteration renewal rebuilds atom set;
According to the reconstructed spectrum signal of minimum value determination kth time iteration wherein, Y represents the measured value of spectrum signal to be reconstructed, represent that the calculation matrix after time iteration renewal of described kth rebuilds atom set, the measured value of the described spectrum signal rebuild after representing described kth time iteration.
2. method according to claim 1, is characterized in that, described basis the reconstructed spectrum signal of minimum value determination kth time iteration afterwards, also comprise:
According to the i that described kth time iteration obtains mk, described in determining i-th mksubsignal z belonging to individual element j; Wherein z j∈ { z c, z 1..., z j, according to upgrade described subsignal z jsub-support set wherein, n represents described subsignal z jlength, J represents the number of described frequency spectrum subsignal to be reconstructed, J and N be more than or equal to 1 integer; z crepresent the common ground of described J frequency spectrum subsignal to be reconstructed, z 1..., z jrepresent the unique portion of described J frequency spectrum subsignal to be reconstructed;
Determine described frequency spectrum subsignal z jsub-support set element number whether reach described signal z jdegree of rarefication K j∈ { K c, K 1..., K j, wherein, degree of rarefication K c, K 1..., K jrepresent subsignal z c, z 1..., z jmiddle nonzero element number; If reached, then kth+1 iteration finds i m (k+1)time do not consider subsignal z jany element.
3. method according to claim 2, is characterized in that, the described i obtained according to described kth time iteration mk, described in determining i-th mksubsignal z belonging to individual element jalso comprise: described in determining i-th mkindividual element is at described subsignal z jin position i', for finding i in kth+1 iteration m (k+1)time do not consider i-th in each subsignal ' individual element.
4. according to the method in claim 2 or 3, it is characterized in that, the described reconstructed spectrum signal obtained according to kth-1 iteration in i-th element and described support set in weighted term corresponding to minimum range determination spectral signal vector i-th element of each element, comprising:
According to the reconstructed spectrum signal that described kth-1 iteration obtains in i-th element and the support set of reconstructed spectrum signal that obtains of described kth-1 iteration in the minimum range of each element, determine the nonzero probability of described reconstructed spectrum signal i-th element;
Weighted term v (the p of described reconstructed spectrum signal i-th element is determined according to the nonzero probability of i-th element in described reconstructed spectrum signal i).
5. method according to claim 4, is characterized in that, the described reconstructed spectrum signal obtained according to described kth-1 iteration in i-th element and the support set of reconstructed spectrum signal that obtains of described kth-1 iteration in the minimum range of each element, determine the nonzero probability of described reconstructed spectrum signal i-th element, comprising:
According to
p i=C·exp(-n i/α)
Determine the nonzero probability of described reconstructed spectrum signal i-th element;
Wherein, p irepresent the nonzero probability of i-th element, C represents a constant, is positive count, n irepresent the reconstructed spectrum signal that kth-1 iteration obtains the minimum value of the distance of each element in the support set of the reconstructed spectrum signal that i-th element and described kth-1 iteration obtain, α represents decay factor, control p ithe speed of decay.
6. method according to claim 4, is characterized in that, the described nonzero probability according to i-th element in described reconstructed spectrum signal determines the weighted term of described reconstructed spectrum signal i-th element, comprising:
According to
v ( p i ) = g 2 ( 2 ( K + 1 - k ) - 1 + 2 ( σ g ) 2 ) l o g p i 1 - p i
Determine the weighted term of described reconstructed spectrum signal i-th element;
Wherein, v (p i) representing the weighted term of described reconstructed spectrum signal i-th element, g is the average of nonzero element in the measured value of described spectrum signal, and σ is the standard deviation of noise, and k represents iterations, and K represents the joint sparse degree of spectrum signal to be reconstructed, wherein, K=K c+ K 1+ ... + K j.
7. method according to claim 6, is characterized in that, describedly determines described subsignal z jsub-support set element number whether reach described signal z jdegree of rarefication K j∈ { K c, K 1..., K j, if reached, then next iteration finds i m (k+1)time do not consider subsignal z jany element after, also comprise:
The residual error of the measured value rebuild after upgrading the measured value of described spectrum signal and kth time iteration and iterations k;
Judge whether described iterations k meets k>K, if met, then stop iteration.
8. method according to claim 7, is characterized in that, describedly judges whether described iterations k meets k>K, if met, then, after stopping iteration, also comprises:
By described reconstructed spectrum signal matrix be split as reconstruction subsignal obtain the reconstructed results of described frequency spectrum subsignal
9. a spectrum signal sensing device, is characterized in that, comprising:
Determination module, for:
According to the reconstructed spectrum signal that kth-1 iteration obtains in i-th element and described support set in the minimum range of each element, determine the weighted term that spectral signal vector i-th element is corresponding; I and k be more than or equal to 1 integer;
Determine to make or reach the value of the i of maximum, in kth time iteration, this value is designated as i mk; Wherein, a irepresent the i-th row of known calculation matrix, r k-1the residual error of the measured value rebuild after representing the measured value of spectrum signal and kth-1 iteration, v (p i) represent the weighted term that spectral signal vector i-th element is corresponding; I and k be more than or equal to 1 integer;
Update module, for: by the i determined in described kth time iteration mkvalue joins in the support set of the reconstructed spectrum signal that kth-1 iteration obtains, form the support set after upgrading in kth time iteration, and the i-th row of calculation matrix described in described kth time iteration are added in the calculation matrix reconstruction atom set obtained in kth-1 iteration, the calculation matrix after formation kth time iteration upgrades rebuilds atom set; The support set of described reconstructed spectrum signal represents the set of nonzero element position in described reconstructed spectrum signal;
Described determination module, also for:
According to the reconstructed spectrum signal of minimum value determination kth time iteration wherein, Y represents the measured value of spectrum signal to be reconstructed, represent that the calculation matrix of kth time iteration rebuilds atom set; the measured value of the described spectrum signal rebuild after representing described kth time iteration.
CN201511021220.3A 2015-12-30 2015-12-30 Spectrum signal cognitive method and device Active CN105471529B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511021220.3A CN105471529B (en) 2015-12-30 2015-12-30 Spectrum signal cognitive method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511021220.3A CN105471529B (en) 2015-12-30 2015-12-30 Spectrum signal cognitive method and device

Publications (2)

Publication Number Publication Date
CN105471529A true CN105471529A (en) 2016-04-06
CN105471529B CN105471529B (en) 2018-02-02

Family

ID=55608881

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511021220.3A Active CN105471529B (en) 2015-12-30 2015-12-30 Spectrum signal cognitive method and device

Country Status (1)

Country Link
CN (1) CN105471529B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108566254A (en) * 2018-03-27 2018-09-21 广东工业大学 A kind of frequency spectrum sensing method and device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101764652A (en) * 2010-01-18 2010-06-30 哈尔滨工业大学 Signal detection method having compression perception process based on orthogonal matching pursuit
CN102984711A (en) * 2012-11-21 2013-03-20 北京邮电大学 Multi-user collaborative spectrum sensing method based on single bit compression sensing technology
CN103139828A (en) * 2012-12-29 2013-06-05 东北大学 Broadband spectrum sensing device and method
CN103532645A (en) * 2013-10-10 2014-01-22 南京邮电大学 Compressive spectrum sensing method for observing matrix optimization
CN103873170A (en) * 2014-03-26 2014-06-18 哈尔滨工业大学 Compressed sensing spectrum detecting method under blind sparse condition
CN103929256A (en) * 2014-04-25 2014-07-16 哈尔滨工业大学 Multi-frame compression perception signal spectrum detection method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101764652A (en) * 2010-01-18 2010-06-30 哈尔滨工业大学 Signal detection method having compression perception process based on orthogonal matching pursuit
CN102984711A (en) * 2012-11-21 2013-03-20 北京邮电大学 Multi-user collaborative spectrum sensing method based on single bit compression sensing technology
CN103139828A (en) * 2012-12-29 2013-06-05 东北大学 Broadband spectrum sensing device and method
CN103532645A (en) * 2013-10-10 2014-01-22 南京邮电大学 Compressive spectrum sensing method for observing matrix optimization
CN103873170A (en) * 2014-03-26 2014-06-18 哈尔滨工业大学 Compressed sensing spectrum detecting method under blind sparse condition
CN103929256A (en) * 2014-04-25 2014-07-16 哈尔滨工业大学 Multi-frame compression perception signal spectrum detection method

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108566254A (en) * 2018-03-27 2018-09-21 广东工业大学 A kind of frequency spectrum sensing method and device

Also Published As

Publication number Publication date
CN105471529B (en) 2018-02-02

Similar Documents

Publication Publication Date Title
Mendlik et al. Selecting climate simulations for impact studies based on multivariate patterns of climate change
US9665669B2 (en) Techniques for estimating compound probability distribution by simulating large empirical samples with scalable parallel and distributed processing
Rodrigues et al. Constraints on galaxy formation models from the galaxy stellar mass function and its evolution
CN105893407A (en) Individual user portraying method and system
CN104519494A (en) Device, method, and system for radio-frequency spectrum management, and secondary user equipment
WO2018161729A1 (en) User path recovery method and device
US10659919B1 (en) System and method for automated commissioning of one or more networks of electronic devices
Risser et al. A probabilistic gridded product for daily precipitation extremes over the United States
CN104506378A (en) Data flow prediction device and method
CN112100841B (en) Method and device for predicting attitude of shield tunneling machine, terminal equipment and storage medium
CN105049141A (en) Analysis method for interference relation among signals based on multidimensional interference space model
Malakar et al. Optimal scheduling of in-situ analysis for large-scale scientific simulations
US20170344284A1 (en) Method and apparatus for storage device latency/bandwidth self monitoring
CA3117872A1 (en) Clustering techniques for machine learning models
CN104618045A (en) Collected data-based wireless channel transmission model establishing method and system
CN104952039A (en) Distributed compressive sensing reconstruction method for images
CN105139282A (en) Power grid index data processing method, device and calculation device
CN105471529A (en) Spectrum signal sensing method and device
Ma et al. Ordinal optimization-based performance model estimation method for HDFS
CN114817846A (en) Dynamic OD matrix estimation method and device under prior time characteristic constraint
CN113342418B (en) Distributed machine learning task unloading method based on block chain
CN109525339A (en) Frequency spectrum sensing method, device, equipment and the storage medium of cognitive radio
CN107124727B (en) PCI optimization method and device
CN102930458A (en) Method and system for electricity price evaluation simulation
CN105071877A (en) Broadband distributed Bayes compression spectrum sensing method

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