CN105141385B - Multiband cooperative cognitive frequency spectrum sensing method - Google Patents

Multiband cooperative cognitive frequency spectrum sensing method Download PDF

Info

Publication number
CN105141385B
CN105141385B CN201510595891.4A CN201510595891A CN105141385B CN 105141385 B CN105141385 B CN 105141385B CN 201510595891 A CN201510595891 A CN 201510595891A CN 105141385 B CN105141385 B CN 105141385B
Authority
CN
China
Prior art keywords
prime
snr
user
primary
noise ratio
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510595891.4A
Other languages
Chinese (zh)
Other versions
CN105141385A (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.)
Yuhang District Artificial Intelligence Industry Intellectual Property Alliance In Hangzhou City
Original Assignee
Ningbo University
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 Ningbo University filed Critical Ningbo University
Priority to CN201510595891.4A priority Critical patent/CN105141385B/en
Publication of CN105141385A publication Critical patent/CN105141385A/en
Application granted granted Critical
Publication of CN105141385B publication Critical patent/CN105141385B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Telephonic Communication Services (AREA)

Abstract

本发明涉及多频段协作认知频谱感知方法,各次用户分别将自身信噪比、置信度、业务需求以及对多个主用户的频谱检测结果发送给频谱感知融合中心,由频谱感知融合中心计算、筛选初选协作次用户,然后分配频段数量给初选协作次用户;计算所有初选次用户信噪比均方根值与各初选次用户的信噪比间的商值,根据各商值与信噪比预设阈值间的关系,选定参与协作的复选次用户,并调整复选次用户的检测概率;结合各复选次用户的信噪比和调整后的平均检测概率,重新返回选择复选次用户的步骤,得到最终参与协作的终选次用户,并以加权的OR准则协作的全局检测概率为频谱感知融合中心的最终检测结果,避免了低信噪比或性能差的次用户对整体检测性能的恶劣影响。

The invention relates to a multi-band cooperative cognitive spectrum sensing method. Each user sends its own signal-to-noise ratio, confidence level, service requirements, and spectrum detection results of multiple primary users to a spectrum sensing fusion center, and the spectrum sensing fusion center calculates 1. Screen the primary cooperative sub-users, and then assign the number of frequency bands to the primary cooperative sub-users; calculate the quotient value between the root mean square value of the signal-to-noise ratio of all primary sub-users and the signal-to-noise ratio of each primary sub-user, according to each quotient value and the preset threshold of signal-to-noise ratio, select the secondary users participating in the collaboration, and adjust the detection probability of secondary users; combine the signal-to-noise ratio of each secondary user and the adjusted average detection probability, Go back to the step of selecting secondary users to get the final selected secondary users who finally participate in the cooperation, and use the weighted OR criterion to cooperate with the global detection probability as the final detection result of the spectrum sensing fusion center, avoiding low signal-to-noise ratio or poor performance Secondary users have a bad influence on the overall detection performance.

Description

多频段协作认知频谱感知方法Multi-band Cognitive Cognitive Spectrum Sensing Method

技术领域technical field

本发明涉及频谱检测领域,尤其涉及一种多频段协作认知频谱感知方法。The invention relates to the field of spectrum detection, in particular to a multi-band cooperative cognitive spectrum sensing method.

背景技术Background technique

随着无线通信技术的不断发展,以LTE、Wi-Fi、卫星通信及协同通信等为标志的新兴技术相继涌现,层出不穷。这些通信技术对无线频谱资源提出了更高的需求,从而令频谱资源变的趋于紧张,认知无线电技术(Cognitive Radio,CR)在此背景下应运而生。With the continuous development of wireless communication technology, emerging technologies marked by LTE, Wi-Fi, satellite communication and collaborative communication have emerged one after another, emerging in endlessly. These communication technologies put forward higher demands on wireless spectrum resources, thus making the spectrum resources more tense, and cognitive radio technology (Cognitive Radio, CR) emerges under this background.

认知无线电的基本途径是,首先次用户(或称感知用户、认知用户)采用频谱感知对周围环境中的已授权频谱资源进行持续检测;然后在保证主用户(又称授权用户)能够优先占用该段频谱且传输性能几乎不受影响的条件下,次用户自适应地调整收发设备,并将收发设备调整至空闲频谱上进行通信。当次用户感知(或称检测)到有主用户信号出现时,次用户则要快速腾出信道供主用户使用,进而避免对主用户的正常通信进行干扰,从而提供频谱资源利用率。The basic approach of cognitive radio is that first, the secondary user (or cognitive user, cognitive user) uses spectrum sensing to continuously detect the authorized spectrum resources in the surrounding environment; Under the condition that the frequency spectrum is occupied and the transmission performance is almost unaffected, the secondary user adaptively adjusts the transceiver equipment, and adjusts the transceiver equipment to the idle frequency spectrum for communication. When the secondary user perceives (or detects) that there is a signal of the primary user, the secondary user must quickly vacate the channel for the primary user to use, thereby avoiding interference to the normal communication of the primary user, thereby improving spectrum resource utilization.

为了减少实际环境中多径衰落、阴影效应和噪声不确定性等诸多因素对检测性能的不利影响,基于多个次用户的协作频谱感知方法被不断提出。通过将每个次用户的检测结果发送给频谱感知融合中心进行融合,以达到对频谱进行感知的目的。In order to reduce the adverse effects of many factors such as multipath fading, shadowing effect and noise uncertainty on the detection performance in the actual environment, cooperative spectrum sensing methods based on multiple secondary users have been proposed continuously. The detection result of each secondary user is sent to the spectrum sensing fusion center for fusion to achieve the purpose of spectrum sensing.

然而,现有的协作频谱感知方法多数只是针对单频段进行感知,为了提高频谱利用率,针对多频段的协作感知方法成为新的研究热点。但是,在多频段的协作感知中,由于各次用户的检测性能、信噪比并非完全是优良的,当所有的次用户均参与到对多频段的协作感知时,低信噪比或检测性能差的次用户会对整体的检测性能造成恶劣影响。However, most of the existing cooperative spectrum sensing methods are only for single-band sensing. In order to improve spectrum utilization, multi-band cooperative sensing methods have become a new research hotspot. However, in multi-band cooperative sensing, since the detection performance and signal-to-noise ratio of each user are not completely excellent, when all secondary users participate in multi-band cooperative sensing, the low signal-to-noise ratio or detection performance Poor secondary users can have a bad impact on the overall detection performance.

发明内容Contents of the invention

本发明所要解决的技术问题是针对上述现有技术提供一种能够避免具有低信噪比或检测性能差的次用户对整体检测性能造成恶劣影响的多频段协作认知频谱感知方法。The technical problem to be solved by the present invention is to provide a multi-band cooperative cognitive spectrum sensing method that can avoid the adverse impact of secondary users with low signal-to-noise ratio or poor detection performance on the overall detection performance.

本发明解决上述技术问题所采用的技术方案为:多频段协作认知频谱感知方法,其特征在于,依次包括如下步骤:The technical solution adopted by the present invention to solve the above-mentioned technical problems is: a multi-band cooperative cognitive spectrum sensing method, which is characterized in that it includes the following steps in sequence:

(1)设定在认知网络中,主用户的数量为N1,次用户的数量为N2,频谱感知融合中心的数量为1,主用户分别独立地占用频谱中的各自频段;N2个次用户分别独立地获取自身信噪比SNRj以及对N1个主用户占用频段的频谱检测结果,并分别将获取的信噪比SNRj、频谱检测结果、检测置信度Pj和业务需求Rj发送至频谱感知融合中心,其中,(1) In the cognitive network, the number of primary users is N 1 , the number of secondary users is N 2 , the number of spectrum sensing fusion centers is 1, and the primary users independently occupy their respective frequency bands in the spectrum; N 2 Each secondary user independently obtains its own SNR j and the spectrum detection results of the frequency bands occupied by N 1 primary users, and respectively obtains the SNR j , spectrum detection results, detection confidence P j and service requirements R j is sent to the spectrum sensing fusion center, where,

主用户记为PUi,次用户记为CRj,频谱感知融合中心记为FC,业务需求Rj∈[0,N1],Pj∈[0,1],检测置信度Pd,ji为次用户CRj对主用户PUi的检测概率;1≤i≤N1,1≤j≤N2,N1≥2,N2≥2;The primary user is denoted as PU i , the secondary user is denoted as CR j , the spectrum sensing fusion center is denoted as FC, business requirements R j ∈ [0,N 1 ], P j ∈ [0,1], detection confidence P d,ji is the detection probability of secondary user CR j to primary user PU i ; 1≤i≤N 1 , 1≤j≤N 2 , N 1 ≥2, N 2 ≥2;

(2)频谱感知融合中心FC接收各次用户CRj发送来的信噪比SNRj、频谱检测结果、检测置信度Pj和业务需求Rj,并判断次用户的信噪比SNRj大于预设的信噪比筛选值SNRchose时,选择此信噪比对应的次用户为参与协作检测的初选次用户,并记初选次用户为CR't,并执行步骤(3),否则,选择具有最高信噪比的次用户所对应的频谱检测结果为频谱感知融合中心FC的最终检测结果;其中,(2) The spectrum sensing fusion center FC receives the signal-to-noise ratio SNR j , spectrum detection results, detection confidence P j and service requirements R j sent by each user CR j , and judges that the signal-to-noise ratio SNR j of the secondary user is greater than the expected When the selected signal-to-noise ratio screening value SNR chosen , select the secondary user corresponding to the signal-to-noise ratio as the primary secondary user participating in the collaborative detection, and record the primary secondary user as CR' t , and perform step (3), otherwise, Select the spectrum detection result corresponding to the secondary user with the highest SNR as the final detection result of the spectrum sensing fusion center FC; where,

