CN113922900B - 一种考虑天线尺寸与天线互耦的无线电信道建模方法 - Google Patents
一种考虑天线尺寸与天线互耦的无线电信道建模方法 Download PDFInfo
- Publication number
- CN113922900B CN113922900B CN202111232635.0A CN202111232635A CN113922900B CN 113922900 B CN113922900 B CN 113922900B CN 202111232635 A CN202111232635 A CN 202111232635A CN 113922900 B CN113922900 B CN 113922900B
- Authority
- CN
- China
- Prior art keywords
- matrix
- antenna
- receiving end
- dimensions
- transmitting
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000008878 coupling Effects 0.000 title claims abstract description 46
- 238000010168 coupling process Methods 0.000 title claims abstract description 46
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 100
- 238000005314 correlation function Methods 0.000 claims abstract description 17
- 230000004044 response Effects 0.000 claims abstract description 10
- 238000012512 characterization method Methods 0.000 claims description 14
- 101000964052 Rattus norvegicus 5-hydroxytryptamine receptor 5B Proteins 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 5
- 230000005855 radiation Effects 0.000 claims description 4
- 238000003491 array Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000003990 capacitor Substances 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 4
- 238000005457 optimization Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000005311 autocorrelation function Methods 0.000 description 9
- 238000004891 communication Methods 0.000 description 9
- 238000004088 simulation Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 230000005672 electromagnetic field Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005562 fading Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004613 tight binding model Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
- H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种考虑天线尺寸与天线互耦的无线电信道建模方法,包括以下步骤:首先通过预设的场景、频段生成无线信道矩阵;然后根据天线阵元尺寸确定天线等效电路的结构以及电路器件的值,并根据天线阵元间距给出描述天线阵元间互耦的矩阵;接着根据电路内部和外部噪声源的统计特性生成噪声电压与电流,并确定发射端和接收端的阻抗匹配网络矩阵、发射端的源电压与源阻抗、接收端的负载阻抗;最后,将上述各电路模块进行级联,获取端到端响应,并由此分析电路系统的空间‑时间‑频率相关函数。本发明以电路多端口网络的形式对无线信道进行建模,加深了无线信道与天线的耦合,有助于天线阵元尺寸、天线阵列排布的设计与优化。
Description
技术领域
本发明属于通信技术领域,尤其涉及一种考虑天线尺寸与天线互耦的无线电信道建模方法。
背景技术
为了实现即将到来的第六代(Sixth Generation,6G)移动通信中的发展愿景与范式转换,研究人员一直在寻找可以提升未来无线通信信息传输速率和和可靠性的天线技术与无线传输技术。无线信道模型是无线通信系统设计、理论分析以及网络部署的基础。在当今的信道建模中,天线与无线信道往往是通过天线远场方向图进行联系的。然而,随着智能反射表面(Intelligent Reflecting Surface,IRS)和全息多输入多输出(Multiple-InputMultiple-Output,MIMO)天线等天线技术的出现,天线与信道的关系愈发紧密。因此,在无线信道建模的过程中需要将天线充分地考虑进来,为系统性能的优化提供更加准确的结果。
面对上述问题,已经有许多研究致力于将无线信道与天线更加紧密地联系起来,例如电磁信息论和通信电路理论。电磁信息论中的场自由度理论通过二元格林函数联系天线源电流分布和观测到的电磁场分布,得出辐射算子并进行特征分析;电磁信息论中的基函数扩展理论采用基函数扩展,提出计算连续电磁空间电磁信道容量的通用框架;电磁信息论中的球形矢量波扩展理论将信道分解为发射天线信道、电磁传播信道、接收天线信道的级联,将电磁场进行球形矢量波扩展。电磁信息论具有如下缺点:1)函数扩展或特征分解需要了解全空间电磁场的完整信息,这在实际应用中是很难实现的;2)在其确定性的模型假设中,电磁场的不确定度或噪声难以度量。通信电路理论为无线通信系统建立了基于电路理论的多端口模型,将无线信道与天线处的电磁场特性以多端口网络的形式描述,然而该理论当前研究中的信道建模多集中于大尺度衰落和视距,且并未全面地考虑MIMO情形和天线等效电路。另外,对通信电路理论中信道统计特性的研究也十分匮乏。
因此,需要一种综合考虑天线等效电路模型和无线信道小尺度衰落的MIMO无线电信道建模方法。
发明内容
本发明目的在于提供一种考虑天线尺寸与天线互耦的无线电信道建模方法,以解决现有无线信道建模方法中天线与信道联系不紧密的技术问题。
为解决上述技术问题,本发明的具体技术方案如下:
一种考虑天线尺寸与天线互耦的无线电信道建模方法,包括以下步骤:
步骤1、通过预设的场景、频段生成无线信道矩阵其中发射端和接收端均配备有垂直极化全向天线组成的均匀线阵,发射端具有MT根天线,接收端具有MR根天线,t为时间,f为频率,Hqp(t,f)为第p根发射天线/>与第q根接收天线/>间的信道传输函数,其中p=1,...,MT,q=1,...,MR;
步骤2、根据垂直极化全向天线的尺寸a确定收发端天线等效电路中电阻、电感、电容的值,并根据发射端天线间距δT和接收端天线间距δR确定发射端天线互耦矩阵CT和接收端天线互耦矩阵CR;
步骤3、根据电路内部和外部噪声源的统计特性生成接收端外部噪声电压源接收端内部噪声电压源vN、接收端内部噪声电流源iN,并确定发射端阻抗匹配网络矩阵ZMT、接收端的阻抗匹配网络矩阵ZMR、发射端的源电压vG,设置收发端各电路端口的源阻抗和负载阻抗均为纯电阻R;
步骤4、将上述步骤1-3中的阻抗匹配网络、噪声源、天线等效电路进级联进行级联,获取发射端源电压vG至接收端负载电压vL的端到端响应,并由此分析电路系统的空间域、时间域和频率域的相关函数。
进一步的,步骤1中与/>间的信道传输函数Hqp(t,f)表示为视距分量/>与非视距分量/>的叠加,由下述公式确定:
其中,e为自然对数的底数,π为圆周率,KR为莱斯因子,Gr为接收端天线增益,Gt为发射端天线增益,fc为载波频率,Nqp(t)为时刻t时/>与/>间散射体簇的总数目,散射体簇序号n遍历n=1,2,...,Nqp(t),Mn为第n个散射体簇内散射体的总数目,散射体序号m遍历m=1,2,...,Mn,/>为时刻t时/>与/>间视距路径的时延,/>和/>分别为在时刻t时/>与/>之间的第n个散射体簇的第m个散射体的功率和时延,/>为在(-π,π]服从均匀分布的非视距路径第n个散射体簇的第m个散射体的初始随机相位。
进一步的,步骤2中,在发射端考虑尺寸大于十倍载频波长的天线,因此发射端天线的等效电路由单位电阻Rp构成,在接收端根据天线尺寸a确定接收端天线/>的等效电路由单位电阻Rq与电感aRq/c并联再同电容a/cRq串联构成,其中c为光速。
进一步的,步骤2中根据发射端天线间距δT和接收端天线间距δR确定托普利兹矩阵形式的发射端天线互耦矩阵CT和接收端天线互耦矩阵CR,计算公式如下:
其中CT为MT×MT维的矩阵,CR为MR×MR维的矩阵,k为波数且k=2πfc/c,函数j0(·)可以定义为j0(x)=sinx/x。
进一步的,步骤3中,接收端的外部噪声电压源由MR×1的向量表征,其统计特性由MR×MR维的协方差矩阵φ控制,由下述公式确定:
其中kb为玻尔兹曼常量,TA为天线等效噪声温度,B为带宽,Rr为天线辐射电阻,表示输入变量的数学期望,(·)H表示矩阵的共轭转置;接收端内部噪声电压源由MR×1的向量vN表征,接收端内部噪声电流源由MR×1的向量iN表征,vN与iN应满足如下条件:
其中β为iN中各元素的方差,RN为噪声电阻,ρ为vN与iN间的复相关系数,IM表示M×M维的单位矩阵;且设置发射端采用功率匹配,接收端采用噪声匹配,则发射端阻抗匹配网络矩阵由2MT×2MT维矩阵ZMT表征,接收端的阻抗匹配网络矩阵由2MR×2MR维矩阵ZMR表征,分别由下式给出:
其中ZMT11、ZMT12、ZMT21、ZMT22均为MT×MT维的矩阵,ZMR11、ZMR12、ZMR21、ZMR22均为MR×MR维的矩阵,Re{·}表示取输入变量的实部,Im{·}表示取输入变量的虚部,OM表示M×M维的全零矩阵,ZAT和ZAR可由如下公式计算:
ZAT=Rp·CT
其中ZAT为MT×MT维的矩阵,ZAR为MR×MR维的矩阵,s为拉普拉斯算子,即s=j2πf;另外,设置发射端的源电压由MT×1的向量vG表征。
进一步的,步骤4中,发射端源电压vG与接收端负载电压vL之间的端到端响应由如下公式计算:
其中系统传输矩阵D为MR×MT维的矩阵,vL和η为MR×1的向量,D和η的计算公式如下:
其中(·)-1表示矩阵的逆矩阵,ZT为MT×MT维的矩阵,ZR为MR×MR维的矩阵,ZRT为MR×MT维的矩阵,可由如下公式确定:
ZR=ZMR11-ZMR12(ZMR22+ZAR)-1ZMR21
ZT=ZMT11-ZMT12(ZMT22+ZAT)-1ZMT21
ZRT=ZMR12(ZMR22+ZAR)-1ZART[ZMT12(ZMT22+ZAT)-1]T
其中(·)T表示矩阵的转置矩阵,ZART为MR×MT维的矩阵,其第q行第p列的元素ZART(q,p)可由如下公式确定:
进一步的,步骤4中根据vG与vL的端到端响应中的系统传输矩阵D,整个系统的空间-时间-频率相关函数由下式计算:
其中(·)*表示输入变量的共轭,Δr为空间间隔,Δt为时间间隔,Δf为频率间隔,Dqp(t,f)为矩阵D第q行第p列的元素在t时刻f频率处的值,为矩阵D第行第/>列的元素在t-Δt时刻f-Δf频率处的值。
本发明的一种考虑天线尺寸与天线互耦的无线电信道建模方法,具有以下优点:
1、本发明在无线信道模型中考虑了天线等效电路、天线互耦以及阻抗匹配网络等与天线有关的效应和电路器件,以电路多端口网络的形式对无线信道进行建模,加深了无线信道与天线的耦合,有助于天线阵元尺寸、天线阵列排布的设计与优化。
2、本发明提出的无线信道建模方法可以通过参数的预设或测量拟合,提供整个通信系统的电路端到端响应,并由此分析系统的空间-时间-频率相关函数。
附图说明
图1为本发明的无线电信道模型框图;
图2为本发明的发射端任意天线阵元与接收端任意天线阵元间的电路端口网络图;
图3为本发明实施例中考虑天线尺寸与天线互耦的无线电信道建模方法与现有不考虑天线互耦的信道建模方法的时间自相关函数随天线阵元数目变化仿真对比图;
图4为本发明实施例中考虑天线尺寸与天线互耦的无线电信道建模方法与现有不考虑天线尺寸的信道建模方法的时间自相关函数随天线尺寸变化仿真对比图;
图5为本发明实施例中虑天线尺寸与天线互耦的无线电信道建模方法与现有不考虑天线互耦的信道建模方法的空间互相关函数随天线阵元位置与天线阵元间距变化仿真对比图;
图6为本发明实施例中考虑天线尺寸与天线互耦的无线电信道建模方法与现有不考虑天线尺寸的信道建模方法的空间互相关函数随天线尺寸变化仿真对比图。
具体实施方式
为了更好地了解本发明的目的、结构及功能,下面结合附图,对本发明一种考虑天线尺寸与天线互耦的无线电信道建模方法做进一步详细的描述。
如图1所示,本发明的无线电信道模型包括了发射端电压源、发射端源阻抗、发射端阻抗匹配网络、考虑与天线尺寸相关的天线等效电路、无线信道以及天线互耦的收发端连接模块、接收端外部噪声源、接收端阻抗匹配网络、接收端内部噪声源、接收端负载,最终得出发射端电压源与接收端负载电压之间的端到端响应。
本发明提出的一种考虑天线尺寸与天线互耦的无线电信道建模方法包括以下步骤:
步骤1、通过预设的场景、频段生成无线信道矩阵其中发射端和接收端均配备有垂直极化全向天线组成的均匀线阵,发射端具有MT根天线,接收端具有MR根天线,t为时间,f为频率,Hqp(t,f)为第p根发射天线/>与第q根接收天线/>间的信道传输函数(p=1,...,MT,q=1,...,MR)。Hqp(t,f)可以表示为视距分量/>与非视距分量的叠加,由下述公式确定:
其中,e为自然对数的底数,π为圆周率,KR为莱斯因子,Gr为接收端天线增益,Gt为发射端天线增益,fc为载波频率(λc为对应波长),Nqp(t)为时刻t时/>与/>间散射体簇的总数目,散射体簇序号n遍历n=1,2,...,Nqp(t),Mn为第n个散射体簇内散射体的总数目,散射体序号m遍历m=1,2,...,Mn,/>为时刻t时/>与/>间视距路径的时延,和/>分别为在时刻t时/>与/>之间的第n个散射体簇的第m个散射体的功率和时延,/>为在(-π,π]服从均匀分布的非视距路径第n个散射体簇的第m个散射体的初始随机相位。
步骤2、确定收发端天线等效电路的结构以及电路器件参数时(参考图2),在发射端考虑尺寸大于十倍载频波长的天线,由此确定发射端天线的等效电路由单位电阻Rp构成,在接收端根据天线尺寸a确定接收端天线/>的等效电路由单位电阻Rq与电感aRq/c并联再同电容a/cRq串联构成,其中c为光速;根据发射端天线间距δT和接收端天线间距δR可以确定托普利兹矩阵(Toeplitz Matrix)形式的发射端天线互耦矩阵CT和接收端天线互耦矩阵CR,计算公式如下:
其中CT为MT×MT维的矩阵,CR为MR×MR维的矩阵,k为波数且k=2πfc/c,函数j0(·)可以定义为j0(x)=sinx/x。
步骤3、通过MR×1的向量表征接收端的外部噪声电压源,其统计特性由MR×MR维的协方差矩阵φ控制,由下述公式确定:
其中kb为玻尔兹曼常量,TA为天线等效噪声温度,B为带宽,Rr为天线辐射电阻,表示输入变量的数学期望,(·)H表示矩阵的共轭转置;通过MR×1的向量vN表征接收端内部噪声电压源,通过MR×1的向量iN表征接收端内部噪声电流源,vN与iN应满足如下条件:
其中β为iN中各元素的方差,RN为噪声电阻,ρ为vN与iN间的复相关系数,IM表示M×M维的单位矩阵;且设置发射端采用功率匹配,接收端采用噪声匹配,则发射端阻抗匹配网络矩阵由2MT×2MT维矩阵ZMT表征,接收端的阻抗匹配网络矩阵由2MR×2MR维矩阵ZMR表征,分别由下式给出:
其中ZMT11、ZMT12、ZMT21、ZMT22均为MT×MT维的矩阵,ZMR11、ZMR12、ZMR21、ZMR22均为MR×MR维的矩阵,Re{·}表示取输入变量的实部,Im{·}表示取输入变量的虚部,OM表示M×M维的全零矩阵,ZAT和ZAR可由如下公式计算:
ZAT=Rp·CT
其中ZAT为MT×MT维的矩阵,ZAR为MR×MR维的矩阵,s为拉普拉斯算子,即s=j2πf;另外,设置发射端的源电压由MT×1的向量vG表征。
步骤4、通过如下公式计算发射端源电压vG与接收端负载电压vL之间的端到端响应:
其中系统传输矩阵D为MR×MT维的矩阵,vL和η为MR×1的向量,D和η的计算公式如下:
其中(·)-1表示矩阵的逆矩阵,ZT为MT×MT维的矩阵,ZR为MR×MR维的矩阵,ZRT为MR×MT维的矩阵,可由如下公式确定:
ZR=ZMR11-ZMR12(ZMR22+ZAR)-1ZMR21
ZT=ZMT11-ZMT12(ZMT22+ZAT)-1ZMT21
ZRT=ZMR12(ZMR22+ZAR)-1ZART[ZMT12(ZMT22+ZAT)-1]T
其中(·)T表示矩阵的转置矩阵,ZART为MR×MT维的矩阵,其第q行第p列的元素ZART(q,p)可由如下公式确定:
根据vG与vL的端到端响应中的系统传输矩阵D,整个系统的空间-时间-频率相关函数可由下式计算:
其中(·)*表示输入变量的共轭,Δr为空间间隔,Δt为时间间隔,Δf为频率间隔,Dqp(t,f)为矩阵D第q行第p列的元素在t时刻f频率处的值,为矩阵D第行第/>列的元素在t-Δt时刻f-Δf频率处的值。
图3是所提考虑天线尺寸与天线互耦的无线电信道建模方法与现有不考虑天线互耦的信道建模方法的时间自相关函数随天线阵元数目变化仿真对比图。其中fc=30GHz,a=0.01λc,δT=λc/2,δR=λc/4,p=1,q=1。可以看出,若采取不考虑天线互耦的信道建模方法,时间自相关函数不会因天线阵列阵元数目变化;若采取本发明的信道建模方法,尤其是在天线阵列阵元数目较多时,所得时间自相关函数会与采取不考虑天线互耦的信道建模方法得出的时间自相关函数有较大的差别。因此,在天线阵元数目较多且存在天线互耦时,本发明的信道建模方法可揭示不考虑天线互耦的信道建模方法难以体现的统计特性。
图4是所提考虑天线尺寸与天线互耦的无线电信道建模方法与现有不考虑天线尺寸的信道建模方法的时间自相关函数随天线尺寸变化仿真对比图。其中fc=30GHz,MR=6,δT=λc/2,δR=λc/4,p=1,q=1。可以看出,采取本发明的信道建模方法与采取不考虑天线尺寸的信道建模方法所得的时间自相关函数有较大的差别,且天线阵元尺寸越大,本发明的信道建模方法得出的时间自相关函数越大。本发明的信道建模方法的正确性可由图中理论值与仿真值的较好拟合验证。
图5是所提考虑天线尺寸与天线互耦的无线电信道建模方法与现有不考虑天线互耦的信道建模方法的空间互相关函数随天线阵元位置与天线阵元间距变化仿真对比图。其中fc=30GHz,MR=12,δT=λc/2,p=1,a=0.05λc。可以看出,本发明的信道建模方法与不考虑天线互耦的信道建模方法所得的空间互相关函数差别较大。本发明的信道建模方法可揭示一定的空域非平稳性,即空间互相关函数随天线阵元的空间位置变化,且天线阵元间距越小,空域非平稳性越不显著。因此,在本发明的信道建模方法可揭示不考虑天线互耦的信道建模方法难以体现的统计特性
图6是所提考虑天线尺寸与天线互耦的无线电信道建模方法与现有不考虑天线尺寸的信道建模方法的空间互相关函数随天线尺寸变化仿真对比图。其中fc=30GHz,MR=23,δT=λc/2,δR=λc/4,p=1,q=1。可以看出,采取本发明的信道建模方法与采取不考虑天线尺寸的信道建模方法所得的空间互相关函数有较大的差别,且本发明的信道建模方法所得的空间互相关函数会随天线尺寸的变化而变化。本发明的信道建模方法的正确性可由图中理论值与仿真值的较好拟合验证。
本发明提出了一种考虑天线尺寸与天线互耦的无线电信道建模方法。与现有的无线信道建模方法相比,本发明所提出的方法考虑了天线等效电路、天线互耦以及阻抗匹配网络等与天线有关的效应和电路器件,以电路多端口网络的形式对无线信道进行建模,加深了无线信道与天线的耦合,有助于天线阵元尺寸、天线阵列排布的设计与优化。本发明提出的无线信道建模方法可以通过参数的预设或测量拟合,提供整个通信系统的电路端到端响应,并由此分析系统的空间-时间-频率相关函数。
可以理解,本发明是通过一些实施例进行描述的,本领域技术人员知悉的,在不脱离本发明的精神和范围的情况下,可以对这些特征和实施例进行各种改变或等效替换。另外,在本发明的教导下,可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此,本发明不受此处所公开的具体实施例的限制,所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。
Claims (1)
1.一种考虑天线尺寸与天线互耦的无线电信道建模方法,其特征在于,包括以下步骤:
步骤1、通过预设的场景、频段生成无线信道矩阵其中发射端和接收端均配备有垂直极化全向天线组成的均匀线阵,发射端具有MT根天线,接收端具有MR根天线,t为时间,f为频率,Hqp(t,f)为第p根发射天线/>与第q根接收天线/>间的信道传输函数,其中p=1,...,MT,q=1,...,MR;
步骤2、根据垂直极化全向天线的尺寸a确定收发端天线等效电路中电阻、电感、电容的值,并根据发射端天线间距δT和接收端天线间距δR确定发射端天线互耦矩阵CT和接收端天线互耦矩阵CR;
步骤3、根据电路内部和外部噪声源的统计特性生成接收端外部噪声电压源接收端内部噪声电压源vN、接收端内部噪声电流源iN,并确定发射端阻抗匹配网络矩阵ZMT、接收端的阻抗匹配网络矩阵ZMR、发射端的源电压vG,设置收发端各电路端口的源阻抗和负载阻抗均为纯电阻R;
步骤4、将上述步骤1-3中的阻抗匹配网络、噪声源、天线等效电路进行级联,获取发射端源电压vG至接收端负载电压vL的端到端响应,并由此分析电路系统的空间域、时间域和频率域的相关函数;
所述步骤1中与/>间的信道传输函数Hqp(t,f)表示为视距分量/>与非视距分量/>的叠加,由下述公式确定:
其中,e为自然对数的底数,π为圆周率,KR为莱斯因子,Gr为接收端天线增益,Gt为发射端天线增益,fc为载波频率,Nqp(t)为时刻t时/>与/>间散射体簇的总数目,散射体簇序号n遍历n=1,2,...,Nqp(t),Mn为第n个散射体簇内散射体的总数目,散射体序号m遍历m=1,2,...,Mn,/>为时刻t时/>与/>间视距路径的时延,/>和/>分别为在时刻t时/>与/>之间的第n个散射体簇的第m个散射体的功率和时延,/>为在(-π,π]服从均匀分布的非视距路径第n个散射体簇的第m个散射体的初始随机相位;
所述步骤2中,在发射端考虑尺寸大于十倍载频波长的天线,因此发射端天线的等效电路由单位电阻Rp构成,在接收端根据天线尺寸a确定接收端天线/>的等效电路由单位电阻Rq与电感aRq/c并联再同电容a/cRq串联构成,其中c为光速;
所述步骤2中根据发射端天线间距δT和接收端天线间距δR确定托普利兹矩阵形式的发射端天线互耦矩阵CT和接收端天线互耦矩阵CR,计算公式如下:
其中CT为MT×MT维的矩阵,CR为MR×MR维的矩阵,k为波数且k=2πfc/c,函数j0(·)定义为j0(x)=sinx/x;
所述步骤3中,接收端的外部噪声电压源由MR×1的向量表征,其统计特性由MR×MR维的协方差矩阵φ控制,由下述公式确定:
其中kb为玻尔兹曼常量,TA为天线等效噪声温度,B为带宽,Rr为天线辐射电阻,E{·}表示输入变量的数学期望,(·)H表示矩阵的共轭转置;接收端内部噪声电压源由MR×1的向量vN表征,接收端内部噪声电流源由MR×1的向量iN表征,vN与iN应满足如下条件:
其中β为iN中各元素的方差,RN为噪声电阻,ρ为vN与iN间的复相关系数,IM表示M×M维的单位矩阵;且设置发射端采用功率匹配,接收端采用噪声匹配,则发射端阻抗匹配网络矩阵由2MT×2MT维矩阵ZMT表征,接收端的阻抗匹配网络矩阵由2MR×2MR维矩阵ZMR表征,分别由下式给出:
其中ZMT11、ZMT12、ZMT21、ZMT22均为MT×MT维的矩阵,ZMR11、ZMR12、ZMR21、ZMR22均为MR×MR维的矩阵,Re{·}表示取输入变量的实部,Im{·}表示取输入变量的虚部,OM表示M×M维的全零矩阵,ZAT和ZAR可由如下公式计算:
ZAT=Rp·CT
其中ZAT为MT×MT维的矩阵,ZAR为MR×MR维的矩阵,s为拉普拉斯算子,即s=j2πf;另外,设置发射端的源电压由MT×1的向量vG表征;
所述步骤4中,发射端源电压vG与接收端负载电压vL之间的端到端响应由如下公式计算:
其中系统传输矩阵D为MR×MT维的矩阵,vL和η为MR×1的向量,D和η的计算公式如下:
其中(·)-1表示矩阵的逆矩阵,ZT为MT×MT维的矩阵,ZR为MR×MR维的矩阵,ZRT为MR×MT维的矩阵,可由如下公式确定:
ZR=ZMR11-ZMR12(ZMR22+ZAR)-1ZMR21
ZT=ZMT11-ZMT12(ZMT22+ZAT)-1ZMT21
ZRT=ZMR12(ZMR22+ZAR)-1ZART[ZMT12(ZMT22+ZAT)-1]T
其中(·)T表示矩阵的转置矩阵,ZART为MR×MT维的矩阵,其第q行第p列的元素ZART(q,p)可由如下公式确定:
所述步骤4中根据vG与vL的端到端响应中的系统传输矩阵D,整个系统的空间-时间-频率相关函数由下式计算:
其中(·)*表示输入变量的共轭,Δr为空间间隔,Δt为时间间隔,Δf为频率间隔,Dqp(t,f)为矩阵D第q行第p列的元素在t时刻f频率处的值,为矩阵D第/>行第列的元素在t-Δt时刻f-Δf频率处的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111232635.0A CN113922900B (zh) | 2021-10-22 | 2021-10-22 | 一种考虑天线尺寸与天线互耦的无线电信道建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111232635.0A CN113922900B (zh) | 2021-10-22 | 2021-10-22 | 一种考虑天线尺寸与天线互耦的无线电信道建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113922900A CN113922900A (zh) | 2022-01-11 |
CN113922900B true CN113922900B (zh) | 2024-01-26 |
Family
ID=79242325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111232635.0A Active CN113922900B (zh) | 2021-10-22 | 2021-10-22 | 一种考虑天线尺寸与天线互耦的无线电信道建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113922900B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114584237B (zh) * | 2022-03-04 | 2024-01-26 | 东南大学 | 大规模天线阵列通信的空域非平稳无线信道容量计算方法 |
CN114598408B (zh) * | 2022-03-10 | 2024-02-27 | 东南大学 | 一种适用于全频段全场景的6g普适信道建模方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107171757A (zh) * | 2017-05-22 | 2017-09-15 | 电子科技大学 | 一种融合智能天线的mimo无线信道建模方法 |
CN109412723A (zh) * | 2017-08-16 | 2019-03-01 | 中兴通讯股份有限公司 | 一种mimo信道分析模型、建模方法及计算机可读存储介质 |
CN109639378A (zh) * | 2019-01-29 | 2019-04-16 | 南京信息工程大学 | 矩形隧道无线传播信道建模方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016174679A2 (en) * | 2015-04-27 | 2016-11-03 | Vayyar Imaging Ltd | System and methods for calibrating an antenna array using targets |
-
2021
- 2021-10-22 CN CN202111232635.0A patent/CN113922900B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107171757A (zh) * | 2017-05-22 | 2017-09-15 | 电子科技大学 | 一种融合智能天线的mimo无线信道建模方法 |
CN109412723A (zh) * | 2017-08-16 | 2019-03-01 | 中兴通讯股份有限公司 | 一种mimo信道分析模型、建模方法及计算机可读存储介质 |
CN109639378A (zh) * | 2019-01-29 | 2019-04-16 | 南京信息工程大学 | 矩形隧道无线传播信道建模方法 |
Non-Patent Citations (1)
Title |
---|
天线互耦对室内MIMO无线信道容量的影响;肖海林;欧阳缮;聂在平;李忻;;电子科技大学学报(04);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113922900A (zh) | 2022-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113922900B (zh) | 一种考虑天线尺寸与天线互耦的无线电信道建模方法 | |
Kailath | Sampling models for linear time-variant filters | |
CN110109051B (zh) | 基于频控阵的互耦阵列doa估计方法 | |
CN109714120B (zh) | 一种模拟耦合多天线室内空间衰落信道传播特性的方法 | |
CN109005133B (zh) | 双稀疏多径信道模型及基于此模型的信道估计方法 | |
Rabault et al. | On the tacit linearity assumption in common cascaded models of RIS-parametrized wireless channels | |
Liang et al. | Experimental characterization of UWB beamformers based on multidimensional beam filters | |
CN114124623A (zh) | 一种无线通信信道估计方法和装置 | |
CN104537171A (zh) | 一种mimo信道空间衰落相关性计算方法及多天线系统 | |
Li et al. | Distributed MIMO radar based on sparse sensing: Analysis and efficient implementation | |
Sanchez et al. | AirNN: Over-the-air computation for neural networks via reconfigurable intelligent surfaces | |
Palacios et al. | Multidimensional orthogonal matching pursuit: theory and application to high accuracy joint localization and communication at mmWave | |
Bai et al. | A mixed-bouncing based non-stationary model for 6G massive MIMO mmWave UAV channels | |
Wang et al. | Joint estimation of TOA and DOA in IR-UWB system using a successive MUSIC algorithm | |
CN103716802B (zh) | 基于双正交极化天线的多输入多输出宽带卫星移动通信信道建模方法 | |
Yu et al. | Joint Nuclear Norm and ℓ 1–2-Regularization Sparse Channel Estimation for mmWave Massive MIMO Systems | |
Li et al. | A novel method for DOA and time delay joint estimation in multipath OFDM environment | |
Wang et al. | A novel geometry-based MIMO channel model for vehicle-to-vehicle communication systems | |
Aquino et al. | A Review of Direction of Arrival Estimation Techniques in Massive MIMO 5G Wireless Communication Systems | |
Sun et al. | A channel prediction scheme with channel matrix doubling and temporal-spatial smoothing | |
Yao et al. | Massive MIMO channel modeling using map-based ray tracing method | |
Van Tien et al. | Hardware Implementation of a MIMO Channel Emulator for high speed WLAN 802.11 ac | |
Kontorovich | An Approach for Reconfigurable Intelligent Surface (RIS) design for MIMO-NOMA Transmission: theoretical study | |
Mai et al. | Codesign of beam pattern and sparse frequency waveforms for MIMO radar | |
Lan et al. | A novel DOA estimation of closely spaced sources using attention mechanism with conformal arrays |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |