CN109547030B - 一种基于脉宽调制的随机解调采样方法 - Google Patents

一种基于脉宽调制的随机解调采样方法 Download PDF

Info

Publication number
CN109547030B
CN109547030B CN201811439185.0A CN201811439185A CN109547030B CN 109547030 B CN109547030 B CN 109547030B CN 201811439185 A CN201811439185 A CN 201811439185A CN 109547030 B CN109547030 B CN 109547030B
Authority
CN
China
Prior art keywords
counter
signal
input signal
comparator
output
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.)
Active
Application number
CN201811439185.0A
Other languages
English (en)
Other versions
CN109547030A (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.)
Fuzhou University
Original Assignee
Fuzhou University
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 Fuzhou University filed Critical Fuzhou University
Priority to CN201811439185.0A priority Critical patent/CN109547030B/zh
Publication of CN109547030A publication Critical patent/CN109547030A/zh
Application granted granted Critical
Publication of CN109547030B publication Critical patent/CN109547030B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M7/00Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
    • H03M7/30Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
    • H03M7/3059Digital compression and data reduction techniques where the original information is represented by a subset or similar information, e.g. lossy compression
    • H03M7/3062Compressive sampling or sensing
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M1/00Analogue/digital conversion; Digital/analogue conversion
    • H03M1/12Analogue/digital converters
    • H03M1/124Sampling or signal conditioning arrangements specially adapted for A/D converters
    • H03M1/1245Details of sampling arrangements or methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明涉及一种基于脉宽调制的随机解调采样方法。采用比较器与计数器,代替传统的ADC,完成观测矩阵所要实现的功能;PN序列控制的锯齿状调制信号,通过比较器调制输入信号,而产生PWM波,其作为计数器的计数信号,控制计数器计数;通过这些步骤,实现了将输入信号的幅度信息转换成了时间信息。与传统的压缩采样系统对比,本系统采用了数字编码的方式,其稳定性大大提升;由于本发明系统稳定性高,功耗低,体积小,在低成本的采集设备,医疗器械及便携式终端等有良好的应用前景。

Description

