CN103428704A - Method and device for sensing frequency spectra - Google Patents

Method and device for sensing frequency spectra Download PDF

Info

Publication number
CN103428704A
CN103428704A CN2013103261174A CN201310326117A CN103428704A CN 103428704 A CN103428704 A CN 103428704A CN 2013103261174 A CN2013103261174 A CN 2013103261174A CN 201310326117 A CN201310326117 A CN 201310326117A CN 103428704 A CN103428704 A CN 103428704A
Authority
CN
China
Prior art keywords
channel
seizure condition
channels
perception
historical
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
CN2013103261174A
Other languages
Chinese (zh)
Other versions
CN103428704B (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.)
Guilin Ceke Communication Equipment Co ltd
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 CN201310326117.4A priority Critical patent/CN103428704B/en
Priority to PCT/CN2013/083253 priority patent/WO2015014005A1/en
Publication of CN103428704A publication Critical patent/CN103428704A/en
Application granted granted Critical
Publication of CN103428704B publication Critical patent/CN103428704B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W16/00Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/02Resource partitioning among network components, e.g. reuse partitioning
    • H04W16/10Dynamic resource partitioning

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Telephonic Communication Services (AREA)

Abstract

The invention discloses a method and a device for sensing frequency spectra, and relates to the technical field of cognitive radio. The method includes clustering all channels in the same business to obtain N channel sets and a representative channel of each channel set; determining a channel set of a channel to be sensed and a representative channel C<n> corresponding to the channel to be sensed, and computing a current occupancy state of the channel to be sensed according to a historical occupancy state of the channel to be sensed and a current occupancy state of the representative channel C<n>. The n is larger than or equal to 1 and is smaller than or equal to N. The method and the device have the advantages that only the current occupancy state of the certain representative channel needs to be detected, so that current occupancy states of all the channels in the corresponding channel sets can be computed and predicted by the method; the channel occupancy state prediction accuracy can be improved by the aid of occupancy dependency among the various channels in the same business.

Description

