CN110212935A - 一种准确评估pcm/fm接收机鉴频器非相干解调性能的方法 - Google Patents

一种准确评估pcm/fm接收机鉴频器非相干解调性能的方法 Download PDF

Info

Publication number
CN110212935A
CN110212935A CN201910375497.8A CN201910375497A CN110212935A CN 110212935 A CN110212935 A CN 110212935A CN 201910375497 A CN201910375497 A CN 201910375497A CN 110212935 A CN110212935 A CN 110212935A
Authority
CN
China
Prior art keywords
noise
frequency discriminator
signal
pcm
bandwidth
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
Application number
CN201910375497.8A
Other languages
English (en)
Other versions
CN110212935B (zh
Inventor
丛波
段庆光
夏利利
王天云
杨建东
余清华
顾福飞
沈小青
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
63686 Troops Of Pla
Original Assignee
63686 Troops Of Pla
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 63686 Troops Of Pla filed Critical 63686 Troops Of Pla
Publication of CN110212935A publication Critical patent/CN110212935A/zh
Application granted granted Critical
Publication of CN110212935B publication Critical patent/CN110212935B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06Receivers
    • H04B1/16Circuits
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B14/00Transmission systems not characterised by the medium used for transmission
    • H04B14/02Transmission systems not characterised by the medium used for transmission characterised by the use of pulse modulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明涉及一种准确评估PCM/FM接收机鉴频器非相干解调性能的方法,基于统计信号处理技术,在深入分析PCM/FM接收机鉴频器非相干解调原理的基础上,并区分两种不同的情况,综合考虑信道滤波器带宽、噪声功率谱密度、解调滤波器带宽等因素,结合工程实际进行合理的近似和简化,提出新的误码率计算公式。然后,通过载波定标的方式,确定指定误码率下的信号电平,对设备的误码率进行统计和分析。本发明能够方便快速准确地评估遥测和通信设备性能,便于系统设计和系统评估工作。本发明提出的方法可以推广到航空遥测、移动通信、无线网络等领域,具有较强的普适性。

Description

一种准确评估PCM/FM接收机鉴频器非相干解调性能的方法
技术领域
本发明涉及一种准确评估PCM/FM接收机鉴频器非相干解调性能的新方法,属于信息与通信工程领域。
背景技术
PCM/FM是一种常用的靶场遥测调制体制,自上世纪70年代至今,已有近50年的历史。其调制体制简单而性能优越,地位至今无法替代。
在一般工程应用中,PCM/FM大都采用基于鉴频器的非相干解调方法。目前,在遥测设备的研制和测试过程中,PCM/FM的理论误码率一般按照FSK非相干解调的理论误码率进行计算,即:
然而,事实上,这个公式并不符合PCM/FM体制的情况。首先,这个公式只适用于FSK匹配滤波器方式的解调过程。而PCM/FM体制基本均采用基于鉴频器的解调方案。二者有一定区别。上述公式的根本问题是未能体现PCM/FM中调制频偏的影响。在FM体制中,调制频偏越大,解调输出的信噪比也就越高,而上述公式完全忽略了这一因素。
其次,采用鉴频器方案的PCM/FM解调的理论误码特性,在很大程度上依赖解调器的具体设计,如前置滤波器、中频带宽、检后低通滤波器等等。而上述公式未能体现这些影响。实际测试结果也表明:上述公式和实际的理论误码率存在着显著的差异。
本方法基于统计信号处理技术,综合考虑信道滤波器带宽、噪声功率谱密度、解调滤波器带宽等因素,结合工程实际进行合理的近似和简化,能够准确地评估通信系统的接收解调性能。
发明内容
本发明所要解决的技术问题是针对上述现有技术提供一种准确评估PCM/FM接收机鉴频器非相干解调性能的新方法,能够方面准确地评估遥测设备的误码率性能,从而便于遥测和通信系统的系统设计和系统评估。
本发明解决上述问题所采用的技术方案为:一种准确评估PCM/FM接收机鉴频器非相干解调性能的方法,该方法基于统计信号处理和随机过程分析,建立基于鉴频器的PCM/FM信号接收解调数学模型,综合考虑信道滤波器带宽、噪声功率谱密度、解调滤波器带宽因素,开发出定量的评估和定标方法,具体包括以下步骤:
步骤一、建立基于鉴频器的PCM/FM信号接收解调数学模型
记中频带通滤波器的带宽为BIF,鉴频器的归一化增益系数为1/2π,低通滤波器的带宽为BL
接收信号的时域表达式为:
xFM(t)=Ac·cos[2πfct+φs(t)]+n(t)
其中:
Ac为接收信号的幅度
fc为接收信号的频率;
φs(t)为数据调制造成的相位动态变化;
m(t)为基带的数据信号,不妨设其幅度为1和-1;
Kf为最大调制频偏;
中频带通滤波器输出的噪声表示为:
n(t)=nI(t)cos2πfct-nQ(t)sin2πfct
=r(t)cos[2πfct+ψ(t)]
其中:
r(t)为噪声的幅度;ψ(t)为噪声的相位,单边带噪声功率谱密度为N0/2,功率为Ni(t)=N0BIF,nI(t)、nQ(t)分别为噪声中的同相分量和正交分量;
步骤二、鉴频器输出噪声功率谱密度的求解
信号和噪声叠加后,总的相位表示为:
鉴频器的输出表示为:
xd(t)=(1/2π)dθ(t)/dt
=(1/2π)dφs(t)/dt+nd(t)
=Kfm(t)+nd(t)
其中,nd(t)是鉴频器输出的噪声,表示为:
当信噪比较大时,认为在[-π,π]间均匀分布,则上式进一步简化为:
由上式可以看出,鉴频器输出的噪声过程是输入噪声过程对时间的微分,则nd(t)的噪声功率谱密度可表示为:
Nd0=(1/Ac)2N0f2
步骤三、根据滤波器带宽选取误码率计算公式
噪声功率的计算分为两种情况:
(1)中频带宽小于或等于两倍低通滤波器带宽
鉴频器输出噪声功率通过对步骤二中获得的噪声功率谱密度积分取得:
则:
每比特码元的能量Eb由载波功率来计算,即
令Rb=1/Tb,h=2Kf/Rb,其中,h为调制指数,Tb为传输数据的码元宽度,Rb为传输数据的码元速率;
则误码率表示为:
其中
(2)中频带宽大于或等于两倍低通滤波器带宽
鉴频器输出噪声功率同样通过对步骤二中获得的噪声功率谱密度的积分来计算,其积分区间由低通滤波器带宽BL决定:
当BL大于或等于1/Tb,满足无码间干扰的要求,则BL=1/Tb,上式重写为:
Nd=(2/3)(1/Ac)2N0(1/Tb)3
则误码率可表示为:
令L=4/3h2,误码率公式可进一步表示为:
步骤四、测试信号定标和测试过程
根据系统测试指定的误码率,由步骤三中的公式计算出所需要的Eb/N0,再根据下列公式计算出载波的信号功率与噪声功率谱密度之比:
S/φ=Eb/N0+10logRb
此处,Rb为遥测数据码速率,
按上式计算得到的S/N0,加上系统的解调性能指标要求Ds,将相应的信号输入至通信系统或遥测系统中。
与现有技术相比,本发明的优点在于:
1)本发明能够方便快速准确地评估遥测和通信设备性能,便于系统设计和系统评估工作。
2)本发明提出的方法可以推广到航空遥测、移动通信、无线网络等领域,具有较强的普适性。
附图说明
图1为PCM/FM体制鉴频解调原理示意图。
图2为三种调制体制理论误码率比较示意图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
本实施例中的一种准确评估PCM/FM接收机鉴频器非相干解调性能的方法,该方法基于统计信号处理和随机过程分析,建立基于鉴频器的PCM/FM信号接收解调数学模型,综合考虑信道滤波器带宽、噪声功率谱密度、解调滤波器带宽等因素,开发出定量的评估和定标方法,提高了通信系统性能评估的准确性。
所述方法具体包括以下步骤:
步骤一、建立基于鉴频器的PCM/FM信号接收解调数学模型
典型的PCM/FM体制接收机的接收解调原理如图1所示。
信号和噪声叠加后,通过一个中频带通滤波器,滤除带外干扰。经过限幅器,克服信道传输过程中的幅度抖动。经过鉴频器,对信号的相位进行微分(实际是检测频率的变化),再通过低通滤波器,获得原始的调制信号。经判决后,得到原始1和0的信息。
记中频带通滤波器的带宽为BIF,鉴频器的归一化增益系数为1/2π,低通滤波器的带宽为BL
接收信号的时域表达式为:
xFM(t)=Ac·cos[2πfct+φs(t)]+n(t)
其中:
Ac为接收信号的幅度
fc为接收信号的频率;
φs(t)为数据调制造成的相位动态变化;
m(t)为基带的数据信号,不妨设其幅度为1和-1;
Kf为最大调制频偏。
中频带通滤波器可以视为以fc为中心的窄带滤波器,从而中频带通滤波器输出的噪声可表示为:
n(t)=nI(t)cos2πfct-nQ(t)sin2πfct
=r(t)cos[2πfct+ψ(t)]
其中:
r(t)为噪声的幅度;ψ(t)为噪声的相位,其单边带噪声功率谱密度为N0/2,功率为Ni(t)=N0BIF。nI(t)、nQ(t)分别为噪声中的同相分量和正交分量。
步骤二、鉴频器输出噪声功率谱密度的求解
信号和噪声叠加后,总的相位可表示为:
这个式子的表示比较复杂,难于运算和推导。鉴于在实际应用中,一般仅考虑信噪比较大的情况。在这种情况下,上式可近似简化为:
鉴频器的作用是对输入信号的相位进行微分,因此鉴频器的输出可表示为:
xd(t)=(1/2π)dθ(t)/dt
=(1/2π)dφs(t)/dt+nd(t)
=Kfm(t)+nd(t)
其中,nd(t)是鉴频器输出的噪声,
上式仍难于求解。为了进一步处理,需进行一定程度的近似。当信噪比较大时,可以认为在[-π,π]间均匀分布,则上式可进一步简化为:
由上式可以看出,鉴频器输出的噪声过程是输入噪声过程对时间的微分。微分运算相当于将输入噪声过程通过一个频率响应为j2πf的滤波器。如果输入是正态过程,输出仍为正态过程。则nd(t)的噪声功率谱密度可表示为:
Nd0=(1/Ac)2N0f2
此处值得注意的是,由于输入噪声经过正负频率轴上的两个中频滤波器(在物理上实为一个),此处输入的噪声功率谱密度应该按N0而不是N0/2计算。
由上式可以看出,鉴频器输出的噪声不再是白噪声,其功率谱密度呈抛物线型。这就造成了PCM/FM非相干检测的理论误码公式要比相干检测方式更为复杂。
步骤三、根据滤波器带宽选取误码率计算公式
鉴频器输出的信号可表示为:
当传输数据为1时,鉴频器输出的信号符合正态分布,其均值为Kf,方差(Nd为噪声功率),其表示式为:
由于受噪声的干扰,当x≤0时,判决器会错误地将1判决为0。错误概率为
其中,Q(·)函数的定义为:
当数据中1和0的概率相等时,解调过程中的误码率为:
Pe=Q(Kf/σ)
然而,噪声功率的计算却要区分两种情况。
噪声功率的计算区分两种情况。
(1)中频带宽小于或等于两倍低通滤波器带宽
中频滤波器等效到基带之后,对噪声的作用等效为一个区间为[-BIF/2,BIF/2]的低通滤波器。如果这个等效低通滤波器的带宽小于鉴频器后固有的硬件的低通滤波器(区间为[-BL,BL]),则鉴频器的输出噪声功率实际上取决于中频滤波器带宽BIF
鉴频器输出噪声功率可通过对步骤二中获得的噪声功率谱密度的积分来计算:
则:
由于鉴频器的输出只能反映输入信号相位的变化,无法反映功率的变化,因此,每比特码元的能量Eb不能由Kf计算,只能由载波功率来计算,即
令Rb=1/Tb,h=2Kf/Rb,其中,h一般被称为调制指数。
则误码率可进一步表示为:
则上式表示为:
考虑到,FSK相干检测方式的误码率为
此处,称L为PCM/FM非相干解调相对于FSK相干解调的恶化因子。为达到某个指定的误码率指标,记PCM/FM非相干解调所需的Eb/N0为Eb/N0|uncoherentCPFSK,记FSK相干解调所需的Eb/N0为Eb/N0|coherentFSK
由于函数较易计算和查表,则在其基础上,将某指定误码率所需的Eb/N0加上10log10L即可得到PCM/FM非相干检测所需的Eb/N0
它们的关系是:
Eb/N0|uncoherentCPFSK=Eb/N0|coherentFSK+10log10L
(2)中频带宽大于或等于两倍低通滤波器带宽
在这种情况下,鉴频器的输出噪声功率实际上取决于低通滤波器带宽BL
鉴频器输出噪声功率同样可通过对步骤二中获得的噪声功率谱密度的积分来计算,但其积分区间由低通滤波器带宽BL决定:
假设传输数据的码元宽度为Tb。一般来说,只要BL大于或等于1/Tb,即可满足无码间干扰的要求。此处假定BL=1/Tb,则上式可重写为:
Nd=(2/3)(1/Ac)2N0(1/Tb)3
则误码率可表示为:
令L=4/3h2,误码率公式可进一步表示为:
步骤四、测试信号定标和测试过程
通信系统或遥测系统的解调性能指标一般定义为:为达到某个指定的误码率,所需要的实际信号功率和理论信号功率的差值,一般用分贝(dB)表示。这个差值越小,说明系统的解调性能越优秀。此处用Ds表示这个指标。
根据系统测试指定的误码率,由步骤三中的公式计算出所需要的Eb/N0,再根据下列公式计算出载波的信号功率与噪声功率谱密度之比:
S/N0=Eb/N0+10logRb
此处,Rb为数据的码速率。
按上式计算得到的S/N0,加上系统的解调性能指标要求Ds,将相应的信号输入至通信系统或遥测系统中。如果系统的误码率低于指定的误码率,则表明该系统满足设计要求,如图2所示。
除上述实施例外,本发明还包括有其他实施方式,凡采用等同变换或者等效替换方式形成的技术方案,均应落入本发明权利要求的保护范围之内。

