CN102664588A - 用于频率调制和幅度调制信号的数字化解调方法及装置 - Google Patents

用于频率调制和幅度调制信号的数字化解调方法及装置 Download PDF

Info

Publication number
CN102664588A
CN102664588A CN2012101499272A CN201210149927A CN102664588A CN 102664588 A CN102664588 A CN 102664588A CN 2012101499272 A CN2012101499272 A CN 2012101499272A CN 201210149927 A CN201210149927 A CN 201210149927A CN 102664588 A CN102664588 A CN 102664588A
Authority
CN
China
Prior art keywords
match
frequency
sequence
signal
sinusoidal
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.)
Pending
Application number
CN2012101499272A
Other languages
English (en)
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.)
Beijing Changcheng Institute of Metrology and Measurement AVIC
Original Assignee
Beijing Changcheng Institute of Metrology and Measurement AVIC
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 Beijing Changcheng Institute of Metrology and Measurement AVIC filed Critical Beijing Changcheng Institute of Metrology and Measurement AVIC
Priority to CN2012101499272A priority Critical patent/CN102664588A/zh
Publication of CN102664588A publication Critical patent/CN102664588A/zh
Pending legal-status Critical Current

Links

Images

Abstract

本发明涉及一种用于频率调制信号和幅度调制信号的数字化解调方法及装置,属于电子测量技术领域。本发明虚拟仪器方式实现频率调制和幅度调制信号的数字化解调测量,利用小于一个波形周期的局部波形获得载波正弦信号序列的幅度、频率等局域模型参数值,将小于一个波形周期的正弦模型在实际测量曲线波形上进行滑动拟合,从而实现调制信号的数字化解调。另外,其具有时间分辨力高、鲁棒性好、解调失真小、绝对收敛、快速解调等特征。本发明时间分辨力高、鲁棒性好、算法绝对收敛、可以方便进行溯源校准。

Description