一种基于脉宽调制的随机解调采样方法
技术领域
本发明涉及一种基于脉宽调制的随机解调采样方法。
背景技术
随着电子信息技术的飞速发展,数字化系统已成为现代电子系统的主流。模数转换器(ADC)承担着模拟信号数字化的主要工作,是现代系统不可缺少的一部分。近年来,随着电子信息技术的广泛发展,各种终端日益往小型化、便携式方向发展。以奈奎斯特采样定理为指导的信号采集方法导致了模字转换器的体积大,功耗高等问题而面临着巨大的挑战。我们知道奈奎斯特采样定理是不失真恢复原始信号的充分条件,而非必要条件。近几年来,数字信号处理领域的专家和学者们提出了一种全新的信息获取方法,称为压缩采样。由压缩采样理论可知:只要信号本身或者在其某一变换域上是稀疏的,可以利用与变换矩阵非相干的观测矩阵将变换系数线性投影为低维观测向量,且投影保持了重建信号所需的信息,通过求解稀疏最优化问题就能够从低维观测向量精确地或高概率精确地重建原始高维信号。
对于稀疏信号,压缩采样理论能够同时完成信号的采样与压缩。因此省掉传统奈奎斯特采样理论中先采样再压缩的低效过程,大大提高了对于稀疏信号的信号采样效率,突破了奈奎斯特采样定理的限制,使ADC的体积和功耗大大降低。
本发明采用了数字编码的方式进行压缩采样,整套采样系统使用了比较器与计数器,代替了ADC模块,避开了使用ADC的弊端,大大降低了整体的功耗与体积。
发明内容
本发明的目的在于提供一种基于脉宽调制的随机解调采样方法,该方法的应用系统稳定性高,功耗低,体积小,在低成本的采集设备,医疗器械及便携式终端等有良好的应用前景。
为实现上述目的,本发明的技术方案是:一种基于脉宽调制的随机解调采样方法,提供有比较器、锯齿波信号发生器、计数器;通过具有信号奈奎斯特频率的伪随机序列控制锯齿波发生器,使锯齿波信号发生器输出的调制信号波形为单个锯齿波与高电平交替出现的波形,将调制信号输入到比较器中,与输入比较器的输入信号进行比较,使得比较器的输出为高电平持续时间不断变化的PWM波信号,并将PWM波信号输入到计数器当中,作为控制计数器计数时长的控制信号,从而实现将输入信号的幅度信息转换成时间信息。
在本发明一实施例中,所述将输入信号的幅度信息转换成时间信息的方式具体如下:
输入信号为Vc(t),三角波的周期为T,对于一个三角波的周期而言,三角波的表达式为:
Figure BDA0001884283110000021
假设这个周期的占空比为δc,在t=δcT时,Vr(t)=Vc(t),即
Figure BDA0001884283110000022
化简(2)可得:
VccT)=Vcmaxδc (3)
由式(3)可知,输入信号的幅值为关于占空比δc的线性函数;即,输入信号的幅值信息反映在占空比的大小中;
假设三角波的频率为f1,计数器的频率为f2,满足f2>>f1;对于一个三角波周期而言,计数器的最高计数次数为
Figure BDA0001884283110000023
次;计数器的输出为C,即计数器计数C次,则计数时间为
Figure BDA0001884283110000024
则可以得到:
Figure BDA0001884283110000025
将式(4)带入式(3),可得:
Figure BDA0001884283110000026
由式(5)可知,输入信号的幅值是关于计数时间t′的线性函数;即,输入信号的幅值信息反映在计数器的计数时长之中,实现了将输入信号的幅值信息转换成了时间信息;
当计数器完成计数工作后,输出一串序列值y[m],m为压缩采样后序列的序号,序列内的每一个值分别对应计数器的每一次输出值C;而后利用子空间基追踪算法进行信号重构。
相较于现有技术,本发明具有以下有益效果:
1、本发明方法利用比较器与计数器设计了压缩采样系统,完成了观测矩阵所要实现的功能;
2、本发明采用数字编码的形式对输入信号进行压缩采样,对比于传统的基于随机解调的压缩采样系统,本发明的稳定性大大地提升;
本发明适用于超低功耗或者超高速,并且不需要ADC模块实现模数转换的应用场景;在低成本的采集设备,医疗器械及便携式终端等有良好的应用前景。
附图说明
图1为PWM波发生器结构。
图2为本发明电路框图。
图3为本发明实际应用电路框图。
图4为本发明实际应用电路各输出波形图。
图5为心电信号。
图6为重构信号。
具体实施方式
下面结合附图,对本发明的技术方案进行具体说明。
本发明提供了一种基于脉宽调制的随机解调采样方法,提供有比较器、锯齿波信号发生器、计数器;通过具有信号奈奎斯特频率的伪随机序列控制锯齿波发生器,使锯齿波信号发生器输出的调制信号波形为单个锯齿波与高电平交替出现的波形,将调制信号输入到比较器中,与输入比较器的输入信号进行比较,使得比较器的输出为高电平持续时间不断变化的PWM波信号,并将PWM波信号输入到计数器当中,作为控制计数器计数时长的控制信号,从而实现将输入信号的幅度信息转换成时间信息。
所述将输入信号的幅度信息转换成时间信息的方式具体如下:
输入信号为Vc(t),三角波的周期为T,对于一个三角波的周期而言,三角波的表达式为:
Figure BDA0001884283110000031
假设这个周期的占空比为δc,在t=δcT时,Vr(t)=Vc(t),即
Figure BDA0001884283110000034
化简(2)可得:
VccT)=Vcmaxδc (3)
由式(3)可知,输入信号的幅值为关于占空比δc的线性函数;即,输入信号的幅值信息反映在占空比的大小中;
假设三角波的频率为f1,计数器的频率为f2,满足f2>>f1;对于一个三角波周期而言,计数器的最高计数次数为
Figure BDA0001884283110000032
次;计数器的输出为C,即计数器计数C次,则计数时间为
Figure BDA0001884283110000033
则可以得到:
Figure BDA0001884283110000041
将式(4)带入式(3),可得:
Figure BDA0001884283110000042
由式(5)可知,输入信号的幅值是关于计数时间t′的线性函数;即,输入信号的幅值信息反映在计数器的计数时长之中,实现了将输入信号的幅值信息转换成了时间信息;
当计数器完成计数工作后,输出一串序列值y[m],m为压缩采样后序列的序号,序列内的每一个值分别对应计数器的每一次输出值C;而后利用子空间基追踪算法进行信号重构。
以下为本发明的具体实现过程。
本发明一种基于脉宽调制的随机解调采样方法,采用了数字编码的方式进行压缩采样,设计了基于脉宽调制的随机解调采样系统,整套采样系统使用了比较器与计数器,代替了ADC模块,避开了使用ADC的弊端,大大降低了整体的功耗与体积。
具体的,本发明方法是基于脉冲宽度调制提供了一种压缩采样系统。脉冲宽度调制(PWM)是一种对模拟信号电平进行数字编码的方法,其结构如图1所示,系统由一个比较器和一个锯齿波信号发生器组成。当输入信号的电平大于锯齿波时,比较器输出高电平,否则输出低电平,因此比较器的输出波形为占空比持续变化的方波,其占空比的大小反映了输入信号的幅度信息。
本发明利用了具有信号奈奎斯特频率的伪随机序列控制锯齿波发生器,使调制信号的波形为单个锯齿波与高电平交替出现的波形。调制信号输入到比较器中,与输入信号进行比较,比较器的输出为高电平持续时间不断变化的PWM波,其输入到计数器当中,作为控制计数器计数时长的控制信号,其结构如图2所示。经过这一步实现了将输入信号的幅度信息转换成了时间信息,具体转换过程如下:
输入信号为Vc(t),三角波的周期为T,对于一个三角波的周期而言,三角波的表达式为:
Figure BDA0001884283110000043
假设这个周期的占空比为δc,在t=δcT时,Vr(t)=Vc(t),即
Figure BDA0001884283110000044
化简(2)可得:
VccT)=Vcmaxδc (3)
由式(3)可知,输入信号的幅值为关于占空比δc的线性函数;即,输入信号的幅值信息反映在占空比的大小中;
假设三角波的频率为f1,计数器的频率为f2,满足f2>>f1;对于一个三角波周期而言,计数器的最高计数次数为
Figure BDA0001884283110000051
次;计数器的输出为C,即计数器计数C次,则计数时间为
Figure BDA0001884283110000052
则可以得到:
Figure BDA0001884283110000053
将式(4)带入式(3),可得:
Figure BDA0001884283110000054
由式(5)可知,输入信号的幅值是关于计数时间t′的线性函数;即,输入信号的幅值信息反映在计数器的计数时长之中,实现了将输入信号的幅值信息转换成了时间信息;
当计数器完成计数工作后,输出一串序列值y[m],m为压缩采样后序列的序号,序列内的每一个值分别对应计数器的每一次输出值C;而后利用子空间基追踪算法进行信号重构。
实施例:
由于PN序列是由1和-1组成的序列,本发明采用的是数字编码的形式,不容易表示负数,因此在进行采样的时候,需要分成两路对输入信号进行采样。其实际应用电路框图如图3所示。输入信号以心电信号为例,电路中各个器件的输出波形如图4所示,调制信号发生器输出两种调制信号W1与W2,分别对应PN序列的1与-1。由于调制信号不同,两个比较器输出的PWM波也是不同的。两个计数器分别输出y1[m]与y2[m],并将其通过串口输入到上位机,在matlab平台上使用子空间基追踪算法完成信号重构。
在实物实现过程中,输入信号x[n]采用从MIT-BIH数据库下载的心电信号,其频率为360Hz,如图5所示。因为输入信号的频率为360Hz,所以PN序列的频率也为360Hz。最后在PC机上的matlab平台,使用子空间基追踪算法,对实物系统输出的y[m]进行重构,重构的结果如图6所示,两者的结构相似度为0.9667。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。

