发明名称 一种基于群时延效应和非线性约束的卫星信道建模方法
摘要 本发明公开了一种卫星信道的建模方法,将高功率放大器的非线性和群时延非线性效应有机结合,解决了现有技术信道模型对真实卫星信道的模拟跟踪性能低下的问题。本发明创新性地将卫星信道视为线性信道、非线性信道、群时延信道1和群时延信道2四个部分构成;每个部分分别由维纳模型的线性部分和非线性部分去逼近,当逼近过程收敛时,保留了维纳模型特性;当逼近过程不收敛时,定义一个控制函数作为非线性信道输入输出的约束函数替代维纳模型的非线性函数,以控制逼近过程收敛。与传统的卫星信道维纳模型相比,本发明提高了收敛速度、降低了均方误差,具有良好的动态跟踪性能,从而保证了卫星通信的效率和质量。
申请公布号 CN103812587B 申请公布日期 2015.09.23
申请号 CN201410048820.8 申请日期 2014.02.12
申请人 南京信息工程大学 发明人 郭业才;徐冉;赵卫娟;黄友锐
分类号 H04B17/391(2015.01)I 主分类号 H04B17/391(2015.01)I
代理机构 南京众联专利代理有限公司 32206 代理人 顾进;叶涓涓
主权项 一种基于群时延效应和非线性约束的卫星信道建模方法,其特征在于,包括如下步骤:步骤A,系统输入信号a(n)经线性信道得到其输出信号b(n);系统输入信号a(n)经群时延信道1得到其输出信号b<sub>1</sub>(n);系统输入信号a(n)经群时延信道2得到其输出信号b<sub>2</sub>(n);其中,n为整数,表示时间序列;步骤B,由步骤A所述的系统输入信号a(n)、线性信道输出信号b(n)和群时延信道1输出信号b<sub>1</sub>(n),经运算器1运算得非线性信道输入信号x(n):x(n)=b(n)b<sub>1</sub>(n)/a(n);步骤C,步骤B所述的非线性信道输入信号x(n)经非线性信道,得其输出信号y(n);步骤D,步骤A所述的系统输入信号a(n)及群时延信道2输出信号b<sub>2</sub>(n)和步骤C所述的非线性信道输出信号y(n),经运算器2运算得卫星信道最终输出信号z(n):z(n)=y(n)b<sub>2</sub>(n)/a(n);所述线性信道由发射滤波器、上行链路通信环境及卫星内部高功率放大器中的无群时延线性部分级联;非线性信道由发射滤波器、上行链路通信环境及卫星内部高功率放大器中的无群时延非线性部分级联;群时延信道1表示发射滤波器、高功率放大器线性部分及上行链路通信环境引起的具有群时延效应的非线性信道;群时延信道2表示发射滤波器、高功率放大器非线性部分引起的具有群时延效应的非线性信道;所述的线性信道由维纳模型的线性部分表示为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>b</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>0</mn></msub><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>0</mn></msub></munderover><msub><mi>c</mi><mrow><mn>0</mn><msub><mi>m</mi><mn>0</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>a</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mi>m</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000758520750000011.GIF" wi="580" he="148" /></maths>式中,m<sub>0</sub>=1,…,M<sub>0</sub>,M<sub>0</sub>为正整数;<img file="FDA0000758520750000012.GIF" wi="153" he="83" />为线性信道C<sub>0</sub>(n)的第m<sub>0</sub>个延时抽头系数;所述的群时延信道1、群时延信道2和非线性信道均由维纳模型的非线性部分表示如下:其中,群时延信道1由维纳模型的非线性部分表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>b</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>1</mn></msub><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><msub><mi>c</mi><mrow><mn>1</mn><msub><mi>m</mi><mn>1</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>a</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>m</mi><mn>1</mn></msub></msup><mo>,</mo></mrow>]]></math><img file="FDA0000758520750000013.GIF" wi="542" he="156" /></maths>式中,m<sub>1</sub>=1,…,M<sub>1</sub>,M<sub>1</sub>为正整数;<img file="FDA0000758520750000014.GIF" wi="145" he="82" />为群时延信道1权向量C<sub>1</sub>(n)的第m<sub>1</sub>个延时抽头系数;群时延信道2由维纳模型的非线性部分表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>b</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>2</mn></msub><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>2</mn></msub></munderover><msub><mi>c</mi><mrow><mn>2</mn><msub><mi>m</mi><mn>2</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>a</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>m</mi><mn>2</mn></msub></msup><mo>,</mo></mrow>]]></math><img file="FDA0000758520750000015.GIF" wi="564" he="155" /></maths>式中,m<sub>2</sub>=1,…,M<sub>2</sub>,M<sub>2</sub>为正整数;<img file="FDA0000758520750000021.GIF" wi="156" he="79" />为群时延信道2权向量C<sub>2</sub>(n)的第m<sub>2</sub>个延时抽头系数;非线性信道调幅‑调幅效应由维纳模型的非线性部分表示为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>G</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>3</mn></msub><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>3</mn></msub></munderover><msub><mi>c</mi><mrow><mn>3</mn><msub><mi>m</mi><mn>3</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>&rho;</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>M</mi><mn>3</mn></msub></msup><mo>,</mo></mrow>]]></math><img file="FDA0000758520750000022.GIF" wi="579" he="156" /></maths>式中,ρ(n)=|x(n)|<sup>2</sup>为非线性信道输入x(n)的幅度平方;m<sub>3</sub>=1,…,M<sub>3</sub>,M<sub>3</sub>为正整数;G(n)表示非线性信道的幅度,<img file="FDA0000758520750000023.GIF" wi="154" he="83" />为非线性信道调幅‑调幅效应权向量C<sub>3</sub>(n)的第m<sub>3</sub>个延时抽头系数;非线性信道调幅‑调相效应由维纳模型的非线性部分表示为:<img file="FDA0000758520750000024.GIF" wi="581" he="154" />式中,<img file="FDA0000758520750000025.GIF" wi="110" he="75" />表示非线性信道的相位;m<sub>4</sub>=1,…,M<sub>4</sub>,M<sub>4</sub>为正整数;<img file="FDA0000758520750000026.GIF" wi="157" he="81" />为非线性信道调幅‑调相效应权向量C<sub>4</sub>(n)的第m<sub>4</sub>个延时抽头系数;非线性信道输入x(n)‑输出y(n)间的关系表示为<img file="FDA0000758520750000027.GIF" wi="1313" he="187" />式中,<img file="FDA0000758520750000028.GIF" wi="179" he="86" />为虚数单位,下同;当线性信道由维纳模型的线性部分,群时延信道1、群时延信道2和非线性信道均由维纳模型的非线性部分表示时,若这种模型在逼近真实卫星信道的过程不收敛时,需对群时延信道1、群时延信道2和非线性信道的输入与输出间的关系进行约束,约束函数定义为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup></mrow><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000758520750000029.GIF" wi="447" he="195" /></maths>式中,α(n)表示输入,f(n)表示输出;对群时延信道1,α(n)表示系统输入a(n),f(n)表示群时延信道1的输出b<sub>1</sub>(n);对群时延信道2,α(n)表示系统输入a(n),f(n)表示群时延信道2的输出b<sub>2</sub>(n);对非线性信道,α(n)表示非线性信道输入x(n),f(n)表示非线性信道输出y(n);使用约束函数后,非线性信道输入x(n)为<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>0</mn></msub><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>0</mn></msub></munderover><msub><mi>c</mi><mrow><mn>0</mn><msub><mi>m</mi><mn>0</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>a</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mi>m</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><msub><mi>c</mi><mrow><mn>1</mn><msub><mi>m</mi><mn>1</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>a</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>m</mi><mn>1</mn></msub></msup></mrow></msup><mo>)</mo></mrow></mrow><mrow><mi>a</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><msub><mi>c</mi><mrow><mn>1</mn><msub><mi>m</mi><mn>1</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>a</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>m</mi><mn>1</mn></msub></msup></mrow></msup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000758520750000031.GIF" wi="856" he="379" /></maths>卫星信道总输出为<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>z</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>0</mn></msub><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>0</mn></msub></munderover><msub><mi>c</mi><mrow><mn>0</mn><msub><mi>m</mi><mn>0</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>a</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mi>m</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><msub><mi>c</mi><mrow><mn>1</mn><msub><mi>m</mi><mn>1</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>a</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>M</mi><mn>1</mn></msub></msup></mrow></msup><mo>)</mo></mrow></mrow><mrow><mi>a</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><msub><mi>c</mi><mrow><mn>1</mn><msub><mi>m</mi><mn>1</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>a</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>M</mi><mn>1</mn></msub></msup></mrow></msup><mo>)</mo></mrow></mrow></mfrac><mo>&CenterDot;</mo><mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>2</mn></msub><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>2</mn></msub></munderover><msub><mi>c</mi><mrow><mn>2</mn><msub><mi>m</mi><mn>2</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>a</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>M</mi><mn>2</mn></msub></msup></mrow></msup><mo>)</mo></mrow><mrow><mi>a</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>2</mn></msub><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>2</mn></msub></munderover><msub><mi>c</mi><mrow><mn>2</mn><msub><mi>m</mi><mn>2</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>a</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>M</mi><mn>2</mn></msub></msup></mrow></msup><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>3</mn></msub><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>3</mn></msub></munderover><msub><mi>c</mi><mrow><mn>3</mn><msub><mi>m</mi><mn>3</mn></msub></mrow></msub><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>&rho;</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>M</mi><mn>3</mn></msub></msup><mo>)</mo><msup><mi>e</mi><mrow><mi>j</mi><munderover><mo>&Sigma;</mo><mrow><msub><mi>m</mi><mn>4</mn></msub><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>4</mn></msub></munderover><msub><mi>c</mi><mrow><mn>4</mn><msub><mi>m</mi><mn>4</mn></msub></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mrow><mo>(</mo><mi>&rho;</mi><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>M</mi><mn>4</mn></msub></msup></mrow></msup></mrow></mtd></mtr></mtable><mo>.</mo></mrow>]]></math><img file="FDA0000758520750000032.GIF" wi="1580" he="694" /></maths>
地址 210044 江苏省南京市宁六路219号