CN103984833A - 一种简化的有向天线极化建模方法 - Google Patents
一种简化的有向天线极化建模方法 Download PDFInfo
- Publication number
- CN103984833A CN103984833A CN201410230924.0A CN201410230924A CN103984833A CN 103984833 A CN103984833 A CN 103984833A CN 201410230924 A CN201410230924 A CN 201410230924A CN 103984833 A CN103984833 A CN 103984833A
- Authority
- CN
- China
- Prior art keywords
- prime
- theta
- phi
- sin
- antenna
- 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
Links
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明提供一种简化的有向天线极化建模方法,当有向天线旋转时,其场波瓣图会发生复杂的形变,计算起来是较为复杂的。但是在实际应用中,有向天线一般应用于基站,基站天线的旋转角度是固定的,且一般下倾角不会超过20°。在这种情况下,本发明提出一种简化的极化天线建模方法来进行基站有向天线的极化建模,其中心思想是在容忍一定误差的基础上,用场波瓣图的直接旋转来代替旋转后形变的场波瓣图,以达到简化计算的目的。
Description
技术领域
本发明属于无线通信系统中多天线(MIMO)信道建模领域,涉及一种简化的有向天线极化建模方法。
背景技术
天线的极化建模就是计算天线的场波瓣图(field pattern)。在信道建模时,尤其是在建立三维的信道模型时,天线的场波瓣图是计算信道系数的重要参数。在已知有向天线在自身坐标系中的辐射波瓣图(radiation pattern)的情况下,可以采用已存在的方法计算得有向天线在全局坐标系中的场波瓣图,计算方法和有向天线自身坐标系相对于全局坐标系的旋转角度有关。当有向天线相对于全局坐标系旋转时,尤其是向下倾斜时,天线的场波瓣图不会简单地进行旋转,而是在旋转的同时发生较为复杂的形变,带来较大的计算复杂度。
发明内容
本发明的目的在于提供一种简化的根据天线旋转角来计算有向天线的场波瓣图的方法。
为达到上述目的,本发明采用了以下技术方案为:
包括以下步骤:
1)、设定全局坐标系GCS,得到第一三维坐标系[x,y,z];
2)、首先,将全局坐标系GCS绕z轴旋转角度α,得到第二三维坐标系然后,第二三维坐标系绕轴旋转角度β,得到第三三维坐标系第三三维坐标系绕轴旋转角度γ,得到的局部坐标系LCS坐标系其坐标轴记为[x′,y′,z′],其中,α为方位角,β为下倾角,γ为倾斜角;
3)、设LCS坐标系中天线的辐射波瓣图为A(θ′,φ′),其中,(θ′,φ′)为LCS球坐标系中的天顶角和方位角,天顶角和方位角定义了LCS球坐标系中的一个方向,设(θ,φ)为该方向在GCS球坐标系中的天顶角和方位角,其中,φ′=φ-α,θ′=θ-β;
4)、根据(θ′,φ′)和(θ,φ),计算天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ),并根据天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ),建立有向天线极化模型,以获得信道系数。
根据(θ′,φ′)和(θ,φ),计算天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ)包括:
2.1、将(φ′,θ′)映射到区间φ′∈(-180°,180°],θ′∈[0°,180°]上;
2.2、使用下式计算Fθ(θ,φ)和Fφ(θ,φ):
式中,
其中,ζ为天线极化倾斜角。
根据(θ′,φ′)和(θ,φ),计算天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ)包括:
3.1、所述将(φ′,θ′)映射到区间φ′∈(-180°,180°],θ′∈[0°,180°]上;
3.2、使用下式计算Fθ(θ,φ)和Fφ(θ,φ):
其中,ζ为天线极化倾斜角。
所述天线的下倾角β≤20°,倾斜角γ=0°。与现有技术比较,本发明的有与现有技术比较,本发明的益效果为:
本发明提供了一种简化的有向天线极化建模方法,通过将定义天线的自身坐标系的局部坐标系LCS转化到定义网络拓扑中基站和移动台的位置的全局坐标系GCS中,最终获得天线在GCS中的垂直极化波瓣图和水平极化波瓣图,使得用天线场波瓣图的直接旋转来代替波瓣图的复杂变换,简化了运算。
附图说明
图1为有向天线的辐射方向图,其中,(a)为水平切面上的方向图,(b)为垂直切面上的方向图;
图2为GCS到LCS的转化方法;
图3为球坐标系以及极化方向的定义;
图4为+45°极化有向天线旋转前的场波瓣图,其中,(a)为垂直极化场波瓣图,(b)为水平极化场波瓣图;
图5为+45°极化有向天线下倾45°后的场波瓣图,其中,(a)为垂直极化场波瓣图,(b)为水平极化场波瓣图,(c)为垂直极化场波瓣图,(d)为水平极化场波瓣图;
图6为+45°极化有向天线下倾20°后的场波瓣图,其中,(a)为垂直极化场波瓣图,(b)为水平极化场波瓣图,(c)为垂直极化场波瓣图,(d)为水平极化场波瓣图;
图7为在β=20°时的场波瓣图的形变差值,其中,(a)为垂直极化场波瓣图,(b)为水平极化场波瓣图;
图8为+45°极化有向天线的平均形变比例,(a)为垂直极化场波瓣图,(b)为水平极化场波瓣图;
图9为+45°极化有向天线的最大形变比例,(a)为垂直极化场波瓣图,(b)为水平极化场波瓣图。
具体实施方式
下面结合附图对本发明做详细描述。
本发明天线极化建模方法的应用条件是天线的下倾角(downtilt)β≤20°,天线的倾斜角γ=0°。图1为有向天线的辐射方向图,其中,(a)为水平切面上的方向图,(b)为垂直切面上的方向图。
考虑信道模型的坐标系统由全局坐标系GCS(Global CoordinateSystem)和局部坐标系LCS(Local Coordinate System)组成,其中GCS定义网络拓扑中基站和移动台的位置,而LCS则定义天线的波瓣图,LCS就是天线的自身坐标系,LCS相对于GCS的旋转即为基站或移动台的天线相对于基站或移动台的旋转,如图2所示,图2为GCS到LCS的转化方法。
本发明提供了一种简化的有向天线极化建模方法,包括以下步骤:
1)、设定全局坐标系GCS,得到第一三维坐标系[x,y,z],如图2所示,
2)首先,将全局坐标系GCS绕z轴旋转角度α,得到第二三维坐标系然后,第二三维坐标系绕轴旋转角度β,得到第三三维坐标系第三三维坐标系绕轴旋转角度γ,得到的局部坐标系LCS坐标系其坐标轴记为[x′,y′,z′],其中,α为方位角,β为下倾角,γ为倾斜角;
3)、设LCS坐标系中天线的辐射波瓣图为A(θ′,φ′),其中,(θ′,φ′)为LCS球坐标系中的天顶角和方位角,天顶角和方位角定义了三维空间(LCS球坐标系)中的一个方向,设(θ,φ)为该方向在GCS球坐标系中的天顶角和方位角,其中,φ′=φ-α,θ′=θ-β,如图3所示,图3为球坐标系以及极化方向的定义;
4)、根据(θ′,φ′)和(θ,φ),计算天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ),并根据天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ),建立有向天线极化模型,以获得信道系数。
根据(θ′,φ′)和(θ,φ),计算天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ)本发明提出了两种模型,在应用时可以选择一个模型进行应用,它们的步骤分别介绍如下:
模型1:
2.1、将(φ′,θ′)映射到区间φ′∈(-180°,180°],θ′∈[0°,180°]上;
2.2、使用下式计算Fθ(θ,φ)和Fφ(θ,φ):
式中, 其中,ζ为天线极化倾斜角。
模型2:
3.1、所述将(φ′,θ′)映射到区间φ′∈(-180°,180°],θ′∈[0°,180°]上;
3.2、使用下式计算Fθ(θ,φ)和Fφ(θ,φ):
其中,ζ为天线极化倾斜角。
原来通用的计算方法如下所示,
现有方法中:
模型1:
第一步:计算
第二步:计算
其中,
第三步:计算
其中,
已有模型2:
第一步:计算
第二步:计算
其中,
可以看出,在β≤20°,γ=0°的条件下,本发明能简化GCS中场波瓣图的运算复杂度。
下面结合附图对本发明的原理作进一步说明。
本发明提出的极化天线简化建模方法的适用条件是天线的下倾角不大于20°,天线的倾斜角为0°。主要的思路是:在满足上述条件时,天线的旋转带来的天线场波瓣图的改变可以忽略,进而可以用天线场波瓣图的直接旋转来代替波瓣图的复杂变换。
具体解释如下:
天线方向的旋转和天线场波瓣图有复杂的对应关系,计算起来较为麻烦,事实上,因为基站天线的物理下倾角一般不会超过20°,且不会进行自身的倾斜,存在一种计算天线场波瓣图的简化近似方法。
首先我们来看天线朝向的旋转对天线场波瓣图的影响,如图4就是,一个朝向为(α=0,β=0)的+45°极化有向天线的场波瓣图,其中,(a)为垂直极化场波瓣图,(b)为水平极化场波瓣图;。
事实上,通过仿真结果可以得知水平旋转角度α只能改变天线场波瓣图的中心方向,改变不了场波瓣图的形状,但是垂直旋转角度β能改变天线场波瓣图的形状,图5的(a)和(b)就是图4中的天线垂直旋转45°后(即β=45°)的场波瓣图,作为对比,图5的(c)和(d)为图4直接垂直旋转45°后的图形。
从图5中我们可以看出,当垂直旋转角β改变时,天线的场波瓣图会发生形变。但是在实际应用中,基站天线的下倾角一般都不会太大,在此我们假设其不会超过20°。取β=20°,按照绘制图5的方法重新绘图得图6,其中,(a)为垂直极化场波瓣图,(b)为水平极化场波瓣图,(c)为垂直极化场波瓣图,(d)为水平极化场波瓣图。
从图6中我们可以看到,在垂直旋转β=20°时,天线的场波瓣图的形状虽然也有变化,但其变化的程度并不大,和原来的形状基本上是相似的。图7是β=20°时的天线场波瓣图相对于直接旋转20°的场波瓣图的差值的绝对值,即形变差值,其中,(a)为垂直极化图,(b)为水平极化图。
从图7可以看出,形变差量并不是很大。在此定义最大形变比例和平均形变比例作为衡量极化响应图形变的一种指标,设垂直旋转角β对应的场波瓣图为Fβ(θ,φ)而原来的场波瓣图直接旋转β角后的场波瓣图为则最大形变比例ΔFβ,max定义为:
而平均形变比例ΔFβ,mean定义为:
对于图6,可以求出,在垂直极化响应上,平均形变比例是0.0444,最大的形变比例是0.2494;在水平极化响应上,平均形变比例是0.0444,最大形变比例是0.25。可以看出,形变比较小的。
当β<20°时,形变会更小,图8为+45°极化天线的平均形变比例随着β变化的曲线,图9为同样条件下最大形变比例随着β变化的曲线,其中,(a)为垂直极化图,(b)为水平极化图。。
在垂直旋转角β=20°的情况下,使用不同极化模型(模型1或模型2)和天线极化倾斜角ζ(0°、+45°或-45°)时,得到的平均形变比例和最大形变比例如表1所示。
表1β=20°时各个场景下的极化响应图形变参数
在表1中,由于在ζ=0°时恒有所以计算不出平均形变比例和最大形变比例。
从表1中我们可以看出,在β=20°时的平均形变比例和最大形变比例都能维持在一个较小的范围内,那么当β<20°时,平均形变比例和最大形变比例将会更小,这说明当垂直旋转角不大于20°时,天线的极化响应图的形变是较小的。
既然形变是较小的,那么在容忍一定误差的条件下,就可以用天线场波瓣图的直接旋转来代替原来复杂的旋转对应关系,进而达到简化计算的目的,这也就是本发明的原理。
Claims (4)
1.一种简化的有向天线极化建模方法,其特征在于,包括以下步骤:
1)、设定全局坐标系GCS,得到第一三维坐标系[x,y,z];
2)、首先,将全局坐标系GCS绕z轴旋转角度α,得到第二三维坐标系然后,第二三维坐标系绕轴旋转角度β,得到第三三维坐标系第三三维坐标系绕轴旋转角度γ,得到的局部坐标系LCS坐标系其坐标轴记为[x′,y′,z′],其中,α为方位角,β为下倾角,γ为倾斜角;
3)、设LCS坐标系中天线的辐射波瓣图为A(θ′,φ′),其中,(θ′,φ′)为LCS球坐标系中的天顶角和方位角,天顶角和方位角定义了LCS球坐标系中的一个方向,设(θ,φ)为该方向在GCS球坐标系中的天顶角和方位角,其中,φ′=φ-α,θ′=θ-β;
4)、根据(θ′,φ′)和(θ,φ),计算天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ),并根据天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ),建立有向天线极化模型,以获得信道系数。
2.根据权利要求1所述的一种简化的有向天线极化建模方法,其特征在于,根据(θ′,φ′)和(θ,φ),计算天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ)包括:
2.1、将(φ′,θ′)映射到区间φ′∈(-180°,180°],θ′∈[0°,180°]上;
2.2、使用下式计算Fθ(θ,φ)和Fφ(θ,φ):
式中,
其中,ζ为天线极化倾斜角。
3.根据权利要求1所述的一种简化的有向天线极化建模方法,其特征在于,根据(θ′,φ′)和(θ,φ),计算天线在GCS中的垂直极化波瓣图Fθ(θ,φ)和水平极化波瓣图Fφ(θ,φ)包括:
3.1、所述将(φ′,θ′)映射到区间φ′∈(-180°,180°],θ′∈[0°,180°]上;
3.2、使用下式计算Fθ(θ,φ)和Fφ(θ,φ):
其中,ζ为天线极化倾斜角。
4.根据权利要求1至3任一项所述的一种简化的有向天线极化建模方法,其特征在于,所述天线的下倾角β≤20°,倾斜角γ=0°。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410230924.0A CN103984833B (zh) | 2014-05-28 | 2014-05-28 | 一种简化的有向天线极化建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410230924.0A CN103984833B (zh) | 2014-05-28 | 2014-05-28 | 一种简化的有向天线极化建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103984833A true CN103984833A (zh) | 2014-08-13 |
CN103984833B CN103984833B (zh) | 2017-04-26 |
Family
ID=51276802
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410230924.0A Active CN103984833B (zh) | 2014-05-28 | 2014-05-28 | 一种简化的有向天线极化建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103984833B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109214021A (zh) * | 2017-07-03 | 2019-01-15 | 上海东峻信息科技有限公司 | 一种电磁仿真中斜口面激励源的生成方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101789812A (zh) * | 2010-01-15 | 2010-07-28 | 重庆邮电大学 | 基于双极化天线的lte mimo通信传播信道建模方法 |
EP2256864A1 (de) * | 2009-05-30 | 2010-12-01 | Delphi Delco Electronics Europe GmbH | Antenne für zirkulare Polarisation mit einer leitenden Grundfläche |
CN103402227A (zh) * | 2013-08-09 | 2013-11-20 | 西安电子科技大学 | 一种winner信道仿真中的三维极化投影算法 |
-
2014
- 2014-05-28 CN CN201410230924.0A patent/CN103984833B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2256864A1 (de) * | 2009-05-30 | 2010-12-01 | Delphi Delco Electronics Europe GmbH | Antenne für zirkulare Polarisation mit einer leitenden Grundfläche |
CN101789812A (zh) * | 2010-01-15 | 2010-07-28 | 重庆邮电大学 | 基于双极化天线的lte mimo通信传播信道建模方法 |
CN103402227A (zh) * | 2013-08-09 | 2013-11-20 | 西安电子科技大学 | 一种winner信道仿真中的三维极化投影算法 |
Non-Patent Citations (2)
Title |
---|
GUO X 等: "A CANDECOMP/PARAFAC perspective on uniqueness of DOA estimation using a vector sensor array", 《SIGNAL PROCESSING, IEEE TRANSACTIONS ON》 * |
束咸荣 等: "斜极化和圆极化天线波瓣图免失真测试研究", 《研究与探讨》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109214021A (zh) * | 2017-07-03 | 2019-01-15 | 上海东峻信息科技有限公司 | 一种电磁仿真中斜口面激励源的生成方法 |
CN109214021B (zh) * | 2017-07-03 | 2023-04-21 | 上海东峻信息科技有限公司 | 一种电磁仿真中斜口面激励源的生成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103984833B (zh) | 2017-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107808395B (zh) | 一种基于slam的室内定位方法 | |
Gao et al. | Toward 5G NR high-precision indoor positioning via channel frequency response: A new paradigm and dataset generation method | |
CN104375115B (zh) | 极化敏感阵列下非圆信号doa与极化参数联合估计方法 | |
CN107015191A (zh) | 一种在多径干扰环境下单偶极子极化敏感阵列降维doa估计方法 | |
US7535425B2 (en) | Method and system for generating three-dimensional antenna radiation patterns | |
CN107888241A (zh) | 一种曲面共形圆极化相控阵天线波束合成算法 | |
CN109444809B (zh) | 一种基于智能天线的无人机测控方法 | |
CN108111208A (zh) | 一种灵活码本配置的毫米波系统角度域信道估计方法 | |
CN106025550A (zh) | 一种以电性能为目标的双反射面天线副面位置调整方法 | |
CN104181513A (zh) | 一种雷达天线阵元位置的校正方法 | |
CN104102844B (zh) | 基于11范数约束的有理函数模型参数求解和全参数优化方法 | |
WO2021027948A1 (zh) | 定位方法及其装置 | |
CN107144814B (zh) | 一种适于定向辐射场景的无源uhf rfid定位精度的评价方法 | |
CN103984836B (zh) | 一种3d‑mimo天线极化增益的联合存储方法 | |
CN103984833A (zh) | 一种简化的有向天线极化建模方法 | |
Pihlajasalo et al. | Absolute positioning with unsupervised multipoint channel charting for 5G networks | |
CN110444903A (zh) | 基于多环阵列的低副瓣涡旋波束产生方法 | |
CN108200547B (zh) | 基于测量距离的刚体定位方法 | |
CN103777179B (zh) | 一种用于雷达立体共形阵列的子阵降维方法 | |
CN103473421B (zh) | 地基天线温度补偿方法 | |
CN107144813B (zh) | 一种构建四阵元立体阵列的方法和装置 | |
CN103780296B (zh) | 一种均匀圆形阵列天线的信号接收及失真矫正方法 | |
CN102509898B (zh) | 一种低剖面椭圆波束卡氏天线的设计方法 | |
CN107466063B (zh) | 一种通信卫星多波束无线测试方法 | |
CN104268406A (zh) | 空间任意位置抛物面天线的拟合方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210104 Address after: 221600 North Side of Xiaohe Road, Peixian Development Zone, Xuzhou City, Jiangsu Province Patentee after: Xuzhou Guolong Electric Power Parts Foundry Co.,Ltd. Address before: 710049 No. 28 West Xianning Road, Shaanxi, Xi'an Patentee before: XI'AN JIAOTONG University |
|
TR01 | Transfer of patent right |