专利名称：一种检测bpsk卫星通信系统误码性能恶化程度的方法
技术领域：
本发明涉及卫星通信系统误码性能检测技术领域，特别是涉及一种对二进制相移键控(BPSK)卫星通信系统在受单频干扰情况下误码性能恶化程度进行检测的方法。
背景技术：
单载干扰对采用BPSK调制体制的卫星通信系统误码性能恶化的影响，一般体现在载波的跟踪性能恶化和检测判决性能恶化两个方面，在单频干扰强度小于信号载波强度的情况下，其影响主要体现在检测判决性能恶化方面。
现有对BPSK卫星通信系统误码性能恶化程度的检测方法主要分为三类第一类是将单频干扰折算为在接收频带内的等效噪声分量，然后对该等效噪声分量进行检测分析；第二类将单频干扰与接收噪声区分开来，但是却将单频干扰相位作为确定量进行检测；第三类虽然将干扰相位作为随机量引入检测分析中，但是却仅仅给出了定性分析结果，没有形成量化的检测方法。
上述三类方法均存在不同程度的缺陷，会直接导致检测结果产生偏差甚至出现错误，影响工程应用中的抗干扰设计分析。首先第一类方法存在理论上的漏洞，BPSK通信系统无干扰时的误码率计算公式是在假定系统噪声满足高斯分布时得到的统计分析结果，将单频干扰等效为噪声分量后，系统噪声并不满足高斯分布，因此再采用原有的误码率计算公式，理论上不严谨。其次，第二类方法仅适用于干扰与接收端距离不变的通信系统，而在卫星通信系统中信号传输距离是不断改变的，将会导致干扰信号相位不断变化，因此将单频信号相位作为确定值的分析方法显然不具有一般性。最后，第三类方法没有给出量化的检测方法，在工程实践中不具备实用性。
鉴于上述现有检测方法的缺陷，本发明人基于从事绕月探测工程数据接收系统研制和设计的经验和心得，并配合通信原理及概率论的知识运用，积极加以研究创新，以期创设一种准确的、具有一般性和实用性的对BPSK卫星通信系统受单频干扰时的误码性能恶化程度进行检测的方法，经过不断的演算、推理分析终于创设出本发明。
发明内容
(一)要解决的技术问题 有鉴于此，本发明的主要目的在于提供一种检测BPSK卫星通信系统误码性能恶化程度的方法，以实现对系统误码性能恶化程度的精确检测。
(二)技术方案 为达到上述目的，本发明提供了一种检测BPSK卫星通信系统误码性能恶化程度的方法，该方法包括 分别计算BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率； 根据该BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率，计算BPSK卫星通信系统受单频干扰时的误码性能恶化量； 比较该误码性能恶化量与预先设定的误码率恶化量门限值，如果该误码性能恶化量大于该门限值，则判定误码性能重度恶化。
上述方案中，所述计算BPSK卫星通信系统受单频干扰时的误码率包括将BPSK卫星通信系统中的载波与噪声频率差Δf、码速率Rb、信道噪声功率谱密度N0、接收信噪比Eb/N0以及载干比(A/Ai)2代入公式
计算BPSK卫星通信系统受单频干扰时的误码率Pe，其中
Δf≠0时
为BPSK调制信号的幅度，Ai为单载波干扰信号的幅度，
T＝1/Rb为码元宽度，Δf＝fc-fi为载波与干扰信号的频率差，Eb表示每个码元周期内的信号能量，Rb为传输码速率，N0为信道噪声功率谱密度。
上述方案中，所述计算BPSK卫星通信系统无单频干扰时的误码率包括将BPSK卫星通信系统中的系统信噪比Eb/N0代入公式
计算BPSK卫星通信系统无单频干扰时的误码率P0，其中
Eb为每个码元周期内的信号能量，N0为信道噪声功率谱密度。
上述方案中，所述根据该BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率，计算BPSK卫星通信系统受单频干扰时的误码性能恶化量包括将该BPSK卫星通信系统受单频干扰时的误码率Pe和该无单频干扰时的误码率P0代入公式H＝lg(Pe/P0)，计算得到BPSK卫星通信系统受单频干扰时的误码性能恶化量H。
上述方案中，所述比较该误码性能恶化量与预先设定的误码率恶化量门限值，进一步包括如果该误码性能恶化量小于或等于该门限值，则判定误码性能没有重度恶化，无需针对此单频干扰对该BPSK卫星通信系统进行抗干扰处理。
上述方案中，所述误码性能恶化程度不能被容忍，该方法进一步包括针对此单频干扰对该BPSK卫星通信系统进行抗干扰处理。
(三)有益效果 从上述技术方案可以看出，本发明具有以下有益效果 1、本发明提供的这种检测BPSK卫星通信系统误码性能恶化程度的方法，完善和改进了现有BPSK卫星通信系统受单频干扰时误码性能恶化检测的方法，提供了一种准确的、具有一般性和实用性的检测方法。
2、本发明提供的这种检测BPSK卫星通信系统误码性能恶化程度的方法，不是简单的将单频干扰折算为接收频带内的等效噪声进行检测分析，而是将单频干扰作为系统接收信号的一部分独立进行考虑，并将其相位作为随机变量，采用最大似然估计，统计在检测判决时刻的误码概率，通过比较无干扰和有干扰时的误码率，实现了对系统误码性能恶化程度的精确检测。该方法分析更为严谨，结果更加准确可信，并且更具实用性。
3、本发明提供的这种检测BPSK卫星通信系统误码性能恶化程度的方法，成功应用在了绕月探测工程中，为我国与日本探月卫星信号频率协调工作以及地面接收站的抗干扰设计提供了依据，对以后深空卫星通信中类似情况的分析处理提供依据。
图1是本发明提供的检测BPSK卫星通信系统误码性能恶化程度的理论分析框图； 图2是本发明提供的检测BPSK卫星通信系统误码性能恶化程度的方法流程图。
具体实施例方式 为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。
首先，介绍本发明的实现原理。考虑系统接收信号由BPSK调制信号、高斯白噪声和单频干扰三部分组成，其中BPSK调制信号载波幅度和频率、噪声功率谱密度、单频干扰幅度和频率为确定量。另外，不仅将信道噪声作为按高斯分布的随机变量，而且将单频干扰信号相位作为均匀分布的随机变量统一进行概率分析，采用最大似然估计，统计在检测判决时刻的错误概率，从而得到系统在不同信噪比、干扰频率和干扰强度时的误码率；然后将该误码率与系统无干扰时的误码率进行比对，得到当前干扰情况下的误码性能恶化量，为系统抗干扰设计提供依据。
请参照图1，图1是本发明提供的检测BPSK卫星通信系统误码性能恶化程度的理论分析框图。先考虑存在单频干扰的情况。不失一般性，将发送信息码元近似为随机序列进行考虑。基带波形为双极性非归零码；BPSK调制后的发送波形用si(t)表示——幅度为A，载频为fc，码速率为Rb的BPSK信号(码元宽度T＝1/Rb)。
其中
Eb表示每个码元周期内si(t)的能量。
将干扰信号表示为幅度为Ai，频率为fi的单载波信号(fi在系统带通内)，用I(t)表示
其中
是在(0，2π)均匀分布的随机变量。
n(t)是信道的加性高斯白噪声，双边功率谱密度为N0/2。
上述三个信号同时被送入BPSK接收机进行相干解调和检测得到发送的信息码元。
如说明书附图1，当发送端发送s1(t)时，接收机前端接收信号为
与本地载波相乘后输出信号为
令Δf＝fc-fi，则 1)、Δf＝0时 在t＝nT时(n是自然数)时刻的采样值为
由于fc和fc+fi很大，上式近似为
其中第一项为所需的信号分量，第二项为单频信号干扰项，第三项为信道高斯白噪声分量。
令随机变量
由概率论可得，X的概率分布函数为 Y满足均值
方差
的正态分布。
因此发送s1(t)时的误码概率P(e/s1)为 令
X、Y的概率密度函数为
由于X、Y独立，则
当y＞b时，P(e/s1)＝0； 当y＜-b时，
其中 当-b≤y≤b时， 最后可得 同理可得，当发送s0(t)时的误码概率为 2)，Δf≠0时，有 在t＝nT时(n是自然数)时刻的采样值为
由于fc，fc+fi很大，上式近似为
令
与1)中推理相同，可得 发送s1(t)时的误码概率为 发送s0(t)时的误码概率为 由于发送s1(t)和s0(t)的统计概率相等，综上所述接收误码率为 其中Δf＝0时，
Δf≠0时，
式(13)即为BPSK卫星通信系统受单频干扰时的误码率计算公式，根据这个公式，通过接收系统和干扰信号的相关参数即可求出当前情况的下的系统接收误码率。
而BPSK通信系统无干扰时的误码率为 因此BPSK卫星通信系统受单频干扰时的误码性能恶化量为 H＝lg(Pe/P0) (15) 若设定系统所能容忍的误码率恶化量为H0，则当H＞H0时，系统需要针对此单频干扰进行抗干扰处理。
基于上述本发明的实现原理，图2示出了本发明提供的检测BPSK卫星通信系统误码性能恶化程度的方法流程图，该方法包括以下步骤 步骤1分别计算BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率； 步骤2根据该BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率，计算BPSK卫星通信系统受单频干扰时的误码性能恶化量； 步骤3比较该误码性能恶化量与预先设定的误码率恶化量门限值，如果该误码性能恶化量大于该门限值，则判定误码性能重度恶化，针对此单频干扰对该BPSK卫星通信系统进行抗干扰处理；如果该误码性能恶化量小于或等于该门限值，则判定误码性能没有重度恶化，无需针对此单频干扰对该BPSK卫星通信系统进行抗干扰处理。
其中，步骤1中所述计算BPSK卫星通信系统受单频干扰时的误码率包括将BPSK卫星通信系统中的载波与噪声频率差Δf、码速率Rb、信道噪声功率谱密度N0、接收信噪比Eb/N0以及载干比(A/Ai)2代入公式
计算BPSK卫星通信系统受单频干扰时的误码率Pe，其中Δf＝0时
Δf≠0时
为BPSK调制信号的幅度，Ai为单载波干扰信号的幅度，
T＝1/Rb为码元宽度，Δf＝fc-fi为载波与干扰信号的频率差，Eb表示每个码元周期内的信号能量，Rb为传输码速率，N0为信道噪声功率谱密度。
步骤1中所述计算BPSK卫星通信系统无单频干扰时的误码率包括将BPSK卫星通信系统中的系统信噪比Eb/N0代入公式
计算BPSK卫星通信系统无单频干扰时的误码率P0，其中
Eb为每个码元周期内的信号能量，N0为信道噪声功率谱密度。
步骤2中所述根据该BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率，计算BPSK卫星通信系统受单频干扰时的误码性能恶化量包括将该BPSK卫星通信系统受单频干扰时的误码率Pe和无单频干扰时的误码率P0代入公式H＝lg(Pe/P0)，计算得到BPSK卫星通信系统受单频干扰时的误码性能恶化量H。
本发明成功应用在了绕月探测工程中，为我国与日本探月卫星信号频率协调工作以及地面接收站的抗干扰处理提供了依据。下面以此为例对本发明进行详细阐述。
近年来，随着世界各国运载火箭、远距离测控以及空间探测设备等技术的日趋成熟，深空探测活动日益频繁。出现了好几个国家的卫星同时对同一个天体进行探测的现象。这就使得各国的卫星与地面站的通信存在相互干扰的可能。
2007年9月，日本研制的SELENE探月卫星(以下简称“日星”)发射成功，SELENE探月卫星由一颗母卫星和两颗子卫星构成，两颗子卫星分别为“中继星”(Rstar)和“干涉测量卫星”(Vstar)。其中Rstar工作在距月球100km×2400km、倾角为90°的椭圆轨道上，执行在母卫星转到月球背面时与地面的数据中继任务。Vstar则工作在距月球100km×800km、倾角为90°的椭圆轨道上，与母卫星一起执行反演月球重力场等任务。
同年10月底，我国研制的绕月探测工程“嫦娥一号(CE-1)”探月卫星顺利升空。它运行在距月球200km的极轨道上，执行月面三维地形测绘、探测月球元素分布、评估月壤厚度及监测空间环境等任务。
CE-1的TD3波束作为探测数据下行的唯一途径，对整个绕月探测工程至关重要。而日星向国际电联上报的VS3波束频率落在了TD3传输频带内，将会对CE-1下行探测数据的接收造成干扰。
因此，我们需要在仔细研究卫星信号和接收系统各项参数的基础上，就日星干扰对CE-1数据接收系统造成的误码率恶化情况进行评估。然后根据其评估结果，与日方进行卫星频率协调或对CE-1数据接收系统进行抗干扰处理。
CE-1的TD3波束采用BPSK调制体制，用于探测数据下行，中心频率2290MHz，带宽6MHz。根据日本向国际电联申报的资料，日星VS3波束采用单载波形式，用于测控下行，中心频率2287.3MHz，带宽150kHz，见下表1，表1是中日卫星信号参数。则该问题实际上等效于评估高斯白噪声背景下的BPSK数据传输系统，在受单载波干扰时的误码率恶化情况。
表1 因此由中日卫星的信号参数可得，两者频率差约为2.7MHz，而在最坏情况下的载干比约为34.5dB。CE-1数据接收系统要求在3Mbps码速率下，接收误码率优于1×10-6，因此取系统的门限信噪比10.53dB。带入(13)、(14)分别计算无干扰和有干扰时的误码率。采用MALAB数字计算工具，其计算程序如下 function[theoryPe，Pe]＝bpskpe(deltaf，SNR，CNI，N0，Rb) ％计算在信噪比SNR时的理论误码率 ％计算在fc-fi＝deltaf，信噪比SNR信干比SNI时的BPSK接收机误码率 ％deltaf表示信号与干扰频率差(Hz)，SNR表示接收系统无干扰时信噪比(dB)，CNI表示载噪比(dB)，N0表示解调机输入端噪声谱密度(mw/Hz)，Rb表示传输码速率 pi＝3.1415； ％圆周率 T＝1/Rb；％信号码元周期 ％信号功率 Eb＝N0*10^(0.1*SNR)；％单位(瓦) A＝sqrt(2*Eb/T)；％信号幅度 ％干扰功率 Ai＝sqrt(A^2/10^(0.1*CNI))； ％变量计算 ceta＝sqrt(Eb*N0/2)；％方差 u＝A*A*T/2； ％均值 ％载波与干扰信号的频率差(fc＝fi) it deltaf＝＝0； d＝A*Ai*T/2； else d＝A*Ai*sin(pi*deltaf*T)/(2*pi*deltaf)； end ％无干扰理论误码率 theoryPe＝0.5*erfc(u/(ceta*1.4142))； ％有干扰时误码率 P1＝0.5*erfc((d+u)/(ceta*1.4142))； fun1＝@(y)(pi-acos(-y/d)).*exp(-0.5*((y-u)/ceta).^2)； P2＝quadl(fun1，-d，d，1e-9)； fun2＝@(y)acos(-y/d).*exp(-0.5*((y+u)/ceta).^2)； P3＝quadl(fun2，-d，d，1e-9)； Pe＝P1+0.5/(pi*ceta*sqrt(2*pi))*(P2+P3)； 可得接收误码率从无干扰时的P0＝9.9932×10-7恶化到Pe＝9.9988×10-7，造成的误码率恶化量为H＝2.4×10-4dB。由于CE-1数据接收系统对可靠性要求较高，因此将系统可容忍的误码性能恶化门限设为H0＝0.1dB。很明显H＜H0，说明日星干扰对CE-1数据接收系统的影响极其微小，不用对系统另外进行抗干扰处理。
综上所述，根据此MATLAB计算程序，输入系统接收信噪比、噪声谱密度、传输码速率、发送载波与干扰信号的频率差以及载干比，可以方便的求出BPSK通信系统无干扰时的理论接收误码率和受单频干扰时的理论接收误码率，从而对系统的误码恶化情况进行检测，具有一定的应用价值。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
权利要求
1.一种检测二进制相移键控卫星通信系统误码性能恶化程度的方法，其特征在于，该方法包括
分别计算二进制相移键控BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率；
根据该BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率，计算BPSK卫星通信系统受单频干扰时的误码性能恶化量；
比较该误码性能恶化量与预先设定的误码率恶化量门限值，如果该误码性能恶化量大于该门限值，则判定误码性能重度恶化。
2.根据权利要求1所述的检测二进制相移键控卫星通信系统误码性能恶化程度的方法，其特征在于，所述计算BPSK卫星通信系统受单频干扰时的误码率包括
将BPSK卫星通信系统中的载波与噪声频率差Δf、码速率Rb、信道噪声功率谱密度N0、接收信噪比Eb/N0以及载干比(A/Ai)2代入公式计算BPSK卫星通信系统受单频干扰时的误码率Pe，其中Δf＝0时
Δf≠0时
为BPSK调制信号的幅度，Ai为单载波干扰信号的幅度，
T＝1/Rb为码元宽度，Δf＝fc-fi为载波与干扰信号的频率差，Eb表示每个码元周期内的信号能量，Rb为传输码速率，N0为信道噪声功率谱密度。
3.根据权利要求1所述的检测二进制相移键控卫星通信系统误码性能恶化程度的方法，其特征在于，所述计算BPSK卫星通信系统无单频干扰时的误码率包括
将BPSK卫星通信系统中的系统信噪比Eb/N0代入公式
计算BPSK卫星通信系统无单频干扰时的误码率P0，其中
Eb为每个码元周期内的信号能量，N0为信道噪声功率谱密度。
4.根据权利要求1所述的检测二进制相移键控卫星通信系统误码性能恶化程度的方法，其特征在于，所述根据该BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率，计算BPSK卫星通信系统受单频干扰时的误码性能恶化量包括
将该BPSK卫星通信系统受单频干扰时的误码率Pe和该无单频干扰时的误码率P0代入公式H＝lg(Pe/P0)，计算得到BPSK卫星通信系统受单频干扰时的误码性能恶化量H。
5.根据权利要求1所述的检测二进制相移键控卫星通信系统误码性能恶化程度的方法，其特征在于，所述比较该误码性能恶化量与预先设定的误码率恶化量门限值，进一步包括
如果该误码性能恶化量小于或等于该门限值，则判定误码性能没有重度恶化，无需针对此单频干扰对该移相键控卫星通信系统进行抗干扰处理。
6.根据权利要求1所述的检测二进制相移键控卫星通信系统误码性能恶化程度的方法，其特征在于，所述判定误码性能重度恶化，该方法进一步包括
针对此单频干扰对该移相键控卫星通信系统进行抗干扰处理。
全文摘要
本发明公开了一种检测BPSK卫星通信系统误码性能恶化程度的方法，包括分别计算BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率；根据该BPSK卫星通信系统受单频干扰时的误码率和无单频干扰时的误码率，计算BPSK卫星通信系统受单频干扰时的误码性能恶化量；比较该误码性能恶化量与预先设定的误码率恶化量门限值，如果该误码性能恶化量大于该门限值，则判定误码性能重度恶化。利用本发明，完善和改进了现有BPSK卫星通信系统受单频干扰时误码性能恶化检测的方法，提供了一种准确的、具有一般性和实用性的检测方法。另外，本发明分析更为严谨，结果更加准确可信，并且更具实用性。
文档编号H04W84/06GK101814941SQ20101014439
公开日2010年8月25日 申请日期2010年4月12日 优先权日2010年4月12日
发明者郑磊, 苏彦, 李春来 申请人:中国科学院国家天文台