发明名称 气象卫星信道概率统计特性仿真模型的建立方法
摘要 本发明提供了一种气象卫星信道概率统计特性仿真模型的建立方法,包括下列步骤:一)计算气象卫星发送信号;二)根据气象卫星发射的信号在不同的天气状态下因多径效应、散射效应和阴影遮蔽效应等,使得到达地面接收站点的接收信号由不同成分组成,推导相应接收信号的包络概率密度分布,根据不同天气状态下的气象卫星地面接收系统接收的信号表达形式建立对应的仿真模型。本发明建立的仿真模型物理意义明确,易于实现,可以用于气象卫星通信系统的仿真,为优化气象卫星通信传输体制,设计气象卫星通信系统和通信技术提供实用、方便和有效的概率统计仿真模型。
申请公布号 CN102255647A 申请公布日期 2011.11.23
申请号 CN201110162630.5 申请日期 2011.06.16
申请人 南京信息工程大学 发明人 张秀再;郭业才
分类号 H04B7/185(2006.01)I;H04L12/24(2006.01)I 主分类号 H04B7/185(2006.01)I
代理机构 南京众联专利代理有限公司 32206 代理人 顾进
主权项 1.气象卫星信道概率统计特性仿真模型的建立方法,其特征在于,包括下列步骤:一)计算气象卫星发送信号,信号调制方式可采用BPSK或QPSK的调制方式,气象卫星发送信号可以表示为s(t)=Acos(ω<sub>c</sub>t+θ)                    (1)式中,A为发射信号的振幅,ω<sub>c</sub>为载波信号的角频率,θ根据调制方式的不同取值也不同;二)根据气象卫星发射的信号s(t)=Acos(ω<sub>c</sub>t+θ<sub>k</sub>),k=4,采用QPSK调制方式,在不同的天气状态下因多径效应、散射效应和阴影遮蔽效应等,使得到达地面接收站点的接收信号由不同成分组成,推导相应的包络概率密度分布,根据不同天气状态下的气象卫星地面接收系统接收的信号表达形式,建立对应的仿真模型:1)阴雨天气时,气象卫星信道物理环境较差,假设地面接收站的接收信号只包含一定强度的多径分量,接收信号的包络和相位都是随机缓变的,则接收信号r(t)可表示为r(t)=a(t)cos[ω<sub>c</sub>t+θ(t)]=μ<sub>c</sub>(t)cosω<sub>c</sub>t-μ<sub>s</sub>(t)sinω<sub>c</sub>t    (2)式中μ<sub>c</sub>(t)=a(t)cos[θ(t)]           (3)μ<sub>s</sub>(t)=a(t)sin[θ(t)]           (4)式中,a(t)和θ(t)分别是r(t)的随机包络和随机相位,ω<sub>c</sub>为载波频率,a(t)≥0,由式(2)知:接收信号r(t)的统计特性可以由μ<sub>c</sub>(t)和μ<sub>s</sub>(t)的统计特性确定,由中心极限定理可知,μ<sub>c</sub>(t)和μ<sub>s</sub>(t)为正态随机过程,其概率密度函数分别为<maths num="0001"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>&sigma;</mi><mi>c</mi></msub><msqrt><mn>2</mn><mi>&pi;</mi></msqrt></mrow></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><msub><mi>&mu;</mi><mi>c</mi></msub><mn>2</mn></msup><msup><msub><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>c</mi></msub><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>&sigma;</mi><mi>s</mi></msub><msqrt><mn>2</mn><mi>&pi;</mi></msqrt></mrow></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><msub><mi>&mu;</mi><mi>s</mi></msub><mn>2</mn></msup><msup><msub><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>s</mi></msub><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>进而推导得包络a(t)的一维概率密度函数为<maths num="0003"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mo>&infin;</mo></mrow><mo>&infin;</mo></msubsup><mi>f</mi><mrow><mo>(</mo><mi>a</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>d&theta;</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></msubsup><mfrac><mi>a</mi><msubsup><mrow><mn>2</mn><mi>&pi;&sigma;</mi></mrow><mi>r</mi><mn>2</mn></msubsup></mfrac><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac><msup><mi>a</mi><mn>2</mn></msup><msubsup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>r</mi><mn>2</mn></msubsup></mfrac><mo>]</mo><mi>d&theta;</mi><mo>=</mo><mfrac><mi>a</mi><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup></mfrac><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac><msup><mi>a</mi><mn>2</mn></msup><msubsup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>r</mi><mn>2</mn></msubsup></mfrac><mo>]</mo><mo>,</mo><mi>a</mi><mo>&GreaterEqual;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>该式表明,接收信号包络a(t)的概率密度服从Rayleigh分布,由式(2)、(3)、(4)、(5)和(6)可以建立图2Rayleigh信道仿真模型;2)多云天气时,因云层较厚和云量较多,云在地面接收站上空一定的范围散开分布形成部分阴影遮挡,地面接收站接收信号如果含有直射信号和一定强度的多径信号,假设直射信号只有时延没有衰减,则接收的混合信号为r(t)=Acos[ω<sub>c</sub>(t-τ)+θ]+r′(t)=Acos[ω<sub>c</sub>t-ω<sub>c</sub>τ+θ]+r′(t)=Acos[ω<sub>c</sub>t+θ′]+r′(t)                (8)式中,τ表示直射信号的时延值,θ′是直射信号经时延到接收站点的相位,r′(t)如式(2)所示。将式(8)展开推导得r(t)的包络υ(t)的概率密度函数为<maths num="0004"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>&upsi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>&upsi;</mi><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><msup><mi>&upsi;</mi><mn>2</mn></msup><mo>+</mo><msup><mi>A</mi><mn>2</mn></msup></mrow><msubsup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>r</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>I</mi><mn>0</mn></msub><mrow><mo>(</mo><mfrac><mi>A&upsi;</mi><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>,</mo><mi>&upsi;</mi><mo>&GreaterEqual;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>即接收信号r(t)的包络υ(t)的概率密度服从Rice分布,由式(8)可以建立图3Rice信道仿真模型;3)多云天气时,因云层较厚和云量较多,在地面接收站上空一定的范围分布较厚的云层形成完全阴影遮挡,云的垂直密度是非均匀分布的,可以在云的垂直方向上划分N层不同密度的层结,而传播信号每经过一个层结都会因阴影遮蔽效应造成不同程度的衰减,且衰减是时变的,假设云的第i层结阴影遮蔽因子是K<sub>i</sub>(t),则整个云层的阴影遮蔽因子<img file="FDA0000068829680000023.GIF" wi="327" he="121" />该阴影遮蔽因子是信道的乘性干扰,假设发射信号到达云顶时没有衰减只有时延τ,地面接收站接收的信号如果只含有经较厚云层阴影遮蔽后的直射信号,则接收信号可表示为r(t)=k(t)·s(t-τ)=k(t)Acos(ω<sub>c</sub>t+θ′)=ξ<sub>c</sub>(t)cosω<sub>c</sub>t-ξ<sub>s</sub>(t)sinω<sub>c</sub>t    (10)式中ξ<sub>c</sub>(t)=k(t)Acosθ′                            (11)ξ<sub>s</sub>(t)=k(t)Asinθ′                            (12)则r(t)的包络和相位分别为<maths num="0005"><![CDATA[<math><mrow><mi>&xi;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msqrt><msubsup><mi>&xi;</mi><mi>c</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&xi;</mi><mi>s</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msqrt><mo>=</mo><mi>Ak</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>&xi;</mi><mo>&GreaterEqual;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths><img file="FDA0000068829680000032.GIF" wi="1313" he="122" />由式(13)推导出<maths num="0006"><![CDATA[<math><mrow><mi>ln</mi><mi>&xi;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>ln</mi><mi>A</mi><mo>+</mo><mi>ln</mi><mi>k</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>ln</mi><mi>A</mi><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>ln</mi><msub><mi>K</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow></math>]]></maths>由中心极限定理可知,<img file="FDA0000068829680000034.GIF" wi="211" he="119" />服从正态分布,进而可以推导r(t)包络的概率密度函数<maths num="0007"><![CDATA[<math><mrow><msub><mi>f</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>&sigma;&xi;</mi><msqrt><mn>2</mn><mi>&pi;</mi></msqrt></mrow></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>ln</mi><mi>&xi;</mi><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><mn>2</mn></msup><msup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow></math>]]></maths>即r(t)包络的概率密度服从Lognormal分布,由式(10)可以建立图4所示Lognormal信道仿真模型;4)晴朗天气时,假设地面接收信号主要包括没有衰减只有时延的直射信号和很弱的多径信号,则接收信号表示为r(t)=s(t-τ)+r″(t)≈s(t-τ)                    (17)式中,r″(t)是很弱的多径信号,即接收信号的信噪比很大,因为与Rice分布时的接收信号的组成分量类似,仅信噪比相对要大很多,则由式(9)可知,当信噪比<img file="FDA0000068829680000036.GIF" wi="190" he="127" />很大时,有<img file="FDA0000068829680000037.GIF" wi="299" he="118" />则式(9)可为<maths num="0008"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>&upsi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>&upsi;</mi><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><msup><mi>&upsi;</mi><mn>2</mn></msup><mo>+</mo><msup><mi>A</mi><mn>2</mn></msup></mrow><msubsup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>r</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><mn>1</mn><msqrt><mn>2</mn><mi>&pi;</mi><mfrac><mi>A&upsi;</mi><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup></mfrac></msqrt></mfrac><mi>exp</mi><mrow><mo>(</mo><mfrac><mi>A&upsi;</mi><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>=</mo><mfrac><mi>&upsi;</mi><mrow><msqrt><mn>2</mn><mi>&pi;A&upsi;</mi></msqrt><mo>&CenterDot;</mo><msub><mi>&sigma;</mi><mi>r</mi></msub></mrow></mfrac><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>&upsi;</mi><mo>-</mo><mi>A</mi><mo>)</mo></mrow><mn>2</mn></msup><msubsup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>r</mi><mn>2</mn></msubsup></mfrac><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow></math>]]></maths>当υ≈A时,f(υ)可近似表示为高斯分布,即接收信号包络的概率密度函数为<maths num="0009"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>&upsi;</mi><mo>)</mo></mrow><mo>&ap;</mo><mfrac><mn>1</mn><mrow><msqrt><mn>2</mn><mi>&pi;</mi></msqrt><mo>&CenterDot;</mo><msub><mi>&sigma;</mi><mi>r</mi></msub></mrow></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>&upsi;</mi><mo>-</mo><mi>A</mi><mo>)</mo></mrow><mn>2</mn></msup><msubsup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>r</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>,</mo><mi>&upsi;</mi><mo>&GreaterEqual;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mrow></math>]]></maths>即r(t)包络的概率密度服从Gauss分布,根据Gauss信道统计特性模型,将Rice信道仿真模型中的多径信号的能量取小就可以实现对Gauss信道的仿真。
地址 210044 江苏省南京市浦口区盘城新街114号