Claims (1)

1.一种基于脉宽调制的随机解调采样方法,其特征在于,提供有比较器、锯齿波信号发生器、计数器;通过具有信号奈奎斯特频率的伪随机序列控制锯齿波发生器,使锯齿波信号发生器输出的调制信号波形为单个锯齿波与高电平交替出现的波形,将调制信号输入到比较器中,与输入比较器的输入信号进行比较,使得比较器的输出为高电平持续时间不断变化的PWM波信号,并将PWM波信号输入到计数器当中,作为控制计数器计数时长的控制信号,从而实现将输入信号的幅度信息转换成时间信息;所述将输入信号的幅度信息转换成时间信息的方式具体如下:
输入信号为Vc(t),三角波的周期为T,对于一个三角波的周期而言,三角波的表达式为:
Figure FDA0003165955050000011
假设这个周期的占空比为δe,在t=δcT时,Vr(t)=Vc(t),即
Figure FDA0003165955050000012
化简(2)可得:
VccT)=Vcmaxδc (3)
由式(3)可知,输入信号的幅值为关于占空比δe的线性函数;即,输入信号的幅值信息反映在占空比的大小中;
假设三角波的频率为f1,计数器的频率为f2,满足f2>>f1;对于一个三角波周期而言,计数器的最高计数次数为
Figure FDA0003165955050000013
次;计数器的输出为C,即计数器计数C次,则计数时间为
Figure FDA0003165955050000014
则可以得到:
Figure FDA0003165955050000015
将式(4)带入式(3),可得:
Figure FDA0003165955050000016
由式(5)可知,输入信号的幅值是关于计数时间t′的线性函数;即,输入信号的幅值信息反映在计数器的计数时长之中,实现了将输入信号的幅值信息转换成了时间信息;
当计数器完成计数工作后,输出一串序列值y[m],m为压缩采样后序列的序号,序列内的每一个值分别对应计数器的每一次输出值C;而后利用子空间基追踪算法进行信号重构。
CN201811439185.0A 2018-11-29 2018-11-29 一种基于脉宽调制的随机解调采样方法 Active CN109547030B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811439185.0A CN109547030B (zh) 2018-11-29 2018-11-29 一种基于脉宽调制的随机解调采样方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811439185.0A CN109547030B (zh) 2018-11-29 2018-11-29 一种基于脉宽调制的随机解调采样方法