初选次用户数量为N'2,初选次用户CR't对应的信噪比为SNR't、检测置信度为P't、业务需求为R't,1≤t≤N'2≤N2;频谱感知融合中心FC接收到的频谱检测结果数量为N1×N2个;The number of primary secondary users is N' 2 , the signal-to-noise ratio corresponding to the primary secondary user CR' t is SNR' t , the detection confidence is P' t , and the service requirement is R' t , 1≤t≤N' 2 ≤ N 2 ; the number of spectrum detection results received by the spectrum sensing fusion center FC is N 1 × N 2 ;

(3)频谱感知融合中心FC根据初选次用户CR't的置信度P't、业务需求R't,对初选次用户CR't分配需要检测的频段数量Ct(3) The spectrum sensing fusion center FC assigns the number of frequency bands C t to be detected to the primary secondary user CR' t according to the confidence P' t of the primary secondary user CR' t and the service demand R' t :

(3-1)根据各初选次用户CR't的置信度P't,分别对每一个初选次用户CR't的置信度P't进行归一化,得到每个初选次用户CR't的归一化置信度值 (3-1) According to the confidence degree P' t of each primary selected secondary user CR' t , respectively normalize the confidence P' t of each primary selected secondary user CR' t to obtain each primary selected secondary user Normalized confidence value of CR't

(3-2)根据步骤(3-1)所得每个初选次用户CR't对应的归一化置信度值计算频谱感知融合中心FC分配初选次用户CR't需要检测的频段数量Ct(3-2) The normalized confidence value corresponding to each primary user CR' t obtained according to step (3-1) Calculate the number of frequency bands C t that need to be detected when the spectrum sensing fusion center FC assigns the primary user CR' t to be detected:

(4)频谱感知融合中心FC根据参与协作检测的初选次用户CR't的信噪比SNR't,计算所有初选次用户的信噪比均方根值并令信噪比SNR't=γt,其中,信噪比均方根值的计算如下:(4) The spectrum sensing fusion center FC calculates the root mean square value of the signal-to-noise ratio of all primary secondary users according to the signal-to-noise ratio SNR' t of primary secondary users CR' t participating in cooperative detection And let the signal-to-noise ratio SNR' t = γ t , where the root mean square value of the signal-to-noise ratio is calculated as follows:

(5)频谱感知融合中心FC分别依次计算所有初选次用户的信噪比均方根值与各初选次用户CR't的信噪比SNR't之间的商值ηt,其中,(5) The spectrum sensing fusion center FC calculates the root mean square value of the signal-to-noise ratio of all primary secondary users in turn The quotient η t between the signal-to-noise ratio SNR' t of each primary secondary user CR' t , where

(6)频谱感知融合中心FC计算、获取信噪比预设阈值λ和信噪比最佳阈值λoptimal,并分别根据各信噪比商值ηt与信噪比预设阈值λ之间的大小关系,选定参与协作的复选次用户CR”k,复选次用户CR”k的信噪比为SNR”k,其中,(6) The spectrum sensing fusion center FC calculates and obtains the SNR preset threshold λ and the SNR optimal threshold λ optimal , and respectively calculates and obtains the SNR quotient value η t and the SNR preset threshold λ according to The size relationship, select the secondary user CR” k to participate in the collaboration, and the signal-to-noise ratio of the secondary user CR” k is SNR” k , where,

