[go: up one dir, main page]

CN101399586B - 信号处理方法 - Google Patents

信号处理方法 Download PDF

Info

Publication number
CN101399586B
CN101399586B CN2007101754311A CN200710175431A CN101399586B CN 101399586 B CN101399586 B CN 101399586B CN 2007101754311 A CN2007101754311 A CN 2007101754311A CN 200710175431 A CN200710175431 A CN 200710175431A CN 101399586 B CN101399586 B CN 101399586B
Authority
CN
China
Prior art keywords
user
mrow
arrival
angle
deviation
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.)
Expired - Fee Related
Application number
CN2007101754311A
Other languages
English (en)
Other versions
CN101399586A (zh
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.)
TD Tech Ltd
Original Assignee
TD Tech Ltd
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 TD Tech Ltd filed Critical TD Tech Ltd
Priority to CN2007101754311A priority Critical patent/CN101399586B/zh
Publication of CN101399586A publication Critical patent/CN101399586A/zh
Application granted granted Critical
Publication of CN101399586B publication Critical patent/CN101399586B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种信号处理方法,包括:根据用户的不同到达角用户信号的方向性矢量、所述用户对应的下行赋形矢量以及接收天线在不同到达角的方向性增益,确定所述用户的不同到达角用户信号的强度;根据所述不同到达角用户信号的强度,确定所述用户的方位估计值。另一种信号处理方法,包括:确定用户的方位估计值;根据设定的衰减门限对应的两个用户信号到达角的差,估算所述用户的方位估计值存在的偏差;利用估算出的偏差,对所述用户的方位估计值进行修正。本发明提供的两种方法,可以减少甚至消除用户方位估计值与真实值之间的偏差。

Description