Publications (2)

Publication Number Publication Date
CN109547030A CN109547030A (zh) 2019-03-29
CN109547030B true CN109547030B (zh) 2021-08-31

Family

ID=65851821

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811439185.0A Active CN109547030B (zh) 2018-11-29 2018-11-29 一种基于脉宽调制的随机解调采样方法

Country Status (1)

Country Link
CN (1) CN109547030B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114978183A (zh) * 2022-06-29 2022-08-30 山东大学 一种基于类锯齿波比较型adc的多通道采样装置及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10222344A (ja) * 1997-02-10 1998-08-21 Nec Corp 波形処埋装置
CN103391099A (zh) * 2013-07-15 2013-11-13 中国科学院微电子研究所 适用于一维缓变信号的随机采样器
CN103532552A (zh) * 2013-10-10 2014-01-22 中国科学院微电子研究所 一种随机采样斜率模数转换器
CN103795422A (zh) * 2014-03-07 2014-05-14 中国科学院微电子研究所 一种基于压缩感知的一维信号随机采样方法
CN205374571U (zh) * 2015-12-28 2016-07-06 上海致远绿色能源股份有限公司 一种隔离采样装置
CN107078745A (zh) * 2014-11-14 2017-08-18 索尼公司 信号处理装置、控制方法、摄像元件及电子设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3295107A (en) * 1961-05-17 1966-12-27 Magnavox Co Electronic time compression device
BE758978A (fr) * 1969-11-20 1971-04-30 Western Electric Co Systeme de transmission pour transmettre des signaux numeriquesa niveaux multiples
CN101304247B (zh) * 2008-04-29 2011-10-12 哈尔滨工业大学深圳研究生院 多周期随机数字脉宽调制电路及方法
CN103346988B (zh) * 2013-06-13 2016-03-23 电子科技大学 一种fsk数字解调器
CN104124976B (zh) * 2014-07-24 2017-07-04 福州大学 有限新息率信号结构化亚奈奎斯特率采样方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10222344A (ja) * 1997-02-10 1998-08-21 Nec Corp 波形処埋装置
CN103391099A (zh) * 2013-07-15 2013-11-13 中国科学院微电子研究所 适用于一维缓变信号的随机采样器
CN103532552A (zh) * 2013-10-10 2014-01-22 中国科学院微电子研究所 一种随机采样斜率模数转换器
CN103795422A (zh) * 2014-03-07 2014-05-14 中国科学院微电子研究所 一种基于压缩感知的一维信号随机采样方法
CN107078745A (zh) * 2014-11-14 2017-08-18 索尼公司 信号处理装置、控制方法、摄像元件及电子设备
CN205374571U (zh) * 2015-12-28 2016-07-06 上海致远绿色能源股份有限公司 一种隔离采样装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
PWM的工作原理;电子工程专辑;《https://www.sohu.com/a/120467759_465219》;20161202;第1-6页 *

