CN101977176A - 一种对不同调制信号的通用解调实现方法 - Google Patents
一种对不同调制信号的通用解调实现方法 Download PDFInfo
- Publication number
- CN101977176A CN101977176A CN2010105045161A CN201010504516A CN101977176A CN 101977176 A CN101977176 A CN 101977176A CN 2010105045161 A CN2010105045161 A CN 2010105045161A CN 201010504516 A CN201010504516 A CN 201010504516A CN 101977176 A CN101977176 A CN 101977176A
- Authority
- CN
- China
- Prior art keywords
- signal
- demodulation module
- universal demodulation
- obtains
- input
- 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
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 54
- 238000001914 filtration Methods 0.000 claims abstract description 26
- 230000011664 signaling Effects 0.000 claims description 46
- 230000015572 biosynthetic process Effects 0.000 claims description 25
- 238000003786 synthesis reaction Methods 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 22
- 238000002156 mixing Methods 0.000 claims description 11
- 238000011156 evaluation Methods 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 abstract description 5
- 238000001308 synthesis method Methods 0.000 abstract 1
- 101150118300 cos gene Proteins 0.000 description 47
- 230000008030 elimination Effects 0.000 description 5
- 238000003379 elimination reaction Methods 0.000 description 5
- 230000006854 communication Effects 0.000 description 4
- 101100113692 Caenorhabditis elegans clk-2 gene Proteins 0.000 description 3
- 102100026758 Serine/threonine-protein kinase 16 Human genes 0.000 description 3
- 101710184778 Serine/threonine-protein kinase 16 Proteins 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 101100234408 Danio rerio kif7 gene Proteins 0.000 description 2
- 101100221620 Drosophila melanogaster cos gene Proteins 0.000 description 2
- 101100398237 Xenopus tropicalis kif11 gene Proteins 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000012467 final product Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000035559 beat frequency Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明公开了一种对不同调制信号的通用解调实现方法。先采用直接数字频率合成方法,通过扩展收敛域的CORDIC算法中圆周坐标下旋转模式,产生两路本地正交载波信号,与接收机端输出信号合成同向支路数据xI(t)和正交支路数据xQ(t),再经过低通滤波器,滤除二次谐波分量,然后将两路信号输入到扩展收敛域CORDIC算法圆周坐标下向量模式的通用解调模块,根据不同的调制方式选择平方根运算值或反正切运算值经过滤波或差分,得到解调信息。本发明将多种解调方式集成在同一硬件中,灵活性强,通用性能突出,同时将复杂的算术运算转化成了移位和加法运算,可以大大降低硬件实现的复杂度,且具有较好的抗载频失配特性。
Description
技术领域
本发明涉及通信信号处理领域,尤其涉及一种对不同调制信号的通用解调实现方法。
背景技术
随着现代通信技术的不断发展,不少接收机已经突破单一的频段和模式,从传统的单信道、单频段、单模式向标准化的多信道、多波段、多模式系统发展,以满足在全部通信频段、各种信道调制方式下的通信需求。目前常用的模拟和数字调制方式有AM、FM、PM、DSB、SSB,ASK、FSK、PSK和QAM等,在常规方法中,每种调制信号都需要一个专门的硬件电路来解调,如果要实现多信号的通用解调,则接收机会极其复杂、体积重量都很大,硬件资源消耗也很大。一种常用的实现通用解调器的技术途径是利用FPGA芯片,采用数字式解调器,用不同的算法实现载波跟踪、码元定时及信号解调。其主要做法是把这些不同算法的设计文件经过编译后形成FPGA的位流文件存入计算机的存储器中,当需要使用某个解调算法时,监控单元把该文件配置到FPGA中,以实现选用的解调方式。本发明所采用的方法,用同一硬件电路对不同调制信号进行解调,且不需要每次更改FPGA配置文件,更加方便快捷。
基于坐标旋转的计算机(CORDIC)算法的基本思想是用一系列与运算基数相关的角度不断偏摆从而逼近所旋转的角度,它可以把圆周旋转、直线旋转和双曲线旋转到同一个CORDIC迭代方程里,为同一硬件可编程器件实现多功能解调奠定基础,而且由于CORDIC算法将复杂的算术运算转化成了移位和加法运算,大大降低了硬件实现的复杂度,在非扩展收敛域下,CORDIC算法所能覆盖的角度范围为[-99.8°,+99.8°],而采用扩展收敛域后,适用范围更广。
发明内容
本发明的目的是提供一种对不同调制信号的通用解调实现方法。
对不同调制信号的通用解调实现方法的步骤如下:
1)从接收机端获取信号载波信息,利用扩展收敛域CORDIC算法的圆周坐标下旋转模式,生成两路正交的本地载波信号;
2)将接收机端输出信号与两路正交的本地载波信号混频,生成两路带有信息的信号;
3)将两路带有信息的信号经过低通滤波器,滤除二次谐波分量,得到两路正交信号;
4)将两路正交信号输入通用解调模块xin输入端和通用解调模块yin输入端,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,在通用解调模块xout输出端得到两路正交信号的平方根运算值,在通用解调模块zout输出端得到两路正交信号的反正切运算值;
5)根据不同调制信号,对平方根运算值滤波或对反正切运算值差分,得到信号的解调信息。
所述的从接收机端获取信号载波信息,利用扩展收敛域CORDIC算法的圆周坐标下旋转模式,生成两路正交的本地载波信号步骤为:
根据接收机端收到的信号载波信息,选取系统时钟频率fclk和频率控制字W输入到相位累加器,生成相位信息θ,并输入到的直接数字频率合成模块zin输入端,经采用扩展收敛域CORDIC算法的圆周坐标下旋转模式迭代后,在直接数字频率合成模块xout输出端得到本地载波余弦信号,在直接数字频率合成模块yout输出端得到本地载波正弦信号。
采用扩展收敛域CORDIC算法的圆周坐标下旋转模式迭代公式为:
xj+1=xj-σjyj(1-2j-2)
j≤0时,yj+1=yj+σjxj(1-2j-2)
zj+1=zj-σjtan-1(1-2j-2)
xj+1=xj-σjyj2-j
j>0时,yj+1=yj+σjxj2-j
zj+1=zj-σjtan-12-j
式中,
其中,σj为本级运算控制字,xj、yj、zj为本级迭代器输入端,xj+1、yj+1、zj+1为下一级迭代器输入端,j=0,0,0,1,2,...N-1为迭代序列。
设初始值为:xin=1/K,yin=0,zin=θ,采用扩展收敛域后,尺度因子K=3.29352,最大收敛角度θmax=3.3141。
迭代输出值为:xout=cosθ,yout=sinθ,zout=0,即θ=2πfct,fc即为本地端载波频率。fclk和fc的关系式为:
式中,fclk是系统时钟频率,W是频率控制字,N是W的位宽,fc是直接数字频率合成的输出频率。
所述的将两路正交信号输入通用解调模块xin输入端和通用解调模块yin输入端,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,在通用解调模块xout输出端得到两路正交信号的平方根运算值,在通用解调模块zout输出端得到两路正交信号的反正切运算值步骤为:
将两路正交信号输入通用解调模块xin输入端和通用解调模块yin输入端,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,在通用解调模块xout输出端得到两路正交信号的平方根运算值,在通用解调模块zout输出端得到两路正交信号的反正切运算值。
采用扩展收敛域CORDIC算法的圆周坐标下向量模式迭代公式为:
xj+1=xj-σjyj(1-2j-2)
对于j≤0 yj+1=yj+σjxj(1-2j-2)
zj+1=zj-σjtan-1(1-2j-2)
xj+1=xj-σjyj2-j
对于j>0 yj+1=yj+σjxj2-j
zj+1=zj-σjtan-12-j
式中,
其中,σj为本级运算控制字,xj、yj、zj为本级迭代器输入端,xj+1、yj+1、zj+1为下一级迭代器输入端,j=0,0,0,1,2,...N-1为迭代序列。
初始值为xin=xI(t)、yin=xQ(t)、zin=0,xI(t)、xQ(t)为两路滤除二次谐波分量后的正交信号,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,最终迭代结果为:
yout=0
其中,通用解调模块xout输出端得到两路正交信号的平方根运算值,通用解调模块zout输出端得到两路正交信号的反正切运算值。综合考虑尺度因子和收敛角度范围的折衷办法,当我们选取迭代序列j=0,0,0,1,2,...时,最大收敛角度为:
尺度因子为:
所述的根据不同调制信号,对平方根运算值滤波或对反正切运算值差分,得到信号的解调信息的步骤为:
根据不同的调制信号,调幅AM信号由通用解调模块xout输出端信号减去平均直流分量得到解调信息,抑制载波的双边带DSB信号由通用解调模块xout输出端直接得到解调信息,抑制载波单边带SSB信号由通用解调模块xout输出端直接得到解调信息,调频FM信号由通用解调模块zout输出端信号经差分后滤除直流分量得到解调信息,M进制振幅键控MASK信号由通用解调模块zout输出端直接得到解调信息,M进制移频键控MFSK信号由通用解调模块zout输出端信号差分后得到解调信息,M进制移相键控MPSK信号由通用解调模块zout输出端直接得到解调信息,M进制正交振幅调制MQAM信号由通用解调模块xout输出端信号滤波和通用解调模块zout输出端信号差分后联合求值得到解调信息。
本发明与现有技术相比具有有益效果
1)直接数字频率合成和通用解调模块共同构成本发明,且两者都采用扩展收敛域的CORDIC算法,综合考虑尺度因子,可使收敛角度扩大到3.3141,根据需要选择合适的迭代序列,收敛角度可以更大,这就使得无需对输入信号进行预先处理,扩展了适用范围。
2)本发明采用正交化线性模型解调,这里仅列出八种常用解调方式,其他可使用正交化方法解调的单载波调制方式都能使用本发明进行解调,且不需要增加额外硬件电路,通用性强。
3)本发明对部分调制信号的频偏不敏感,即使接收端本地载波存在频偏,仍能有效地进行解调,具有良好的抗载频失配特性。
4)本发明解调精度跟迭代次数有关,增加迭代次数可以提高运算精度,当数据位宽为带符号12比特时,可表示为-4096~4095,当数据位宽为带符号16比特时,可表示为-32768~32767。而采用CORDIC算法流水线结构可进一步提高运算速度。
5)由于CORDIC算法将复杂的算术运算转化成了移位和加法运算,可以大大降低硬件实现的复杂度且易于实现,同时本发明将多种解调方式集成在同一个硬件模块中,灵活性强,通用性能突出。
附图说明
图1是对不同调制信号的通用解调实现方法原理图
图2是直接数字频率合成结构框图
图3是采用CORDIC算法实现直接数字频率合成原理图
图4是CORDIC算法一级迭代结构图
图5是采用扩展收敛域CORDIC算法实现直接数字频率合成迭代结构图
图6是采用扩展收敛域CORDIC算法实现通用解调模块迭代结构图
具体实施方式
对不同调制信号的通用解调实现方法的步骤如下:
1)从接收机端获取信号载波信息,利用扩展收敛域CORDIC算法的圆周坐标下旋转模式,生成两路正交的本地载波信号;
2)将接收机端输出信号与两路正交的本地载波信号混频,生成两路带有信息的信号;
3)将两路带有信息的信号经过低通滤波器,滤除二次谐波分量,得到两路正交信号;
4)将两路正交信号输入通用解调模块xin输入端和通用解调模块yin输入端,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,在通用解调模块xout输出端得到两路正交信号的平方根运算值,在通用解调模块zout输出端得到两路正交信号的反正切运算值;
5)根据不同调制信号,对平方根运算值滤波或对反正切运算值差分,得到信号的解调信息。
所述的从接收机端获取信号载波信息,利用扩展收敛域CORDIC算法的圆周坐标下旋转模式,生成两路正交的本地载波信号步骤为:
根据接收机端收到的信号载波信息,选取系统时钟频率fclk和频率控制字W输入到相位累加器,生成相位信息θ,并输入到的直接数字频率合成模块zin输入端,经采用扩展收敛域CORDIC算法的圆周坐标下旋转模式迭代后,在直接数字频率合成模块xout输出端得到本地载波余弦信号,在直接数字频率合成模块yout输出端得到本地载波正弦信号。
采用扩展收敛域CORDIC算法的圆周坐标下旋转模式迭代公式为:
xj+1=xj-σjyj(1-2j-2)
j≤0时,yj+1=yj+σjxj(1-2j-2)
zj+1=zj-σjtan-1(1-2j-2)
xj+1=xj-σjyj2-j
j>0时,yj+1=yj+σjxj2-j
zj+1=zj-σjtan-12-j
式中,
其中,σj为本级运算控制字,xj、yj、zj为本级迭代器输入端,xj+1、yj+1、zj+1为下一级迭代器输入端,j=0,0,0,1,2,...N-1为迭代序列。
设初始值为:xin=1/K,yin=0,zin=θ,采用扩展收敛域后,尺度因子K=3.29352,最大收敛角度θmax=3.3141。
迭代输出值为:xout=cosθ,yout=sinθ,zout=0,即θ=2πfct,fc即为本地端载波频率。fclk和fc的关系式为:
式中,fclk是系统时钟频率,W是频率控制字,N是W的位宽,fc是直接数字频率合成的输出频率。
所述的将两路正交信号输入通用解调模块xin输入端和通用解调模块yin输入端,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,在通用解调模块xout输出端得到两路正交信号的平方根运算值,在通用解调模块zout输出端得到两路正交信号的反正切运算值步骤为:
将两路正交信号输入通用解调模块xin输入端和通用解调模块yin输入端,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,在通用解调模块xout输出端得到两路正交信号的平方根运算值,在通用解调模块zout输出端得到两路正交信号的反正切运算值。
采用扩展收敛域CORDIC算法的圆周坐标下向量模式迭代公式为:
xj+1=xj-σjyj(1-2j-2)
对于j≤0 yj+1=yj+σjxj(1-2j-2)
zj+1=zj-σjtan-1(1-2j-2)
xj+1=xj-σjyj2-j
对于j>0 yj+1=yj+σjxj2-j
zj+1=zj-σjtan-12-j
式中,
其中,σj为本级运算控制字,xj、yj、zj为本级迭代器输入端,xj+1、yj+1、zj+1为下一级迭代器输入端,j=0,0,0,1,2,...N-1为迭代序列。
初始值为xin=xI(t)、yin=xQ(t)、zin=0,xI(t)、xQ(t)为两路滤除二次谐波分量后的正交信号,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,最终迭代结果为:
yout=0
其中,通用解调模块xout输出端得到两路正交信号的平方根运算值,通用解调模块zout输出端得到两路正交信号的反正切运算值。综合考虑尺度因子和收敛角度范围的折衷办法,当我们选取迭代序列j=0,0,0,1,2,...时,最大收敛角度为:
尺度因子为:
所述的根据不同调制信号,对平方根运算值滤波或对反正切运算值差分,得到信号的解调信息的步骤为:
根据不同的调制信号,调幅AM信号由通用解调模块xout输出端信号减去平均直流分量得到解调信息,抑制载波的双边带DSB信号由通用解调模块xout输出端直接得到解调信息,抑制载波单边带SSB信号由通用解调模块xout输出端直接得到解调信息,调频FM信号由通用解调模块zout输出端信号经差分后滤除直流分量得到解调信息,M进制振幅键控MASK信号由通用解调模块zout输出端直接得到解调信息,M进制移频键控MFSK信号由通用解调模块zout输出端信号差分后得到解调信息,M进制移相键控MPSK信号由通用解调模块zout输出端直接得到解调信息,M进制正交振幅调制MQAM信号由通用解调模块xout输出端信号滤波和通用解调模块zout输出端信号差分后联合求值得到解调信息。具体分析为:
(1)调幅AM信号解调
设AM调制信号为s(t)=[1+Qm(t)]cos2πfct
式中Q为调幅系数,m(t)为调制信息,fc为载波频率。假设接收端本地载波有一个频差为fc+Δf,经过直接频率合成后两路本地正交载频与接收信号混频后可得
滤除二次谐波并整理后可得
I=[1+Qm(t)]cos2πΔft,
Q=[1+Qm(t)]sin2πΔft,
经通用解调模块做平方根运算得
调制信号m(t)由1+Qm(t)减去平均直流分量1或直接将1+Qm(t)滤去直流得到。
(2)抑制载波的双边带DSB信号解调
设DSB调制信号的时域表达式为
sDBS(t)=m(t)cosωct
其中m(t)是调制信息,ωc是载波频率。假设在接收机端本地载波有一个频差,与直接频率合成两路本地载波正交信号混频后,可得
滤出二次谐波并整理后可得
即可直接得到调制信号。
(3)抑制载波单边带SSB信号解调
同样对信号分别乘以cos(ωct)、sin(ωct)两路本地载波正交信号可得
整理后可得
并滤去二次谐波,预处理后的I路信号为单边带信号。
(4)调频FM、调相PM信号解调
FM信号的时域表达式为SFM(t)=Acos[2πfct+KF∫m(t)dt],
式中,m(t)为调制信号,KF为调频灵敏度,fc为载波频率,在接收端信号经与本地载波混频预处理后可以得到两路信号I、Q,设载波频偏为Δf,
将I、Q两路信号做反正切运算,可以得到幅角信息,即实现信号的鉴相
对上式做微分运算后就可以求得调制信号:
去除直流分量后在进行适当的幅度调整就可恢复出调制信号。
相位信号为FM信号的积分,在离散运算时,对FM解调后的信号Δφ=φn-φn-1求累加和即为相位。
Δφn=φn-φn-1,Δφn+1=φn+1-φn
Δφn+1+Δφn=φn+1-φn+φn-φn-1=φn+1-φn-1
相位
(5)MASK信号解调
MASK信号的时域表达式为
eMASK(t)=s(t)cosωct
式中,s(t)为M进制单极性非归零码,同样经过DDS两路信号后I、Q中包含二进制码信息,设接收端载波为ωc+Δω,Δω为频偏,可得
滤去二次谐波,并整理后可得
经通用解调模块做平方根运算即可得调制信息
(6)M进制移频键控MFSK信号解调
MFSK信号时域表达式为
eMFSK(t)=Acos(ωc±iΔω)t,i=1,3,…,M/2
式中,A为MFSK信号的包络,Δω=2πΔf是对应中心频率的最小偏移量,2Δω是相邻传输在载频间的间隔。
接收机端信号与本地正交载波混频后可得
滤除二次谐波后得
(7)MPSK信号解调
MPSK信号时域表达式为
式中,Ts为M进制信号的码元宽度,将式子右边展开得
滤除二次谐波分量可得
(8)M进制正交振幅调制(MQAM)
MQAM信号的时域表达式为
eMQAM(t)=I(t)cosωct+Q(t)sinωct
式中,I(t)、Q(t)为双极性多进制信号,ωc为载波频率,同样在接收机端与两路本地载波混频后可得
滤除二次谐波后为
由于I(t)、Q(t)为双极性多进制信号,I、Q信号经扩展收敛域后可得幅度信息及幅角信息:
如图1所示,接收机端信号与直接数字频率合成的两路正交本地载波信号混频,生成两路带有调制信息的信号,经低通滤波器滤去二次谐波后输入通用解调模块,经过扩展收敛域CORDIC算法即可得到不同调制信号的解调信息。
如图2所示,直接数字频率合成(DDS)的思路是按一定的时钟节拍从存放有正弦函数表的ROM中读出这些离散的代表正弦幅值的二进制数,然后经过D/A变换并滤波,得到一个模拟的正弦波,改变读数的节拍频率或者取点的个数,就可以改变正弦波的频率,图2中fclk和fc的关系式为:
如图3所示,相位累加器和扩展收敛域的CORDIC算法相连,输入信号fclk分别输入到相位累加器和CORDIC算法模块中,频率控制字W与输入频率信号经相位累加器生成相位信息θ,经扩展收敛域的CORDIC算法后生成两路相互正交的本地载波余弦信号cosθ和正弦信号sinθ,θ=2πfct。
如图4所示,CORDIC算法一级迭代结构为:初始值通过一路进入到寄存器中,在时钟到来时,x、y寄存器中的数据分成两路,一路直接送到加减法器x、y,另一路分别送到两个移位器,经本级移位控制字j移位后送到与本数据相异的y、x加减法器,在本级运算控制字σj控制下生成数据存储到下一级x、y寄存器中,z寄存器中的数据直接送到独立加减法器z,与运算控制字σj和本级旋转角度αj生成数据存储到下一级z寄存器,y、z加减法器数据经选择器生成下一级运算控制字σj+1。
如图5所示,采用扩展收敛域CORDIC算法DDS迭代结构为:前一级x、y、z寄存器分别直接与后一级x、y、z寄存器相连,且前一级运算控制字σ直接控制下一级运算,并由本级迭代序列j和旋转角度αj控制本级的移位器和加减法器,扩展收敛域后迭代序列为j=0,0,0,1,2,...N-1,迭代公式为
xj+1=xj-σjyj(1-2j-2)
j≤0时,yj+1=yj+σjxj(1-2j-2)
zj+1=zj-σjtan-1(1-2j-2)
xj+1=xj-σjyj2-j
j>0时,yj+1=yj+σjxj2-j
zj+1=zj-σjtan-12-j
式中,
设初始值为:xin=1/K,yin=0,zin=θ,采用扩展收敛域后,尺度因子K=3.29352,最大收敛角度θmax=3.3141(189.9°),
迭代输出值为:xout=cosθ,yout=sinθ,zout=0,即θ=2πfct,fc即为本地端载波频率。
如图6所示,采用扩展收敛域CORDIC算法通用解调模块迭代结构为:前一级x、y、z寄存器分别直接与后一级x、y、z寄存器相连,且前一级运算控制字σ直接控制下一级运算,并由本级迭代序列j和旋转角度αj控制本级的移位器和加减法器,扩展收敛域后迭代序列为j=0,0,0,1,2,...N-1,迭代公式为
xj+1=xj-σjyj(1-2j-2)
对于j≤0 yj+1=yj+σjxj(1-2j-2)
zj+1=zj-σjtan-1(1-2j-2)
xj+1=xj-σjyj2-j
对于j>0 yj+1=yj+σjxj2-j
zj+1=zj-σjtan-12-j
式中,
初始值为xin=xI(t)、yin=xQ(t)、zin=0,采用扩展收敛域的CORDICS算法,综合考虑尺度因子和收敛角度范围的折衷办法,当我们选取迭代序列j=0,0,0,1,2,...时,最大收敛角度为:
尺度因子为:
最终迭代结果为:
yf=0
即经过扩展收敛域CORDIC算法迭代后,X输出端为两路正交信号幅度信息,Z输出端为两路正交信号幅角信息。
实施例
对调幅AM信号的解调实现方法的步骤如下:
设AM调制信号为s(t)=[1+Qm(t)]cos2πfct,式中Q为调幅系数,m(t)为调制信息,fc为载波频率。
1)根据接收机端收到的信号载波信息,选取系统时钟频率fclk和频率控制字W输入到相位累加器,生成相位信息θ,并输入到的直接数字频率合成模块zin输入端,经采用扩展收敛域CORDIC算法的圆周坐标下旋转模式迭代后,假设接收端本地载波有一个频差为fc+Δf,在直接数字频率合成模块xout输出端得到本地载波余弦信号cos(2π(fc+Δf)t),在直接数字频率合成模块yout输出端得到本地载波正弦信号sin(2π(fc+Δf)t)。
2)将接收机端输出信号与两路正交的本地载波信号混频,生成两路带有信息的信号:
3)将两路带有信息的信号经过低通滤波器,滤除二次谐波分量,得到两路正交信号;
I=[1+Qm(t)]cos2πΔft,
Q=[1+Qm(t)]sin2πΔft,
4)将两路正交信号输入通用解调模块xin输入端和通用解调模块yin输入端,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,在通用解调模块xout输出端得到两路正交信号的平方根运算值
5)调制信号m(t)由1+Qm(t)减去平均直流分量1或直接将1+Qm(t)滤去直流得到。
Claims (4)
1.一种对不同调制信号的通用解调实现方法,其特征在于它的步骤如下:
1)从接收机端获取信号载波信息,利用扩展收敛域CORDIC算法的圆周坐标下旋转模式,生成两路正交的本地载波信号;
2)将接收机端输出信号与两路正交的本地载波信号混频,生成两路带有信息的信号;
3)将两路带有信息的信号经过低通滤波器,滤除二次谐波分量,得到两路正交信号;
4)将两路正交信号输入通用解调模块xin输入端和通用解调模块yin输入端,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,在通用解调模块xout输出端得到两路正交信号的平方根运算值,在通用解调模块zout输出端得到两路正交信号的反正切运算值;
5)根据不同调制信号,对平方根运算值滤波或对反正切运算值差分,得到信号的解调信息。
2.根据权利要求1所述的一种对不同调制信号的通用解调实现方法,其特征在于所述的从接收机端获取信号载波信息,利用扩展收敛域CORDIC算法的圆周坐标下旋转模式,生成两路正交的本地载波信号步骤为:
根据接收机端收到的信号载波信息,选取系统时钟频率fclk和频率控制字W输入到相位累加器,生成相位信息θ,并输入到的直接数字频率合成模块zin输入端,经采用扩展收敛域CORDIC算法的圆周坐标下旋转模式迭代后,在直接数字频率合成模块xout输出端得到本地载波余弦信号,在直接数字频率合成模块yout输出端得到本地载波正弦信号。
采用扩展收敛域CORDIC算法的圆周坐标下旋转模式迭代公式为:
xj+1=xj-σjyj(1-2j-2)
j≤0时,yj+1=yj+σjxj(1-2j-2)
zj+1=zj-σjtan-1(1-2j-2)
xj+1=xj-σjyj2-j
j>0时,yj+1=yj+σjxj2-j
zj+1=zj-σjtan-12-j
式中,
其中,σj为本级运算控制字,xj、yj、zj为本级迭代器输入端,xj+1、yj+1、 zj+1为下一级迭代器输入端,j=0,0,0,1,2,...N-1为迭代序列。
设初始值为:xin=1/K,yin=0,zin=θ,采用扩展收敛域后,尺度因子K=3.29352,最大收敛角度θmax=3.3141。
迭代输出值为:xout=cosθ,yout=sinθ,zout=0,即θ=2πfct,fc即为本地端载波频率。fclk和fc的关系式为:
式中,fclk是系统时钟频率,W是频率控制字,N是W的位宽,fc是直接数字频率合成的输出频率。
3.根据权利要求1所述的一种对不同调制信号的通用解调实现方法,其特征在于所述的将两路正交信号输入通用解调模块xin输入端和通用解调模块yin输入端,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,在通用解调模块xout输出端得到两路正交信号的平方根运算值,在通用解调模块zout输出端得到两路正交信号的反正切运算值步骤为:
将两路正交信号输入通用解调模块xin输入端和通用解调模块yin输入端,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,在通用解调模块xout输出端得到两路正交信号的平方根运算值,在通用解调模块zout输出端得到两路正交信号的反正切运算值。
采用扩展收敛域CORDIC算法的圆周坐标下向量模式迭代公式为:
xj+1=xj-σjyj(1-2j-2)
对于j≤0 yj+1=yj+σjxj(1-2j-2)
zj+1=zj-σjtan-1(1-2j-2)
xj+1=xj-σjyj2-j
对于j>0 yj+1=yj+σjxj 2-j
zj+1=zj-σjtan-12-j
式中,
其中,σj为本级运算控制字,xj、yj、zj为本级迭代器输入端,xj+1、yj+1、zj+1为下一级迭代器输入端,j=0,0,0,1,2,...N-1为迭代序列。
初始值为xin=xI(t)、yin=xQ(t)、zin=0,xI(t)、xQ(t)为两路滤除二次谐波分量后的正交信号,采用扩展收敛域CORDIC算法的圆周坐标下向量模式,最终迭代 结果为:
yout=0
其中,通用解调模块xout输出端得到两路正交信号的平方根运算值,通用解调模块zout输出端得到两路正交信号的反正切运算值。综合考虑尺度因子和收敛角度范围的折衷办法,当我们选取迭代序列j=0,0,0,1,2,...时,最大收敛角度为:
4.根据权利要求1所述的一种对不同调制信号的通用解调实现方法,其特征在于所述的根据不同调制信号,对平方根运算值滤波或对反正切运算值差分,得到信号的解调信息的步骤为:
根据不同的调制信号,调幅AM信号由通用解调模块xout输出端信号减去平均直流分量得到解调信息,抑制载波的双边带DSB信号由通用解调模块xout输出端直接得到解调信息,抑制载波单边带SSB信号由通用解调模块xout输出端直接得到解调信息,调频FM信号由通用解调模块zout输出端信号经差分后滤除直流分量得到解调信息,M进制振幅键控MASK信号由通用解调模块zout输出端直接得到解调信息,M进制移频键控MFSK信号由通用解调模块zout输出端信号差分后得到解调信息,M进制移相键控MPSK信号由通用解调模块zout输出端直接得到解调信息,M进制正交振幅调制MQAM信号由通用解调模块xout输出端信号滤波和通用解调模块zout输出端信号差分后联合求值得到解调信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010504516 CN101977176B (zh) | 2010-10-12 | 2010-10-12 | 一种对不同调制信号的通用解调实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010504516 CN101977176B (zh) | 2010-10-12 | 2010-10-12 | 一种对不同调制信号的通用解调实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101977176A true CN101977176A (zh) | 2011-02-16 |
CN101977176B CN101977176B (zh) | 2013-02-27 |
Family
ID=43577021
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010504516 Expired - Fee Related CN101977176B (zh) | 2010-10-12 | 2010-10-12 | 一种对不同调制信号的通用解调实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101977176B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102200550A (zh) * | 2011-05-31 | 2011-09-28 | 中国航空无线电电子研究所 | 一种用于高精度测量相位差的延迟正交数字中频鉴相方法 |
CN102437900A (zh) * | 2011-12-28 | 2012-05-02 | 大唐微电子技术有限公司 | 一种数字对讲机以及用于降低对讲机误码率的方法 |
CN102685066A (zh) * | 2012-05-22 | 2012-09-19 | 杭州电子科技大学 | 大频偏条件短时猝发高扩频通信接收系统的载波恢复方法 |
CN102946370A (zh) * | 2012-12-05 | 2013-02-27 | 天津光电通信技术有限公司 | 一种基于fpga实现fm调频和解调数字逻辑电路的方法 |
CN104483554A (zh) * | 2014-12-31 | 2015-04-01 | 思澜科技(成都)有限公司 | 用于生物电阻抗测量的数字相位解调方法及其系统 |
CN105703773A (zh) * | 2014-11-28 | 2016-06-22 | 中国科学院沈阳自动化研究所 | 一种旋转变压器解码方法 |
CN106506116A (zh) * | 2016-10-25 | 2017-03-15 | 武汉大学 | 一种基于fpga的ais基带信号处理系统及方法 |
CN108011616A (zh) * | 2017-12-26 | 2018-05-08 | 中山大学花都产业科技研究院 | 低复杂度iir数字移频移相器及方法 |
CN109067679A (zh) * | 2018-08-31 | 2018-12-21 | 北京北广科技股份有限公司 | 一种抵消频偏的dsb解调方法 |
CN110430154A (zh) * | 2019-08-05 | 2019-11-08 | 上海无线电设备研究所 | 基于fpga的多路副载波数字鉴频与解调方法 |
CN110429985A (zh) * | 2018-05-15 | 2019-11-08 | 东南大学 | 一种全集成低成本高速高精度异相调制器 |
CN110474859A (zh) * | 2019-08-28 | 2019-11-19 | 哈尔滨工程大学 | 一种零中频基带2fsk信号相干解调方法 |
CN111865862A (zh) * | 2020-07-17 | 2020-10-30 | 电子科技大学 | 一种基于cordic算法的ofdm基带信号生成及解调方法 |
CN111901272A (zh) * | 2020-07-03 | 2020-11-06 | 航天恒星科技有限公司 | 多进制双极性正交波形调制信号的接收处理方法及装置 |
CN112671682A (zh) * | 2020-12-14 | 2021-04-16 | 天津光电通信技术有限公司 | 一种一体化数字解调方法及装置 |
CN115189991A (zh) * | 2022-05-20 | 2022-10-14 | 南昌大学 | 一种用于相位调制同步解调的装置及方法 |
CN117348685A (zh) * | 2023-10-25 | 2024-01-05 | 灿芯半导体(成都)有限公司 | 一种高精度dds发生器的信号处理方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101499775A (zh) * | 2009-02-20 | 2009-08-05 | 武汉大学 | 利用cordic对线性调频信号进行脉冲压缩处理的方法和装置 |
-
2010
- 2010-10-12 CN CN 201010504516 patent/CN101977176B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101499775A (zh) * | 2009-02-20 | 2009-08-05 | 武汉大学 | 利用cordic对线性调频信号进行脉冲压缩处理的方法和装置 |
Non-Patent Citations (2)
Title |
---|
《集成电路应用》 20100131 何晓华等 "基于扩展收敛域CORDIC算法的FM 信号数字化解调实现" 第54-55页 1 , 第1期 * |
何晓华等: ""基于扩展收敛域CORDIC算法的FM 信号数字化解调实现"", 《集成电路应用》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102200550A (zh) * | 2011-05-31 | 2011-09-28 | 中国航空无线电电子研究所 | 一种用于高精度测量相位差的延迟正交数字中频鉴相方法 |
CN102200550B (zh) * | 2011-05-31 | 2013-03-27 | 中国航空无线电电子研究所 | 一种用于高精度测量相位差的延迟正交数字中频鉴相方法 |
CN102437900A (zh) * | 2011-12-28 | 2012-05-02 | 大唐微电子技术有限公司 | 一种数字对讲机以及用于降低对讲机误码率的方法 |
CN102437900B (zh) * | 2011-12-28 | 2014-08-20 | 大唐微电子技术有限公司 | 一种数字对讲机以及用于降低对讲机误码率的方法 |
CN102685066B (zh) * | 2012-05-22 | 2014-11-05 | 杭州电子科技大学 | 大频偏条件短时猝发高扩频通信接收系统的载波恢复方法 |
CN102685066A (zh) * | 2012-05-22 | 2012-09-19 | 杭州电子科技大学 | 大频偏条件短时猝发高扩频通信接收系统的载波恢复方法 |
CN102946370A (zh) * | 2012-12-05 | 2013-02-27 | 天津光电通信技术有限公司 | 一种基于fpga实现fm调频和解调数字逻辑电路的方法 |
CN102946370B (zh) * | 2012-12-05 | 2015-04-22 | 天津光电通信技术有限公司 | 一种基于fpga实现fm调频和解调数字逻辑电路的方法 |
CN105703773A (zh) * | 2014-11-28 | 2016-06-22 | 中国科学院沈阳自动化研究所 | 一种旋转变压器解码方法 |
CN104483554A (zh) * | 2014-12-31 | 2015-04-01 | 思澜科技(成都)有限公司 | 用于生物电阻抗测量的数字相位解调方法及其系统 |
CN104483554B (zh) * | 2014-12-31 | 2023-07-04 | 思澜科技(成都)有限公司 | 用于生物电阻抗测量的数字相位解调方法及其系统 |
CN106506116A (zh) * | 2016-10-25 | 2017-03-15 | 武汉大学 | 一种基于fpga的ais基带信号处理系统及方法 |
CN108011616A (zh) * | 2017-12-26 | 2018-05-08 | 中山大学花都产业科技研究院 | 低复杂度iir数字移频移相器及方法 |
CN108011616B (zh) * | 2017-12-26 | 2021-06-22 | 中山大学花都产业科技研究院 | 低复杂度iir数字移频移相器及方法 |
CN110429985A (zh) * | 2018-05-15 | 2019-11-08 | 东南大学 | 一种全集成低成本高速高精度异相调制器 |
CN109067679A (zh) * | 2018-08-31 | 2018-12-21 | 北京北广科技股份有限公司 | 一种抵消频偏的dsb解调方法 |
CN110430154B (zh) * | 2019-08-05 | 2021-11-02 | 上海无线电设备研究所 | 基于fpga的多路副载波数字鉴频与解调方法 |
CN110430154A (zh) * | 2019-08-05 | 2019-11-08 | 上海无线电设备研究所 | 基于fpga的多路副载波数字鉴频与解调方法 |
CN110474859A (zh) * | 2019-08-28 | 2019-11-19 | 哈尔滨工程大学 | 一种零中频基带2fsk信号相干解调方法 |
CN111901272A (zh) * | 2020-07-03 | 2020-11-06 | 航天恒星科技有限公司 | 多进制双极性正交波形调制信号的接收处理方法及装置 |
CN111901272B (zh) * | 2020-07-03 | 2023-06-06 | 航天恒星科技有限公司 | 多进制双极性正交波形调制信号的接收处理方法及装置 |
CN111865862A (zh) * | 2020-07-17 | 2020-10-30 | 电子科技大学 | 一种基于cordic算法的ofdm基带信号生成及解调方法 |
CN112671682A (zh) * | 2020-12-14 | 2021-04-16 | 天津光电通信技术有限公司 | 一种一体化数字解调方法及装置 |
CN115189991A (zh) * | 2022-05-20 | 2022-10-14 | 南昌大学 | 一种用于相位调制同步解调的装置及方法 |
CN115189991B (zh) * | 2022-05-20 | 2024-06-11 | 南昌大学 | 一种用于相位调制同步解调的装置及方法 |
CN117348685A (zh) * | 2023-10-25 | 2024-01-05 | 灿芯半导体(成都)有限公司 | 一种高精度dds发生器的信号处理方法 |
CN117348685B (zh) * | 2023-10-25 | 2024-05-10 | 灿芯半导体(成都)有限公司 | 一种dds发生器的信号处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101977176B (zh) | 2013-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101977176B (zh) | 一种对不同调制信号的通用解调实现方法 | |
US8107551B2 (en) | Systems and methods for signal modulation and demodulation using phase | |
Farzana et al. | High Performance VLSI Architecture for Advanced QPSK Modems | |
US7400904B2 (en) | Dual mode AM-FM receiver with I-Q decoding, such as for utility data collection | |
CN100571230C (zh) | 统一的正交二元偏移键控调制和解调方法 | |
CN101223701A (zh) | 用于功率和带宽高效通信的多载波恒包络信号方案 | |
CN105403865B (zh) | 多载波信号恒定包络调制方法 | |
Song et al. | Design and implement of QPSK modem based on FPGA | |
CN102209056A (zh) | 一种导航信号调制方法 | |
CN112350970B (zh) | 一种多相位频移键控调制、解调方法及设备 | |
CN107135174A (zh) | 信号发送装置、载波相位恢复装置及方法 | |
US8965290B2 (en) | Amplitude enhanced frequency modulation | |
JP3185867B2 (ja) | 誤差検出方法および装置、信号復調方法および装置 | |
CN102655491A (zh) | 相干解调频移键控调制信号的频率偏移估计方法及系统 | |
CN101834818B (zh) | Gmsk调制装置及方法 | |
CN102946370B (zh) | 一种基于fpga实现fm调频和解调数字逻辑电路的方法 | |
CN102014092A (zh) | 一种基于级联模式的四进制msk调制方法及装置 | |
CN105915481B (zh) | 多路高速宽带信号模拟相位调制方法 | |
Babu et al. | GMSK modulator for GSM system, an economical implementation on FPGA | |
Al Safi et al. | Toward digital transmitters with amplitude shift keying and quadrature amplitude modulators implementation examples | |
US8005161B2 (en) | Method, hardware product, and computer program product for performing high data rate wireless transmission | |
CN110011687A (zh) | Oqpsk数据的msk收发器 | |
CN100438520C (zh) | 1比特采样的差分四相相移键控的解调电路 | |
CN203014854U (zh) | 偏移正交相移键控信号发射机 | |
CN106059708A (zh) | 一种多码率数据无线传输系统 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130227 Termination date: 20141012 |
|
EXPY | Termination of patent right or utility model |