信号处理方法
技术领域
本发明涉及无线通信技术,特别涉及信号处理方法。
背景技术
基于用户位置的移动定位业务(LBS)的实现依赖于用户方位的准确估计。现有的用户方位估计方法,通常根据用户的多个来波方向的用户信号强度估计该用户的方位。具体的实现方法包括:
根据基站端天线阵列接收到的某一用户的信号估计信道冲击响应,计算下行赋形矢量;根据计算所得的下行赋形矢量,确定基站端天线阵列接收到的该用户在不同到达角的用户信号强度;将用户信号强度最大的到达角作为该用户的方位估计值。
其中,基站端天线阵列接收到用户的不同到达角信号的强度与对应天线阵列的方向性增益有关。在使用扇区天线的系统中,来自不同方向的用户信号在接收天线处的方向性增益不同,越靠近天线的法线方向,天线的方向性增益越大;相反越远离天线的法线方向,天线的方向性增益越小。上述这种方向性增益的不同导致不同到达角的用户信号的增益不均衡,也就是说,当来自用户的不同信号的强度相同时,如果这两个信号的到达角不同,则基站端天线阵列接收到的这两个信号的强度会由于方向性增益的不同而不同,进而造成用户方位估计的偏差。
比如,当用户发出的信号通过两个不同方向到达基站端天线阵列,这两个信号的强度相同,用户的真实方位靠近两个到达方向的中间角度。用户信号的一个到达方向在天线的法线方向,另一个到达方向在偏离法线方向约90度的方向,在这种情况下,因为靠近法线的方向性增益大,而远离法线的方向性增益小,因此基站端天线阵列接收到的两个信号的强度将差别很大,导致最后估计出的用户方位靠近法线方向,偏离真实值。
总之,由于在扇区天线系统中,存在天线的方向性增益差别,因此用户信号到达角的扩散范围越大,天线的方向性增益对用户方位估计值的准确性影响就越大。如果用户方位估计不准确,将直接影响基于用户位置的LBS业务的实现。
发明内容
有鉴于此,本发明实施例提供了两种信号处理方法,减少用户方位估计值与真实值之间的偏差。
本发明实施例提供的一种信号处理方法,包括:
根据用户的不同到达角用户信号的方向性矢量、所述用户对应的下行赋形矢量以及接收天线在不同到达角的方向性增益,确定所述用户的不同到达角用户信号的强度;
根据所述不同到达角用户信号的强度,确定所述用户的方位估计值。
本发明实施例提供的另一种信号处理方法,包括:
确定用户的方位估计值;
根据设定的衰减门限对应的两个用户信号到达角的差,估算所述用户的方位估计值存在的偏差;
利用估算出的偏差,对所述用户的方位估计值进行修正。
由以上技术方案的描述可知,本发明提供的第一种信号处理方法,在确定基站端天线阵列接收到的不同到达角的用户信号强度gi时,考虑接收天线的方向性增益在不同到达角上的差异,之后再根据确定出的gi确定用户方位估计值,进而可以减少用户方位的估计偏差。本发明提供的另一种信号处理方法,首先采用现有技术中的用户方位估计方法确定出用户方位估计值,然后再根据预设的衰减门限对应的两个用户信号到达角的差,估算用户方位估计值存在的偏差;利用估算出的偏差对确定出的用户方位估计值进行修正,从而可以减少用户方位估计值的偏差。
附图说明
图1为本发明信号处理方法实施例一的流程图;
图2为本发明实施例与现有技术相比gi与到达角的曲线关系比较图;
图3为接收天线方向性增益函数曲线的具体实例;
图4为本发明信号处理方法实施例二的流程图。
具体实施方式
为使本发明的目的、技术方案和有益效果更加清楚明白,下面结合实施例和附图,对本发明做进一步地详细说明。
本发明实施例提供的第一种信号处理方法,首先采用现有技术中的用户方位估计方法确定出用户方位估计值,然后再根据预设的衰减门限对应的两个用户信号到达角的差,估算用户方位估计值存在的偏差;利用估算出的偏差对确定出的用户方位估计值进行修正。
图1为本发明信号处理方法实施例一的流程图,该流程包括:
步骤101,确定用户方位估计值。
具体可以采用如下两种方式:
一种方式包括:
选择角度集合,如θi=0~360°,i=0,1,...,预先定义所选定角度的方向性矢量(steering vector)si,且满足‖si‖=1,i=0,1,.....;
根据基站端天线阵列接收到的某一用户的信号估计信道冲击响应,计算下行赋形矢量w;
采用公式(1)计算同一用户对应的不同到达角的信号的方向性矢量si与下行赋形矢量w的内积绝对值gi,gi表征基站端天线阵列接收到的不同到达角的用户信号的强度,但并不表示不同到达角信号到达基站天线阵列的实际强度。
gi=‖wHsi‖;    (1)
其中,wH为下行赋形矢量w的共轭转置,si为第i个用户信号对应的方向矢量。
选择使gi取得最大值的角度
Figure S2007101754311D00041
作为用户方位的估计值,如公式(2)所示:
θ ^ = arg θ i { max ( g i ) } ; - - - ( 2 )
另一种方式包括:
采用与第一种方式相同的方法确定出同一用户对应的全部gi,之后确定出大于设定门限值的gi对应的所有角度,如公式(3)所示:
j;gj≥max(gi)-Δ}    (3)
然后,根据不同的权重,对确定出的所有角度进行加权平均,得出用户的方位估计值,如公式(4)所示:
θ ^ = Σ θ j α j θ j ; g j ≥ max ( g i ) - Δ - - - ( 4 )
其中,加权系数αj和门限值Δ可以根据仿真结果确定,由系统预先设定,比如可以取Δ=3~5dB。
步骤102,根据设定的衰减门限对应的两个用户信号到达角的差,估算用户方位估计值存在的偏差。
以图2中虚线所示的gi与到达角的曲线关系为例,说明估算用户方位估计值所存在偏差的方法:
图2中虚线表示采用步骤101所述方式确定的gi与到达角的曲线关系,从图中可见,虚线中最大gi对应的角度为-40度,假设预先设定的衰减门限为3dB,则用最大gi减去3dB,确定该值对应的两个到达角的差Ω,Ω为图2中虚线所示曲线在设定衰减门限对应的波束宽度,Ω越大,表示用户到达角的扩散程度越大,则接收天线的方向性增益的不平衡对用户方位估计的影响就越大。
将确定出的Ω代入公式(5),估算用户方位估计值所存在偏差;
{ G ( θ ^ - a ( Ω ) ) - G ( θ ^ + b ( Ω ) ) } × c , - - - ( 5 )
其中,
Figure S2007101754311D00052
表示当到达角为
Figure S2007101754311D00053
时,接收天线的方向性增益;
Figure S2007101754311D00054
表示当到达角为
Figure S2007101754311D00055
时,接收天线的方向性增益。
Figure S2007101754311D00056
为步骤101确定出的用户方位估计值;a、b和c为预先设定的常数,这些常数随接收天线参数的变化而变化,其具体数值可以根据仿真结果确定。
需要说明的是,接收天线的方向性增益函数
Figure S2007101754311D00057
表征接收天线的方向性增益随到达角的变化,通常到达角以天线的法线为0度。如图3所示为接收天线方向性增益函数曲线的具体实例。
另外,在实际实现中,即使用户信号的到达角扩展为0,其对应的波束宽度Ω也并非为0,而是趋向一个确定的最小值,为了进一步修正用户方位估计值,可以采用公式(6),估算用户方位估计值所存在偏差;
{ G ( θ ^ - a ( Ω - Ω 0 ) ) - G ( θ ^ + b ( Ω - Ω 0 ) ) } × c , - - - ( 6 )
其中Ω0为用户信号的到达角扩展为0时,对应的波束宽度。Ω为设定的衰减门限对应的两个用户信号到达角的差;Ω0为用户信号的到达角扩展为0时,对应的波束宽度,Ω0为常数;
Figure S2007101754311D00059
为步骤101确定出的用户方位估计值;
Figure S2007101754311D000510
表示当到达角为
Figure S2007101754311D000511
时,接收天线的方向性增益;
Figure S2007101754311D000512
表示当到达角为
Figure S2007101754311D000513
时,接收天线的方向性增益;a、b和c为设定的常数。a、b、c以及Ω0随接收天线参数的变化而变化,其具体数值可以根据仿真结果确定。
步骤103,利用估算出的用户方位估计值存在的偏差,对步骤101确定出的方位估计值进行修正。
具体的实现方法可以为:利用公式(7)确定修正后的用户方位估计值值。
Figure S2007101754311D000514
本发明实施例提供的第二种信号处理方法,在确定基站端天线阵列接收到的不同到达角的用户信号强度gi时,考虑接收天线的方向性增益在不同到达角上的差异,之后采用步骤101中所述的两种方式中的任一种,根据确定出的gi确定用户方位估计值。
当然,为了进一步提高用户方位估计值的精确度,还可以采用步骤102和步骤103所述的方法对确定的用户方位估计值进行修正。
图4为本发明信号处理方法实施例二的流程图,该流程包括:
步骤401,定义选定角度的方向性矢量si,使其满足‖si‖=1,i=0,1,.....。
步骤402,根据基站端天线阵列接收到的某一用户的信号估计信道冲击响应,计算该用户对应的下行赋形矢量w。
计算下行赋形矢量的方法为现有技术,这里不再赘述。
步骤403,根据该用户的不同到达角用户信号的方向性矢量、该用户对应的下行赋形矢量以及接收天线在不同到达角的方向性增益,确定该用户的不同到达角用户信号的强度gi
一种较佳的具体实现方式可以为:采用公式(8)确定基站端天线阵列接收到的某一用户的不同到达角用户信号的强度gi
g i = | | W H S i | | / | G ( θ i ) | - - - ( 8 )
其中,设定接收天线的方向性增益函数为G(θ),|G(θi)|为到达角为θi时对应的接收天线方向性增益的绝对值。采用公式(8)计算出的gi可以在一定程度上补偿由于接收天线的方向性增益差别造成的用户方位估计偏差。
当然,公式(8)给出的计算方式只是一种较佳方式,在实际实现中还可以采用其它方式考虑接收天线的方向性增益在不同到达角上的差异,使得确定出的用户信号强度gi更接近真实值,进而减小用户方位估计的偏差。
步骤404,根据计算出的所有gi,确定用户方位估计值。
具体的实现方式可以为:
选择使gi取得最大值的角度
Figure S2007101754311D00071
作为用户方位的估计值,如公式(2)所示
也可以为:
首先确定出大于设定门限值的gi对应的所有角度,然后,根据设定的权重,对确定出的所有角度进行加权平均,得出用户的方位估计值,如公式(3)和(4)所示。
利用该实施例的方法实现简单,可以在一定程度上减小用户方位估计值的偏差。图2中实线所示曲线即为通过该方法确定的gi与用户信号到达角的关系,图2所示曲线对应的用户的实际位置在-45度,利用现有技术方法确定的用户方位估计值如虚线所示为-40度,利用该实施例确定的用户方位估计值为-43度。由此可见该实施例提供的方法可以减小用户方位估计值的偏差。
该实施例实现简单,但其没有考虑用户信号到达角的分布和用户方位对用户方位估计值的影响。为了进一步提高用户方位估计值的精确度,还可以在步骤404确定出用户方位估计值后,进一步采用步骤102和103所述的方法对确定出的用户方位估计值进行修正。
以图2实线所示gi与到达角的曲线关系为例,说明估算用户方位估计值所存在偏差的方法:
图2中实线表示采用步骤301-303所述方法确定的gi与到达角的曲线关系,从图中可见,实线曲线中最大gi对应的角度为-43度,假设预先设定的衰减门限为3dB,则用最大gi减去3dB,确定该值对应的两个到达角的差Ω,Ω为图2中实线所示曲线在设定衰减门限对应的波束宽度,Ω越大,表示用户到达角的扩散程度越大,则扇区天线的方向性增益的不平衡对用户方位估计的影响就越大。
将确定出的Ω代入公式(5)或公式(6),估算用户方位估计值所存在偏差;之后再利用估算出的用户方位估计值存在的偏差,对步骤404确定出的方位估计值进行修正,具体可以利用公式(7)确定修正后的用户方位估计值。通过公式(7)修正后的用户的方位估计偏差几乎为0,因此,利用本发明实施例提供的信号处理方法可以显著提高用户方位估计值的精确度,为基于用户位置的LBS业务的实现提供了保证。
总之,以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种信号处理方法,其特征在于,包括:
根据用户的不同到达角用户信号的方向性矢量、所述用户对应的下行赋形矢量以及接收天线在不同到达角的方向性增益,确定所述用户的不同到达角用户信号的强度;
根据所述不同到达角用户信号的强度,确定所述用户的方位估计值。
2.如权利要求1所述的方法,其特征在于,所述确定用户的不同到达角用户信号的强度,具体为:
所述用户的一个到达角用户信号的强度,等于所述用户的下行赋形矢量的共轭转置与所述用户的该到达角用户信号的方向性矢量的内积绝对值,再除以该到达角接收天线的方向性增益的绝对值。
3.如权利要求1所述的方法,其特征在于,所述确定用户的方位估计值,具体为:
选择所述不同到达角用户信号的强度中,最大值对应的到达角;
或者,首先确定出大于设定门限值的所述用户信号的强度对应的角度,然后,根据设定的权重,对确定出的所有角度进行加权平均。
4.如权利要求1所述的方法,其特征在于,确定出所述用户的方位估计值后,该方法进一步包括:
根据设定的衰减门限对应的两个用户信号到达角的差,估算所述用户的方位估计值存在的偏差;
利用估算出的偏差,对所述用户的方位估计值进行修正。
5.如权利要求4所述的方法,其特征在于,估算所述用户的方位估计值存在的偏差,具体为:
采用下述公式计算所述偏差,
{ G ( θ ^ - a ( Ω ) ) - G ( θ ^ + b ( Ω ) ) } × c ,
其中,Ω为设定的衰减门限对应的两个用户信号到达角的差;
Figure FSB00000703385700012
为所述用户的方位估计值;
Figure FSB00000703385700021
表示当到达角为
Figure FSB00000703385700022
时,接收天线的方向性增益;
Figure FSB00000703385700023
表示当到达角为
Figure FSB00000703385700024
时,接收天线的方向性增益;a、b和c为设定的常数。
6.如权利要求4所述的方法,其特征在于,估算所述用户的方位估计值存在的偏差,具体为:
采用下述公式计算所述偏差,
{ G ( θ ^ - a ( Ω - Ω 0 ) ) - G ( θ ^ + b ( Ω - Ω 0 ) ) } × c ,
其中,Ω为设定的衰减门限对应的两个用户信号到达角的差;Ω0为所述用户的用户信号的到达角扩展为0时,对应的波束宽度;
Figure FSB00000703385700026
为所述用户的方位估计值;表示当到达角为时,接收天线的方向性增益;
Figure FSB00000703385700029
表示当到达角为
Figure FSB000007033857000210
时,接收天线的方向性增益;a、b和c为设定的常数。
CN2007101754311A 2007-09-29 2007-09-29 信号处理方法 Expired - Fee Related CN101399586B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2007101754311A CN101399586B (zh) 2007-09-29 2007-09-29 信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2007101754311A CN101399586B (zh) 2007-09-29 2007-09-29 信号处理方法

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN201210010306.6A Division CN103780288A (zh) 2007-09-29 2007-09-29 信号处理方法