Claims (1)

1.一种准确评估PCM/FM接收机鉴频器非相干解调性能的方法,其特征在于:该方法基于统计信号处理和随机过程分析,建立基于鉴频器的PCM/FM信号接收解调数学模型,综合考虑信道滤波器带宽、噪声功率谱密度、解调滤波器带宽因素,开发出定量的评估和定标方法,具体包括以下步骤:
步骤一、建立基于鉴频器的PCM/FM信号接收解调数学模型
记中频带通滤波器的带宽为BIF,鉴频器的归一化增益系数为1/2π,低通滤波器的带宽为BL,
接收信号的时域表达式为:
xFM(t)=Ac·cos[2πfct+φs(t)]+n(t)
其中:
Ac为接收信号的幅度
fc为接收信号的频率;
φs(t)为数据调制造成的相位动态变化;
m(t)为基带的数据信号,不妨设其幅度为1和-1;
Kf为最大调制频偏;
中频带通滤波器输出的噪声表示为:
n(t)=nI(t)cos2πfct-nQ(t)sin2πfct
=r(t)cos[2πfct+ψ(t)]
其中:
r(t)为噪声的幅度;ψ(t)为噪声的相位,单边带噪声功率谱密度为N0/2,功率为Ni(t)=N0BIF,nI(t)、nQ(t)分别为噪声中的同相分量和正交分量;
步骤二、鉴频器输出噪声功率谱密度的求解
信号和噪声叠加后,总的相位表示为:
鉴频器的输出表示为:
xd(t)=(1/2π)dθ(t)/dt
=(1/2π)dφs(t)/dt+nd(t)
=Kfm(t)+nd(t)
其中,nd(t)是鉴频器输出的噪声,表示为:
当信噪比较大时,认为在[-π,π]间均匀分布,则上式进一步简化为:
由上式可以看出,鉴频器输出的噪声过程是输入噪声过程对时间的微分,则nd(t)的噪声功率谱密度可表示为:
Nd0=(1/Ac)2N0f2
步骤三、根据滤波器带宽选取误码率计算公式
噪声功率的计算分为两种情况:
(1)中频带宽小于或等于两倍低通滤波器带宽
鉴频器输出噪声功率通过对步骤二中获得的噪声功率谱密度积分取得:
则:
每比特码元的能量Eb由载波功率来计算,即
令Rb=1/Tb,h=2Kf/Rb,其中,h为调制指数,Tb为传输数据的码元宽度,Rb为传输数据的码元速率;
则误码率表示为:
其中
(2)中频带宽大于或等于两倍低通滤波器带宽
鉴频器输出噪声功率同样通过对步骤二中获得的噪声功率谱密度的积分来计算,其积分区间由低通滤波器带宽BL决定:
当BL大于或等于1/Tb,满足无码间干扰的要求,则BL=1/Tb,上式重写为:
Nd=(2/3)(1/Ac)2N0(1/Tb)3
则误码率可表示为:
令L=4/3h2,误码率公式可进一步表示为:
步骤四、测试信号定标和测试过程
根据系统测试指定的误码率,由步骤三中的公式计算出所需要的Eb/N0,再根据下列公式计算出载波的信号功率与噪声功率谱密度之比:
S/φ=Eb/N0+10logRb
此处,Rb为遥测数据码速率,
按上式计算得到的S/N0,加上系统的解调性能指标要求Ds,将相应的信号输入至通信系统或遥测系统中。
CN201910375497.8A 2018-11-06 2019-05-07 一种准确评估pcm/fm接收机鉴频器非相干解调性能的方法 Expired - Fee Related CN110212935B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2018113145168 2018-11-06
CN201811314516 2018-11-06