用于频率调制和幅度调制信号的数字化解调方法及装置
技术领域
本发明涉及一种用于频率调制信号和幅度调制信号的数字化解调方法及装置,属于电子测量技术领域。
背景技术
调制与解调技术是无线电基本技术之一,广泛应用于广播、电视、通信、导航、测量、测试、仪器、仪表等行业领域中。
调幅信号解调基本原理是使用峰值检波硬件实现,调频信号的解调原理之一即是使用鉴频器将频率变化转换成幅度变化,然后再使用峰值检波手段完成。特点是解调实时性好、速度快,不足是解调电路复杂、误差环节因素较多、解调失真较大、解调灵敏度受硬件条件限制不易达到很高。在如彩色电视信号等复合调制信号的解调中过程和原理更加复杂,需要分别进行处理和分别解调。
计算机和数字化采样技术的进展,产生并推动了软件无线电技术的发展,使得人们可以使用数字化方式进行调制信号解调,不仅可以达到较低的解调失真度,还能方便地进行量值溯源,解决调制解调仪器设备的计量校准问题。在数字化解调过程中,最大的困难在于解调的时间分辨力问题,一些解调方法如过零点检测法解调频率调制信号,会遇到分辨力不能低于载波周期并且会随时间变化的问题,其本质原因是算法在少于一个周期波形的情况下无法有效获得其载波的模型参数。解调精度较低,失真较大。
发明内容
本发明的目的是为了克服现有技术的局陷,提出一种用于频率调制和幅度调制信号的数字化解调方法及装置,使用不足一个波形周期的正弦波形(例如十分之一个波形周期)获取已调信号局部波形序列的载波幅度和频率模型值,装置借助于高速数据采集和量化技术,以数字化方式同时实现幅度调制信号和频率调制信号的精确解调,并最终解决幅度调制信号和频率调制信号的溯源问题。
本发明是通过以下技术方案实现的。
本发明的一种用于频率调制和幅度调制信号的数字化解调方法,具体步骤为:
1)对待测频率调制信号或者幅度调制信号进行信号调理,利用模数转换电路对调理后的信号进行采样,将采样得到的波形数据序列yi,i=1,2,...,n送入计算机,其中i表示同步采样序列中采样点序号;
模数转换电路采集波形的序列长度n以及采样速率v根据被测调制信号的载波频率设定,设定原则为保证每个载波波形周期内要多于20个采样点;
2)首先,在波形采集序列前缘截取少于一个载波周期的包含m1个采样点的波形段数据序列yi,i=1,2,...,m1;假设截取的波形数据序列近似为正弦波,对该正弦波信号进行数学描述,其波形为
Figure BSA00000717805400021
其中,A1为正弦波幅值,f1为正弦波频率值,
Figure BSA00000717805400022
为正弦波初始相位值;
3)用计算机对步骤2)截取的波形数据序列进行正弦拟合,获得拟合正弦波的幅度A1、频率f1,具体过程为:
3.1 在波形数据序列前缘截取少于一个载波周期的包含m1个采样点的波形段数据序列yi,i=1,2,...,m1,则
yi=y(ti)=y((i-1)·Δτ),i=1,2,...,m1
其中,采样时间间隔Δτ=1/v,v为采样速率;
3.2 用计算机对截取的数据序列yi(i=1,2,...,m1)进行正弦波形四参数拟合,获得拟合信号:
Figure BSA00000717805400023
其中,
Figure BSA00000717805400024
为正弦波幅度的拟合值,
Figure BSA00000717805400025
为正弦波角频率的拟合值,
Figure BSA00000717805400026
为正弦波初始相位的拟合值,
Figure BSA00000717805400027
为正弦波直流分量值的拟合值;
频率拟合值
Figure BSA00000717805400028
f ^ 1 = ω ^ 1 · v 2 π ;
4)存储输出正弦模型幅度A1的估计值
Figure BSA000007178054000210
正弦模型频率f1的估计值
5)以小于
Figure BSA000007178054000212
的周期长度对应的序列时间长度截取下一段拟合序列长度m2,该拟合序列的中心位置从m1/2处向后移动一个采样点,在新截取的拟合序列上执行步骤3)、4)所述过程,获得输出正弦模型的幅度正弦模型的频率
Figure BSA00000717805400032
6)按照步骤5)的方式重复向后推动截取拟合过程,直至数据序列达到终点,获得幅度解调波形序列:
Figure BSA00000717805400033
为幅度调制信号的解调波形输出,获得频率解调波形序列:
Figure BSA00000717805400034
为频率调制信号的解调波形输出。
所述正弦波形四参数拟合具体过程为:
(1)设定初始参数,具体包括:待测的正弦信号的频率预估值为f0,角频率预估值ω0=2πf0/v,波形采集序列所含信号不足一个周期且个数为p,波形采集序列占用时间长度为τ,则f0≤1/τ,选取另一个足够小的正数因子q=1×10-5,使得f0>q/τ,此时f0∈[q/τ,2/τ];
(2)设定拟合迭代停止条件为一个接近于0的足够小的正数he=1×10-20
(3)从已知时刻t1,t2,...,tn的正弦波采集样本y1,y2,...yn,使用计点法获得信号波形占用时间长度为τ=(n-1)/v,确定目标频率f0的存在区间[q/τ,2/τ];
(4)确定迭代左边界频率fL=q/τ和迭代左边界角频率ωL=2πfL/v以及迭代右边界频率.fR=2/τ和迭代右边界角频率ωR=2πfR/v,令中值角频率ωM=(ωRL)/2;
(5)在迭代左边界角频率ωL、迭代右边界角频率ωR和中值角频率ωM上分别利用频率已知的三参数正弦波形拟合公式计算各自的拟合残差ρ(ωL)、ρ(ωR)和ρ(ωM);
(6)若ρ(ωL)<η·ρ(ωM),则令迭代右边界角频率ωR=ωM,迭代左边界角频率ωL不变,重复执行步骤(5)~步骤(6),其中η为判据因子且η取值范围为1~1.5;
(7)若ρ(ωL)≥η·ρ(ωM),则必有ωR<2ω0,确定左边界角频率为ωL,右边界角频率ωR,按照优选法原则,选取两个中值角频率分别为
ωM=ωL+0.618×(ωRL)
ωT=ωR-0.618×(ωRL);
(8)在左边界角频率ωL上执行频率已知的三参数正弦曲线拟合获得拟合幅度AL、拟合相位
Figure BSA00000717805400035
拟合直流分量DL、拟合残差ρL,在右边界角频率ωR上执行频率已知的三参数正弦曲线拟合获得拟合幅度AR、拟合相位
Figure BSA00000717805400036
拟合直流分量DR、拟合残差ρR,在中值角频率ωM上执行频率已知的三参数正弦曲线拟合获得拟合幅度AM、拟合相位
Figure BSA00000717805400041
拟合直流分量DM、拟合残差ρM,在中值角频率ωT上执行频率已知的三参数正弦曲线拟合获得拟合幅度AT、拟合相位
Figure BSA00000717805400042
拟合直流分量DT、拟合残差ρT
(9)若ρM<ρT,则最小拟合残差ρ=ρM,有ω0∈[ωT,ωR],参量更新为ωL=ωT,ωT=ωM,ωM=ωL+0.618×(ωRL);若ρM≥ρT,则最小拟合残差ρ=ρT,有ω0∈[ωL,ωM],参量更新为ωR=ωM,ωM=ωT,ωT=ωR-0.618×(ωRL);
(10)判定若|(ρMT)/ρT|≥he,则重复步骤(5)~步骤(10);若|ρMT)/ρT|<he,则停止迭代,此时又分为两种情况:若最小拟合残差ρ=ρT,获得四参数拟合正弦曲线参数为
Figure BSA00000717805400043
拟合过程结束;若最小拟合残差ρ=ρM,获得四参数拟合正弦曲线参数为 拟合过程结束。
本发明的一种用于频率调制和幅度调制信号的数字化解调装置,用于频率调制信号或者幅度调制信号的解调,包括信号调理电路、模数转换电路、数据寄存器、接口电路、计算机、控制逻辑电路、时钟电路和输入输出电路;
被测调制信号经信号调理电路进行滤波、放大后进入模数转换电路,模数转换电路对进入的信号进行采样和模数转换,将得到的数字化波形数据序列送入数据寄存器进行寄存;接口电路将数据寄存器中的波形数据序列依次读入计算机,计算机对获取的数据序列进行存储、处理、数字化解调运算,将运算结果经过输入输出电路以虚拟仪器的方式提供给操作者并实现人机交互;
计算机通过控制逻辑电路为模数转换电路、数据寄存器和接口电路提供统一的时序和逻辑控制信号,时钟电路为逻辑控制电路和计算机提供统一的时钟信号;
上述接口电路为并行可编程输入接口,采用FPGA编程实现或者采用可编程I/O接口芯片实现。
有益效果
采用借助于通用数据采集技术和波形测量技术,利用虚拟仪器方式,用于进行滑动模型方式的调制信号解调,仅仅使用不足一个载波周期的部分波形获得调制波形参数值,例如在仅有十分之一个波形周期下即可以实现正弦载波模型参数值测量辨识,且能同时完成FM信号解调和AM信号解调,因而具有更高的模型解调分辨力和更高的解调效率,本发明时间分辨力高、鲁棒性好、算法绝对收敛、可以方便进行溯源校准。本发明可用于制作调制度分析仪、测量接收机、AM、FM信号的数字化解调分析测量系统等,并可以用于AM、FM信号参数的计量校准。
附图说明
图1为本发明中数字化解调装置的结构示意图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明。
实施例
一种用于频率调制和幅度调制信号的数字化解调方法,具体步骤为:
1)对待测频率调制信号或者幅度调制信号进行信号调理,利用模数转换电路对调理后的信号进行采样,将采样得到的波形数据序列yi,i=1,2,...,n送入计算机,其中i表示同步采样序列中采样点序号;
模数转换电路采集波形的序列长度n以及采样速率v根据被测调制信号的载波频率设定,设定原则为保证每个载波波形周期内要多于20个采样点;
2)首先,在波形采集序列前缘截取少于一个载波周期的包含m1个采样点的波形段数据序列yi,i=1,2,...,m1;假设截取的波形数据序列近似为正弦波,对该正弦波信号进行数学描述,其波形为
Figure BSA00000717805400051
其中,A1为正弦波幅值,f1为正弦波频率值,
Figure BSA00000717805400052
为正弦波初始相位值;
3)用计算机对步骤2)截取的波形数据序列进行正弦拟合,获得拟合正弦波的幅度A1、频率f1,具体过程为:
3.1在波形数据序列前缘截取少于一个载波周期的包含m1个采样点的波形段数据序列yi(i=1,2,...,m1),则
yi=y(ti)=y((i-1)·Δτ),i=1,2,...,m1
其中,采样时间间隔Δτ=1/v,v为采样速率;
3.2用计算机对截取的数据序列yi(i=1,2,...,m1)进行正弦波形四参数拟合,获得拟合信号:
Figure BSA00000717805400061
其中,
Figure BSA00000717805400062
为正弦波幅度的拟合值,
Figure BSA00000717805400063
为正弦波角频率的拟合值,
Figure BSA00000717805400064
为正弦波初始相位的拟合值,
Figure BSA00000717805400065
为正弦波直流分量值的拟合值;
频率的拟合值
Figure BSA00000717805400066
f ^ 1 = ω ^ 1 · v 2 π
所述正弦波形四参数拟合具体过程为:
(1)设定初始参数,具体包括:待测的正弦信号的频率预估值为f0,角频率预估值ω0=2πf0/v,波形采集序列所含信号不足一个周期且个数为ρ,波形采集序列占用时间长度为τ,则f0≤1/τ,选取另一个足够小的正数因子q=1×10-5,使得f0>q/τ,此时f0∈[q/τ,2/τ];
(2)设定拟合迭代停止条件判据为一个接近于0的足够小的正数he
(3)从已知时刻t1,t2,...,tn的正弦波采集样本y1,y2,...yn,使用计点法获得信号波形占用时间长度为τ=(n-1)/v,确定目标频率f0的存在区间[q/τ,2/τ];
(4)确定迭代左边界频率fL=q/τ和迭代左边界角频率ωL=2πfL/v以及迭代右边界频率fR=2/τ和迭代右边界角频率ωR=2πfR/v,令中值角频率ωM=(ωRL)/2;
(5)在迭代左边界角频率ωL、迭代右边界角频率ωR和中值角频率ωM上分别利用频率已知的三参数正弦波形拟合公式计算各自的拟合残差ρ(ωL)、ρ(ωR)和ρ(ωM);
(6)若ρ(ωL)<η·ρ(ωM),则令迭代右边界角频率ωR=ωM,迭代左边界角频率ωL不变,重复执行步骤(5)~步骤(6),其中η为判据因子且η取值范围为1~1.5;
(7)若ρ(ωL)≥η·ρ(ωM),则必有ωR<2ω0,确定左边界角频率为ωL,右边界角频率ωR,按照优选法原则,选取两个中值角频率分别为
ωM=ωL+0.618×(ωRL)
ωT=ωR-0.618×(ωRL);
(8)在左边界角频率ωL上执行频率已知的三参数正弦曲线拟合获得拟合幅度AL、拟合相位
Figure BSA00000717805400068
拟合直流分量DL、拟合残差ρL,在右边界角频率ωR上执行频率已知的三参数正弦曲线拟合获得拟合幅度AR、拟合相位
Figure BSA00000717805400069
拟合直流分量DR、拟合残差ρR,在中值角频率ωM上执行频率已知的三参数正弦曲线拟合获得拟合幅度AM、拟合相位
Figure BSA00000717805400071
拟合直流分量DM、拟合残差ρM,在中值角频率ωT上执行频率已知的三参数正弦曲线拟合获得拟合幅度AT、拟合相位
Figure BSA00000717805400072
拟合直流分量DT、拟合残差ρT
(9)若ρM<ρT,则最小拟合残差ρ=ρM,有ω0∈[ωT,ωR],参量更新为ωL=ωT,ωT=ωM,ωM=ωL+0.618×(ωRL);若ρM≥ρT,则最小拟合残差ρ=ρT,有ω0∈[ωL,ωM],参量更新为ωR=ωM,ωM=ωT,ωT=ωR-0.618×(ωRL);
(10)判定若|(ρMT)/ρT|≥he,则重复步骤(5)~步骤(10);若|(ρMT)/ρT|<he,则停止迭代,此时又分为两种情况:若最小拟合残差ρ=ρT,获得四参数拟合正弦曲线参数为
Figure BSA00000717805400073
拟合过程结束;若最小拟合残差ρ=ρM,获得四参数拟合正弦曲线参数为
Figure BSA00000717805400074
Figure BSA00000717805400075
拟合过程结束;
4)存储输出正弦模型幅度A1的估计值
Figure BSA00000717805400076
正弦模型频率f1的估计值
5)以小于
Figure BSA00000717805400078
的周期长度(例如
Figure BSA00000717805400079
对应的序列时间长度截取下一段拟合序列长度m2,该拟合序列的中心位置从m1/2处向后移动一个采样点,在新截取的拟合序列上执行步骤3)、4)所述过程,获得输出正弦模型的幅度
Figure BSA000007178054000710
正弦模型的频率
Figure BSA000007178054000711
6)按照步骤5)的方式重复向后推动截取拟合过程,直至数据序列达到终点,获得幅度解调波形序列:
Figure BSA000007178054000712
为幅度调制信号的解调波形输出,获得频率解调波形序列:
Figure BSA000007178054000713
为频率调制信号的解调波形输出。
一种用于频率调制和幅度调制信号的数字化解调装置,用于频率调制信号或者幅度调制信号的解调,其结构如图1所示,包括信号调理电路、模数转换电路、数据寄存器、接口电路、计算机、控制逻辑电路、时钟电路和输入输出电路;
被测调制信号经信号调理电路进行滤波、放大后进入模数转换电路,模数转换电路对进入的信号进行采样和模数转换,将得到的数字化波形数据序列送入数据寄存器进行寄存;接口电路将数据寄存器中的波形数据序列依次读入计算机,计算机对获取的数据序列进行存储、处理、数字化解调运算,将运算结果经过输入输出电路以虚拟仪器的方式提供给操作者并实现人机交互;
上述计算机进行数字化解调运算过程中,获得的正弦模型幅度随时间变化的曲线是幅度调制信号波形解调结果曲线,获得的正弦模型频率随时间变化的曲线是频率调制信号波形解调结果曲线;
计算机通过控制逻辑电路为模数转换电路、数据寄存器和接口电路提供统一的时序和逻辑控制信号,时钟电路为逻辑控制电路和计算机提供统一的时钟信号;
上述接口电路为并行可编程输入接口,采用可编程I/O接口芯片8255实现。
目前,射频信号发生器和任意波发生器均具有频率调制信号和幅度调制信号输出功能,令其分别输出相应的频率调制或幅度调制信号,使用通用数字存储示波器执行波形测量,使用本发明所述方法与过程,可以获得其频率调制或幅度调制信号的解调波形及参数,可以获得优于0.1%的参数准确度以及优于0.01%解调参数分辨力。
以上所述为本发明的较佳实施例而已,本发明不应该局限于该实施例和附图所公开的内容。凡是不脱离本发明所公开的精神下完成的等效或修改,都落入本发明保护的范围。

