CN1238622A - 用于数字解调调频信号的方法和设备 - Google Patents
用于数字解调调频信号的方法和设备 Download PDFInfo
- Publication number
- CN1238622A CN1238622A CN99102580A CN99102580A CN1238622A CN 1238622 A CN1238622 A CN 1238622A CN 99102580 A CN99102580 A CN 99102580A CN 99102580 A CN99102580 A CN 99102580A CN 1238622 A CN1238622 A CN 1238622A
- Authority
- CN
- China
- Prior art keywords
- phase place
- sample value
- expansion phase
- expansion
- grad
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/18—Phase-modulated carrier systems, i.e. using phase-shift keying
- H04L27/22—Demodulator circuits; Receiver circuits
- H04L27/233—Demodulator circuits; Receiver circuits using non-coherent demodulation
- H04L27/2332—Demodulator circuits; Receiver circuits using non-coherent demodulation using a non-coherent carrier
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Stereo-Broadcasting Methods (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Superheterodyne Receivers (AREA)
- Noise Elimination (AREA)
Abstract
本发明代表在数字解调例如用于fm音频广播的调频信号的领域中新的解决方案,本发明算法改善载波噪声比在fm门限以下的解调信号的信噪比。根据本发明,检测和消除对于载波噪声比在fm门限以下的信号一般出现的解调信号中的尖峰,利用调频信号的展开相位Φu的内插与其他准则进行检测,将展开相位Φu计算为不限于范围[-π,…,π]的绝对相位。
Description
本发明涉及用于数字解调例如立体声广播信号的调频信号的方法和设备,特别涉及数字调频(fm)接收机的fm门限的扩展(extension)。
对调制信号s(t)进行频率调制:
其中WC表示载波频率,ΔW表示频率偏移。
图10表示利用CORDIC算法的本领域数字fm解调器。将所接收的fm调制信号afm(t)利用模拟前端模块1下变换为IF(中频)信号b(t),将这个IF信号b(t)由A/D变换器2利用抽样周期T进行抽样,
并且由IQ生成装置3变换为复合基带信号
存在几种用于数字fm解调的算法。图11所示的fm解调器使用CORDIC算法。fm解调在原理上通过微分复合基带信号CT(KT)的角φC(KT)来完成:
利用近似式:
得到:
从中能看出:利用CORDIC模块4和微分模块7执行的fm解调通过CORDIC模块4中的复合基带信号CT(KT)的两个样值的角判定和微分模块7中的减法来完成。
CORDIC算法利用从笛卡坐标到极坐标的变换计算复合基带信号CT(KT)的模2π除的相位φ(KT)。
这个CORDIC算法是众所周知的。图11和下表1给出此算法的简短介绍。图11表示具有复合基带信号CT(KT)的5个样值点的复平面。CORDIC算法的幅角范围限制为[-π,…,π]。因此,利用CORDIC算法计算的相位φ(KT)进行模2π除。下表1表示复合基带信号CT(KT)的绝对相位φc(KT)和利用CORDIC算法计算的模2π除的相位φ(KT)。
K | 1 | 2 | 3 | 4 | 5 |
φ(KT) | 2π/5 | 4π/5 | -4π/5 | -π/5 | π/5 |
φc(KT) | 2π/5 | 4π/5 | 6π/5 | 9π/5 | 11π/5 |
表1
使用CORDIC算法的本领域解调器利用式(4)中所示的减法来解调fm调制的基带信号CT(KT)。由于CORDIC算法如上所述计算模2π除的相位φ(KT),所以随即进行减法并得到下面的调制信号S(KT):
下表2表示在这个操作中计算的模2π除的相位φ(KT)和绝对相位φc(KT)的相位差:
K | 2 | 3 | 4 | 5 |
φ(KT)-φ((K-1)T) | 2π/5 | -8π/5=-2π+2π/5 | 3π/5 | 2π/5 |
φc(KT)-φc((K-1)T) | 2π/5 | 2π/5 | 3π/5 | 2π/5 |
表2
能看出:在利用CORDIC算法计算的模2π除的相位φ(KT)的相位差的抽样值K=3时,出现-2π的误差。在本领域解调器中,这个误差利用通过固定点实施中的溢出完成的模2π操作来消除。因此,相位差φ(KT)-φ((K-1)T)的有效范围在固定点实施中限制为[-π,…,π],于是值-8π/5=-2π+2π/5引起溢出,并因此错误值-8π/5进行模2π除而变为正确值2π/5。
本发明的目的是改善载波噪声比在fm门限以下的解调信号的信噪比。
根据本发明,这个目的通过展开在从笛卡儿坐标至极坐标的变换期间计算的模2π除的相位φ(KT)并检测和消除展开相位中的相位跃变以生成随后在微分级中使用的校正相位来实现。
根据本发明的方法定义在独立权利要求1中,从属权利要求2至19表示本发明的方法的最佳实施例,能完成本发明方法的设备定义在权利要求20和21中。
从下面对现有技术的解调器的分析和结合附图仅作为示例的本发明的最佳实施例的详细描述,将会更好地理解本发明的其他目的、特性和优点,其中:
图1表示根据本发明一个实施例的fm解调器的原理方框图;
图2表示绝对相位φc(KT),即具有无限载波噪声比的理想相位,和具有载波噪声比10dB的频率调制信号的在展开状态中通过从笛卡儿坐标变换为极坐标计算的相位φu(KT),这里,调制信号是如同在fm广播中使用的具有音乐信息的立体声多路复用信号;
图3表示放大比例的具有载波噪声比10dB的频率调制信号的另一个展开的相位信号φu(KT),这里,调制信号还是具有音乐信息的立体声多路复用信号;
图4表示展开的相位φu(KT)生成的流程图;
图5表示图2与3中所示的展开相位φu(KT)失真的检测与消除的原理概念流程图;
图6表示具有载波噪声比10dB的展开相位φu(KT)和根据本发明校正的相位φe(KT);
图7表示根据本发明第二实施例的展开相位、失真的检测与消除以及微分的方框流程图;
图8表示根据本发明第二实施例的展开相位、失真的检测与消除以微分的直接实现模块;
图9表示现有技术fm解调器和根据本发明的fm解调器的解调器输出上的信噪比;
图10表示使用CORDIC算法的现有技术fm解调器的方框图;和
图11表示具有复合基带信号CT(KT)的5个样值点的复平面。
图1表示根据本发明一个实施例的用于数字解调频率调制信号afm(t)的设备的原理框图。具有与图10中所示的并结合图10中所示的现有技术设备描述的功能相同的功能的方框利用相同的标号来表示。现有技术设备的微分模块7产生本发明不允许的溢出,所以根据本发明的设备包括不产生溢出的微分模块6。除了现有技术设备之外,根据本发明的解调器还包括数字门限扩展单元5,此单元从CORDIC模块4中接收相位φ(KT),CORDIC模块4计算随后输入给微分模块6的校正相位φe(KT)。附加的数字IF滤波器能放置在CORDIC模块4的前面,IF滤波器对于相邻信道抑制是必需的。
根据本发明的实施例,数字门限扩展单元5包括从CORDIC模块4的相位输出φ(KT)中计算展开相位φu(KT)的相位展开单元5a和校正单元5b,校正单元5b检测和消除展开相位φu(KT)中的2π相位跃变并输出随后输入微分模块6的校正相位φe(KT)。展开相位φu(KT)计算为不限于范围[-π,…,π]的绝对相位。
利用展开相位φu(KT)而不利用校正单元5b进行的fm解调得到与上述的现有技术fm解调器完全相同的值,这个情况与fm信号的载波噪声比(CNR)无关。
图2表示具有利用无限载波噪声比计算的递减斜率的调频立体声多路复用信号的绝对相位φe(KT),这个绝对相位利用虚线表示。图2也表示具有载波噪声比10dB的相应的展开相位φu(KT),这个展开相位φu(KT)利用实线表示,并且能清楚地看出:相对于绝对相位,从值K=0开始,此相位φu(KT)包括2π的相位跃变。
在样值K=0上的展开相位信号φu(KT)的2π相位跃变对于在fm门限以下的载波噪声比是典型的。在微分之后,这个相位跃变在解调信号中产生噪声峰值,这样的噪声峰值也在执行现有技术fm解调器的模2π除的相位φ(KT)的微分时出现。频率解调信号S(KT)中的这些峰值是动态降低信噪比在fm门限以下的原因。
展开相位φu(KT)的这些+2π相位跃变仅在φu(KT)的递减斜率上出现在fm接收机操作的正常模式中,具有值-2π的相位跃变仅出现在上升斜率上。
根据本发明,检测和消除这样的相位跃变,以获得具有解调信号中显著减少的噪声峰值数的fm解调器,并因此获得比现有技术fm解调器更好的fm门限之下的信噪比(SNR)。
图3表示具有载波噪声比10dB的调频立体声多路复用信号的展开相位φu(KT)的放大视图,这次具有上升斜率,这个展开相位信号φu(KT)在样值K=0时具有-2π的相位跃变。
图4表示根据本发明如何从模2π除的相位φ(KT)中计算展开的相位φu(KT)的原理流程图。在步骤S1中开始之后,在步骤S2检查样值K上的相位φ(KT)与在样值K的前一样值上的相位φ((K-1)T)的差是否超过π。如果不是这种情况,则在步骤S3中检查此差是否低于-π。如果也不是这种情况,即如果所述差在范围[-π,…,π]内,则将展开相位φu(KT)设置为相位φ(KT)。随后,此算法利用步骤S2对下一样值K=K+1继续进行。
另一方面,如何在步骤S2中确定样值K上的相位φ(KT)与样值K的前一样值的相位φ((K-1)T)之间的所述差超过π,则在利用步骤S6继续此算法之前将在步骤S4中从实际样值的相位φ(KT)和所有后面样值的相位中减去2π。
如果在步骤实际样值的相位φ(KT)与前一样值的相位φ((K-1)T)之间的差落在-π以下,则在利用步骤S6继续此算法之前,将在步骤S5中以相同的方式将2π加到实际样值的相位φ(KT)和所有后面样值的相位上。
利用图1所示的校正单元5b检测和消除展开相位φu(KT)的相位跃变,此校正单元5b用作检测与消除模块,在可能的相位跃变之前与之后结合内插检查下降斜率上展开相位信号φu(KT)的正相位跃变和上升斜率上负的相位跃变。如果检测到相位跃变,此模块通过将2π加到随后样值的展开相位信号φu(KT)上或从随后样值的展开相位信号φu(KT)中减去2π来消除相位跃变。
检测相位跃变的第一准则是上升或下降斜率。在要检查相位跃变出现的测试样值之前和之后结合内插确定此准则,下面解释此内插。
此内插结合高抽样率的展开相位信号φu(KT)使用调频信号的频谱分布。在fm立体声广播中,立体声多路复用信号S(t)主要包含低频成分。由于基带信号CT(KT)的带宽,解调器的抽样率得选择得高,例如,1/T>300KHz。因此,相位跃变具有比绝对相位φc(KT)更强的梯度并且能检测到。
如上所述,图3表示在样值K=0时具有-2π的相位跃变的fm解调信号的展开相位φu(KT)。对于样值K的内插,得计算K之前的梯度,即gradbefφu(KT)和K之后的梯度,即下面称为gradaftφu(KT)。这应该通过平均多于1个的样值的梯度来完成,这是因为展开相位信号φu(KT)被噪声干扰。而且,直接相邻相位跃变的样值也不应用于内插,因为这些样值被相位跃变干扰,如能从图3中所看到的。因此,利用下式计算K之前和K之后的梯度。
其中p和q表示用于平均的样值数,并且o和r表示用于计算梯度的样值距K的最小距离。例如,如果从距K的距离o=r=4的样值开始通过平均p=q=8个样值计算K之前与之后的梯度,则利用下式计算梯度:
梯度(7)现在用于内插。利用下式利用K之前的梯度在样值K上计算展开相位φu(KT)的内插误差:
φerrbefd.e(kT)=φU((k-d-1)T)+(d+e+1)·gradbef(kT)-φU((k+e)T) (8a)
其中d和e表示内插样值距实际样值K的距离。在下面,d和e假定为相同的值(d=e),但这不是必须的。
由于直接相邻样值K上的相位跃变的样值被干扰并因此不应用于展开相位的内插,所以应选择d和e≥2。
为了提高内插误差的可靠性,应对于多于1个的样值进行内插。下式表示d=e=2、3、4、5的内插误差φerrbefd,,e(KT):
φerrbef2.2(kT)=φU((k-3)T)+5·graabef(kT)-φU((k+2)T) (8b)
φerrbef3.3(kT)=φU((k-4)T)+7·gradbef(kT)-φU((k+3)T)
φerrbef4.4(kT)=φU((k-5)T)+9·gradbef(kT)-φU((k+4)T)
φerrbef5.5(kT)=φU((k-6)T)+11·gradbef(kT)-φU((k+5)T)
作为批准内插可靠性的另一个点,也应利用相位跃变之后的梯度来完成。随后,将利用下面的内插计算相位误差φerrafef,,g(KT):
φerraftf.g(kT)=φU((k-f-1)T)+(f+g+1)·gradaft(kT)-φU((k+g)T) (9a)
其中f和g表示内插样值距K的距离,再次假设f和g是相同的值(f=g),即使这不是必须的。
由于在K上相位跃变相邻样值的不可靠性,内插样值距K的最小距离应再次进行选择,f,g≥2。下面的示例表示f=g=2、3、4、5时内插误差:φerrafef·g(KT):φerraft.2,2(kT)=φU((k-3)T)+5·gradaft(kT)-φU((k+2)T) (9b)φerraft.3.3(kT)=φU((k-4)T)+7·gradaft(kT)-φU((k+3)T)φerraft.4.4(kT)=φU((k-5)T)+9·gradaft(kT)-φU((k+4)T)φerraft.5.5(kT)=φU((k-6)T)+11·gradaft(kT)-φU((k+5)T)
所有内插值相加得到总的内插误差:
其中do和fo表示内插样值距K的最小距离,而dmax与fmax表示内插样值距K的最大距离。在上面的示例中,do=fo=2并且dmax=fmax=5,在利用(dmax-do+fmax-fo+2)=8个样值进行内插时,所有内插值相加得到总的内插误差:
φerr(kT)=2·(φU((k-3)T)+φU((k-4)T)+φU((k-5)T)+φU((k-6)T)) (10b)
-2·(φU((k+2)T)+φU((k+3)T)+φU((k+4)T)+φU((k+5)T))
+32·(gradbef(kT)+gradaft(kT))
=2·(φU((k-3)T)+φU((k-4)T)+φU((k-5)T)+φU((k-6)T))
-2·(φU((k+2)T)+φU((k+3)T)+φU((k+4)T)+φU((k+5)T))
+4·(φU((k-4)T)-φU((k-12)T)+φU((k+11)T)-φU((k+3)T))
在理想的相位跃变情况中,能利用下式计算总的内插误差:
φerrideal=2π·(dmax-do+fmax-fo+2) (11a)
在上面示例中,do=fo=2并且dmax=fmax=5,则总的内插误差为:
φerrideal=2π·8=50,27 (11b)
另一方面,实际经验已证明,φerr.thr=1/2φerrideal(12)的内插误差门限在约410KHz的抽样率1/T上得到好的性能。对于上面的示例,选择φerr.ehr=26。
如前所述,得检查展开相位φu(KT)的相位跃变。由于值为+2π的相位跃变仅出现在展开相位信号φu(KT)的下降斜率上,所以能以下面的方式检测相位跃变,能以展开相位的梯度完成下降斜率的检测。下面选择KTgradbef(KT)之前的梯度来检查此条件(gradbef(KT)<-gradthr),也可能选择KT(gradaft(KT)<-gradthr)之后的梯度或(gradbef(KT)+gradaft(KT)<-gradthr)的组合,其中gradthr是最好等于零或靠近零的小值的门限值。
(φU(kT)-φU((k-1)T>φthr1(kT)) (13)
AND(gradbef(kT)<-gradthr)
AND(φerr(kT)<-φerr.thr)
相位跃变不总是能从图2与3中所看见的尖跃变。一些相位跃变延伸到一个以上的样值上,以致根据本发明下面条件包含在相位跃变的检查中:
φU(kT)-φU((k-2)T)>φthr2(kT) (14)
利用式(13)与(14),随后是检测展开相位信号φu(KT)下降斜率上的相位跃变的条件:
[(φU(kT)-φU((k-1)T)>φthr1(kT))OR(φU(kT)-φU((k-2)T)>φthr2(kT))] (15)
AND(gradbef(kT)<-gradthr)
AND(φerr(kT)<-φerr.thr)
能以类似方式描述上升斜率φu(KT)的条件:[(φU(kT)-φU((k-1)T)<-φthr1(kT))OR(φU(kT)-φU((k-2)T)<-φthr2(kT))] (16)AND(gradbef(kT)>gradthr)AND(φerr(kT)>φerr.thr
上升斜率gradbef(KT)>gradthr的检查能再次利用gradaft(KT)>gradthr或(gradbef(KT)+gradaft(KT))>gradthr来代替。
实践经验已证明:上面公式(15)与(16)应如下以约410KHZ的抽样率1/T进行选择:
φthr1=1,5
φthr1=2 (17)
gradthr=0
检测与消除的流程图表示在图5中,这个流程图包括上面讨论的结合相位跃变之前与之后的内插φu(KT)下降斜率上正相位跃变和φu(KT)上升斜率上负相位跃变的检查。如果检测到相位跃变,此模块通过给展开相位信号φu(KT)的实际样值和随后样值加上或减去2π来消除相位跃变。
在步骤S10中开始之后,将在步骤S11中检查条件(16)。如果此条件不满足,将在步骤S12中检查条件(15)。如果此条件也不满足,即如果未检测到相位跃变,将在步骤S15中将校正相位φe(KT)设置为展开相位φu(KT)的值。此后,利用下一样值K=K+1从步骤S11开始重复此算法。
另一方面,如果在步骤S11中确定满足条件(16),则将2π加到展开相位φu(KT)的实际样值上,并对步骤S15之前所有后续样值执行相同操作。同样,如果在步骤S12中确定满足条件(15),则将从展开相位φu(KT)的实际样值中减去2π,并对步骤S15之前所有后续样值执行同样操作。
图6表示校正单元5b之后的均具有载波噪声比10dB的展开相位φu(KT)和已消除相位跃变的相位φe(KT),从中能清楚地看出:已消除相位跃变的相位φe(KT)更好地对应图2所示的具有无限载波噪声比的绝对相位φc(KT),则展开相位φu(KT)没有2π相位跃变的检测与消息。
已消除相位跃变φe(KT)的相位现在能进行微分而没有fm门限引起的噪声。这表示:能得到在fm门限以下的解调信号的信噪比的高增益。
在根据本发明的第二实施例中,考虑不是很容易加上或从相位或展开相位的实际样值与所有随后样值中减去2π,即不可能以容易的方式例如来解答图4的步骤S4中所示的方程式:对于i=k至∞ (18)
φ(iT)=φ(iT)-2π
结束
因此,根据本发明第二实施例,计算展开相位φu(KT)的窗口限制为此算法所必需的值。在上面的示例中,窗口限制为范围[K-O-P…K+q+r-1],例如[k-12…K+11],展开相位φu(KT)的计算利用辅助变量W来完成,此变量W包含有关给模2π除的相位φ(KT)加上或减去2π的信息:
利用这个辅助变量W(KT)能以下面方法从样值(K-N)T开始在窗口中计算展开相位φu(KT):
φu((k-N)T)=φ((k-N)T)(20a)
φu((k-N+1)T)=φ((k-N+1)T)+2πw((k-N+1)T)…
其中N表示用于展开相位φ(KT)的窗口的开始,V表示样值(K-N+V)T距用于展开相位的窗口开始(K-N)T的距离。
利用N=o+p=12,则上面示例为:φu((k-12)T)=φ((k-12)T) (20b)φu((k-11)T)=φ((k-11)T)+2πw((k-11)T)φu((k-10)T)=φ((k-10)T)+2π(w((k-11)T)+w((k-10)T))…
利用这个等式,等式(7b)的示例为:
将式10(b)变换为另一表示式得到:
-(φu((k+2)T)+φu((k+3)T)+φu((k+4)T)+φu((k+5)T))
+16·gradaft(kT)+16·gradbef(kT)
=φu((k-3)T)+φu((k-4)T)+φu((k-5)T)+φu((k-6)T)
-(φu((k+2)T)+φu((k+3)T)+φu((k+4)T)+φu((k+5)T))
+2·(φu((k+11)T)-φu((k+3)T))+16·gradbef(kT)
利用N=6的用于展开相位计算的另一窗口定义:
φu((k-6)T)=φ((k-6)T) (23)
φu((k-5)T)=φ((k-5)T)+2πw(k-5)
…
得到:
-φ((k+2)T)-3φ((k+3)T)-φ((k+4)T)-φ((k+5)T)
+2π·(-w((k-5)T)-2w((k-4)T)-3w((k-3)T)-4w((k-2)T)
-4w((k-1)T)-4w(kT)-4w((k+1)T)-4w((k+2)T)
-3w((k+3)T)+w((k+5)T)+2w((k+6)T)+2w((k+7)T)
+2w((k+8)T)+2w((k+9)T)+2w((k+10)T)+2w((k+11)T)}
+16·gradbef(kT)
式(15)、(16)与图5中所示的条件等式能利用下式替代:
φu(kT)-φu((k-1)T)=φ(kT)-φ((k-1)T)+2πw(kT) (25)
φu(kT)-φu((k-2)T)=φ(kT)-φ((k-2)T)+2π(w(kT)+w((k-1)T))
图5所示的上面示例的流程图利用式(21)、(24)与(25)在利用下面条件替代在步骤S11中确定的条件时能简化为如下:
AND gradbef>gradthr AND φerr>φerr,thr
在步骤S12中确定的条件利用下面的条件替代:
AND gradbef<-gradthr ANDφerr<-φerr,thr
并且用于生成频率解调信号S(KT)的微分在步骤S15中包括在下式中:
s(kT)=φ(kT)-φ((k-1)T)+2πw(kT) (28)
k=k+1
图7表示对应根据本发明的第二实施例中执行的程序的流程图。在步骤S20中程序开始之后,在步骤S21中确定条件(26)。如果这个条件不满足,则在步骤S22中检查条件(27)。如果这个条件也不满足,则在从步骤S21开始利用下一样值K=K+1重复程序之前,根据等式(28)在步骤S25中执行用于生成频率解调信号S(KT)的微分。
另一方面,如果在步骤S21中确定条件(26)不满足,则在步骤S23中将变量W(KT)设置为W(KT)=W(KT)+1。对应于此,当在执行步骤S25之前在步骤S22中确定条件(27)满足时,在步骤S24中将W(KT)设置为W(KT)=W(KT)-1。
图8表示根据本发明第二实施例的数字门限扩展单元5与微分单元6的直接实现,这个方框图利用级为T的延迟单元、乘法器、加法器、减法器和两个比较单元来实现。此方框图可直接从上面等式与条件(19)至(25)中导出。
图9表示根据本发明的具有(实线)和没有(虚线)门限扩展的频率解调音乐信号的信噪比。利用用于fm音频广播的参数进行模拟,以使频偏ΔW=2π-75KHz并且去加重滤波器包括在解调器输出上,能发现fm门限近似为CNR=17dB。由门限扩展引起的信噪比的增益约为9dB,SNR的增益取决于调制信号的频偏,频偏ΔW越高,SNR的增益越高。
Claims (21)
1.一种用于数字解调频率调制信号(afm(t))的方法,包括以下步骤:
a)计算接收的频率调制信号(afm(t))的数字复合基带信号CT(KT)
的相位(φ);
b)计算相位(φ)的展开相位(φu);
c)检测和消除展开相位(φu)中的2π相位跃变以生成校正的相位
(φe);和
d)微分校正的相位(φe)。
2.根据权利要求1的方法,其特征在于,所述步骤b)对于每一个输入样值的相位(φ)包括下面步骤:
b1)如果实际样值的相位(φu(KT))与前一样值相位(φ((K-1)T))之间的差超过π,则从实际样值相位(φ(KT))和预定数量(I)的随后样值的相位(φ((K+1)T)…(φ((K+I)T))中减去2π;和
b2)如果实际样值相位(φ(KT))与前一样值相位(φ((K-1)T))之间的差低于-π,则将2π加到实际样值相位(φ(KT))与预定数量(I)的随后样值相位(φ((K+1)T))…(φ((K+I)T))上。
3.根据权利要求1或2的方法,其特征在于,所述步骤c)对于每个样值的展开相位(φu)包括下列步骤:
c1)如果实际样值的展开相位(φu(KT))与前一样值的展开相位(φu((K-1)T))之间的差低于负的第一门限值(-φthri)并且实际样值的展开相位(φu(KT))之前的展开相位的梯度(gradbef(KT))或实际样值的展开相位(φu(KT))之后的展开相位的梯度(gradaft(KT))或二者的组合高于梯度门限(gradthr),则将2π加到实际样值的展开相位(φu(KT))与预定数量(I)的随后样值的展开相位(φu((K+1)T)…φu((K+I)T))上;和
c2)如果实际样值的展开相位(φu(KT))与前一样值的展开相位(φu((K-1)T))之间的差超过正的第一门限值(φthr1)并且实际样值的展开相位(φu(KT))之前的展开相位的梯度(gradbef(KT))或实际样值的展开相位(φu(KT))之后的展开相位的梯度(gradaft(KT))或二者的组合低于负的梯度门限(-gradthr),则从实际相值的展开相位(φu(KT))与预定数量(I)的随后样值的展开相位(φu((K+1)T)…φu((K+I)T))中减去2π。
4.根据权利要求1-3中任何一个权利要求的方法,其特征在于,所述步骤c)对于每个样值的展开相位(φu)包括以下步骤:
cⅠ)如果实际样值的展开相位(φu(KT))与前一样值的展开相位(φu((K-1)T))之间的差不低于负的第一门限值(-φthr1),但实际样值的展开相位(φu(KT))与前一样值之前的样值的展开相位(φu((K-2)T))之间的差低于负的第二门限值(-φthr2),并且实际样值的展开相位(φu(KT))之前的展开相位梯度(gradbef(KT))或实际样值的展开相位(φu(KT))之后的展开相位梯度(gradaft(KT))或二者的组合高于梯度门限(gradthr(KT)),则将2π加到实际样值的展开相位(φu(KT))与预定数量(I)的随后样值的展开相位(φu((K+1)T)…φ((K+I)T))上;和
cⅡ)如果实际样值的展开相位(φu(KT))与前一样值的展开相位(φu((K-1)T))之间的差不超过正的第一门限值(φthr1),但实际样值的展开相位(φu(KT))与前一样值之前的样值的展开相位(φu((K-2)T))之间的差超过正的第二门限值(φthr2),并且实际样值的展开相位(φu(KT))之前的展开相位梯度(gradbef(KT))或实际样值的展开相位(φu(KT))之后的展开相位梯度(gradaft(KT))或二者的组合低于负的梯度门限(-gradthr),则从实际样值的展开相位(φu(KT))与预定数量(I)的随后样值的展开相位(φu((K+1)T)…φu((K+I)T))中减去2π。
5.根据权利要求2-4中任何一个权利要求的方法,其特征在于,所述预定数量(I)的随后样值是所有的随后样值。
6.根据权利要求1的方法,其特征在于,所述步骤b)对于每一个输入样值的相位(φ)包括以下步骤:
b1)如果实际样值的相位(φu(KT))与前一样值的相位(φ(K-1)T))之间的差超过π,设置辅助变量(W(KT))为-1;
b2)如果实际样值的相位(φ(KT))与前一样值的相位(φ((K-1)T))之间的差低于-π,则将辅助变量(W(KT)设置为1;和
b3)如果实际相位的相位(φ(KT))与前一样值的相位(φ((K-1)T))之间的差位于-π与π之间,则设置辅助变量(W(KT))为0。
7.根据权利要求6的方法,其特征在于,所述步骤c)对于每个样值的展开相位(φu)包括以下步骤:
c1)如果实际样值的展开相位(φu(KT))与前一样值的展开相位(φu((K-1)T))之间的差低于负的第一门限值(-φthr1),并且实际样值的展开相位(φu(KT))之前的展开相位梯度(gradbef(KT))或实际样值的展开相位(φu(KT))之后的展开相位梯度(gradaft(KT))或二者的组合高于梯度门限(gradthr),则将所述辅助变量(W(KT))加1;和
c2)如果实际样值的展开相位(φu(KT))与前一样值的展开相位(φu((K-1)T))之间的差超过正的第一门限值(φthr1),并且实际样值的展开相位(φu(KT))之前的展开相位梯度(gradbef(KT))或实际样值的展开相位(φu(KT))之后的展开相位梯度(gradaft(KT))或二者的组合低于负的梯度门限(-gradthr),则从所述辅助变量(W(KT))中减去1。
8.根据权利要求6或7的方法,其特征在于,所述步骤c)对于每个样值的展开相位(φu)包括以下步骤:
(Ⅰ)如果实际样值的展开相位(φu(KT))与前一样值的展开相位
(φu((K-1)T))之间的差不低于负的第一门限值(-φthr1),但
实际样值的展开相位(φu(KT))与前一样值之前的样值的展开
相位(φu((K-2)T))之间的差低于负的第二门限值(-φthr2),
并且实际样值的展开相位(φu(KT))之前的展开相位梯度
(gradbef(KT))或实际样值的展开相位(φu(KT))之后的展开相
位梯度(gradaft(KT))或二者的组合高于梯度门限(gradthr),
则将所述辅助变量(W(KT))加1;和
(Ⅱ)如果实际样值的展开相位(φu(KT))与前一样值的展开相位
(φu((K-1)T))之间的差不超过正的第一门限值(φthr1),但实
际样值的展开相位(φu(KT))与前一样值之前的样值的展开相
位(φu((K-2)T))之间的差超过正的第二门限值(φthr2),并
且实际样值的展开相位(φu(KT))之前的展开相位的梯度
(gradbef(KT))或实际样值的展开相位(φu(KT))之后的展开相
位梯度(gradatf(KT))或二者的组合低于负的梯度门限(-
gradthr),从所述辅助变量(W(KT))中减1。
9.根据权利要求6-8中任何一个权利要求的方法,其特征在于,从样值(K-N)T开始的下面方法在窗口中计算展开相位(φu(KT)):
φu((k-N)T)=φ((k-N)T)
φu((k-N+1)T)=φ((k-N+1)T)+2πw((k-N+1)T)…
其中N表示用于展开相位(φu(KT))的窗口的开始,而V表示样值(K-N+V)T距用于展开相位的窗口的开始(K-N)T的距离。
10.根据权利要求3-5或7-9中任何权利要求的方法,其特征在于,如果实际样值的展开相位(φu(KT))上的展开相位的相位误差(φerr(KT))超过正的相位误差门限值(φerr,thr),则只执行所述步骤C1)或/和CI),并且如果实际样值的展开相位(φu(KT))上的展开相位的相位误差(φerr(KT))低于负的相位误差门限值(-φerr,thr),则只执行所述步骤c2)或/和CII)。
11.根据权利要求3-5或7-10中任何一个权利要求的方法,其特征在于,所述第一门限值(φthr1)、所述第二门限值(φhtr2)和所述梯度门限(gradthr)根据抽样速率(1/T)利用实验来确定。
12.根据权利要求11的方法,其特征在于,所述第一门限值(φthres1)等于1.5,所述第二门限值(φthres2)等于2,并且所述梯度门限(gradthr)等于0,当所述抽样速率(1/T)在410KHz范围中时。
13.根据权利要求3-5或7-12中任何一个权利要求的方法,其特征在于,实际样值的展开相位(φu(KT))上的展开相位的所述相位误差(φerr(KT))根据下式进行计算:
其中gradbef(KT)对应实际样值的展开相位(φu(KT))之前的展开相位梯度,而gradaft(KT)对应实际样值的展开相位(φu(KT))之后的展开相位的梯度,d0与f0表示内插样值距K的最小距离,dmax与fmax表示内插样值距K的最大距离,并且d、e、f与g表示用于内插的样值距K的距离。
14.根据权利要求13的方法,其特征在于,do=fo=2,dmax=fmax=5,并且e=d,g=f。
15.根据权利要求3-5或7-14中任何一个权利要求的方法,其特征在于,实际样值的展开相位(φu(KT))之前的展开相位的所述梯度(gradbef(KT))通过根据下式平均不与实际样值展开相位(φu(KT))直接相邻的P个样值的展开相位(φu)来计算,距实际样值K的最小距离为0,并且实际样值的展开相位(φu(KT))之后的展开相位的所述梯度(gradaft(KT))通过根据下式平均不与实际样值的展开相位(φu(KT))直接相邻的q个样值的展开相位(φu)来计算,距实际样值K的最小距离为r:
16.根据权利要求15的方法,其特征在于,o=r=4,并且p=q=8。
17.根据权利要求10-16中任何一个权利要求的方法,其特征在于,根据抽样速率(1/T)利用实验确定所述相位误差门限值(φerr,thr),并且此实验好的开始值是:
_φerr,thr≈1/2 φerrideal=1/2·2π(dmax-do+fmax-fo+2)。
18.根据权利要求17的方法,其特征在于,当所述抽样速率(1/T)在410KHz范围中时,相位误差门限值(φerr,thr)等于26。
19.根据权利要求1-18中任何一个权利要求的方法,其特征在于,利用CORDIC算法执行数字复合基带信号(CT(KT))相位的所述计算。
20.一种用于数字解调频率调制信号(afm(t))的设备,包括:
a)变换单元(4),计算接收的调频信号(afm(t))的数字复合基
带信号CT(KT)的相位(φ);
b)相位展开单元(5a),计算变换单元(4)的相位(4)输出的
展开相位(φu);
c)校正单元(5b),检测与消除展开相位(φu)中的2π相位跃变并
输出校正的相位(φe);和
d)微分单元(6),微分校正的相位(φe)。
21.根据权利要求20的设备,其特征在于,它根据上面的权利要求1-19中任何一个权利要求中定义的方法工作。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP98103735A EP0940958A1 (en) | 1998-03-03 | 1998-03-03 | Method and device for digitally demodulating a frequency modulated signal |
EP98103735.1 | 1998-03-03 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1238622A true CN1238622A (zh) | 1999-12-15 |
CN1115833C CN1115833C (zh) | 2003-07-23 |
Family
ID=8231518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN99102580A Expired - Fee Related CN1115833C (zh) | 1998-03-03 | 1999-03-03 | 用于数字解调调频信号的方法和设备 |
Country Status (5)
Country | Link |
---|---|
US (1) | US6075410A (zh) |
EP (1) | EP0940958A1 (zh) |
JP (1) | JP2000059253A (zh) |
CN (1) | CN1115833C (zh) |
CA (1) | CA2261442C (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114785434A (zh) * | 2022-03-23 | 2022-07-22 | 北京睿信丰科技有限公司 | S波段频道可用性检测方法、装置、设备、系统及介质 |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19860402A1 (de) * | 1998-12-28 | 2000-06-29 | Bosch Gmbh Robert | Verfahren und Schaltungsanordnung zur Demodulation eines digitalen frequenzmodulierten Signals |
US7388934B2 (en) | 2000-06-27 | 2008-06-17 | Koninklijke Philips Electronics N.V. | Detection and correction of phase jumps in a phase sequence |
US6969984B2 (en) * | 2000-12-21 | 2005-11-29 | Tropian, Inc. | Direct phase and frequency demodulation |
EP1303060A1 (en) * | 2001-10-11 | 2003-04-16 | Sony International (Europe) GmbH | Receive beam forming with shared phase estimation |
US7356091B2 (en) | 2003-12-09 | 2008-04-08 | M/A-Com, Inc. | Apparatus, methods and articles of manufacture for signal propagation using unwrapped phase |
CN101253681B (zh) * | 2005-07-04 | 2011-05-25 | Nxp股份有限公司 | 解调器以及解调方法 |
US8355752B2 (en) * | 2009-07-09 | 2013-01-15 | Qualcomm Incorporated | Using LO shifting to prevent a local transceiver from interfering with an FM radio |
US9391625B1 (en) * | 2015-03-24 | 2016-07-12 | Innophase Inc. | Wideband direct modulation with two-point injection in digital phase locked loops |
US9621205B2 (en) * | 2015-03-26 | 2017-04-11 | Vt Idirect, Inc. | Apparatus and method for frequency estimation ambiguity removal of a burst signal |
US10826738B2 (en) | 2019-01-07 | 2020-11-03 | Innophase Inc. | Systems and methods for maximizing power efficiency of a digital power amplifier in a polar transmitter |
US10651876B1 (en) | 2019-06-12 | 2020-05-12 | Innophase Inc. | System and method for dividing the carrier center frequency of an RF modulated signal by a non-integer divisor |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3700457C1 (de) * | 1987-01-09 | 1988-06-23 | Ant Nachrichtentech | Verfahren und Anordnung zur Synchronisation eines Empfaengers bei digitalen UEbertragungssystemen |
FR2633471B1 (fr) * | 1988-06-28 | 1990-10-05 | Trt Telecom Radio Electr | Procede de demodulation coherente d'un signal module numeriquement en phase continue et a enveloppe constante |
US4887050A (en) * | 1989-03-31 | 1989-12-12 | Motorola, Inc. | Frequency control apparatus and method for a digital radio receiver |
GB2234411A (en) * | 1989-07-03 | 1991-01-30 | Marconi Instruments Ltd | Integrated circuit for digital demodulation |
SG44793A1 (en) * | 1991-12-07 | 1997-12-19 | Philips Electronics Nv | Mobile radio receiver for a radio transmission system |
FR2711028B1 (fr) * | 1993-10-08 | 1995-11-17 | Thomson Csf | Procédé et dispositif de démodulation cohérente par blocs de signaux complexes modulés en phase. |
US5970086A (en) * | 1997-05-22 | 1999-10-19 | Itt Manufacturing Enterprises | Radio remote interface for modulating/demodulating data in a digital communication system |
-
1998
- 1998-03-03 EP EP98103735A patent/EP0940958A1/en not_active Withdrawn
-
1999
- 1999-02-11 CA CA002261442A patent/CA2261442C/en not_active Expired - Fee Related
- 1999-03-01 US US09/259,230 patent/US6075410A/en not_active Expired - Fee Related
- 1999-03-03 CN CN99102580A patent/CN1115833C/zh not_active Expired - Fee Related
- 1999-03-03 JP JP11056181A patent/JP2000059253A/ja not_active Withdrawn
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114785434A (zh) * | 2022-03-23 | 2022-07-22 | 北京睿信丰科技有限公司 | S波段频道可用性检测方法、装置、设备、系统及介质 |
CN114785434B (zh) * | 2022-03-23 | 2022-12-02 | 北京睿信丰科技有限公司 | S波段频道可用性检测方法、装置、设备、系统及介质 |
Also Published As
Publication number | Publication date |
---|---|
EP0940958A1 (en) | 1999-09-08 |
CA2261442A1 (en) | 1999-09-03 |
CA2261442C (en) | 2006-09-05 |
JP2000059253A (ja) | 2000-02-25 |
US6075410A (en) | 2000-06-13 |
CN1115833C (zh) | 2003-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1115833C (zh) | 用于数字解调调频信号的方法和设备 | |
CN1093341C (zh) | 多级干扰消除器及该干扰消除器中的实验判定方法 | |
CN1310447C (zh) | 广播接收器 | |
CN1259785C (zh) | 获取时分同步cdma(td-scdma)用户终端的载波频偏的方法 | |
CN1079619C (zh) | 相位误差检测方法及相位跟踪环路 | |
CN1098569C (zh) | 双模式调频/码分多址发射机和接收机 | |
US9106472B1 (en) | Channel state information (CSI) estimation and applications for in-band on-channel radio receivers | |
CN1754308A (zh) | 用于干扰抑制的、使用加权lsb/usb总和的相干am解调器 | |
US20080298440A1 (en) | Method and Apparatus for Implementing Seek and Scan Functions for an FM Digital Radio Signal | |
CN1585396A (zh) | 多载波调制信号的脉冲噪声降低的方法及其装置 | |
CN101047688A (zh) | 一种估计信噪比的方法及装置 | |
CN1224229C (zh) | 用于计算软决策值的64元qam解调装置和方法 | |
CN1652472A (zh) | 信号分离装置 | |
CN1801797A (zh) | 正交频分多路复用解调器及方法 | |
CN101053016A (zh) | 用于感知语音分析的频率补偿 | |
CN100340122C (zh) | 移动通信系统中的传播路径推断方法 | |
CN106998213B (zh) | 接收器电路、方法及系统 | |
CN1201508C (zh) | Cdma移动通信系统中的cdma接收设备和接收信号功率测量设备 | |
CN1275843A (zh) | 零码元位置检测方法,零码元位置检测设备和接收机 | |
GB2342546A (en) | Measuring channel state from a received signal and discriminating digital values from a received signal, suitable for use in COFDM | |
CN1897585A (zh) | I信道和q信道之间的失配校正 | |
KR20080068915A (ko) | Am 인-밴드 온 채널 라디오 수신기의 캐리어 트래킹 방법및 장치 | |
CN1370359A (zh) | 正交频分复用通信装置和传播路径估计方法 | |
JP2004336715A (ja) | 特性の異なる信号の干渉を除去するための装置およびその除去方法 | |
JP2005321270A (ja) | 電波方探装置とその方探処理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: GR Ref document number: 1050038 Country of ref document: HK |
|
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20030723 Termination date: 20100303 |