Publications (2)

Publication Number Publication Date
CN110212935A true CN110212935A (zh) 2019-09-06
CN110212935B CN110212935B (zh) 2021-05-11

Family

ID=67785582

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910375497.8A Expired - Fee Related CN110212935B (zh) 2018-11-06 2019-05-07 一种准确评估pcm/fm接收机鉴频器非相干解调性能的方法

Country Status (1)

Country Link
CN (1) CN110212935B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5711009A (en) * 1996-04-01 1998-01-20 The United States Of America As Represented By The Secretary Of The Navy Method and apparatus for minimizing the effects of frequency selective fading on a PCM/FM data signal
US20070178866A1 (en) * 2006-01-31 2007-08-02 Ibm Corporation Receiver and integrated AM-FM/IQ demodulators for gigabit-rate data detection
CN101640654A (zh) * 2009-07-27 2010-02-03 北京航空航天大学 一种用于卫星通信系统的超低码速率psk解调器
CN101741786A (zh) * 2009-12-18 2010-06-16 中国人民解放军理工大学 数字通信系统超宽带接收机及其信号处理方法
CN103501203A (zh) * 2013-09-26 2014-01-08 北京空间飞行器总体设计部 一种适用于遥感卫星激光星地通信链路测试系统
CN104486288A (zh) * 2014-12-01 2015-04-01 北京理工大学 一种适用于pcm/fm遥测接收机的载波频偏抑制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5711009A (en) * 1996-04-01 1998-01-20 The United States Of America As Represented By The Secretary Of The Navy Method and apparatus for minimizing the effects of frequency selective fading on a PCM/FM data signal
US20070178866A1 (en) * 2006-01-31 2007-08-02 Ibm Corporation Receiver and integrated AM-FM/IQ demodulators for gigabit-rate data detection
CN101640654A (zh) * 2009-07-27 2010-02-03 北京航空航天大学 一种用于卫星通信系统的超低码速率psk解调器
CN101741786A (zh) * 2009-12-18 2010-06-16 中国人民解放军理工大学 数字通信系统超宽带接收机及其信号处理方法
CN103501203A (zh) * 2013-09-26 2014-01-08 北京空间飞行器总体设计部 一种适用于遥感卫星激光星地通信链路测试系统
CN104486288A (zh) * 2014-12-01 2015-04-01 北京理工大学 一种适用于pcm/fm遥测接收机的载波频偏抑制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
QUANG-CHINH THAM,ETC: ""Optimal and Suboptimal Perfomance of a PCM/FM Communication System"", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 *
杨俊: ""遥测中调制解调及信噪比估计技术研究"", 《中国优秀硕士学位论文全文数据库(信息科技辑)》 *