Claims (6)

1.一种用于频率调制和幅度调制信号的数字化解调方法,具体步骤为:
1)对待测频率调制信号或者幅度调制信号进行信号调理,利用模数转换电路对调理后的信号进行采样,将采样得到的波形数据序列yi,i=1,2,...,n送入计算机,其中i表示同步采样序列中采样点序号;
2)首先,在波形采集序列前缘截取少于一个载波周期的包含m1个采样点的波形段数据序列yi,i=1,2,...,m1;假设截取的波形数据序列近似为正弦波,对该正弦波信号进行数学描述,其波形为
Figure FSA00000717805300011
其中,A1为正弦波幅值,f1为正弦波频率值,
Figure FSA00000717805300012
为正弦波初始相位值;
3)用计算机对步骤2)截取的波形数据序列进行正弦拟合,获得拟合正弦波的幅度A1、频率f1,具体过程为:
3.1 在波形数据序列前缘截取少于一个载波周期的包含m1个采样点的波形段数据序列yi,i=1,2,...,m1,则
yi=y(ti)=y((i-1)·Δτ),i=1,2,...,m1
其中,采样时间间隔ΔτΩ1/v,v为采样速率;
3.2 用计算机对截取的数据序列yi(i=1,2,...,m1)进行正弦波形四参数拟合,获得拟合信号:
Figure FSA00000717805300013
其中,
Figure FSA00000717805300014
为正弦波幅度的拟合值,
Figure FSA00000717805300015
为正弦波角频率的拟合值,
Figure FSA00000717805300016
为正弦波初始相位的拟合值,
Figure FSA00000717805300017
为正弦波直流分量值的拟合值;
频率拟合值
Figure FSA00000717805300018
f ^ 1 = ω ^ 1 · v 2 π ;
4)存储输出正弦模型幅度A1的估计值正弦模型频率f1的估计值
Figure FSA000007178053000111
5)以小于
Figure FSA000007178053000112
的周期长度对应的序列时间长度截取下一段拟合序列长度m2,该拟合序列的中心位置从m1/2处向后移动一个采样点,在新截取的拟合序列上执行步骤3)、4)所述过程,获得输出正弦模型的幅度
Figure FSA000007178053000113
正弦模型的频率
Figure FSA000007178053000114
6)按照步骤5)的方式重复向后推动截取拟合过程,直至数据序列达到终点,获得幅度解调波形序列:
Figure FSA00000717805300021
为幅度调制信号的解调波形输出,获得频率解调波形序列:
Figure FSA00000717805300022
为频率调制信号的解调波形输出。
2.根据权利要求1所述的一种用于频率调制和幅度调制信号的数字化解调方法,其特征在于,所述步骤2.2中进行正弦波形四参数拟合的过程为:
(1)设定初始参数,具体包括:待测的正弦信号的频率预估值为f0,角频率预估值ω0=2πf0/v,波形采集序列所含信号不足一个周期且个数为p,波形采集序列占用时间长度为τ,则f0≤1/τ,选取另一个足够小的正数因子q=1×10-5,使得f0>q/τ,此时f0∈[q/τ,2/τ];
(2)设定拟合迭代停止条件判据为一个接近于0的足够小的正数he
(3)从已知时刻t1,t2,...,tn的正弦波采集样本y1,y2,...yn,使用计点法获得信号波形占用时间长度为τ=(n-1)/v,确定目标频率f0的存在区间[q/τ,2/τ];
(4)确定迭代左边界频率fL=q/τ和迭代左边界角频率ωL=2πfL/v以及迭代右边界频率fR=2/τ和迭代右边界角频率ωR=2πfR/v,令中值角频率ωM=(ωRL)/2;
(5)在迭代左边界角频率ωL、迭代右边界角频率ωR和中值角频率ωM上分别利用频率已知的三参数正弦波形拟合公式计算各自的拟合残差ρ(ωL)、ρ(ωR)和ρ(ωM);
(6)若ρ(ωL)<η·ρ(ωM),则令迭代右边界角频率ωR=ωM,迭代左边界角频率ωL不变,重复执行步骤(5)~步骤(6),其中η为判据因子且η取值范围为1~1.5;
(7)若ρ(ωL)≥η·ρ(ωM),则必有ωR<2ω0,确定左边界角频率为ωL,右边界角频率ωR,按照优选法原则,选取两个中值角频率分别为
ωM=ωL+0.618×(ωRL)
ωT=ωR-0.618×(ωRL);
(8)在左边界角频率ωL上执行频率已知的三参数正弦曲线拟合获得拟合幅度AL、拟合相位
Figure FSA00000717805300023
拟合直流分量DL、拟合残差ρL,在右边界角频率ωR上执行频率已知的三参数正弦曲线拟合获得拟合幅度AR、拟合相位
Figure FSA00000717805300024
拟合直流分量DR、拟合残差ρR,在中值角频率ωM上执行频率已知的三参数正弦曲线拟合获得拟合幅度AM、拟合相位
Figure FSA00000717805300031
拟合直流分量DM、拟合残差ρM,在中值角频率ωT上执行频率已知的三参数正弦曲线拟合获得拟合幅度AT、拟合相位
Figure FSA00000717805300032
拟合直流分量DT、拟合残差ρT
(9)若ρM<ρT,则最小拟合残差ρ=ρM,有ω0∈[ωT,ωR],参量更新为ωL=ωT,ωT=ωM,ωM=ωL+0.618×(ωRL);若ρM≥ρT,则最小拟合残差ρ=ρT,有ω0∈[ωL,ωM],参量更新为ωR=ωM,ωM=ωT,ωT=ωR-0.618×(ωRL);
(10)判定若|(ρMT)/ρT|≥he,则重复步骤(5)~步骤(10);若|(ρMT)/ρT|<he,则停止迭代,此时又分为两种情况:若最小拟合残差ρ=ρT,获得四参数拟合正弦曲线参数为
Figure FSA00000717805300033
拟合过程结束;若最小拟合残差ρ=ρM,获得四参数拟合正弦曲线参数为
Figure FSA00000717805300034
Figure FSA00000717805300035
拟合过程结束。
3.根据权利要求1所述的一种用于频率调制和幅度调制信号的数字化解调方法,其特征在于,所述步骤2)中模数转换电路采集波形的序列长度n以及采样速率v根据被测调制信号的载波频率设定,设定原则为保证每个载波波形周期内要多于20个采样点。
4.根据权利要求2所述的一种用于频率调制和幅度调制信号的数字化解调方法,其特征在于,所述步骤(2)中he的优选取值范围为1×10-40~1×10-20
5.根据权利要求2所述的一种用于频率调制和幅度调制信号的数字化解调方法,其特征在于,所述步骤(6)中η的优选取值范围为1~1.5。
6.一种用于频率调制和幅度调制信号的数字化解调装置,其特征在于,包括信号调理电路、模数转换电路、数据寄存器、接口电路、计算机、控制逻辑电路、时钟电路和输入输出电路;
被测调制信号经信号调理电路进行滤波、放大后进入模数转换电路,模数转换电路对进入的信号进行采样和模数转换,将得到的数字化波形数据序列送入数据寄存器进行寄存;接口电路将数据寄存器中的波形数据序列依次读入计算机,计算机对获取的数据序列进行存储、处理、数字化解调运算,将运算结果经过输入输出电路以虚拟仪器的方式提供给操作者并实现人机交互;
计算机通过控制逻辑电路为模数转换电路、数据寄存器和接口电路提供统一的时序和逻辑控制信号,时钟电路为逻辑控制电路和计算机提供统一的时钟信号;
上述接口电路为并行可编程输入接口,采用FPGA编程实现或者采用可编程I/O接口芯片实现。
CN2012101499272A 2012-05-14 2012-05-14 用于频率调制和幅度调制信号的数字化解调方法及装置 Pending CN102664588A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012101499272A CN102664588A (zh) 2012-05-14 2012-05-14 用于频率调制和幅度调制信号的数字化解调方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012101499272A CN102664588A (zh) 2012-05-14 2012-05-14 用于频率调制和幅度调制信号的数字化解调方法及装置