A kind of frequency spectrum sensing method and device
Technical field
The present invention relates to the cognitive radio technology field, particularly a kind of frequency spectrum sensing method and device.
Background technology
Along with the rapid growth of radio communication service amount, frequency spectrum problem in short supply is more and more serious.Frequency spectrum detecting result by the research countries in the world is found, authorizes the utilance of frequency range generally lower.Yet fixing spectrum allocation may policy makes the mandate frequency range of poor efficiency can't be used by unauthorized user at one's leisure, causes the great wasting of resources.The solution that appears as this problem of cognitive radio technology provides good approach and scheme.The Related Work group (as 802.22,802.16) of Institute of Electrical and Electric Engineers (being called for short IEEE) and International Telecommunications Union's (being called for short ITU) successively formulate or are actively formulating series of standards to promote the development of this technology under various application scenarioss.In cognitive radio, secondary user's can be in the situation that interfere with primary users communication insertion authority frequency range dynamically not, thereby improves the utilance of channel.In order not affect primary user's communication, secondary user's wanted the perception of advanced line frequency spectrum to determine that whether channel is occupied before the access channel, just secondary user's only can access when channel idle.So the frequency spectrum perception technology is to realize the key of cognitive radio technology.Frequency spectrum detection is the usual way of frequency spectrum perception technology.
Traditional frequency spectrum detection is taked the mode that channel is detected one by one, and the expense on time and energy is all very large.In the situation that huge, a plurality of secondary user's of channel number exists, this expense is especially obvious.The effective scheme addressed this problem is exactly spectrum prediction.Spectrum prediction is to predict next seizure condition of each channel constantly according to the account of the history of each channel occupancy.Secondary user's is pressed idle probability sequence detection from high to low to channel, thereby strengthens the probability of access, saves detection time and energy.The Chinese patent that application number is " 201010204687.2 ", disclose a kind of underwater multiplex communication means based on frequency spectrum perception and prediction.This invention utilizes the history of each channel is carried out to next acquistion probability constantly of predicted channel with statistics and the modeling of state, thereby judges its seizure condition.Its implementation procedure is: sensing node by predetermined time gap periods each channel is detected, obtain sufficient historical data; Historical occupied information by each channel is used Markov model to carry out modeling, obtains acquistion probability and the state transition probability of each channel; Then, utilize the Prediction of Markov model established at any one time the acquistion probability of each channel to be predicted.Replace the mode of direct-detection by the predicted channel seizure condition, can effectively improve message transmission rate, and save energy.
Yet, owing to carrying out channel estimating or while estimating, prior art is only considered the effect of the historic state information of channel to prediction, the influencing factor of consideration is too single, is only calculated according to the historical seizure condition of each channel, is difficult to guarantee the accuracy of prediction; In addition, only utilize the historical information of channel to be predicted the channel occupancy state, along with the time growth of prediction, can produce the phenomenon of error propagation.Therefore, concerning the spectrum prediction technology, the accuracy that how to improve as much as possible prediction is a problem of needing solution badly.
Summary of the invention
(1) technical problem to be solved
The object of the present invention is to provide a kind of frequency spectrum sensing method and device, to improve the accuracy of channel occupancy status predication, and improve frequency spectrum perception efficiency.
(2) technical scheme
In order to solve the problems of the technologies described above, the present invention proposes a kind of frequency spectrum sensing method, described method comprises:
All channels under same business are carried out to cluster, obtain the channel that represents of N channel set and each channel set, N is positive integer;
Determine the channel set treat under channel perception and the corresponding channel C that represents n, 1≤n≤N, treat the historical seizure condition of channel perception and represent channel C according to described nCurrent seizure condition, calculate the described current seizure condition for the treatment of channel perception.
Optionally, all channels under same business being carried out to cluster specifically comprises:
According to the historical seizure condition of all channels, calculate the coefficient correlation of every two interchannels, and according to coefficient correlation, all channels are carried out to cluster.
Optionally, the computing formula of the correlation coefficient ρ of every two interchannels is:
&rho; = E E + D ,
Wherein, the identical number of corresponding bit place value in the historical seizure condition sequence that E is two channels, the different number of corresponding bit place value in the historical seizure condition sequence that D is two channels.
Optionally, according to coefficient correlation, all channels being carried out to cluster specifically comprises:
Construct correlation matrix A according to coefficient correlation;
According to default coefficient correlation thresholding Th ρ, correlation matrix A is converted into to relevant judgement matrix B;
Convert respectively every a line of relevant judgement matrix B to the correlated channels set;
Adopt greedy algorithm, obtain the channel that represents of a described N channel set and each channel set according to described correlated channels set.
Optionally, the element of correlation matrix A is designated as ρ Ij, ρ IjFor channel c iWith channel c jBetween coefficient correlation, 1≤i≤W, 1≤j≤W, the number that W is all channels under described same business;
The element of relevant judgement matrix B is designated as b Ij, and have:
b ij = 1 &rho; ij &GreaterEqual; Th &rho; 0 &rho; ij < Th &rho; ;
Described correlated channels set is designated as s i, and s i={ c j| b Ij=1}.
Optionally, described greedy algorithm specifically comprises:
The set of all channels under described same business is designated as to U;
Calculate and make | s iThe correlated channels set s of ∩ U| maximum i, and by U ∩ s iBe designated as the first channel set, s iCorresponding channel c iBe designated as the channel that represents of described the first channel set;
Make U=U-s i, repeat previous step, obtain the second channel set and represent channel;
Repeat previous step, until
Figure BDA00003593426000033
Thereby, obtain the channel that represents of a described N channel set and each channel set.
Optionally, calculate and describedly treat that the current seizure condition of channel perception specifically comprises:
Obtain representing channel C by frequency spectrum detection nCurrent seizure condition i c
Adopt Markov model, according to described, treat the historical seizure condition sequence of channel perception and represent channel C nHistorical seizure condition sequence, calculate the described state transition probability p (m that treats channel perception c| m h) and describedly treat channel perception and represent channel C nBetween state transition probability p (i c| m cm h), wherein, m cFor the described current seizure condition for the treatment of channel perception, m hFor the described historical seizure condition sequence for the treatment of channel perception;
Calculate joint probability p (i according to condition probability formula cm cm h):
p(i cm cm h)=p(i c|m cm h)p(m c|m h)p(m h),
Wherein, p (m h) treat that for described the historical seizure condition sequence of channel perception is m hProbability;
Adopt the MAP criterion, according to joint probability p (i cm cm h) obtain m cEstimated value
Figure BDA00003593426000041
Thereby obtain the described current seizure condition for the treatment of channel perception
Optionally,
Figure BDA00003593426000042
Computing formula be:
m ^ c = arg max m c p ( i c m c m h ) .
The present invention has proposed a kind of frequency spectrum sensing device simultaneously, and described device comprises cluster module and estimation module, wherein:
Described cluster module, carry out cluster for all channels by under same business, obtains the channel that represents of N channel set and each channel set, and N is positive integer;
Described estimation module, for determining the channel set treat under channel perception and the corresponding channel C that represents n, 1≤n≤N, treat the historical seizure condition of channel perception and represent channel C according to described nCurrent seizure condition, calculate the described current seizure condition for the treatment of channel perception.
Optionally, described device further comprises memory module and detection module, wherein:
Described memory module is for storing the historical seizure condition of all channels under described same business;
Described detection module is for obtaining representing channel C by frequency spectrum detection nCurrent seizure condition;
Described cluster module is further used for obtaining from described memory module the historical seizure condition of all channels described same business, and carries out cluster according to the historical seizure condition of all channels;
Described estimation module is further used for obtaining the described historical seizure condition for the treatment of channel perception from described memory module, from described detection module, obtains and represents channel C nCurrent seizure condition, and calculate the described current seizure condition for the treatment of channel perception.
(3) beneficial effect
The correlation that the frequency spectrum sensing method that the present invention proposes and device have utilized each interchannel under same business, carry out cluster by all channels under same business, and one of selection represents channel from every class channel, carry out seizure condition when prediction treating channel perception, in conjunction with the historical seizure condition information for the treatment of channel perception self and the current seizure condition that represents channel, predicted the outcome.Adopt the method, only need to detect a certain current seizure condition that represents channel, just can calculate and predict the current seizure condition of whole channels in corresponding channel set.Because the method has not only been considered the historic state information of channel, also considered the correlation that takies of interchannel, therefore, compared with prior art, the accuracy of prediction is higher, and the success rate of access is larger, and efficiency is higher.In addition, because the state information that represents channel is used in the estimation for the treatment of the channel perception seizure condition, so each all will be detected the channel that represents in each channel set constantly, can role of correcting be arranged to mistake like this, thereby error propagation is had to obvious inhibitory action.
The accompanying drawing explanation
Fig. 1 is the realization flow figure of one embodiment of the invention intermediate frequency spectrum cognitive method.
Fig. 2 is the schematic diagram of channel cluster result in one embodiment of the invention.
Embodiment
Below in conjunction with drawings and Examples, the specific embodiment of the present invention is described in further detail.
In existing scheme when carrying out the channel occupancy status predication, all only just from each channel angle independently, take situation for the history of each channel respectively and carry out model and forecast, but ignored this part important information of relevance of each interchannel objective reality under same business.In fact, by actual measurement, find, in real communication environment, have very high relevance between the seizure condition of the different channels under same business.If can be used to the high relevance between different channels in this same business, can play very large castering action to forecasting accuracy and the detection efficiency of channel.
In view of this, the embodiment of the present invention provides a kind of frequency spectrum sensing method, and as shown in Figure 1, described method comprises:
All channels under same business are carried out to cluster, obtain the channel that represents of N channel set and each channel set, N is positive integer;
Determine the channel set treat under channel perception and the corresponding channel C that represents n, 1≤n≤N, treat the historical seizure condition of channel perception and represent channel C according to described nCurrent seizure condition, calculate the described current seizure condition for the treatment of channel perception.
The frequency spectrum sensing method that the present invention proposes has utilized the correlation of each interchannel under same business, carry out cluster by all channels under same business, and one of selection represents channel from every class channel, carry out seizure condition when prediction treating channel perception, in conjunction with the historical seizure condition information for the treatment of channel perception self and the current seizure condition that represents channel, predicted the outcome.Adopt the method, only need to detect a certain current seizure condition that represents channel, just can calculate and predict the current seizure condition of whole channels in corresponding channel set.Because the method has not only been considered the historic state information of channel, also considered the correlation that takies of interchannel, therefore, compared with prior art, the accuracy of prediction is higher, and the success rate of access is larger, and efficiency is higher.In addition, because the state information that represents channel is used in the estimation for the treatment of the channel perception seizure condition, so each all will be detected the channel that represents in each channel set constantly, can role of correcting be arranged to mistake like this, thereby error propagation is had to obvious inhibitory action.
The precondition that the technical scheme that the present invention proposes is carried out is to have carried out long detection on each channel of target service, has enough history to detect the deposit of data, in order to future channel status predicted and estimated.That is, at first need to utilize the historical historical seizure condition information that data obtain each channel under same business that detects.Usually channel seizure condition at a time means with bit 0 or 1, and 0 means idlely, and 1 expression takies, the historical seizure condition information of channel be exactly in fact one by 0,1 sequence formed.
The frequency spectrum sensing method that the present invention proposes comprises in channel cluster and class estimates two steps.
When the channel cluster, need to be according to the judgement of each inter-channel correlation under same business, utilize the channel clustering algorithm that all channels are polymerized to several classes (being channel set), and select one represent channel in each channel set, this represent channel be in class with the maximally related channel of other channels.
The channel clustering algorithm is mainly to carry out channel packet by the degree of correlation of interchannel seizure condition in same business, and what the degree of correlation of channel occupancy state was large is divided into one group.
Concrete, according to the historical seizure condition of all channels, calculate the coefficient correlation of every two interchannels, and according to coefficient correlation, all channels are carried out to cluster.
Preferably, the computing formula of the correlation coefficient ρ of every two interchannels is:
&rho; = E E + D ,
Wherein, the identical number of corresponding bit place value in the historical seizure condition sequence that E is two channels, the different number of corresponding bit place value in the historical seizure condition sequence that D is two channels.For example, the historical seizure condition sequence of two channels is 00001111000 and 00010010110, E=5 here, D=6.What in fact the correlation coefficient ρ here meaned is the correlation of the seizure condition of channel, and the correlation of the interchannel of the present invention's research all refers to the correlation of seizure condition.
Correlation coefficient ρ is larger, and the correlation of two channels is larger, illustrates that the seizure condition sequence of two interchannels is more similar.When the ρ of two interchannels is enough large, just can be more exactly the seizure condition of state information estimation one other channel by a channel.Extreme case, if ρ=1, can be by detecting a channel, the seizure condition of one other channel is estimated on zero defect ground.
The purpose of channel clustering algorithm is to guarantee ρ>=Th ρ(Th ρExpression coefficient correlation thresholding) in situation, all channels under same business are gathered for several channel sets, and select the channel that represents of a needs detection in each channel set, make packet count minimum.That is to say, in the situation that guarantee that the predictablity rate realization detects least number of times.
Preferably, the detailed process of channel cluster is as follows:
Construct correlation matrix A according to coefficient correlation;
According to default coefficient correlation thresholding Th ρ, correlation matrix A is converted into to relevant judgement matrix B;
Convert respectively every a line of relevant judgement matrix B to the correlated channels set;
Adopt greedy algorithm, obtain the channel that represents of a described N channel set and each channel set according to described correlated channels set.
For instance, for a certain business, suppose that it contains W channel altogether, is designated as C={c 1, c 2..., c W.
The element of correlation matrix A is designated as ρ Ij, ρ IjFor channel c iWith channel c jBetween coefficient correlation, 1≤i≤W, 1≤j≤W, matrix A is expressed as follows:
Figure BDA00003593426000081
The element of relevant judgement matrix B is designated as b Ij, and have:
b ij = 1 &rho; ij &GreaterEqual; Th &rho; 0 &rho; ij < Th &rho;
Matrix B has recorded the interchannel coefficient correlation and has been greater than default thresholding Th ρSituation.Known according to the above, in matrix B, be that 1 two channels corresponding to element meet very high correlation.Utilize matrix B to carry out channel packet and to select the most reasonably representing channel by algorithm for design.This problem and set covering problem are similar, will guarantee to meet Th exactly ρSituation under, cover C={c with minimum set number 1, c 2..., c W.
For this problem is converted into to set covering problem, establish S set={ s 1, s 2... be a selection set corresponding to matrix B, each element in S is a correlated channels set, can be designated as s i, s iEvery a line by matrix B converts respectively, and s i={ c j| b Ij=1}.
The example that S builds is as follows:
1 1 0 1 1 1 0 1 1 c 1 c 2 c 3 &RightArrow; B &RightArrow; S { { c 1 , c 2 } , { c 1 , c 2 , c 3 } , { c 2 . c 3 } }
The maximum s of containing element in each S iMean when detecting channel c iThe time, the classification that estimable channel is maximum.
Adopt didactic greedy algorithm to be solved this problem, described greedy algorithm specifically comprises:
The set of all channels under described same business is designated as to U;
Calculate and make | s iThe correlated channels set s of ∩ U| maximum i, and by U ∩ s iBe designated as the first channel set, s iCorresponding channel c iBe designated as the channel that represents of described the first channel set;
Make U=U-s i, repeat previous step, obtain the second channel set and represent channel;
Repeat previous step, until
Figure BDA00003593426000091
Thereby, obtain the channel that represents of a described N channel set and each channel set.
Said process also can be expressed as follows:
1 arranges U ← { c 1, c 2..., c W,
Figure BDA00003593426000092
2 while
Figure BDA00003593426000093
3 select s iMake | s i∩ U| maximum
4 g←U∩ is
5 U←U- is
6 G←G∪{g}
7 Return G
Here g is a channel set (also can say a class), s iCorresponding channel c iBe the channel that represents of g.Set G has comprised all channel sets.Fig. 2 is the schematic diagram of channel cluster result in one embodiment of the invention, and in figure, each square frame represents a channel set, and the white round dot in each square frame represents the channel that represents of this channel set, and the black round dot represents the channel perception for the treatment of in this channel set.Owing to adopting greedy formula algorithm, select maximum group number to be divided into groups, so the group separated after number of channel ratio in the group first separated is large at every turn.
After the channel cluster, need to be by detecting the current seizure condition that represents channel in a channel set, accurately estimate the busy-idle condition of all channels in this channel set in conjunction with self historical seizure condition information for the treatment of channel perception.
Below algorithm for estimating in class is described in detail.
At first determine the channel set treat under channel perception and the corresponding channel C that represents n, then to described, treat that the current seizure condition of channel perception carries out estimating in class.
The process of estimating in class is specific as follows:
Obtain representing channel C by frequency spectrum detection nCurrent seizure condition i c
Adopt Markov model, according to described, treat the historical seizure condition sequence of channel perception and represent channel C nHistorical seizure condition sequence, calculate the described state transition probability p (m that treats channel perception c| m h) and describedly treat channel perception and represent channel C nBetween state transition probability p (i c| m cm h), wherein, m cFor the described current seizure condition for the treatment of channel perception, m hFor the described historical seizure condition sequence for the treatment of channel perception;
Calculate joint probability p (i according to condition probability formula cm cm h):
p(i cm cm h)=p(i c|m cm h)p(m c|m h)p(m h),
Wherein, p (m h) treat that for described the historical seizure condition sequence of channel perception is m hProbability;
Adopt the MAP criterion, according to joint probability p (i cm cm h) obtain m cEstimated value
Figure BDA00003593426000101
Thereby, obtain the described current seizure condition for the treatment of channel perception.
For instance, suppose the channel that represents that the I channel is certain channel set, the M channel is that in this channel set treats channel perception.Consider to mean the impact of M channel history state information on current state with the steady Markov model in n rank, with p (m c| m h) mean that the state transition probability of M channel has:
p(m c=X t|m h=X t-1X t-2...X t-n...)=p(m c=X t|m h=X t-1X t-2...X t-n)
Wherein, m cThe state that means current time M channel, m hMean its historic state sequence, X mean the seizure condition value be 0 or 1, t mean corresponding constantly.
Due to I channel and M channel height correlation, so the state transition probability p (i of definition interchannel c| m cm h), m cm hThe status switch that can mean the M channel, comprise current time state and historic state sequence, i cIt is I channel current time state.As now selected the steady Markov model in n+1 rank to describe the state transitions process of interchannel, have following formula to set up:
p(i c=Y t|m cm h=X tX t-1X t-2...X t-n...)=p(i c=Y t|m cm h=X tX t-1X t-2...X t-n)
Because markoff process is stably, state-transition matrix time to time change not, can come out according to historical data.
In conjunction with the relevant information of historic state information and interchannel, according to condition probability formula, calculate joint probability p (i cm cm h):
p(i cm cm h)=p(i c|m cm h)p(m c|m h)p(m h)
P(i cm cm h) expression i cm cm hThe joint probability of state.
In conjunction with above-mentioned several new probability formula, adopt the MAP criterion to estimate m cValue.Because the historic state sequence of I channel current time state and M channel is known, according to MAP criterion, m cEstimated value
Figure BDA00003593426000111
The most reasonably be chosen as and make p (i cm cm h) that m that value is larger cValue, be shown below:
m ^ c = arg max p m c ( i c m c m h )
So far, just estimate to have obtained the current seizure condition of M channel.
The present invention has proposed a kind of frequency spectrum sensing device simultaneously, and described device comprises cluster module and estimation module, wherein:
Described cluster module, carry out cluster for all channels by under same business, obtains the channel that represents of N channel set and each channel set, and N is positive integer;
Described estimation module, for determining the channel set treat under channel perception and the corresponding channel C that represents n, 1≤n≤N, treat the historical seizure condition of channel perception and represent channel C according to described nCurrent seizure condition, calculate the described current seizure condition for the treatment of channel perception.
Preferably, described device further comprises memory module and detection module, wherein:
Described memory module is for storing the historical seizure condition of all channels under described same business;
Described detection module is for obtaining representing channel C by frequency spectrum detection nCurrent seizure condition;
Described cluster module is further used for obtaining from described memory module the historical seizure condition of all channels described same business, and carries out cluster according to the historical seizure condition of all channels;
Described estimation module is further used for obtaining the described historical seizure condition for the treatment of channel perception from described memory module, from described detection module, obtains and represents channel C nCurrent seizure condition, and calculate the described current seizure condition for the treatment of channel perception.
The correlation that the frequency spectrum sensing method that the present invention proposes and device have utilized each interchannel under same business, carry out cluster by all channels under same business, and one of selection represents channel from every class channel, carry out seizure condition when prediction treating channel perception, in conjunction with the historical seizure condition information for the treatment of channel perception self and the current seizure condition that represents channel, predicted the outcome.Adopt the method, only need to detect a certain current seizure condition that represents channel, just can calculate and predict the current seizure condition of whole channels in corresponding channel set.Because the method has not only been considered the historic state information of channel, also considered the correlation that takies of interchannel, therefore, compared with prior art, the accuracy of prediction is higher, and the success rate of access is larger, and efficiency is higher.In addition, because the state information that represents channel is used in the estimation for the treatment of the channel perception seizure condition, so each all will be detected the channel that represents in each channel set constantly, can role of correcting be arranged to mistake like this, thereby error propagation is had to obvious inhibitory action.
Simultaneously, the technical scheme that the present invention proposes also can be secondary user's and provides convenience.In existing scheme, if, in the situation that not communication between secondary user's all can may be selected access in idle channel at same group, be easy to cause the competition conflict.Under the cognitive system existed in a plurality of secondary user's, technical scheme of the present invention can arrange secondary user's to carry out dynamic access in different channel sets, the conflict of just can not competing between secondary user's like this.
The above is only the preferred embodiment of the present invention; it should be pointed out that for the person of ordinary skill of the art, under the prerequisite that does not break away from the technology of the present invention principle; can also make some improvement and replacement, these improvement and replacement also should be considered as protection scope of the present invention.

Claims (10)

1. a frequency spectrum sensing method, is characterized in that, described method comprises:
All channels under same business are carried out to cluster, obtain the channel that represents of N channel set and each channel set, N is positive integer;
Determine the channel set treat under channel perception and the corresponding channel C that represents n, 1≤n≤N, treat the historical seizure condition of channel perception and represent channel C according to described nCurrent seizure condition, calculate the described current seizure condition for the treatment of channel perception.
2. the method for claim 1, is characterized in that, all channels under same business carried out to cluster and specifically comprise:
According to the historical seizure condition of all channels, calculate the coefficient correlation of every two interchannels, and according to coefficient correlation, all channels are carried out to cluster.
3. method as claimed in claim 2, is characterized in that, the computing formula of the correlation coefficient ρ of every two interchannels is:
&rho; = E E + D ,
Wherein, the identical number of corresponding bit place value in the historical seizure condition sequence that E is two channels, the different number of corresponding bit place value in the historical seizure condition sequence that D is two channels.
4. method as claimed in claim 2 or claim 3, is characterized in that, according to coefficient correlation, all channels carried out to cluster and specifically comprise:
Construct correlation matrix A according to coefficient correlation;
According to default coefficient correlation thresholding Th ρ, correlation matrix A is converted into to relevant judgement matrix B;
Convert respectively every a line of relevant judgement matrix B to the correlated channels set;
Adopt greedy algorithm, obtain the channel that represents of a described N channel set and each channel set according to described correlated channels set.
5. method as claimed in claim 4, is characterized in that, the element of correlation matrix A is designated as ρ Ij, ρ IjFor channel c iWith channel c jBetween coefficient correlation, 1≤i≤W, 1≤j≤W, the number that W is all channels under described same business;
The element of relevant judgement matrix B is designated as b Ij, and have:
b ij = 1 &rho; ij &GreaterEqual; Th &rho; 0 &rho; ij < Th &rho; ;
Described correlated channels set is designated as s i, and s i={ c j| b Ij=1}.
6. method as claimed in claim 5, is characterized in that, described greedy algorithm specifically comprises:
The set of all channels under described same business is designated as to U;
Calculate and make | s iThe correlated channels set s of ∩ U| maximum i, and by U ∩ s iBe designated as the first channel set, s iCorresponding channel c iBe designated as the channel that represents of described the first channel set;
Make U=U-s i, repeat previous step, obtain the second channel set and represent channel;
Repeat previous step, until
Figure FDA00003593425900025
Thereby, obtain the channel that represents of a described N channel set and each channel set.
7. the method for claim 1, is characterized in that, calculates describedly to treat that the current seizure condition of channel perception specifically comprises:
Obtain representing channel C by frequency spectrum detection nCurrent seizure condition i c
Adopt Markov model, according to described, treat the historical seizure condition sequence of channel perception and represent channel C nHistorical seizure condition sequence, calculate the described state transition probability p (m that treats channel perception c| m h) and describedly treat channel perception and represent channel C nBetween state transition probability p (i c| m cm h), wherein, m cFor the described current seizure condition for the treatment of channel perception, m hFor the described historical seizure condition sequence for the treatment of channel perception;
Calculate joint probability p (i according to condition probability formula cm cm h):
p(i cm cm h)=p(i c|m cm h)p(m c|m h)p(m h),
Wherein, p (m h) treat that for described the historical seizure condition sequence of channel perception is m hProbability;
Adopt the MAP criterion, according to joint probability p (i cm cm h) obtain m cEstimated value Thereby obtain the described current seizure condition for the treatment of channel perception.
8. method as claimed in claim 7, is characterized in that,
Figure FDA00003593425900023
Computing formula be:
m ^ c = arg max m c p ( i c m c m h ) .
9. a frequency spectrum sensing device, is characterized in that, described device comprises cluster module and estimation module, wherein:
Described cluster module, carry out cluster for all channels by under same business, obtains the channel that represents of N channel set and each channel set, and N is positive integer;
Described estimation module, for determining the channel set treat under channel perception and the corresponding channel C that represents n, 1≤n≤N, treat the historical seizure condition of channel perception and represent channel C according to described nCurrent seizure condition, calculate the described current seizure condition for the treatment of channel perception.
10. device as claimed in claim 9, is characterized in that, described device further comprises memory module and detection module, wherein:
Described memory module is for storing the historical seizure condition of all channels under described same business;
Described detection module is for obtaining representing channel C by frequency spectrum detection nCurrent seizure condition;
Described cluster module is further used for obtaining from described memory module the historical seizure condition of all channels described same business, and carries out cluster according to the historical seizure condition of all channels;
Described estimation module is further used for obtaining the described historical seizure condition for the treatment of channel perception from described memory module, from described detection module, obtains and represents channel C nCurrent seizure condition, and calculate the described current seizure condition for the treatment of channel perception.
CN201310326117.4A 2013-07-30 2013-07-30 A kind of frequency spectrum sensing method and device Expired - Fee Related CN103428704B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201310326117.4A CN103428704B (en) 2013-07-30 2013-07-30 A kind of frequency spectrum sensing method and device
PCT/CN2013/083253 WO2015014005A1 (en) 2013-07-30 2013-09-10 Method and apparatus for sensing spectrum

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310326117.4A CN103428704B (en) 2013-07-30 2013-07-30 A kind of frequency spectrum sensing method and device

Publications (2)

Publication Number Publication Date
CN103428704A true CN103428704A (en) 2013-12-04
CN103428704B CN103428704B (en) 2016-08-10

Family

ID=49652725

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310326117.4A Expired - Fee Related CN103428704B (en) 2013-07-30 2013-07-30 A kind of frequency spectrum sensing method and device

Country Status (2)

Country Link
CN (1) CN103428704B (en)
WO (1) WO2015014005A1 (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103957537A (en) * 2014-03-26 2014-07-30 北京邮电大学 Frequency spectrum perception method and apparatus based on channel cluster
CN106231603A (en) * 2016-07-27 2016-12-14 北京邮电大学 A kind of spectrum prediction method and device
CN107592174A (en) * 2017-09-12 2018-01-16 河南工业大学 A kind of high-efficiency frequency spectrum cognitive method in intelligent grid communication
CN108010282A (en) * 2017-11-30 2018-05-08 中冶南方工程技术有限公司 The method for drafting of figure and the rescue skills in toxic gas place when toxic gas is distributed
CN108616318A (en) * 2016-08-29 2018-10-02 电子科技大学 A kind of safe frequency spectrum sensing method
CN110610192A (en) * 2019-08-07 2019-12-24 电子科技大学 Spectrum space channel clustering method
CN111770504A (en) * 2020-07-06 2020-10-13 北京邮电大学 Method and device for determining occupancy rate of downlink resources of WCDMA (wideband code division multiple Access) base station

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106685549B (en) * 2016-11-17 2020-06-12 北京邮电大学 Primary user spectrum sensing method and device

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101370282A (en) * 2007-08-16 2009-02-18 华为技术有限公司 Method and apparatus for signal channel perception by cognition system
US20090323835A1 (en) * 2008-06-30 2009-12-31 Motorola, Inc . Method and apparatus for detection of orthogonal frequency division multiplexing (ofdm) signals by cognitive radios
CN101711032A (en) * 2009-11-23 2010-05-19 哈尔滨工业大学 Cognitive radio electric dynamic smart frequency spectrum access method for unknown environmental model characteristics
CN101867420A (en) * 2010-06-18 2010-10-20 华南理工大学 Underwater multi-path communication method based on frequency spectrum perception and prediction

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101447571B1 (en) * 2008-01-21 2014-10-07 삼성전자주식회사 Cognitve radio communication method for controlling sensing operation and cognitve radio device enabling the method
CN102638802B (en) * 2012-03-26 2014-09-03 哈尔滨工业大学 Hierarchical cooperative combined spectrum sensing algorithm
CN103139874B (en) * 2012-12-17 2019-04-12 中国人民解放军陆军工程大学 Based on the channel selecting method of time series forecasting in cognitive radio

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101370282A (en) * 2007-08-16 2009-02-18 华为技术有限公司 Method and apparatus for signal channel perception by cognition system
US20090323835A1 (en) * 2008-06-30 2009-12-31 Motorola, Inc . Method and apparatus for detection of orthogonal frequency division multiplexing (ofdm) signals by cognitive radios
US8451917B2 (en) * 2008-06-30 2013-05-28 Motorola Solutions, Inc. Method and apparatus for detection of orthogonal frequency division multiplexing (OFDM) signals by cognitive radios
CN101711032A (en) * 2009-11-23 2010-05-19 哈尔滨工业大学 Cognitive radio electric dynamic smart frequency spectrum access method for unknown environmental model characteristics
CN101867420A (en) * 2010-06-18 2010-10-20 华南理工大学 Underwater multi-path communication method based on frequency spectrum perception and prediction

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103957537A (en) * 2014-03-26 2014-07-30 北京邮电大学 Frequency spectrum perception method and apparatus based on channel cluster
CN103957537B (en) * 2014-03-26 2017-11-14 北京邮电大学 Frequency spectrum sensing method and device based on channel cluster
CN106231603A (en) * 2016-07-27 2016-12-14 北京邮电大学 A kind of spectrum prediction method and device
CN106231603B (en) * 2016-07-27 2019-05-10 北京邮电大学 A kind of spectrum prediction method and device
CN108616318A (en) * 2016-08-29 2018-10-02 电子科技大学 A kind of safe frequency spectrum sensing method
CN108616318B (en) * 2016-08-29 2020-05-12 电子科技大学 Secure spectrum sensing method
CN107592174A (en) * 2017-09-12 2018-01-16 河南工业大学 A kind of high-efficiency frequency spectrum cognitive method in intelligent grid communication
CN108010282A (en) * 2017-11-30 2018-05-08 中冶南方工程技术有限公司 The method for drafting of figure and the rescue skills in toxic gas place when toxic gas is distributed
CN110610192A (en) * 2019-08-07 2019-12-24 电子科技大学 Spectrum space channel clustering method
CN110610192B (en) * 2019-08-07 2022-03-15 电子科技大学 Spectrum space channel clustering method
CN111770504A (en) * 2020-07-06 2020-10-13 北京邮电大学 Method and device for determining occupancy rate of downlink resources of WCDMA (wideband code division multiple Access) base station

Also Published As

Publication number Publication date
WO2015014005A1 (en) 2015-02-05
CN103428704B (en) 2016-08-10

Similar Documents

Publication Publication Date Title
CN103428704A (en) Method and device for sensing frequency spectra
JP6741810B2 (en) Electronic device, method and computer-readable medium
CN103987056A (en) Wireless network telephone traffic prediction method based on big-data statistical model
CN107070569A (en) Multipoint cooperative frequency spectrum sensing method based on HMM model
CN105188124A (en) Robustness gaming power control method under imperfect CSI for multi-user OFDMA relay system
CN109039505B (en) Channel state transition probability prediction method in cognitive radio network
CN104796899A (en) Self-adapted cooperative frequency spectrum sensing method and system based on optimal relaying
CN103916969A (en) Combined authorized user perception and link state estimation method and device
Subramanian Large deviations of max-weight scheduling policies on convex rate regions
CN103249050B (en) Based on the multiple dimensioned frequency spectrum access method of business demand
Karmokar et al. Physical layer‐optimal and cross‐layer channel access policies for hybrid overlay–underlay cognitive radio networks
CN103117821B (en) Weighing cooperation spectrum sensing method based on rayleigh quotient
CN105119669A (en) Clustering cooperative spectrum sensing method for cognitive radio network
CN102857305A (en) Multi-node combination spectrum sensing method and system
CN106972900A (en) Based on broad sense T2The blind frequency spectrum sensing method of statistic
CN103281701B (en) A kind of collaboration frequency spectrum Forecasting Methodology of cognitive radio networks
Yifei et al. QoS Provisioning energy saving dynamic access policy for overlay cognitive radio networks with hidden Markov channels
Tong et al. Cooperative spectrum sensing based on a modified shuffled frog leaping algorithm in 5G network
CN105007585A (en) Power distribution method based on maximum outage probability energy efficiency
Zhang et al. Estimation of channel state transition probabilities based on Markov chains in cognitive radio.
CN104469811A (en) Clustering cooperative spectrum sensing hard fusion method for cognitive wireless sensor network
CN104994511B (en) Authorized user&#39;s method for parameter estimation in cognitive radio system based on cooperative sensing
Li et al. Joint spectrum sensing and transmission strategy for energy-efficient cognitive radio networks
Yao et al. Cognitive radio spectrum decision based on channel usage prediction
CN101860884A (en) Shared-frequency neighboring cell channel estimation 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
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20171108

Address after: 541004 Guangxi Guilin Qixing District Chaoyang Road Information Industry Park D-04 thinking

Patentee after: Guilin CEKE Communication Equipment Co.,Ltd.

Address before: 100876 Beijing city Haidian District Xitucheng Road No. 10

Patentee before: Beijing University of Posts and Telecommunications

TR01 Transfer of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160810

CF01 Termination of patent right due to non-payment of annual fee