(6-1)频谱感知融合中心FC根据接收的N'2个初选从用户对应的信噪比集合{SNR't},获取初选从用户信噪比集合{SNR't}中的信噪比最大值,记该信噪比最大值为SNR'max;(6-1) The spectrum sensing fusion center FC acquires the signal-to-noise ratio set {SNR' t } of the primary-selected slave users according to the received SNR sets {SNR' t } corresponding to the N' 2 primary-selected slave users. The maximum value of the noise ratio, record the maximum value of the signal-to-noise ratio as SNR'max;

(6-2)以获取的信噪比最大值SNR'max为参考,并将N'2个初选从用户CR't的信噪比SNR't分别与信噪比最大值SNR'max作商处理,计算得到各初选从用户信噪比SNR't所对应的初始阈值λt,其中,(6-2) Take the obtained SNR maximum value SNR'max as a reference, and make the SNR' t of the N' 2 primary users CR' t with the SNR maximum value SNR'max respectively quotient processing, and calculate the initial threshold λ t corresponding to each primary user SNR' t , where,

λt=|SNR't/SNR'max|,t=1,2,…,N'2,N'2≤N2λ t = |SNR' t /SNR' max |, t=1,2,...,N' 2 , N' 2 ≤ N 2 ;

(6-3)根据各初选从用户CR't的归一化置信度值和信噪比商值ηt,计算各初选从用户CR't的联合筛选参数值ξt,并根据联合筛选参数值ξt,选取参与协作的复选从用户CR”k,其中,复选从用户CR”k的数量为M,t=1,2,…,N'2,k=1,2,…,M,M≤N'2(6-3) According to the normalized confidence value of each primary selection from the user CR' t and the signal-to-noise ratio quotient η t , calculate the joint screening parameter value ξ t of each primary secondary user CR' t , and select the secondary user CR" k participating in the collaboration according to the joint screening parameter value ξ t , where the complex The number of selected users CR” k is M, t=1,2,...,N' 2 , k=1,2,...,M, M≤N' 2 :

若联合筛选参数值ξt位于预设数值区间范围[ξab]内,即ξa≤ξt≤ξb时,则选取该联合筛选参数值ξt对应的初选从用户为复选从用户,并参与协作检测;否则,该初选从用户不予选取;If the joint screening parameter value ξ t is within the preset value range [ξ a , ξ b ], that is, when ξ a ≤ ξ t ≤ ξ b , then select the primary user corresponding to the joint screening parameter value ξ t as complex Select a secondary user and participate in collaborative detection; otherwise, the primary secondary user will not be selected;

(6-4)根据步骤(6-3)中的信噪比预设阈值λ,获取M个复选从用户CR”k分别在OR准则和AND准则下的协作检测性能曲线,其中,(6-4) According to the signal-to-noise ratio preset threshold λ in step (6-3), obtain the cooperative detection performance curves of M re-selected slave users CR" k under the OR criterion and the AND criterion respectively, wherein,

OR准则: OR criterion:

AND准则:k=1,2,…,M,M≤N'2≤N2AND criteria: k=1,2,...,M, M≤N' 2 ≤N 2 ;

其中,Pd,k为第k个复选从用户CR”k的平均检测概率,Pfa,k为第k个复选从用户CR”k的平均虚警概率;Pd,ks为复选从用户CR”k对其所分配第s个频段的检测概率,Pfa,ks为复选从用户CR”k对其所分配第s个频段的虚警概率;Qd为协作检测后的全局检测概率,Qfa为协作检测后的全局虚警概率;ωk表示信噪比CR”k的权重系数,SNR”k是第k个复选从用户CR”k的信噪比,SNR”max表示M个复选从用户的信噪比最大值,SNR”min表示M个复选从用户的信噪比最小值;Among them, P d, k is the average detection probability of the k-th re-selected slave user CR” k , P fa, k is the average false alarm probability of the k-th re-selected slave user CR” k ; P d, ks is the re-selected The detection probability of the s-th frequency band allocated from user CR” k , P fa,ks is the false alarm probability of the s-th frequency band allocated from user CR” k to it; Q d is the global Detection probability, Q fa is the global false alarm probability after collaborative detection; ω k represents the weight coefficient of signal-to-noise ratio CR" k , SNR" k is the signal-to-noise ratio of the kth re-selected user CR" k , SNR" max Indicates the maximum value of the SNR of the M secondary users, and SNR" min represents the minimum value of the SNR of the M secondary users;

(6-5)根据OR准则和AND准则下的协作检测性能曲线,分别得到在OR准则和AND准则下的最大检测概率Q(OR,d)-max、Q(AND,d)-max,得到Q(OR,d)-max和Q(AND,d)-max的最大值Qd-max,并以该最佳检测性能值Qd-max所对应的信噪比预设阈值为信噪比最佳阈值,记信噪比最佳阈值为λoptimal;其中,Qd-max=max(Q(OR,d)-max,Q(AND,d)-max);(6-5) According to the collaborative detection performance curves under the OR criterion and the AND criterion, the maximum detection probabilities Q (OR, d)-max and Q (AND, d)-max under the OR criterion and the AND criterion are respectively obtained, and The maximum value of Q (OR, d)-max and Q (AND, d)-max Q d-max , and the signal-to-noise ratio preset threshold corresponding to the best detection performance value Q d-max is the signal-to-noise ratio The optimal threshold, the optimal threshold of SNR is λ optimal ; where, Q d-max = max(Q (OR, d)-max , Q (AND, d)-max );

(7)根据获取的信噪比最佳阈值λoptimal,得到该信噪比最佳阈值λoptimal对应的复选从用户CR”,获取该复选从用户CR”的调整因子α以及其他M-1个复选从用户CR”k的调整因子αk,并分别根据调整因子α、αk对应调整复选从用户CR”、CR”k的平均虚警概率,复选从用户CR”调整后的平均虚警概率记为Pfa,复选从用户CR”k调整后的平均虚警概率记为Pfa,k;其中,(7) According to the obtained SNR optimal threshold λ optimal , obtain the multiple-selection slave user CR” corresponding to the SNR optimal threshold λ optimal , obtain the adjustment factor α of the multiple-selection slave user CR” and other M- The adjustment factor α k of one reselected slave user CR” k , and adjust the average false alarm probability of reselected slave users CR” and CR” k according to the adjustment factors α and α k respectively. After the reselected slave user CR” is adjusted The average false alarm probability is denoted as P fa , and the average false alarm probability adjusted from user CR” k is denoted as P fa,k ; where,

Pfa,k=αk·Pfa,k=1,2,…,M-1;P fa,kk P fa ,k=1,2,...,M-1;

其中,αk为复选从用户CR”k的调整因子,用来根据复选从用户CR”k自身的信噪比SNR”k实现对其平均虚警概率大小的调整;SNR”k为第k个复选从用户CR”k的信噪比;Among them, α k is the adjustment factor of the secondary user CR" k , which is used to adjust the average false alarm probability according to the signal-to-noise ratio SNR" k of the secondary user CR" k itself; SNR of k multiple choices from user CR"k;

(8)根据步骤(7)中获取的M个复选从用户的调整因子αk以及对应调整后的平均虚警概率Pfa,k,计算复选从用户CR”k调整后的判决阀值λ”k和平均检测概率Pd,k,其中,(8) According to the adjustment factor α k of the M secondary users obtained in step (7) and the corresponding adjusted average false alarm probability P fa,k , calculate the adjusted decision threshold of the secondary user CR” k λ” k and average detection probability P d,k , where,

其中,k=1,2,…,M,M≤N'2;n为采样点数;in, k=1,2,...,M, M≤N'2; n is the number of sampling points;

(9)根据步骤(8)中M个复选从用户的信噪比SNR”k以及得到的调整后的平均检测概率Pd,k,返回步骤(6),重新在M个复选从用户中选择,得到参与协作的T个终选从用户CR”'t,并以加权的OR准则协作后的全局检测概率为频谱感知融合中心FC的最终检测结果,其中1≤t≤T≤M≤N'2(9) According to the signal-to-noise ratio SNR” k of the M reselected slave users in step (8) and the adjusted average detection probability P d,k obtained, return to step (6), and re-select the M reselected slave users Select among T final users CR"' t participating in the collaboration, and use the weighted OR criterion to obtain the global detection probability after collaboration as the final detection result of the spectrum sensing fusion center FC, where 1≤t≤T≤M≤ N'2 .

进一步地,所述步骤(9)中加权的OR准则如下:Further, the weighted OR criterion in the step (9) is as follows:

其中,P'd,ts为终选从用户CR”'t对其所分配第s个频段的检测概率,P'fa,ts为终选从用户CR”'t对其所分配第s个频段的虚警概率;P'd,t为第t个重新选择的终选从用户CR”'t的平均检测概率,P'fa,t为第t个重新选择的终选从用户CR”'t的平均虚警概率;Q'd为协作检测后的全局检测概率,Q'fa为协作检测后的全局虚警概率;M'为重新选择的终选从用户的数目;ω't为重新选择的终选从用户CR”'t的加权系数。Among them, P' d, ts is the detection probability of the sth frequency band allocated by the final selected slave user CR"' t , and P' fa, ts is the sth frequency band allocated by the final selected slave user CR"' t False alarm probability; P' d,t is the average detection probability of the tth reselected final user CR"' t , P' fa,t is the tth reselected final user CR"' t The average false alarm probability of ; Q' d is the global detection probability after cooperative detection, Q' fa is the global false alarm probability after cooperative detection; M' is the number of reselected final slave users; ω' t is reselected The weighting coefficient of the final selection from the user CR”' t .

与现有技术相比,本发明的优点在于:各次用户分别将自身信噪比、置信度、业务需求以及对多个主用户的频谱检测结果发送给频谱感知融合中心,由频谱感知融合中心计算、筛选初选协作次用户,删除检测性能差和低信噪比的次用户,然后分配频段数量给初选协作次用户;计算所有初选次用户信噪比均方根值与各初选次用户的信噪比间的商值,根据各商值与信噪比预设阈值间的关系,选定参与协作的复选次用户,并调整复选次用户的检测概率;结合各复选次用户的信噪比和调整后的平均检测概率,重新返回选择复选次用户的步骤,得到最终参与协作的终选次用户,并以加权的OR准则协作的全局检测概率为频谱感知融合中心的最终检测结果,避免了低信噪比或性能差的次用户对整体检测性能的恶劣影响。Compared with the prior art, the present invention has the advantage that each user sends its own signal-to-noise ratio, confidence level, service requirements and spectrum detection results of multiple primary users to the spectrum sensing fusion center, and the spectrum sensing fusion center Calculate and screen primary cooperative secondary users, delete secondary users with poor detection performance and low signal-to-noise ratio, and then allocate the number of frequency bands to primary cooperative secondary users; calculate the root mean square value of the signal-to-noise ratio of all primary primary secondary users and The quotient value between the signal-to-noise ratios of the secondary users, according to the relationship between each quotient value and the preset threshold value of the signal-to-noise ratio, select the reselected secondary users participating in the collaboration, and adjust the detection probability of the reselected secondary users; The signal-to-noise ratio of the secondary user and the adjusted average detection probability return to the step of selecting multiple secondary users to obtain the final selected secondary users participating in the collaboration, and use the global detection probability of the weighted OR criterion cooperation as the spectrum sensing fusion center The final detection results can avoid the bad impact of low signal-to-noise ratio or poor performance secondary users on the overall detection performance.

附图说明Description of drawings

图1为本发明实施例中的认知网络结构示意图;FIG. 1 is a schematic diagram of a cognitive network structure in an embodiment of the present invention;

图2为本发明实施例的多频段协作认知频谱感知方法流程示意图;FIG. 2 is a schematic flow chart of a multi-band cooperative cognitive spectrum sensing method according to an embodiment of the present invention;

图3为本发明实施例的多频段协作认知频谱感知方法与传统AND准则协作频谱感知的仿真结果示意图。FIG. 3 is a schematic diagram of a simulation result of a multi-band cooperative cognitive spectrum sensing method and a traditional AND criterion cooperative spectrum sensing according to an embodiment of the present invention.

具体实施方式detailed description

以下结合附图实施例对本发明作进一步详细描述。The present invention will be further described in detail below in conjunction with the accompanying drawings and embodiments.

如图1所示,设定在认知网络中,主用户的数量为N1,次用户的数量为N2,频谱感知融合中心的数量为1,其中,主用户记为PUi,次用户记为CRj,频谱感知融合中心记为FC;1≤i≤N1,1≤j≤N2,N1≥2,N2≥2。As shown in Figure 1, in the cognitive network, the number of primary users is N 1 , the number of secondary users is N 2 , and the number of spectrum sensing fusion centers is 1, where the primary user is denoted as PU i , and the secondary user It is denoted as CR j , and the spectrum sensing fusion center is denoted as FC; 1≤i≤N 1 , 1≤j≤N 2 , N 1 ≥2, N 2 ≥2.

以下结合图1和图2,对本实施例中多频段协作认知频谱感知方法作出说明。其中,该多频段协作认知频谱感知方法,依次包括如下步骤:The following describes the multi-band cooperative cognitive spectrum sensing method in this embodiment with reference to FIG. 1 and FIG. 2 . Wherein, the multi-band cooperative cognitive spectrum sensing method includes the following steps in sequence:

(1)设定在认知网络中,主用户的数量为N1,次用户的数量为N2,频谱感知融合中心的数量为1,主用户分别独立地占用频谱中的各自频段;N2个次用户分别独立地获取自身信噪比SNRj以及对N1个主用户占用频段的频谱检测结果,并分别将获取的信噪比SNRj、频谱检测结果、检测置信度Pj和业务需求Rj发送至频谱感知融合中心,其中,(1) In the cognitive network, the number of primary users is N 1 , the number of secondary users is N 2 , the number of spectrum sensing fusion centers is 1, and the primary users independently occupy their respective frequency bands in the spectrum; N 2 Each secondary user independently obtains its own SNR j and the spectrum detection results of the frequency bands occupied by N 1 primary users, and respectively obtains the SNR j , spectrum detection results, detection confidence P j and service requirements R j is sent to the spectrum sensing fusion center, where,

主用户记为PUi,次用户记为CRj,频谱感知融合中心记为FC,业务需求Rj∈[0,N1],Pj∈[0,1],检测置信度Pd,ji为次用户CRj对主用户PUi的检测概率;1≤i≤N1,1≤j≤N2,N1≥2,N2≥2;The primary user is denoted as PU i , the secondary user is denoted as CR j , the spectrum sensing fusion center is denoted as FC, business requirements R j ∈ [0,N 1 ], P j ∈ [0,1], detection confidence P d,ji is the detection probability of secondary user CR j to primary user PU i ; 1≤i≤N 1 , 1≤j≤N 2 , N 1 ≥2, N 2 ≥2;

(2)频谱感知融合中心FC接收各次用户CRj发送来的信噪比SNRj、频谱检测结果、检测置信度Pj和业务需求Rj,并判断次用户的信噪比SNRj大于预设的信噪比筛选值SNRchose时,选择此信噪比对应的次用户为参与协作检测的初选次用户,并记初选次用户为CR't,并执行步骤(3),否则,选择具有最高信噪比的次用户所对应的频谱检测结果为频谱感知融合中心FC的最终检测结果;其中,(2) The spectrum sensing fusion center FC receives the signal-to-noise ratio SNR j , spectrum detection results, detection confidence P j and service requirements R j sent by each user CR j , and judges that the signal-to-noise ratio SNR j of the secondary user is greater than the expected When the selected signal-to-noise ratio screening value SNR chosen , select the secondary user corresponding to the signal-to-noise ratio as the primary secondary user participating in the collaborative detection, and record the primary secondary user as CR' t , and perform step (3), otherwise, Select the spectrum detection result corresponding to the secondary user with the highest SNR as the final detection result of the spectrum sensing fusion center FC; where,

初选次用户数量为N'2,初选次用户CR't对应的信噪比为SNR't、检测置信度为P't、业务需求为R't,1≤t≤N'2≤N2;频谱感知融合中心FC接收到的频谱检测结果数量为N1×N2个;The number of primary secondary users is N' 2 , the signal-to-noise ratio corresponding to the primary secondary user CR' t is SNR' t , the detection confidence is P' t , and the service requirement is R' t , 1≤t≤N' 2 ≤ N 2 ; the number of spectrum detection results received by the spectrum sensing fusion center FC is N 1 × N 2 ;

(3)频谱感知融合中心FC根据初选次用户CR't的置信度P't、业务需求R't,对初选次用户CR't分配需要检测的频段数量Ct(3) The spectrum sensing fusion center FC assigns the number of frequency bands C t to be detected to the primary secondary user CR' t according to the confidence P' t of the primary secondary user CR' t and the service demand R' t :

(3-1)根据各初选次用户CR't的置信度P't,分别对每一个初选次用户CR't的置信度P't进行归一化,得到每个初选次用户CR't的归一化置信度值 (3-1) According to the confidence degree P' t of each primary selected secondary user CR' t , respectively normalize the confidence P' t of each primary selected secondary user CR' t to obtain each primary selected secondary user Normalized confidence value of CR't

(3-2)根据步骤(3-1)所得每个初选次用户CR't对应的归一化置信度值计算频谱感知融合中心FC分配初选次用户CR't需要检测的频段数量Ct(3-2) The normalized confidence value corresponding to each primary user CR' t obtained according to step (3-1) Calculate the number of frequency bands C t that need to be detected when the spectrum sensing fusion center FC assigns the primary user CR' t to be detected:

(4)频谱感知融合中心FC根据参与协作检测的初选次用户CR't的信噪比SNR't,计算所有初选次用户的信噪比均方根值并令信噪比SNR't=γt,其中,信噪比均方根值的计算如下:(4) The spectrum sensing fusion center FC calculates the root mean square value of the signal-to-noise ratio of all primary secondary users according to the signal-to-noise ratio SNR' t of primary secondary users CR' t participating in cooperative detection And let the signal-to-noise ratio SNR' t = γ t , where the root mean square value of the signal-to-noise ratio is calculated as follows:

(5)频谱感知融合中心FC分别依次计算所有初选次用户的信噪比均方根值与各初选次用户CR't的信噪比SNR't之间的商值ηt,其中,(5) The spectrum sensing fusion center FC calculates the root mean square value of the signal-to-noise ratio of all primary secondary users in turn The quotient η t between the signal-to-noise ratio SNR' t of each primary secondary user CR' t , where

(6)频谱感知融合中心FC计算、获取信噪比预设阈值λ和信噪比最佳阈值λoptimal,并分别根据各信噪比商值ηt与信噪比预设阈值λ之间的大小关系,选定参与协作的复选次用户CR”k,复选次用户CR”k的信噪比为SNR”k,其中,(6) The spectrum sensing fusion center FC calculates and obtains the SNR preset threshold λ and the SNR optimal threshold λ optimal , and respectively calculates and obtains the SNR quotient value η t and the SNR preset threshold λ according to The size relationship, select the secondary user CR” k to participate in the collaboration, and the signal-to-noise ratio of the secondary user CR” k is SNR” k , where,

(6-1)频谱感知融合中心FC根据接收的N'2个初选从用户对应的信噪比集合{SNR't},获取初选从用户信噪比集合{SNR't}中的信噪比最大值,记该信噪比最大值为SNR'max;(6-1) The spectrum sensing fusion center FC acquires the signal-to-noise ratio set {SNR' t } of the primary-selected slave users according to the received SNR sets {SNR' t } corresponding to the N' 2 primary-selected slave users. The maximum value of the noise ratio, record the maximum value of the signal-to-noise ratio as SNR'max;

(6-2)以获取的信噪比最大值SNR'max为参考,并将N'2个初选从用户CR't的信噪比SNR't分别与信噪比最大值SNR'max作商处理,计算得到各初选从用户信噪比SNR't所对应的初始阈值λt,其中,(6-2) Take the obtained SNR maximum value SNR'max as a reference, and make the SNR' t of the N' 2 primary users CR' t with the SNR maximum value SNR'max respectively quotient processing, and calculate the initial threshold λ t corresponding to each primary user SNR' t , where,

λt=|SNR't/SNR'max|,t=1,2,…,N'2,N'2≤N2λ t = |SNR' t /SNR' max |, t=1,2,...,N' 2 , N' 2 ≤ N 2 ;

(6-3)根据各初选从用户CR't的归一化置信度值和信噪比商值ηt,计算各初选从用户CR't的联合筛选参数值ξt,并根据联合筛选参数值ξt,选取参与协作的复选从用户CR”k,其中,复选从用户CR”k的数量为M,t=1,2,…,N'2,k=1,2,…,M,M≤N'2(6-3) According to the normalized confidence value of each primary selection from the user CR' t and the signal-to-noise ratio quotient η t , calculate the joint screening parameter value ξ t of each primary secondary user CR' t , and select the secondary user CR" k participating in the collaboration according to the joint screening parameter value ξ t , where the complex The number of selected users CR” k is M, t=1,2,...,N' 2 , k=1,2,...,M, M≤N' 2 :

若联合筛选参数值ξt位于预设数值区间范围[ξab]内,即ξa≤ξt≤ξb时,则选取该联合筛选参数值ξt对应的初选从用户为复选从用户,并参与协作检测;否则,该初选从用户不予选取;If the joint screening parameter value ξ t is within the preset value range [ξ a , ξ b ], that is, when ξ a ≤ ξ t ≤ ξ b , then select the primary user corresponding to the joint screening parameter value ξ t as complex Select a secondary user and participate in collaborative detection; otherwise, the primary secondary user will not be selected;

(6-4)根据步骤(6-3)中的信噪比预设阈值λ,获取M个复选从用户CR”k分别在OR准则和AND准则下的协作检测性能曲线,其中,(6-4) According to the signal-to-noise ratio preset threshold λ in step (6-3), obtain the cooperative detection performance curves of M re-selected slave users CR" k under the OR criterion and the AND criterion respectively, wherein,

OR准则: OR criterion:

AND准则:k=1,2,…,M,M≤N'2≤N2AND criteria: k=1,2,...,M, M≤N' 2 ≤N 2 ;

其中,Pd,k为第k个复选从用户CR”k的平均检测概率,Pfa,k为第k个复选从用户CR”k的平均虚警概率;Pd,ks为复选从用户CR”k对其所分配第s个频段的检测概率,Pfa,ks为复选从用户CR”k对其所分配第s个频段的虚警概率;Qd为协作检测后的全局检测概率,Qfa为协作检测后的全局虚警概率;ωk表示信噪比CR”k的权重系数,SNR”k是第k个复选从用户CR”k的信噪比,SNR”max表示M个复选从用户的信噪比最大值,SNR”min表示M个复选从用户的信噪比最小值;Among them, P d, k is the average detection probability of the k-th re-selected slave user CR” k , P fa, k is the average false alarm probability of the k-th re-selected slave user CR” k ; P d, ks is the re-selected The detection probability of the s-th frequency band allocated from user CR” k , P fa,ks is the false alarm probability of the s-th frequency band allocated from user CR” k to it; Q d is the global Detection probability, Q fa is the global false alarm probability after collaborative detection; ω k represents the weight coefficient of signal-to-noise ratio CR" k , SNR" k is the signal-to-noise ratio of the kth re-selected user CR" k , SNR" max Indicates the maximum value of the SNR of the M secondary users, and SNR" min represents the minimum value of the SNR of the M secondary users;

(6-5)根据OR准则和AND准则下的协作检测性能曲线,分别得到在OR准则和AND准则下的最大检测概率Q(OR,d)-max、Q(AND,d)-max,得到Q(OR,d)-max和Q(AND,d)-max的最大值Qd-max,并以该最佳检测性能值Qd-max所对应的信噪比预设阈值为信噪比最佳阈值,记信噪比最佳阈值为λoptimal;其中,Qd-max=max(Q(OR,d)-max,Q(AND,d)-max);(6-5) According to the collaborative detection performance curves under the OR criterion and the AND criterion, the maximum detection probabilities Q (OR, d)-max and Q (AND, d)-max under the OR criterion and the AND criterion are respectively obtained, and The maximum value of Q (OR, d)-max and Q (AND, d)-max Q d-max , and the signal-to-noise ratio preset threshold corresponding to the best detection performance value Q d-max is the signal-to-noise ratio The optimal threshold, the optimal threshold of SNR is λ optimal ; where, Q d-max = max(Q (OR, d)-max , Q (AND, d)-max );

(7)根据获取的信噪比最佳阈值λoptimal,得到该信噪比最佳阈值λoptimal对应的复选从用户CR”,获取该复选从用户CR”的调整因子α以及其他M-1个复选从用户CR”k的调整因子αk,并分别根据调整因子α、αk对应调整复选从用户CR”、CR”k的平均虚警概率,复选从用户CR”调整后的平均虚警概率记为Pfa,复选从用户CR”k调整后的平均虚警概率记为Pfa,k;其中,(7) According to the obtained SNR optimal threshold λ optimal , obtain the multiple-selection slave user CR” corresponding to the SNR optimal threshold λ optimal , obtain the adjustment factor α of the multiple-selection slave user CR” and other M- The adjustment factor α k of one reselected slave user CR” k , and adjust the average false alarm probability of reselected slave users CR” and CR” k according to the adjustment factors α and α k respectively. After the reselected slave user CR” is adjusted The average false alarm probability is denoted as P fa , and the average false alarm probability adjusted from user CR” k is denoted as P fa,k ; where,

Pfa,k=αk·Pfa,k=1,2,…,M-1;P fa,kk P fa ,k=1,2,...,M-1;

其中,αk为复选从用户CR”k的调整因子,用来根据复选从用户CR”k自身的信噪比SNR”k实现对其平均虚警概率大小的调整;SNR”k为第k个复选从用户CR”k的信噪比;Among them, α k is the adjustment factor of the secondary user CR" k , which is used to adjust the average false alarm probability according to the signal-to-noise ratio SNR" k of the secondary user CR" k itself; SNR of k multiple choices from user CR"k;

(8)根据步骤(7)中获取的M个复选从用户的调整因子αk以及对应调整后的平均虚警概率Pfa,k,计算复选从用户CR”k调整后的判决阀值λ”k和平均检测概率Pd,k,其中,(8) According to the adjustment factor α k of the M secondary users obtained in step (7) and the corresponding adjusted average false alarm probability P fa,k , calculate the adjusted decision threshold of the secondary user CR” k λ” k and average detection probability P d,k , where,

其中,k=1,2,…,M,M≤N'2;n为采样点数;in, k=1,2,...,M, M≤N'2; n is the number of sampling points;

(9)根据步骤(8)中M个复选从用户的信噪比SNR”k以及得到的调整后的平均检测概率Pd,k,返回步骤(6),重新在M个复选从用户中选择,得到参与协作的T个终选从用户CR”'t,并以加权的OR准则协作后的全局检测概率为频谱感知融合中心FC的最终检测结果,其中,加权的OR准则如下:(9) According to the signal-to-noise ratio SNR” k of the M reselected slave users in step (8) and the adjusted average detection probability P d,k obtained, return to step (6), and re-select the M reselected slave users Select among the T final users CR”' t participating in the collaboration, and use the weighted OR criterion to obtain the global detection probability after collaboration as the final detection result of the spectrum sensing fusion center FC, where the weighted OR criterion is as follows:

其中,P'd,ts为终选从用户CR”'t对其所分配第s个频段的检测概率,P'fa,ts为终选从用户CR”'t对其所分配第s个频段的虚警概率;P'd,t为第t个重新选择的终选从用户CR”'t的平均检测概率,P'fa,t为第t个重新选择的终选从用户CR”'t的平均虚警概率;Q'd为协作检测后的全局检测概率,Q'fa为协作检测后的全局虚警概率;M'为重新选择的终选从用户的数目;ω't为重新选择的终选从用户CR”'t的加权系数,1≤t≤T≤M≤N'2Among them, P' d, ts is the detection probability of the sth frequency band allocated by the final selected slave user CR"' t , and P' fa, ts is the sth frequency band allocated by the final selected slave user CR"' t False alarm probability; P' d,t is the average detection probability of the tth reselected final user CR"' t , P' fa,t is the tth reselected final user CR"' t The average false alarm probability of ; Q' d is the global detection probability after cooperative detection, Q' fa is the global false alarm probability after cooperative detection; M' is the number of reselected final slave users; ω' t is reselected The weighting coefficient of the final selection from the user CR"' t , 1≤t≤T≤M≤N' 2 .

图3给出了本发明中多频段协作认知频谱感知方法的仿真结果示意图,并同时对传统AND协作检测方法进行了仿真。其中,在传统的AND协作检测方法中,所有的次用户均参与协作频谱感知检测。仿真条件如下:Fig. 3 shows a schematic diagram of the simulation results of the multi-band cooperative cognitive spectrum sensing method in the present invention, and simulates the traditional AND cooperative detection method at the same time. Wherein, in the traditional AND cooperative detection method, all secondary users participate in cooperative spectrum sensing detection. The simulation conditions are as follows:

设定在认知无线电网络中,主用户PU的数量N1分别为2、3、5和7,从用户CR的数量N2逐渐从2个增加至11个;从用户的信噪比分别为-5dB、-8dB、-10dB、-13dB、-16dB、-17dB、-19dB、-23dB、-25dB和-27dB,从用户均采用能量检测。It is assumed that in the cognitive radio network, the number N 1 of primary users PU is 2, 3, 5 and 7 respectively, and the number N 2 of slave users CR gradually increases from 2 to 11; the signal-to-noise ratios of slave users are respectively -5dB, -8dB, -10dB, -13dB, -16dB, -17dB, -19dB, -23dB, -25dB and -27dB, all adopt energy detection from the user.

由图3可以看出,在主用户数量和从用户数量一定的条件下,本发明中的全局检测概率要大于传统AND协作检测的检测概率,这表明本发明中的协作频谱感知方法具有更好的检测性能;同时,也可以发现,在协作检测方法和主用户数量一定的条件下,随着从用户数量的增加,协作检测后的全局检测概率逐渐增大;在协作检测方法和从用户数量一定的条件下,随着主用户数量的增加,协作检测后的全局检测概率则逐渐减小。可知,相对传统的AND协作检测方法,本发明实施例中的多频段认知协作频谱感知方法因避免了低信噪比或检测性能差的次用户的恶劣影响,整体的协作检测性能得到很大的提高,具有更好的检测性能。As can be seen from Figure 3, under the condition that the number of master users and the number of slave users are certain, the global detection probability in the present invention is greater than the detection probability of traditional AND cooperative detection, which shows that the cooperative spectrum sensing method in the present invention has better At the same time, it can also be found that under the conditions of the cooperative detection method and the number of master users, with the increase of the number of slave users, the global detection probability after cooperative detection gradually increases; in the cooperative detection method and the number of slave users Under certain conditions, with the increase of the number of primary users, the global detection probability after cooperative detection gradually decreases. It can be seen that, compared with the traditional AND cooperative detection method, the multi-band cognitive cooperative spectrum sensing method in the embodiment of the present invention avoids the bad influence of secondary users with low signal-to-noise ratio or poor detection performance, and the overall cooperative detection performance is greatly improved. The improvement has better detection performance.

Claims (1)

1.多频段协作认知频谱感知方法,其特征在于,依次包括如下步骤:1. The multi-band cooperative cognitive spectrum sensing method is characterized in that, comprising the following steps in sequence: (1)设定在认知网络中,主用户的数量为N1,次用户的数量为N2,频谱感知融合中心的数量为1,主用户分别独立地占用频谱中的各自频段;N2个次用户分别独立地获取自身信噪比SNRj以及对N1个主用户占用频段的频谱检测结果,并分别将获取的信噪比SNRj、频谱检测结果、检测置信度Pj和业务需求Rj发送至频谱感知融合中心,其中,(1) In the cognitive network, the number of primary users is N 1 , the number of secondary users is N 2 , the number of spectrum sensing fusion centers is 1, and the primary users independently occupy their respective frequency bands in the spectrum; N 2 Each secondary user independently obtains its own SNR j and the spectrum detection results of the frequency bands occupied by N 1 primary users, and respectively obtains the SNR j , spectrum detection results, detection confidence P j and service requirements R j is sent to the spectrum sensing fusion center, where, 主用户记为PUi,次用户记为CRj,频谱感知融合中心记为FC,业务需求Rj∈[0,N1],Pj∈[0,1],检测置信度Pd,ji为次用户CRj对主用户PUi的检测概率;1≤i≤N1,1≤j≤N2,N1≥2,N2≥2;The primary user is denoted as PU i , the secondary user is denoted as CR j , the spectrum sensing fusion center is denoted as FC, business requirements R j ∈ [0,N 1 ], P j ∈ [0,1], detection confidence P d,ji is the detection probability of secondary user CR j to primary user PU i ; 1≤i≤N 1 , 1≤j≤N 2 , N 1 ≥2, N 2 ≥2; (2)频谱感知融合中心FC接收各次用户CRj发送来的信噪比SNRj、频谱检测结果、检测置信度Pj和业务需求Rj,并判断次用户的信噪比SNRj大于预设的信噪比筛选值SNRchose时,选择此信噪比对应的次用户为参与协作检测的初选次用户,并记初选次用户为CR't,并执行步骤(3),否则,选择具有最高信噪比的次用户所对应的频谱检测结果为频谱感知融合中心FC的最终检测结果;其中,(2) The spectrum sensing fusion center FC receives the signal-to-noise ratio SNR j , spectrum detection results, detection confidence P j and service requirements R j sent by each user CR j , and judges that the signal-to-noise ratio SNR j of the secondary user is greater than the expected When the selected signal-to-noise ratio screening value SNR chosen , select the secondary user corresponding to the signal-to-noise ratio as the primary secondary user participating in the collaborative detection, and record the primary secondary user as CR' t , and perform step (3), otherwise, Select the spectrum detection result corresponding to the secondary user with the highest SNR as the final detection result of the spectrum sensing fusion center FC; where, 初选次用户数量为N'2,初选次用户CR't对应的信噪比为SNR't、检测置信度为P't、业务需求为R't,1≤t≤N'2≤N2;频谱感知融合中心FC接收到的频谱检测结果数量为N1×N2个;The number of primary secondary users is N' 2 , the signal-to-noise ratio corresponding to the primary secondary user CR' t is SNR' t , the detection confidence is P' t , and the service requirement is R' t , 1≤t≤N' 2 ≤ N 2 ; the number of spectrum detection results received by the spectrum sensing fusion center FC is N 1 × N 2 ; (3)频谱感知融合中心FC根据初选次用户CR't的置信度P't、业务需求R't,对初选次用户CR't分配需要检测的频段数量Ct(3) The spectrum sensing fusion center FC assigns the number of frequency bands C t to be detected to the primary secondary user CR' t according to the confidence P' t of the primary secondary user CR' t and the service demand R' t : (3-1)根据各初选次用户CR't的置信度P't,分别对每一个初选次用户CR't的置信度P't进行归一化,得到每个初选次用户CR't的归一化置信度值 (3-1) According to the confidence degree P' t of each primary selected secondary user CR' t , respectively normalize the confidence P' t of each primary selected secondary user CR' t to obtain each primary selected secondary user Normalized confidence value of CR't PP ′′ tt ‾‾ == PP ′′ tt ΣΣ tt == 11 NN ′′ 22 PP ′′ tt ,, 11 ≤≤ tt ≤≤ NN ′′ 22 ;; (3-2)根据步骤(3-1)所得每个初选次用户CR't对应的归一化置信度值计算频谱感知融合中心FC分配初选次用户CR't需要检测的频段数量Ct(3-2) The normalized confidence value corresponding to each primary user CR' t obtained according to step (3-1) Calculate the number of frequency bands C t that need to be detected when the spectrum sensing fusion center FC assigns the primary user CR' t to be detected: CC tt == PP ′′ tt ‾‾ ·&Center Dot; NN ′′ 22 ,, 11 ≤≤ tt ≤≤ NN ′′ 22 ;; (4)频谱感知融合中心FC根据参与协作检测的初选次用户CR't的信噪比SNR't,计算所有初选次用户的信噪比均方根值并令信噪比SNR't=γt,其中,信噪比均方根值的计算如下:(4) The spectrum sensing fusion center FC calculates the root mean square value of the signal-to-noise ratio of all primary secondary users according to the signal-to-noise ratio SNR' t of primary secondary users CR' t participating in cooperative detection And let the signal-to-noise ratio SNR' t = γ t , where the root mean square value of the signal-to-noise ratio is calculated as follows: γγ ‾‾ == 11 NN ′′ 22 ΣΣ tt == 11 NN ′′ 22 (( SNRSNR ′′ tt )) 22 ,, NN ′′ 22 ≤≤ NN 22 ;; (5)频谱感知融合中心FC分别依次计算所有初选次用户的信噪比均方根值与各初选次用户CR't的信噪比SNR't之间的商值ηt,其中,(5) The spectrum sensing fusion center FC calculates the root mean square value of the signal-to-noise ratio of all primary secondary users in turn The quotient η t between the signal-to-noise ratio SNR' t of each primary secondary user CR' t , where ηη tt == || γγ ‾‾ // γγ tt || ,, tt == 11 ,, 22 ,, ...... ,, NN ′′ 22 ,, NN ′′ 22 ≤≤ NN 22 ;; (6)频谱感知融合中心FC计算、获取信噪比预设阈值λ和信噪比最佳阈值λoptimal,并分别根据各信噪比商值ηt与信噪比预设阈值λ之间的大小关系,选定参与协作的复选次用户CR”k,复选次用户CR”k的信噪比为SNR”k,其中,(6) The spectrum sensing fusion center FC calculates and obtains the SNR preset threshold λ and the SNR optimal threshold λ optimal , and respectively calculates and obtains the SNR quotient value η t and the SNR preset threshold λ according to The size relationship, select the secondary user CR” k to participate in the collaboration, and the signal-to-noise ratio of the secondary user CR” k is SNR” k , where, (6-1)频谱感知融合中心FC根据接收的N'2个初选从用户对应的信噪比集合{SNR't},获取初选从用户信噪比集合{SNR't}中的信噪比最大值,记该信噪比最大值为SNR'max;(6-1) The spectrum sensing fusion center FC acquires the signal-to-noise ratio set {SNR' t } of the primary-selected slave users according to the received SNR sets {SNR' t } corresponding to the N' 2 primary-selected slave users. The maximum value of the noise ratio, record the maximum value of the signal-to-noise ratio as SNR'max; (6-2)以获取的信噪比最大值SNR'max为参考,并将N'2个初选从用户CR't的信噪比SNR't分别与信噪比最大值SNR'max作商处理,计算得到各初选从用户信噪比SNR't所对应的初始阈值λt,其中,(6-2) Take the obtained SNR maximum value SNR'max as a reference, and make the SNR' t of the N' 2 primary users CR' t with the SNR maximum value SNR'max respectively quotient processing, and calculate the initial threshold λ t corresponding to each primary user SNR' t , where, λt=|SNR't/SNR'max|,t=1,2,…,N'2,N'2≤N2λ t = |SNR' t /SNR' max |, t=1,2,...,N' 2 , N' 2 ≤ N 2 ; (6-3)根据各初选从用户CR't的归一化置信度值和信噪比商值ηt,计算各初选从用户CR't的联合筛选参数值ξt,并根据联合筛选参数值ξt,选取参与协作的复选从用户CR”k,其中,复选从用户CR”k的数量为M,t=1,2,…,N′2,k=1,2,…,M,M≤N'2(6-3) According to the normalized confidence value of each primary selection from the user CR' t and the signal-to-noise ratio quotient η t , calculate the joint screening parameter value ξ t of each primary secondary user CR' t , and select the secondary user CR" k participating in the collaboration according to the joint screening parameter value ξ t , where the complex The number of selected users CR” k is M, t=1, 2,..., N' 2 , k=1, 2,..., M, M≤N' 2 : 若联合筛选参数值ξt位于预设数值区间范围[ξab]内,即ξa≤ξt≤ξb时,则选取该联合筛选参数值ξt对应的初选从用户为复选从用户,并参与协作检测;否则,该初选从用户不予选取;If the joint screening parameter value ξ t is within the preset value range [ξ a , ξ b ], that is, when ξ a ≤ ξ t ≤ ξ b , then select the primary user corresponding to the joint screening parameter value ξ t as complex Select a secondary user and participate in collaborative detection; otherwise, the primary secondary user will not be selected; (6-4)根据步骤(6-3)中的信噪比预设阈值λ,获取M个复选从用户CR”k分别在OR准则和AND准则下的协作检测性能曲线,其中,(6-4) According to the signal-to-noise ratio preset threshold λ in step (6-3), obtain the cooperative detection performance curves of M re-selected slave users CR" k under the OR criterion and the AND criterion respectively, wherein, OR准则: OR criterion: AND准则: AND criteria: PP dd ,, kk == ΣΣ kk NN ′′ 22 ΣΣ sthe s == 11 CC tt PP dd ,, kk sthe s NN ′′ 22 ·&Center Dot; CC tt ,, PP ff aa ,, kk == ΣΣ kk NN ′′ 22 ΣΣ sthe s == 11 CC tt PP ff aa ,, kk sthe s NN ′′ 22 ·&Center Dot; CC tt ;; 其中,Pd,k为第k个复选从用户CR”k的平均检测概率,Pfa,k为第k个复选从用户CR”k的平均虚警概率;Pd,ks为复选从用户CR”k对其所分配第s个频段的检测概率,Pfa,ks为复选从用户CR”k对其所分配第s个频段的虚警概率;Qd为协作检测后的全局检测概率,Qfa为协作检测后的全局虚警概率;ωk表示信噪比CR”k的权重系数,SNR”k是第k个复选从用户CR”k的信噪比,SNR”max表示M个复选从用户的信噪比最大值,SNR”min表示M个复选从用户的信噪比最小值;Among them, P d, k is the average detection probability of the k-th re-selected slave user CR” k , P fa, k is the average false alarm probability of the k-th re-selected slave user CR” k ; P d, ks is the re-selected The detection probability of the s-th frequency band allocated from user CR” k , P fa,ks is the false alarm probability of the s-th frequency band allocated from user CR” k to it; Q d is the global Detection probability, Q fa is the global false alarm probability after collaborative detection; ω k represents the weight coefficient of signal-to-noise ratio CR" k , SNR" k is the signal-to-noise ratio of the kth re-selected user CR" k , SNR" max Indicates the maximum value of the SNR of the M secondary users, and SNR" min represents the minimum value of the SNR of the M secondary users; (6-5)根据OR准则和AND准则下的协作检测性能曲线,分别得到在OR准则和AND准则下的最大检测概率Q(OR,d)-max、Q(AND,d)-max,得到Q(OR,d)-max和Q(AND,d)-max的最大值Qd-max,并以该最佳检测性能值Qd-max所对应的信噪比预设阈值为信噪比最佳阈值,记信噪比最佳阈值为λoptimal;其中,Qd-max=max(Q(OR,d)-max,Q(AND,d)-max);(6-5) According to the collaborative detection performance curves under the OR criterion and the AND criterion, the maximum detection probabilities Q (OR, d)-max and Q (AND, d)-max under the OR criterion and the AND criterion are respectively obtained, and The maximum value of Q (OR, d)-max and Q (AND, d)-max Q d-max , and the signal-to-noise ratio preset threshold corresponding to the best detection performance value Q d-max is the signal-to-noise ratio The optimal threshold, the optimal threshold of SNR is λ optimal ; where, Q d-max = max(Q (OR, d)-max , Q (AND, d)-max ); (7)根据获取的信噪比最佳阈值λoptimal,得到该信噪比最佳阈值λoptimal对应的复选从用户CR”,获取该复选从用户CR”的调整因子α以及其他M-1个复选从用户CR”k的调整因子αk,并分别根据调整因子α、αk对应调整复选从用户CR”、CR”k的平均虚警概率,复选从用户CR”调整后的平均虚警概率记为Pfa,复选从用户CR”k调整后的平均虚警概率记为Pfa,k;其中,(7) According to the obtained SNR optimal threshold λ optimal , obtain the multiple-selection slave user CR” corresponding to the SNR optimal threshold λ optimal , obtain the adjustment factor α of the multiple-selection slave user CR” and other M- The adjustment factor α k of one reselected slave user CR” k , and adjust the average false alarm probability of reselected slave users CR” and CR” k according to the adjustment factors α and α k respectively. After the reselected slave user CR” is adjusted The average false alarm probability is denoted as P fa , and the average false alarm probability adjusted from user CR” k is denoted as P fa,k ; where, Pfa,k=αk·Pfa,k=1,2,…,M-1;P fa,kk P fa ,k=1,2,...,M-1; αα kk == 11 ++ SNRSNR ′′ ‾‾ -- SNRSNR ′′ ′′ kk SNRSNR ′′ ′′ ‾‾ ,, kk == 11 ,, 22 ,, ...... ,, Mm -- 11 ;; SNRSNR ′′ ′′ ‾‾ == ΣΣ kk == 11 Mm (( SNRSNR ′′ ′′ kk )) 22 Mm ,, Mm ≤≤ NN ′′ 22 ;; 其中,αk为复选从用户CR”k的调整因子,用来根据复选从用户CR”k自身的信噪比SNR”k实现对其平均虚警概率大小的调整;SNR”k为第k个复选从用户CR”k的信噪比;Among them, α k is the adjustment factor of the secondary user CR" k , which is used to adjust the average false alarm probability according to the signal-to-noise ratio SNR" k of the secondary user CR" k itself; SNR of k multiple choices from user CR"k; (8)根据步骤(7)中获取的M个复选从用户的调整因子αk以及对应调整后的平均虚警概率Pfa,k,计算复选从用户CR”k调整后的判决阀值λ”k和平均检测概率Pd,k,其中,(8) According to the adjustment factor α k of the M secondary users obtained in step (7) and the corresponding adjusted average false alarm probability P fa,k , calculate the adjusted decision threshold of the secondary user CR” k λ” k and average detection probability P d,k , where, λλ == σσ ww 22 [[ 22 nno QQ -- 11 (( PP ff aa ,, kk )) ++ nno ]] == σσ ww 22 [[ 22 nno QQ -- 11 (( δδ ·· PP ff aa )) ++ nno ]] == σσ ww 22 [[ 22 nno QQ -- 11 (( (( 11 ++ SNRSNR ′′ ′′ ‾‾ -- SNRSNR ′′ ′′ kk SNRSNR ′′ ′′ ‾‾ )) ·· PP ff aa )) ++ nno ]] ;; PP dd ,, kk == QQ [[ QQ -- 11 (( PP ff aa ,, kk )) -- nno ·· SNRSNR ′′ ′′ kk ]] ;; nno == 22 [[ QQ -- 11 (( PP ff aa ,, kk )) -- QQ -- 11 (( PP ff aa )) 11 ++ 22 SNRSNR ′′ ′′ kk ]] 22 ·· (( SNRSNR ′′ ′′ kk )) -- 22 ;; 其中,n为采样点数;in, n is the number of sampling points; (9)根据步骤(8)中M个复选从用户的信噪比SNR”k以及得到的调整后的平均检测概率Pd,k,返回步骤(6),重新在M个复选从用户中选择,得到参与协作的T个终选从用户CR”'t,并以加权的OR准则协作后的全局检测概率为频谱感知融合中心FC的最终检测结果,其中1≤t≤T≤M≤N'2;其中:(9) According to the signal-to-noise ratio SNR” k of the M reselected slave users in step (8) and the adjusted average detection probability P d,k obtained, return to step (6), and re-select the M reselected slave users Select among T final users CR"' t participating in the collaboration, and use the weighted OR criterion to obtain the global detection probability after collaboration as the final detection result of the spectrum sensing fusion center FC, where 1≤t≤T≤M≤ N'2; where: 所述加权的OR准则如下:The weighted OR criterion is as follows: QQ ′′ dd == 11 -- ΠΠ tt == 11 Mm ′′ ωω ′′ tt (( 11 -- PP ′′ dd ,, tt )) ,, QQ ′′ ff aa == 11 -- ΠΠ tt == 11 Mm ′′ ωω ′′ tt (( 11 -- PP ′′ ff ,, tt )) ;; ωω ′′ tt == PP ′′ dd ,, tt ΣΣ tt == 11 Mm ′′ PP ′′ dd ,, tt ,, PP ′′ dd ,, tt == ΣΣ tt Mm ′′ ΣΣ sthe s == 11 CC tt PP ′′ dd ,, tt sthe s Mm ′′ ·&Center Dot; CC tt ,, PP ′′ ff aa ,, tt == ΣΣ tt Mm ′′ ΣΣ sthe s == 11 CC tt PP ′′ ff aa ,, tt sthe s Mm ′′ ·· CC tt ,, tt == 11 ,, 22 ,, ...... ,, Mm ′′ ,, Mm ′′ ≤≤ Mm ;; 其中,P'd,ts为终选从用户CR”'t对其所分配第s个频段的检测概率,P'fa,ts为终选从用户CR”'t对其所分配第s个频段的虚警概率;P'd,t为第t个重新选择的终选从用户CR”'t的平均检测概率,P'fa,t为第t个重新选择的终选从用户CR”'t的平均虚警概率;Q'd为协作检测后的全局检测概率,Q'fa为协作检测后的全局虚警概率;M'为重新选择的终选从用户的数目;ω't为重新选择的终选从用户CR”'t的加权系数。Among them, P' d, ts is the detection probability of the sth frequency band allocated by the final selected slave user CR"' t , and P' fa, ts is the sth frequency band allocated by the final selected slave user CR"' t False alarm probability; P' d,t is the average detection probability of the tth reselected final user CR"' t , P' fa,t is the tth reselected final user CR"' t The average false alarm probability of ; Q' d is the global detection probability after cooperative detection, Q' fa is the global false alarm probability after cooperative detection; M' is the number of reselected final slave users; ω' t is reselected The weighting coefficient of the final selection from the user CR”' t .
CN201510595891.4A 2015-09-17 2015-09-17 Multiband cooperative cognitive frequency spectrum sensing method Active CN105141385B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510595891.4A CN105141385B (en) 2015-09-17 2015-09-17 Multiband cooperative cognitive frequency spectrum sensing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510595891.4A CN105141385B (en) 2015-09-17 2015-09-17 Multiband cooperative cognitive frequency spectrum sensing method

Publications (2)

Publication Number Publication Date
CN105141385A CN105141385A (en) 2015-12-09
CN105141385B true CN105141385B (en) 2017-07-11

Family

ID=54726614

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510595891.4A Active CN105141385B (en) 2015-09-17 2015-09-17 Multiband cooperative cognitive frequency spectrum sensing method

Country Status (1)

Country Link
CN (1) CN105141385B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105578483B (en) * 2016-02-15 2022-04-29 天津大学 High-energy-efficiency cooperative spectrum sensing method based on wireless sensor network
CN106341791B (en) * 2016-09-23 2019-11-12 宁波大学 A Dynamic Decision Method for Cognitive Radio Spectrum Sensing

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101908934A (en) * 2009-06-08 2010-12-08 上海贝尔股份有限公司 Frequency spectrum sensing method, device and equipment of frequency spectrum sharing system
CN103036622A (en) * 2011-09-29 2013-04-10 北京邮电大学 Cognitive radio spectrum detection method and device based on self-adaption double thresholds
US20130090145A1 (en) * 2010-06-09 2013-04-11 Industry-Academic Cooperation Foundation, Yonsei University Method for controlling random access for the efficient sensing of the cooperative spectrum in a cognitive radio-based frequency resource sharing system
CN103384174A (en) * 2013-05-10 2013-11-06 江苏科技大学 Method based on cooperation of multiple users and multiple antennas for optimizing spectrum sensing detection probability
CN104378169A (en) * 2014-12-08 2015-02-25 电子科技大学 Collaborative spectrum sensing system and collaborative spectrum sensing method based on power set confidence coefficients
CN104394543A (en) * 2014-12-08 2015-03-04 电子科技大学 Joint frequency spectrum sensing method based on Adaboost algorithm
US20150109947A1 (en) * 2013-10-18 2015-04-23 Research & Business Foundation Sungkyunkwan University Spectrum sensing apparatus and method for cooperative cognitive radio network in non-gaussian noise environment, and fusion center apparatus and cooperative cognitive radio system using the same
CN104780007A (en) * 2015-04-03 2015-07-15 南京邮电大学 Coalitional game based multi-user collaborative spectrum sensing method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101908934A (en) * 2009-06-08 2010-12-08 上海贝尔股份有限公司 Frequency spectrum sensing method, device and equipment of frequency spectrum sharing system
US20130090145A1 (en) * 2010-06-09 2013-04-11 Industry-Academic Cooperation Foundation, Yonsei University Method for controlling random access for the efficient sensing of the cooperative spectrum in a cognitive radio-based frequency resource sharing system
CN103036622A (en) * 2011-09-29 2013-04-10 北京邮电大学 Cognitive radio spectrum detection method and device based on self-adaption double thresholds
CN103384174A (en) * 2013-05-10 2013-11-06 江苏科技大学 Method based on cooperation of multiple users and multiple antennas for optimizing spectrum sensing detection probability
US20150109947A1 (en) * 2013-10-18 2015-04-23 Research & Business Foundation Sungkyunkwan University Spectrum sensing apparatus and method for cooperative cognitive radio network in non-gaussian noise environment, and fusion center apparatus and cooperative cognitive radio system using the same
CN104378169A (en) * 2014-12-08 2015-02-25 电子科技大学 Collaborative spectrum sensing system and collaborative spectrum sensing method based on power set confidence coefficients
CN104394543A (en) * 2014-12-08 2015-03-04 电子科技大学 Joint frequency spectrum sensing method based on Adaboost algorithm
CN104780007A (en) * 2015-04-03 2015-07-15 南京邮电大学 Coalitional game based multi-user collaborative spectrum sensing method

Also Published As

Publication number Publication date
CN105141385A (en) 2015-12-09

Similar Documents

Publication Publication Date Title
CN103220052B (en) A kind of method detecting frequency spectrum cavity-pocket in cognitive radio
RU2508615C2 (en) Improved method and apparatus for combined transmission dispatching in wireless network
JP2010035170A (en) Spectrum sensing method, energy detecting method and apparatus
CN103338458B (en) A kind of cooperative frequency spectrum sensing method for cognitive radio system
CN101951307A (en) Method for selecting cell cooperation set under CoMP
CN104980389A (en) Differential new user selection access method applicable to non-orthogonal multi-access system
CN103260166A (en) Distributed cognition wireless network spectrum allocation method based on coalition games
US9179484B2 (en) Method and device for allocating group identifier
CN101729164B (en) Wireless resource allocation method and cognitive radio user equipment
CN101632275A (en) In 3GPP TD-CDMA system, use and non-ly lead the improved channel equalization that sign indicating number distributes in public
CN104821852B (en) A Spectrum Sensing Method Based on Multi-Antenna Instantaneous Power
CN105141385B (en) Multiband cooperative cognitive frequency spectrum sensing method
CN104202826B (en) A kind of CoMP multi-user dispatching methods based on QoE
CN104918207B (en) More D2D communication resource allocation methods based on frequency spectrum resource distribution in heterogeneous network
CN105375998B (en) The multiband cooperative frequency spectrum sensing method optimized based on sub-clustering
CN105141383B (en) Cognition wireless network primary user's signal detecting method
CN105491572A (en) Joint spectrum sensing method based on decision threshold optimization
CN105515698B (en) Multiband cooperative frequency spectrum sensing method based on dynamic adjustment
CN107734697A (en) Low-complexity multi-user group technology based on NOMA systems
CN105490761A (en) A vehicle network cognitive radio adaptive spectrum sensing method
CN105141384B (en) A kind of cognitive radio cooperative frequency spectrum sensing method
CN104780610A (en) A resource allocating method based on terminal direct communication
CN102413517A (en) A method for selecting resident subbands
CN101772172A (en) Multi-band based method for allocating frequency spectrum resources, device and system thereof
CN108834158B (en) Interference management method for ultra-dense networking

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210202

Address after: Room 1101, building 2, 258 Xiqin street, Wuchang Street, Yuhang District, Hangzhou City, Zhejiang Province 310000

Patentee after: Hangzhou Lechi Xinchuang Artificial Intelligence Technology Service Co.,Ltd.

Address before: 315211, Fenghua Road, Jiangbei District, Zhejiang, Ningbo 818

Patentee before: Ningbo University

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20151209

Assignee: HANGZHOU HONGHAIER TECHNOLOGY Co.,Ltd.

Assignor: Hangzhou Lechi Xinchuang Artificial Intelligence Technology Service Co.,Ltd.

Contract record no.: X2021330000665

Denomination of invention: Multi band cooperative cognitive spectrum sensing method

Granted publication date: 20170711

License type: Common License

Record date: 20211103

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20151209

Assignee: Hangzhou Baowen Network Media Technology Co.,Ltd.

Assignor: Hangzhou Lechi Xinchuang Artificial Intelligence Technology Service Co.,Ltd.

Contract record no.: X2021330000685

Denomination of invention: Multi band cooperative cognitive spectrum sensing method

Granted publication date: 20170711

License type: Common License

Record date: 20211104

Application publication date: 20151209

Assignee: Hangzhou intellectual property operation management Co.,Ltd.

Assignor: Hangzhou Lechi Xinchuang Artificial Intelligence Technology Service Co.,Ltd.

Contract record no.: X2021330000686

Denomination of invention: Multi band cooperative cognitive spectrum sensing method

Granted publication date: 20170711

License type: Common License

Record date: 20211104

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20151209

Assignee: Hangzhou yunzhichuang Technology Co.,Ltd.

Assignor: Hangzhou Lechi Xinchuang Artificial Intelligence Technology Service Co.,Ltd.

Contract record no.: X2021330000767

Denomination of invention: Multi band cooperative cognitive spectrum sensing method

Granted publication date: 20170711

License type: Common License

Record date: 20211130

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20221031

Address after: 101100 705, Floor 7, Building 3, Yard 2, Beihuangmuchang Street, Tongzhou District, Beijing

Patentee after: Beijing Yierbei Health Technology Co.,Ltd.

Address before: Room 1101, building 2, 258 Xiqin street, Wuchang Street, Yuhang District, Hangzhou City, Zhejiang Province 310000

Patentee before: Hangzhou Lechi Xinchuang Artificial Intelligence Technology Service Co.,Ltd.

EC01 Cancellation of recordation of patent licensing contract
EC01 Cancellation of recordation of patent licensing contract

Assignee: Hangzhou yunzhichuang Technology Co.,Ltd.

Assignor: Hangzhou Lechi Xinchuang Artificial Intelligence Technology Service Co.,Ltd.

Contract record no.: X2021330000767

Date of cancellation: 20230105

Assignee: Hangzhou Baowen Network Media Technology Co.,Ltd.

Assignor: Hangzhou Lechi Xinchuang Artificial Intelligence Technology Service Co.,Ltd.

Contract record no.: X2021330000685

Date of cancellation: 20230105

Assignee: Hangzhou intellectual property operation management Co.,Ltd.

Assignor: Hangzhou Lechi Xinchuang Artificial Intelligence Technology Service Co.,Ltd.

Contract record no.: X2021330000686

Date of cancellation: 20230105

Assignee: HANGZHOU HONGHAIER TECHNOLOGY Co.,Ltd.

Assignor: Hangzhou Lechi Xinchuang Artificial Intelligence Technology Service Co.,Ltd.

Contract record no.: X2021330000665

Date of cancellation: 20230105

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20241122

Address after: Room 903, Building A, Intellectual Property Innovation Industrial Park, Yuhang District, Hangzhou City, Zhejiang Province, 311100

Patentee after: Yuhang District Artificial Intelligence Industry Intellectual Property Alliance in Hangzhou City

Country or region after: China

Address before: 101100 705, Floor 7, Building 3, Yard 2, Beihuangmuchang Street, Tongzhou District, Beijing

Patentee before: Beijing Yierbei Health Technology Co.,Ltd.

Country or region before: China