Also Published As

Publication number Publication date
CN109547030A (zh) 2019-03-29

Similar Documents

Publication Publication Date Title
Schreier An empirical study of high-order single-bit delta-sigma modulators
CN105007079A (zh) 逐次逼近型模数转换器的全差分增量采样方法
CN109547030B (zh) 一种基于脉宽调制的随机解调采样方法
CN103326728B (zh) 一种噪声抑制增强σδ调制器结构
US20230393533A1 (en) Time-to-Digital Converter Circuitry
CN103634005A (zh) 一种模数转换器中量化噪声随机化的方法
CN101562454B (zh) 动态元件匹配的数字/模拟转换系统及积分三角调制装置
CN106155177B (zh) 基于数字处理器实现spwm波形的方法
CN109901382B (zh) 一种数字控制系统的规则采样pwm优化方法
CN103490783B (zh) 一种将模拟信号转换为数字信息的方法
CN103178783B (zh) 求和调制信号发生器
CN102098073A (zh) 一种基于重构滤波的时限带通椭圆球面波函数信号产生方法
CN104883192B (zh) 基于混合压缩采样的非合作超宽带信号的重构方法
Zhou et al. Digital-assisted asynchronous compressive sensing front-end
CN104718704A (zh) 转换器
CN114430274A (zh) 一种信号处理方法及信号处理电路
CN102868404A (zh) 基于余弦算法和格雷编码的模数转换方法
CN110868215B (zh) 一种自适应控制的高精度电流/频率转换电路
CN207427126U (zh) 一种提高模数转换器转换速度的电路
Kościelnik et al. Natural compression and expansion characteristics of asynchronous sigma-delta adc
CN112737394A (zh) 一种专用供电装置控制方法
CN103916143B (zh) 甚低频长波发射机非等高阶梯波信号处理方法和装置
Jia et al. Design and analysis of a hardware-efficient compressed sensing architecture for data compression in power quality data acquisition
KR102610404B1 (ko) 슈퍼 레졸루션을 이용한 샘플링 방법
RU89259U1 (ru) Преобразователь напряжения в частоту импульсов

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