Publications (2)

Publication Number Publication Date
CN101399586A CN101399586A (zh) 2009-04-01
CN101399586B true CN101399586B (zh) 2012-04-18

Family

ID=40517879

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2007101754311A Expired - Fee Related CN101399586B (zh) 2007-09-29 2007-09-29 信号处理方法

Country Status (1)

Country Link
CN (1) CN101399586B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106372387B (zh) * 2016-08-24 2019-04-02 重庆大学 一种天线信号的获取方法及装置
JP2020112354A (ja) * 2017-05-11 2020-07-27 シャープ株式会社 端末装置、ロケーションサーバおよび端末装置における方法
CN108226858A (zh) * 2017-12-14 2018-06-29 厦门市美亚柏科信息股份有限公司 定位探测方法及计算机可读存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1668936A (zh) * 2002-09-06 2005-09-14 诺基亚公司 用于估计移动设备的位置的方法和系统
CN1864344A (zh) * 2003-08-06 2006-11-15 香港应用科技研究院有限公司 无线网络中的定位

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1668936A (zh) * 2002-09-06 2005-09-14 诺基亚公司 用于估计移动设备的位置的方法和系统
CN1864344A (zh) * 2003-08-06 2006-11-15 香港应用科技研究院有限公司 无线网络中的定位