Publications (1)

Publication Number Publication Date
CN102664588A true CN102664588A (zh) 2012-09-12

Family

ID=46774023

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012101499272A Pending CN102664588A (zh) 2012-05-14 2012-05-14 用于频率调制和幅度调制信号的数字化解调方法及装置

Country Status (1)

Country Link
CN (1) CN102664588A (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105277801A (zh) * 2015-08-10 2016-01-27 工业和信息化部电信研究院 一种测量调制脉冲参数的方法
CN106124032A (zh) * 2016-08-25 2016-11-16 中国航空工业集团公司北京长城计量测试技术研究所 一种光频调制器调制延迟的数字化测量方法及装置
CN106199608A (zh) * 2016-09-22 2016-12-07 湖北工业大学 一种基于移动正弦拟合的快速解包络算法
CN108333420A (zh) * 2018-04-02 2018-07-27 三峡大学 一种低频正弦信号峰值检测装置和峰值检测方法
CN110138398A (zh) * 2019-04-23 2019-08-16 深圳市华讯方舟微电子科技有限公司 超外差式接收装置调制方法和系统
CN110906875A (zh) * 2019-11-26 2020-03-24 湖北工业大学 一种孔径测量的视觉处理方法
CN111664111A (zh) * 2020-07-09 2020-09-15 苏州格里德医学传感技术有限公司 一种基于电极感应电场变化的电扇摇摆周期测量装置
CN111736014A (zh) * 2020-07-20 2020-10-02 中国空气动力研究与发展中心低速空气动力研究所 一种频率信号采集和数据处理方法及装置
CN111737883A (zh) * 2020-07-30 2020-10-02 哈尔滨工业大学 一种具有输出时滞的非线性双率电路系统鲁棒辨识方法
CN113341219A (zh) * 2021-05-19 2021-09-03 北京航空航天大学 恒频交流供电系统频率调制幅度测量方法及装置
CN113392610A (zh) * 2021-08-18 2021-09-14 苏州贝克微电子有限公司 一种用于大规模集成电路瞬态分析的仿真优化方法及装置
CN114253343A (zh) * 2021-12-23 2022-03-29 中国航空工业集团公司西安航空计算技术研究所 一种任意调幅组件

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005030978A (ja) * 2003-07-09 2005-02-03 Sony Corp 位相差測定装置、位相差測定方法および試験装置
CN102253284A (zh) * 2011-04-18 2011-11-23 中国航空工业集团公司北京长城计量测试技术研究所 一种用于多路正弦波信号的相位差同步测量装置及方法
CN102253254A (zh) * 2011-04-18 2011-11-23 中国航空工业集团公司北京长城计量测试技术研究所 一种多路正弦波相位标准信号的产生装置及方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005030978A (ja) * 2003-07-09 2005-02-03 Sony Corp 位相差測定装置、位相差測定方法および試験装置
CN102253284A (zh) * 2011-04-18 2011-11-23 中国航空工业集团公司北京长城计量测试技术研究所 一种用于多路正弦波信号的相位差同步测量装置及方法
CN102253254A (zh) * 2011-04-18 2011-11-23 中国航空工业集团公司北京长城计量测试技术研究所 一种多路正弦波相位标准信号的产生装置及方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
梁志国 等: "四参数正弦曲线拟合的一种收敛算法", 《仪器仪表学报》, vol. 27, no. 11, 30 November 2006 (2006-11-30), pages 1513 - 1519 *
梁志国 等: "四参数正弦波曲线拟合的快速算法", 《理论与实践》, vol. 26, no. 1, 31 January 2006 (2006-01-31), pages 4 - 7 *
梁志国 等: "正弦波形参数拟合方法述评", 《测量技术学报》, vol. 24, no. 1, 31 January 2010 (2010-01-31), pages 1 - 7 *
梁志国 等: "残周期正弦波形的四参数拟合", 《计量学报》, vol. 30, no. 3, 31 May 2009 (2009-05-31), pages 245 - 249 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105277801A (zh) * 2015-08-10 2016-01-27 工业和信息化部电信研究院 一种测量调制脉冲参数的方法
CN105277801B (zh) * 2015-08-10 2018-08-07 工业和信息化部电信研究院 一种测量调制脉冲参数的方法
CN106124032A (zh) * 2016-08-25 2016-11-16 中国航空工业集团公司北京长城计量测试技术研究所 一种光频调制器调制延迟的数字化测量方法及装置
CN106124032B (zh) * 2016-08-25 2021-11-30 中国航空工业集团公司北京长城计量测试技术研究所 一种光频调制器调制延迟的数字化测量方法及装置
CN106199608A (zh) * 2016-09-22 2016-12-07 湖北工业大学 一种基于移动正弦拟合的快速解包络算法
CN106199608B (zh) * 2016-09-22 2018-05-11 湖北工业大学 一种基于移动正弦拟合的快速解包络算法
CN108333420A (zh) * 2018-04-02 2018-07-27 三峡大学 一种低频正弦信号峰值检测装置和峰值检测方法
CN110138398A (zh) * 2019-04-23 2019-08-16 深圳市华讯方舟微电子科技有限公司 超外差式接收装置调制方法和系统
CN110906875B (zh) * 2019-11-26 2021-07-06 湖北工业大学 一种孔径测量的视觉处理方法
CN110906875A (zh) * 2019-11-26 2020-03-24 湖北工业大学 一种孔径测量的视觉处理方法
CN111664111A (zh) * 2020-07-09 2020-09-15 苏州格里德医学传感技术有限公司 一种基于电极感应电场变化的电扇摇摆周期测量装置
CN111736014A (zh) * 2020-07-20 2020-10-02 中国空气动力研究与发展中心低速空气动力研究所 一种频率信号采集和数据处理方法及装置
CN111736014B (zh) * 2020-07-20 2020-11-24 中国空气动力研究与发展中心低速空气动力研究所 一种频率信号采集和数据处理方法及装置
CN111737883A (zh) * 2020-07-30 2020-10-02 哈尔滨工业大学 一种具有输出时滞的非线性双率电路系统鲁棒辨识方法
CN113341219A (zh) * 2021-05-19 2021-09-03 北京航空航天大学 恒频交流供电系统频率调制幅度测量方法及装置
CN113392610A (zh) * 2021-08-18 2021-09-14 苏州贝克微电子有限公司 一种用于大规模集成电路瞬态分析的仿真优化方法及装置
CN114253343A (zh) * 2021-12-23 2022-03-29 中国航空工业集团公司西安航空计算技术研究所 一种任意调幅组件

Similar Documents

Publication Publication Date Title
CN102664588A (zh) 用于频率调制和幅度调制信号的数字化解调方法及装置
CN101806832B (zh) 一种低频率信号的频率测量方法
CN101825660B (zh) 欠采样下的正弦信号频率的高效测量方法及实施装置
CN100489555C (zh) 通用电力互感器校验装置
CN102879639A (zh) 一种电力系统中频率的实时测量方法
CN102645585A (zh) 一种超低频正弦信号相位差的同步快速测量方法及装置
CN102116798A (zh) 一种电网频率测量方法及装置
CN102809687B (zh) 一种交流电频率的数字化测量方法
CN103983849B (zh) 一种实时高精度的电力谐波分析方法
CN104181577A (zh) 一种基于全数字化技术的束流位置和相位测量系统及方法
CN104597320B (zh) 一种适用于多个频率交流信号计算的方法
CN106841778A (zh) 基于pmu实现的次同步和超同步谐波参数的处理方法
CN203054093U (zh) 利用同步采样在混合信号中提取直流信号的装置
CN107505497B (zh) 一种无源磁电转速传感器信号峰峰值时域测量方法
CN102221639A (zh) 正负序电流实时检测的方法
CN108333426A (zh) 基于傅氏算法的电力系统频率测量方法
CN104022761B (zh) 一种脉冲调制信号的测试方法
CN103969508B (zh) 一种实时高精密的电力谐波分析方法及装置
CN104569581B (zh) 一种电网频率测量的多水平集单周期估计方法
CN103575979B (zh) 一种数字化测量交流电频率的方法
CN102401856B (zh) 一种基于模糊聚类方法的工频电压信号实时频率测量方法
CN105548711B (zh) 一种多频信息滤波递推解调方法
CN103412209A (zh) 一种调谐滤波器的失谐度检测方法
CN103575981A (zh) 一种交流电频率的精确测量方法
CN105866576A (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
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120912