Also Published As

Publication number Publication date
CN110212935B (zh) 2021-05-11

Similar Documents

Publication Publication Date Title
US4381546A (en) System for the quantitative measurement of impairments in the communication channel of a quadrature amplitude modulation data communication system
US6212246B1 (en) Symbol-quality evaluation in a digital communications receiver
AU643727B2 (en) Error rate monitor
CN104301056B (zh) 一种基于信号特征分析的频谱监测方法
CN111800359B (zh) 一种通信信号调制方式的识别方法、装置、设备及介质
US20060072679A1 (en) Method and device for modulation recognition of digitally modulated signals with multi-level magnitudes
EP0682433B1 (en) Method and apparatus for symbol clock phase recovery
CN112737992A (zh) 一种水声信号调制方式自适应类内识别方法
Wong et al. Emitter identification using CNN IQ imbalance estimators
Kadushkin et al. Determination of Occurrence Conditions for Periodically Non-Stationary Random Process with Establishing PSK-n Phase at the Output of the Linear Selective System
CN110212935A (zh) 一种准确评估pcm/fm接收机鉴频器非相干解调性能的方法
Jajoo et al. Blind signal modulation recognition through density spread of constellation signature
KR20190081747A (ko) 직교성을 갖는 동기 워드를 이용한 독립적인 패킷 검출 방법 및 그 수신 장치
Chuangzhan et al. Modulation recognition method of communication signals based on correlation characteristics
JP4969518B2 (ja) 変調方式識別回路および受信装置
US20020094783A1 (en) Space shift keying modulation
JP4361514B2 (ja) 変調型式識別方法、変調型式識別回路および復調装置
Song et al. Quantifying the accuracy of high SNR BER approximation of MPSK and MDPSK in fading channels
JP2003244263A (ja) 信号処理装置
KR102420971B1 (ko) 전력 스펙트럼 밀도 형상을 사용하는 신호 프로세싱을 위한 시스템들 및 방법들
Kuba et al. Automatic communication standard recognition in wireless smart home networks
Dvornikov et al. Recovery of signals of frequency manipulation according to the results of the analysis of edge effects
CN102946368A (zh) 多径衰落信道下含有频偏和相偏的数字调制信号识别方法
CN112202698B (zh) 心跳控制包非相干解调与检测的方法、装置及电子设备
Listova et al. The communications channels models in wireless sensor networks, based on the structural-energetic interaction between signals and interferences

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210511