Also Published As

Publication number Publication date
CN101399586A (zh) 2009-04-01

Similar Documents

Publication Publication Date Title
EP2098875B1 (en) Received field intensity estimating device, and received field intensity estimating program
US8145138B2 (en) Apparatus and method for removing interference in multiple antenna system
EP1014599A1 (en) Radio communication device and transmitting power control method
JP5235560B2 (ja) 端末位置推定システム及び方法ならびに位置推定装置
KR101375425B1 (ko) 배열 안테나에서 전파의 도래각 추정 방법
KR20030069096A (ko) 안테나의 지향각을 제어하는 무선 수신 장치 및 방법
WO2006114045A1 (fr) Procede de formation de faisceau permettant de supprimer les parasites
CN101674140A (zh) 天线校准方法及装置
US20100008268A1 (en) Frequency transformation based transmit beamforming in a communication system
US6501943B1 (en) Adaptive directivity transmission device and method
US7924957B2 (en) Method for creating beamformed multiple-input-multiple-output channels with partial nulling
EP3120417B1 (en) Method for finding signal direction using modal antenna
US20150163073A1 (en) Massive mimo channel estimation
CN101399586B (zh) 信号处理方法
US8577388B2 (en) Apparatus and method for determining position
EP1571761B1 (en) A method for forming a beam of an array antenna and apparatus
US20030021354A1 (en) Transmitter, the method of the same and communication system
JP2007312205A (ja) 基地局装置およびアレイアンテナ制御方法
JP3746983B2 (ja) 移動通信端末、及びそのアレーアンテナ指向性制御方法
JP2003018053A (ja) 通信装置
CN101754362B (zh) 一种双极化智能天线系统中的doa估计方法
JP2007504729A (ja) 無線通信システムにおける信号の伝送方法、ならびに相応する送信局及び受信局
JP4628985B2 (ja) アレーアンテナシステムおよびビーム制御方法
CN103780288A (zh) 信号处理方法
JP4660562B2 (ja) 移動局方向推定方法及び装置

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120418

Termination date: 20210929