超声波流速测量系统

1引言1.1课题研究背景及意义研究利用超声波测量流体流量已经有数十年的历史了。1928年，法国人o.Rutten研制成功了世界上第一台超声流量计。但为了使超声波流量计有一定的精度，时差法超声波流量计要求对时间的测量有相当高的测量精度，这在当时是很难做到的。1955年，应用声循环法的MA*SON流量计在美国研制成功，用于航空燃料油流量的测量。50年代末期，超声波流量计由理论研究阶段进入工业应用时期。但由于电子线路太复杂而未占有结实的地位。到了80年代中后期，单片机技术的应用使超声流量计向高性能、智能化的方向开展。由于使用了单片机作中央处理单元，系统不仅可以进展复杂的数学运算和数据处理、进一步提高了超声波流量计的测量精度，而月.还能设计出友好的人机界面，使系统具有参数设置、自动检错排错功能以及其他一些辅助功能，大大方便了用户的操作和使用。单片机在超声流量计中的应用，是超声流量计开场真正进入工业测量领域。1.2超声波流量计的现状近10年来，基于高速数字信号的处理技术与微处理器技术的进步，基于新型探头材料与工艺的研究，基于声道配置及流动力学的研究，超声流量测量技术取得了长足的进步，显示了它强劲的技术优势，形成了迅猛开展的势头，其潜在的巨大的生命力是显而易见的。在国外，以美国Controlotron公司和Ploysonics公司为代表的产品较多的采用数字信号处理技术，如“同步调制〞和FFT技术，他们广泛的采用以DSP为核心的数字处理电路，从而能够更快实时的处理超声信号，同时能够实现一些复杂的算法，如Ploysonics公司的DDF3088型是该公司的新一代全数字化便携式多普勒流量计，它采用了数字滤波和数字频谱分析技术，能自动识别多普勒信号与噪声信号，抗干扰能力强，采用了高分别率的液晶显示，可以现场对信号进展多普勒分析。在测量方法方面有的利用改良的时差法，用以消除速度受温度的影响，还有将时差法与多普勒法的组合，如Controlotron公司研制的480型超声流量计，这样使得产品的实用*围更广。但由于国外的产品的价格比拟高〔每台约为5～10万元左右〕，所以限制了在国内的大量使用。目前国内的厂家生产的超声波流量计虽然价格比拟廉价，但总体性能差，主要用于测量比拟容易的大管径管道中流体的测量。他们多采用的方法是时差法，以单片机为系统核心进展信号处理。但由于受单片机速度的限制，所以只能进展一些比拟简单的算法，如通过计数来实现流速的计算。同时由于一般的时差法受温度的影响比拟大，所以在精度上面就不会很高。如**亚楠公司生产的ZLC系列流量计是国内率先采用多脉冲发射和接收宽带技术，并带有微机控制的新型超声波时差式流量计。近几年来由于DSP和超大规模集成电路技术的开展，以及以基于DSP为核心的数字信号处理电路技术日益成熟和广泛使用，使得可以采用以DSP为处理核心的超声波流量计来广泛取代国内的以单片机为核心的超声波流量计，从而可利用数字信号处理的一些技术来改善产品的测量精度，如FFT、数字滤波、插值运算、相关运算，从而来提高系统测量的精度。通过这些改良使国内流量计的性能已到达国外的水平，从而能更好的为国民经济效劳。1.2.1超声波流量计的开展趋势在石油和天然气领域快速普及最近几年超声波流量计市场的一切增长，几乎都是由于这种产品在石油和天然气领域销量的增加而带来的。超声波流量计在这一领域的销量比起以前将近翻了一倍。由于AGA9监护运输标准的推广，天然气的监护运输市场已经根本成形。在2003年到2005年间，监护运输量增加了三倍多，而这一领域占超声波流量计销售的5%。只有液体碳氢化合物监护运输的API和OIML标准在行业内得到推广，超声波流量计才能在这个领域普及。超声波流量计具有高准确度和低总体拥有本钱不管是从技术上还是从经济上看，超声波测量仪器都是流量测量的理想选择。通过多光束和数字信号处理，超声波测量仪可以实现很高的测量准确度。与传统的涡轮式仪表不同，它没有移动的元件，因此几乎不需要维修。而且，它也不会阻挡或者减慢管道中气体或者液体的流动。它能够准确地测量液态石油气产品的宽频，而不需像机械型技术那样得到验证。高灵敏度使其可以检测到管道中的任何泄漏，并可以测量和补充各种会影响监护运输领域中的测量准确度的变量。亚洲和中东市场增长最大亚洲和中东地区超声波流量计的增长将会是最大的。中国和印度将会在根底设施和新工厂上大举投资。中国能源缺乏，为了寻找推动经济快速开展的能源，将会对其石油和天然气根底设施进展改造，并建立连接俄罗斯及其它地区供油商的管道。而由于中东地区在石油和天然气生产中的地位，这一地区将仍然是超声波流量计供给商的沃土。这一地区亦将在数个大型发电及海水淡化厂上进展投资。相比之下，北美市场的增长则显得相对平淡，但由于这一地区在石油和天然气根底设施及工业自动化上的投资，北美市场将仍然具有可观的增长。1.3超声波流量计的市场前景随着石油和天然气工业的强势增长，以及行业对于现场设备技术的逐步承受，超声波流量计的全球市场总额将在今后5年内以9.6%的复合年增长率(CAGR)增长。ARCAdvisoryGroup最近的一份报告指出，2005年该市场总值为2.75亿美元，到2010年将超过4.34亿美元。亚洲和中东市场增长最大，亚洲和中东地区超声波流量计的增长将会是最大的。中国和印度将会在根底设施和新工厂上大举投资。中国能源缺乏，为了寻找推动经济快速开展的能源，将会对其石油和天然气根底设施进展改造，并建立连接俄罗斯及其它地区供油商的管道。而由于中东地区在石油和天然气生产中的地位，这一地区将仍然是超声波流量计供给商的沃土。这一地区亦将在数个大型发电及海水淡化厂上进展投资。相比之下，北美市场的增长则显得相对平淡，但由于这一地区在石油和天然气根底设施及工业自动化上的投资，北美市场将仍然具有可观的增长。1.4主要研究内容本课题的任务是应用时差法测量流量。确立了系统的整体方案，在此根底上进展试验、调试，验证和改良设计方案。主要内容如下：〔1〕确立了系统的整体方案，这主要是在查阅国外流量计产品根本性能以及分析超声波流量计根本原理和工程指标的根底上，确立了系统以DSP为信号处理核心，并采用相关算法来计算流速的改良时差法方案；〔2〕硬件电路设计：CPLD采样控制系统；〔3〕程序设计，应用CCS软件，在DSP最小系统上实现相关算法程序设计，进展实验测试。2超声波流量计设计方案及分析2.1超声流量计的分类超声波技术应用于流量测量主要依据是：当超声波入射到流体后，在流体中传播的超声波就会载有流体流速的信息。超声波流量计对信号的发生、传播及检测有着各种不同的设置方法，从而构成了不同原理的超声流量计，其大致可分为传播速度差法〔包括：时差法、相位差法、频差法〕，多普勒法，相关法等等。2.2超声流量测量的根本原理超声流量测量技术的根本原理是利用超声波在流体中传播时所载流体的流速信息来测量流体流量的。超声波流量计一般由超声波换能器、电子线路及流量显示和累积系统三局部组成。超声波换能器将电能转换为超声波能量，将其发射并穿过被测流体，接收换能器接收到超声波信号，经电子线路放大并转换为代表流量的电信号，供显示和积算，这样就实现了流量的检测显示。2.3传播速度差法传播速度差法是根据超声波在流动的流体中，顺流及逆流传播时的速度之差与被测流体流速之间的关系来求流速或流量的方法。按其所测的物理量的不同，传播速度差法又可分为(直接)时差法、相位差法和频差法三种。时差法就是直接测量超声波顺流和逆流传播的时间差，它适用于大、中口径管道及明渠流量的测量，但时差法流量计受温度的影响较大，流体温度变化对这种系统测量准确度的影响是不能无视的，实验说明，在使用在机玻璃声楔的条件下，流体温度每变化10℃，就会给测量增加1%左右的误差，并且声速的温度系数不是常数;另外，当流体的组成或密度变化时，也将引起声速的变化，从而影响测量精度。其开展方向是提高计时精度和设法降低温度对测量精度的影响。频差法是通过测量一定时间内两组闭路循环系统中的循环频率之差来测得流量。频差法精度高、受温度影响较小，但受环境影响较大，工作不稳定。相位差法通过把时间差转换为相位差，防止了测量微小的时间差，可提高测量的精度，但受温度的影响依然很大。2.4时差法测量原理及主要特点本论文研究的是基于时差法的超声波液体流量计，其管段内部构造简化图如图1所示，换能器1，2相对于管道轴线的安装角是θ，管径为D，两个换能器之间的距离为L，流体流动方向如下图：图1时差法测量流量原理图当流体以速度u流动时，超声波的实际传播速度c是声速c0和流体在声道方向上的速度分量ucosθ的叠加，即：超声波信号在流体中顺流和逆流的传播时间分别为：t顺=t2→1==〔2.1〕t逆=t1→2==〔2.2〕式中：D为一次装置管内径，L为两个换能器之间的距离，C0为超声波在静止介质中的传播速度，故：△t=t顺—t逆=〔2.3〕一般情况下，声波在液体中的传播速度C0在1000m/s以上，而多数工业系统中的流速远小于声速，即u2"c02，所以时间差可以近似简化为：△t=t顺—t逆=〔2.4〕即：〔2.5〕体积流量表达式：〔2.6〕式中：A为管道的横截面积，u为流体在管道内的流动速度。由上式可以看出，由于一次装置内径D，换能器与轴线夹角θ可以通过实际情况测量得到，所以当声速一定时，只要测的时间差△t，就可以按式子〔2.5〕求得流速u，进而由式子〔2.6〕求得流体的体积流量。时差法测量流量具有以下几个特点：a).由于时差法超声波流量计采用了精细的时差法检测，并有信号的自动跟踪，温度的自动补偿等先进的技术，因此它具有运行稳定，计量准确可靠，仪表的运算和显示精度较高的特点。b).对被测介质来说，超声波流量计属于非接触测量，在安装过程中不损坏管路，可以满足工厂工艺生产用水不断流的要求，安装极为方便。c).在测量时，管道中没有节流器件，流体不存在压力损失，非常有利于节约能源。d).电子线路的集成度高，几乎不需要维护与修理，编程灵活简单，用户输入的数据为常用的原始数据，不需要经过人工的计算，省时省力。3.超声波信号延迟时间的相关估计超声波信号的延迟时间TOF(TimeofFlight)ADDINEN.CITEADDINEN.CITE.DATA[1,2]或TDOA(TimeDelayofArrival)ADDINEN.CITEADDINEN.CITE.DATA[3]估计方法中，相关估计是一类准确的方法ADDINEN.CITEADDINEN.CITE.DATA[4-7]。尽管有一些超声波测量应用了小波变换、HHT以及FrFT等分析理论与技术ADDINEN.CITEADDINEN.CITE.DATA[8-10],相关估计仍是TOF测量的根底。但TOF较长时，即使采用FFT计算相关函数，较高的计算复杂性仍然阻碍其实时应用。其实，利用TOF估计中鼓励信号的一些特性可以简化相关函数计算。常用的鼓励信号有线性/非线性调频信号(LFM/NLFM)ADDINEN.CITE<EndNote><Cite><Author>Pollakowski</Author><Year>1994</Year><Reum>117</Reum><DisplayTe*t>[11]</DisplayTe*t><record><rec-number>117</rec-number><foreign-keys><keyapp="EN"db-id="2zeddv0thv**syeevvi*v55t552d55vd05e2">117</key></foreign-keys><ref-typename="JournalArticle">17</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">Pollakowski,M.</style></author><author><styleface="normal"font="default"charset="134"size="100%">Ermert,H.</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">ChirpSignalMatchingandSignalPowerOptimizationinPulse-EchoModeUltrasonicNondestructiveTesting</style></title><secondary-title><styleface="normal"font="default"charset="134"size="100%">Ultrasonics,FerroelectricsandFrequencyControl,IEEETransactionson</style></secondary-title></titles><periodical><full-title>Ultrasonics,FerroelectricsandFrequencyControl,IEEETransactionson</full-title></periodical><pages><styleface="normal"font="default"charset="134"size="100%">655-670</style></pages><volume><styleface="normal"font="default"charset="134"size="100%">41</style></volume><number><styleface="normal"font="default"charset="134"size="100%">5</style></number><section><styleface="normal"font="default"charset="134"size="100%">655</style></section><dates><year><styleface="normal"font="default"charset="134"size="100%">1994</style></year></dates><urls></urls></record></Cite></EndNote>[11]、PRBSADDINEN.CITEADDINEN.CITE.DATA[6,12]及混沌调制信号ADDINEN.CITEADDINEN.CITE.DATA[13,14]等。PRBS及混沌信号的相关函数主副瓣差异明显，易于峰值搜索，但超声传感器是一种窄带器件，需采用调制方式进展传输ADDINEN.CITEADDINEN.CITE.DATA[14,15]，且每一位码元必须包含足够多的载波周期，以免起振与拖尾及信号群延迟的影响ADDINEN.CITEADDINEN.CITE.DATA[16,17]。但解调后信号的相关函数主瓣将被展宽。和PRBS及混沌信号相比，LFM信号虽然不易克制多传感器应用中的Cross-talkingADDINEN.CITEADDINEN.CITE.DATA[13,18]，但可以基带传输，更重要的是其相关函数主瓣窄ADDINEN.CITE<EndNote><Cite><Author>Pollakowski</Author><Year>1994</Year><Reum>117</Reum><DisplayTe*t>[11]</DisplayTe*t><record><rec-number>117</rec-number><foreign-keys><keyapp="EN"db-id="2zeddv0thv**syeevvi*v55t552d55vd05e2">117</key></foreign-keys><ref-typename="JournalArticle">17</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">Pollakowski,M.</style></author><author><styleface="normal"font="default"charset="134"size="100%">Ermert,H.</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">ChirpSignalMatchingandSignalPowerOptimizationinPulse-EchoModeUltrasonicNondestructiveTesting</style></title><secondary-title><styleface="normal"font="default"charset="134"size="100%">Ultrasonics,FerroelectricsandFrequencyControl,IEEETransactionson</style></secondary-title></titles><periodical><full-title>Ultrasonics,FerroelectricsandFrequencyControl,IEEETransactionson</full-title></periodical><pages><styleface="normal"font="default"charset="134"size="100%">655-670</style></pages><volume><styleface="normal"font="default"charset="134"size="100%">41</style></volume><number><styleface="normal"font="default"charset="134"size="100%">5</style></number><section><styleface="normal"font="default"charset="134"size="100%">655</style></section><dates><year><styleface="normal"font="default"charset="134"size="100%">1994</style></year></dates><urls></urls></record></Cite></EndNote>[11]，一般为两个基带信号周期ADDINEN.CITE<EndNote><Cite><Author>Zhao</Author><Year>2009</Year><Reum>59</Reum><DisplayTe*t>[19]</DisplayTe*t><record><rec-number>59</rec-number><foreign-keys><keyapp="EN"db-id="2zeddv0thv**syeevvi*v55t552d55vd05e2">59</key></foreign-keys><ref-typename="ConferenceProceedings">10</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">Zhao,*iaochuan</style></author><author><styleface="normal"font="default"charset="134"size="100%">Luo,Qingsheng</style></author><author><styleface="normal"font="default"charset="134"size="100%">Han,Baoling</style></author><author><styleface="normal"font="default"charset="134"size="100%">Li,*iyu</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">Anovelultrasonicrangingsystembasedontheself-correlationofpseudo-randomsequence</style></title><secondary-title><styleface="normal"font="default"charset="134"size="100%">InformationandAutomation,2009.ICIA'09.InternationalConferenceon</style></secondary-title><alt-title><styleface="normal"font="default"charset="134"size="100%">InformationandAutomation,2009.ICIA'09.InternationalConferenceon</style></alt-title></titles><pages><styleface="normal"font="default"charset="134"size="100%">1124-1128</style></pages><keywords><keyword>fastFouriertransforms</keyword><keyword>randomsequences</keyword><keyword>ultrasonictransducers</keyword><keyword>FFT</keyword><keyword>fastFouriertransform</keyword><keyword>modulardesign</keyword><keyword>pseudo-randomsequenceself-correlation</keyword><keyword>ultrasonicpropagationvelocitymeasurement</keyword><keyword>ultrasonicrangingsystem</keyword></keywords><dates><year><styleface="normal"font="default"charset="134"size="100%">2009</style></year></dates><urls></urls></record></Cite></EndNote>[19]。就TOF估计精度及鲁棒性而言，相关函数主副瓣高度的差距只是稳定可靠测量的一个因素，主瓣宽度对测量的可靠性及鲁棒性也有重要影响。本课题提出一种采用离散方波LFM(DRLFM)信号鼓励，用加减运算实现互相关函数的递推算法，以实现TOF的快速估计。文中首先根据超声波传播模型，分析DRLFM作为鼓励信号进展TOF估计的可行性，然后介绍互相关函数递推算法并分析其复杂性，最后用实验与仿真对算法进展验证。3.1TOF相关估计这里以对射式超声波发收组件为例分析TOF的相关估计。图1是信号产生与传输的过程ADDINEN.CITEADDINEN.CITE.DATA[20]，、及为超声波发生、传输及传感环节的传递函数，是鼓励信号，是换能器输出信号，为传感器输出信号，表示多换能器同时工作时的Cross-talkingADDINEN.CITEADDINEN.CITE.DATA[13]，为其它噪声。上述各信号的频域形式分别为及。显然，。图1TOF估计中超声波信号的传输多数介质中，超声波的传输可以看成无畸变传输，即,其中是增益，是信号传输时间，所以可以写成，MACROBUTTONAuroraSupport.PasteReferenceOrEditStyle(IF0=SEQEqChapter\c\*arabic0"""SEQEqChapter\c\*arabic0."IF0=SEQEqSection\c\*arabic0"""SEQEqSection\c\*arabic0."SEQEq\*arabic1)ADDIN定义的互相关函数为，MACROBUTTONAuroraSupport.PasteReferenceOrEditStyle(IF0=SEQEqChapter\c\*arabic0"""SEQEqChapter\c\*arabic0."IF0=SEQEqSection\c\*arabic0"""SEQEqSection\c\*arabic0."SEQEq\*arabic2)ADDIN如果随机且和不相关，待测TOF就是峰值点所对应的时间。由于与基频相位差较小，且可以用稳定的系统误差来描述，因此可以用代替来估计TOF。可以导出，,MACROBUTTONAuroraSupport.PasteReferenceOrEditStyle(IF0=SEQEqChapter\c\*arabic0"""SEQEqChapter\c\*arabic0."IF0=SEQEqSection\c\*arabic0"""SEQEqSection\c\*arabic0."SEQEq\*arabic3)ADDIN这里是的自相关函数，表示单位冲激响应的自相关函数，由超声波换能器的特性决定。可以看出，对于具体的超声波换能器来说，的时域波形取决于鼓励信号的自相关函数。3.2鼓励信号的自相关函数及频谱合理选择鼓励信号有助于稳定可靠地估计TOF。图2是63位幅度为1、周期为127的PRBS信号及其自相关函数、PRBS调制信号及其频谱。载波频率为1MHz，每位PRBS包含8个载波周期，驱动信号帧长约为0.5ms。可以看出，PRBS自相关函数主瓣幅度很高，宽度为0.016ms。图3是一帧0.05ms的LFM信号、自相关函数及其频谱。信号频率从0.83MHz线性增加到1.2MHz,中心频率为1MHz。图4是与图3相对应的DRLFM信号及其相关函数与频谱，采样频率为20Msps。从图3(b)及图4(b)可以看出，相关函数的主瓣很窄，第一副瓣高度约为主瓣高度的46%，而其宽度却比PRBS调制信号窄得多，驱动信号的帧长也比PRBS调制信号短得多。从频谱图上看，LFM信号带宽略宽于PRBS调制信号，但仍在其带宽内，可以高效率传输。由于LFM及DRLFM是连续渐变信号，起振与拖尾现象ADDINEN.CITE<EndNote><Cite><Author>童峰</Author><Year>1998</Year><Reum>5</Reum><DisplayTe*t>[16]</DisplayTe*t><record><rec-number>5</rec-number><foreign-keys><keyapp="EN"db-id="2zeddv0thv**syeevvi*v55t552d55vd05e2">5</key></foreign-keys><ref-typename="JournalArticle">17</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">童峰</style></author><author><styleface="normal"font="default"charset="134"size="100%">许水源</style></author><author><styleface="normal"font="default"charset="134"size="100%">许天增</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">一种高精度超声波测距处理方法</style></title><secondary-title><styleface="normal"font="default"charset="134"size="100%">**大学学报</style></secondary-title></titles><periodical><full-title>**大学学报</full-title></periodical><pages><styleface="normal"font="default"charset="134"size="100%">507-513</style></pages><volume><styleface="normal"font="default"charset="134"size="100%">37</style></volume><number><styleface="normal"font="default"charset="134"size="100%">4</style></number><section><styleface="normal"font="default"charset="134"size="100%">507</style></section><dates><year><styleface="normal"font="default"charset="134"size="100%">1998</style></year></dates><urls></urls></record></Cite></EndNote>[16]的影响不大。图2单位PRBS的自相关函数及其调制信号频谱,PRBS周期为127，载波频率是1MHz(a)63位PRBS；(b)PRBS自相关函数；(c)PRBS调制信号；(d)PRBS调制信号的频谱图3LFM的自相关函数及其频谱(a)LFM信号；(b)自相关函数；(c)LFM的Fourier变换幅度图4DRLFM的自相关函数及其Fourier变换模3.3DRLFM鼓励的超声波TOF递推估计3.3.1计算方法给定离散鼓励信号，是传感器的接收信号。对的互相关函数为，MACROBUTTONAuroraSupport.PasteReferenceOrEditStyle(IF0=SEQEqChapter\c\*arabic0"""SEQEqChapter\c\*arabic0."IF0=SEQEqSection\c\*arabic0"""SEQEqSection\c\*arabic0."SEQEq\*arabic4)ADDIN可以导出，，MACROBUTTONAuroraSupport.PasteReferenceOrEditStyle(IF0=SEQEqChapter\c\*arabic0"""SEQEqChapter\c\*arabic0."IF0=SEQEqSection\c\*arabic0"""SEQEqSection\c\*arabic0."SEQEq\*arabic5)ADDIN其中，MACROBUTTONAuroraSupport.PasteReferenceOrEditStyle(IF0=SEQEqChapter\c\*arabic0"""SEQEqChapter\c\*arabic0."IF0=SEQEqSection\c\*arabic0"""SEQEqSection\c\*arabic0."SEQEq\*arabic6)ADDIN而，图5是一帧幅度为1的DRLFM信号及其，可以看出中大局部为零,其它是2或-2。因此,式REFR_eq_7F13840E77C744B38206E682617CA595(6)中,只需要计算非零的情况。所以有：，这里分别是为+2及-2的序号，和分别为+2及-2的总数。图5及波形图算法复杂性分析对于中心频率为、采样频率的N点DRLFM信号及kN点的,计算约需定点加法运算。而采用FFT计算相关函数,需要计算个数据点的FFT两次，逆FFT一次。为保证FFT精度，层的FFT一般采用浮点运算，其复杂性为浮点运算，需要实浮点乘法次及实浮点加法次。不考虑硬件实现的复杂性，假设浮点乘法、浮点加法运算时间是定点加法运算时间的两倍，FFT基的计算复杂性与递推算法复杂性之比为，MACROBUTTONAuroraSupport.PasteReferenceOrEditStyle(IF0=SEQEqChapter\c\*arabic0"""SEQEqChapter\c\*arabic0."IF0=SEQEqSection\c\*arabic0"""SEQEqSection\c\*arabic0."SEQEq\*arabic7)ADDIN图6是不同的采样倍数,不同的数据长度,据式REFR_eq_69A25455B9154E1AB03858C7E5819955(7)计算出的随鼓励信号时长变化的曲线。可以看出，递推算法比FFT基的估计算法快。图6，不同时随鼓励信号帧长的变化,是变化的3.3.3算法特点1．递推算法只需简单的加减法运算，比基于FFT的算法快，而且可以并行运行，适合基于FPGA或CPLD的硬件实现；2．及的存储空间外，不需要其它存储空间，更易于FPGA或CPLD的硬件实现。3．递推算法是无乘积的单层次运算，防止了FFT算法中多层次运算带来的误差。4仿真与试验为了验证递推算法的有效性，我们构建了TOF估计的实验系统，对递推估计算法进展验证与仿真分析。4.1实验系统如图7所示，实验系统包括DRLFM信号产生及驱动、超声换能器及接收器、信号调理与采集局部以及控制与RS232接口局部。其中驱动信号发生器DRLFM、信号采集控制器Controller由CPLD器件构成。超声波发生器UT与接收器US置于玻璃水槽中，间距可调，最大为50cm。测量时，PC机通过RS232接口向CPU发出测试命令，CPU控制基于DDFS的DLFM发生器，产生一帧定长的DRLFM脉冲序列，与此同时启动采样电路。DRLFM序列经母线电压为12V的桥式驱动电路Driver驱动，输入到超声换能器UT，US接收信号经Conditioning调理后输入ADC，Controller控制ADC转换，并将转换结果写入RAM。采样周期完成后，CPU将存放在RAM中的数据通过RS232口上传到PC机，PC机读取采集数据后对其进展处理。图7实验系统原理图系统中ADC为10bit/20Msps的TLC876，CPLD采用EPF1270T144C5，有源晶振*的频率为80MHz。传感器是中心频率MHz、带宽为=0.2MHz的防水型超声波发收组件。实验中室温度保持27℃。DRLFM为，MACROBUTTONAuroraSupport.PasteReferenceOrEditStyle(IF0=SEQEqChapter\c\*arabic0"""SEQEqChapter\c\*arabic0."IF0=SEQEqSection\c\*arabic0"""SEQEqSection\c\*arabic0."SEQEq\*arabic8)ADDIN，，这里频率变化率。需要说明的是相位截尾产生的杂散谱ADDINEN.CITE<EndNote><Cite><Author>赵正敏</Author><Year>2000</Year><Reum>201</Reum><DisplayTe*t>[21]</DisplayTe*t><record><rec-number>201</rec-number><foreign-keys><keyapp="EN"db-id="2zeddv0thv**syeevvi*v55t552d55vd05e2">201</key></foreign-keys><ref-typename="JournalArticle">17</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">赵正敏</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">相位截尾对DDS输出频谱影响的分析与计算</style></title><secondary-title><styleface="normal"font="default"charset="134"size="100%">东南大学学报:自然科学版</style></secondary-title></titles><periodical><full-title>东南大学学报:自然科学版</full-title></periodical><pages><styleface="normal"font="default"charset="134"size="100%">141-145</style></pages><volume><styleface="normal"font="default"charset="134"size="100%">30</style></volume><number><styleface="normal"font="default"charset="134"size="100%">001</style></number><dates><year><styleface="normal"font="default"charset="134"size="100%">2000</style></year></dates><urls></urls></record></Cite></EndNote>[21]，对相关分析是有益的。4.2短TOF的估计图8是一次短距离测量结果，传感器平面相距约5cm，驱动脉冲帧长40us,Msps。图8短TOF实验结果，=20Msps(a)归一化的鼓励信号，(b)采集信号，(c)相关函数，(d)相关函数(局部)及其拟合结果可以看出，的峰值对应时间为3.55us,测量结果的读数分辨率为0.05us。为了进一步提高分辨率，我们对峰值两边的10个点进展拟合，并求其拟合函数的峰值,(d)是计算结果，时间分辨率准确到ns时，结果为3.558us。图9是一样条件下，重复5000次实验的统计图。统计说明，测量结果呈正态分布，算法具有较好的稳定性。当TOF短时平稳时，采用连续屡次平均的方法可以进一步提高测量的精度。图95000次重复TOF估计的结果分布4．3长TOF的估计TOF较长时，的长度kN很大，信噪比变低。本实验主要考察验对射式淤泥密度测量系统中，算法的精度。为防止淤泥挂搭，探头距离一般在10-40cm。对于距离更远、TOF更长的应用，可以采用粗测与精测两步进展，粗测大致确定TOF的*围。为了提高长距离测量时的信噪比，可以加强鼓励信号幅度，所用的超声信号发生器，最大驱动电压为300Vpp.图10较长TOF的估计实验结果实验中，传感器平面相距约50cm，驱动脉冲帧长为40us，采样数据长度预计为。图10是一组测量结果，(a)是驱动信号，(b)是采集的数据，(c)及(d)是其相关函数，测量结果为TOF=33.55us，分辨率为50ns。此时，各有46个+2和-2，每一点相关函数只需计算加减法92次，总共约需做147200次定点加法运算。采用FFT算法，则需要161841次实数浮点乘法及379000次浮点加法运算。不仅时间长，且需要额外的存储空间。如果采用粗测与精测相结合的方法，计算量还可以进一步降低。5.系统的硬件采用DRLFM鼓励，并采用递推相关函数估计算法超声波TOF测量中，只需要加法减法运算，降低了计算复杂性，不需要额外的存储空间，并且可以并行实现，适合基于FPGA或CPLD的硬件实现。信号采用基带传输，不需要调制解调，相关函数主瓣较窄，不仅提高了信号传输效率，也提高了算法的精度及鲁棒性。这种递推算法为高频、高精度超声测量提供了一种有效的方法。5.1系统的硬件方案图11流量监测系统原理图前面已经简单的给出了系统方案，这里再详细的给出系统的硬件方案，如下图，系统分成三个局部：以模拟电路为主的前端电路，以DSP为核心信号处理电路，以MCU为核心的后端效劳电路。下面将详细的介绍各局部的硬件组成和功能。5.2传感器的选择图.12超声传感器的阻抗特性。(a)阻抗特性图。(b)相位特性图。采用收发一体超声波传感器，中心频率，带宽为0.4MHz,特性如图12所示传感器选择条件：首先要防水，这样不管是用在管道上，还是用在液体中都可以；超声波在水中传输时，中心频1MHz左右时。传输效率较好；收发一体，便于双向交替发送于接收。5.3驱动信号产生图*是本文采用可编程器件实现的DRLFM信号发生器框图，系统主时钟频率为。在每个脉冲串周期中，频率变化率和频率控制字累加器一次，实现频率的线性增加，频率控制字的初值为，相位累加器PA与N位频率控制字累加一次，实现相位的调频变化，即，每个时钟内，N位相位累加器的高M位作为LUT的地址，LUT(Look-Up-Table)中存放的D位的正弦信号波形数据在M位地址信号的寻址后输出送数字比拟器p和给定的参考信号Dref进展比拟,产生方波信号。其调谐方程为，分辨率为；根据超声换能器的中心频、带宽和信号持续时间，可以确定DRLFM的起始频率控制,调频速度。改变FW可以改变DRLFM的中心频率；设定可以设定调频速度，当时输出为频率稳定的方波信号；改变参考电压数据，其方波的占空比可调。图.*基于DDFS的DRLFM发生器5.4参数选择与计算例本文以中心频率为1MHz的防水型超声波发与收对为例，选择并计算DRLFM的相关参数。N=40，M=10，D=8，相位累加器采用流水线构造ADDINEN.CITE<EndNote><Cite><Author>赵正敏</Author><Year>2001</Year><Reum>219</Reum><DisplayTe*t>[9]</DisplayTe*t><record><rec-number>219</rec-number><foreign-keys><keyapp="EN"db-id="2r0*52adfsvfa6e5s0f5a5r3s9f0fzss5rta">219</key></foreign-keys><ref-typename="JournalArticle">17</ref-type><contributors><authors><author>赵正敏</author></authors></contributors><titles><title>DDS流水线构造的改良</title><secondary-title>APPLICATIONOFELECTRONICTECHNIQUE</secondary-title></titles><periodical><full-title>APPLICATIONOFELECTRONICTECHNIQUE</full-title></periodical><volume>27</volume><number>3</number><dates><year>2001</year></dates><publisher>万方数据资源系统</publisher><urls></urls></record></Cite></EndNote>[9]，LUT中的四分之一相限的正弦波形数据为8位256字节，可以在FPGA或CPLD内部实现。p为八位比拟器。当MHz，Hz，超声换能器的中心频MHz、带宽MHz，信号持续时间ms，可以确定DRLFM的起始频率控制；,相对误差为小于0.01%。和其它的DDFS电路一样，当频率控制字为非2的整数幂时，输出信号中会产生轻微的相位截尾误差，但当特别是当很高时，这种相位抖动远小于模拟比拟器及模拟信号中的干扰ADDINEN.CITE<EndNote><Cite><Author>赵正敏</Author><Year>2000</Year><Reum>201</Reum><DisplayTe*t>[10]</DisplayTe*t><record><rec-number>201</rec-number><foreign-keys><keyapp="EN"db-id="2zeddv0thv**syeevvi*v55t552d55vd05e2">201</key></foreign-keys><ref-typename="JournalArticle">17</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">赵正敏</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">相位截尾对DDS输出频谱影响的分析与计算</style></title><secondary-title><styleface="normal"font="default"charset="134"size="100%">东南大学学报:自然科学版</style></secondary-title></titles><periodical><full-title>东南大学学报:自然科学版</full-title></periodical><pages><styleface="normal"font="default"charset="134"size="100%">141-145</style></pages><volume><styleface="normal"font="default"charset="134"size="100%">30</style></volume><number><styleface="normal"font="default"charset="134"size="100%">001</style></number><dates><year><styleface="normal"font="default"charset="134"size="100%">2000</style></year></dates><urls></urls></record></Cite></EndNote>[10]。本系统中采用的MA*II，器件最高工作频率可达300MHz。此外，由于采用相关估计算法估计TDOA，适当的相位抖动对相关估计的精度是有益的。5.5切换单元电路设计这一局部的作用是用来切换两个探头和发射、接收电路之间的连接的。作为对模拟信号的切换，可以有以下三个方案可以选择。用继电器进展切换。当信号接通后，由于继电器实际上就是导线，所以不存在信号失真的现象。但继电器的开关频率有限，而且有一定的总开关次数限制，一般在100万次。采用模拟开关。起初作者采用了这种方案，模拟开关开关频率高，对信号的影响比拟小，但在实用的时候发现，由于切换的双方是接收到的微弱的超声信号和用于发射的高压信号，难以找到既能够承受高压又能使得传输的微弱信号失真较小的芯片，所以最终放弃了此方案。采用分立元器件，利用二极管的开关特性来控制开关。这种方法能够有很高的开关频率，能承受高压，但作者在实验的过程中发现信号失真太大，估计跟发射信号的功率有关。鉴于以上原因，作者暂时选用了继电器方案，这样可以简化工程的难度，而且控制继电器开关所需要的12V电平，发射电路中的CD4504正好能够提供。5.6放大电路设计系统中接收到的超声波信号有以下特点：从超声探头接收到的信号的幅值的*围大概为0.1mV～10mV，而一般ADC需要采样的信号的最大幅值为5V，所以得要放大54dB～94dB，即放大电路的增益为74dB±20dB；接收到的超声信号是一个以2.5MHz为中心的窄带信号；由于后续信号处理采用的是相关算法，所以对放大电路的抗干扰要求不高；鉴于以上特点，作者提出了一种可控谐振三级放大电路的方案。其中前两节是以MA*435/436为核心的固定谐振放大电路，第三级是以AD603为核心的可程控增益放大电路。考虑ADC603，可以实现±20dB的变换。这样MA*435、MA*436两级得要放大74dB，每一级需要37dB即可。如果信号的*围不在0.1mv～10mv之间，则需要通过手工调节可变电阻来实现。由于超声是一个窄带信号，且由图4.5看出在谐振频率时，根本满足每一级37dB的要求，即使实际中不满足，可以通过调节RSET或者RL来到达目的。所以利用MA*435/436谐振放大器的功能来实现超声信号的前两级的固定增益放大的方案是可行的。采用MA*435/436构成放大电路的原理图。为了防止过压，在输入端应该加一对二极管用以电路的保护。该电路能很好的对1MHz的窄带信号进展谐振放大作用。图12放大电路原理图5.7采样电路设计考虑到采样频率等因素，采用TLC876(ADC器件,采样频率fs=20Msps/12bit)的高速信号采集电路，用于对超声信号的回波进展采集。本文设计的采样频率为1MHz，综合考虑整个系统的本钱和精度要求，作出的设计图如下：在使用时作者将ADC芯片的SEL管脚接地，这样使得其输入*围到达最大为0～5V。由于前级电路AD603受MA*512控制后输出为±2.5V，这样得要有一个±2.5V到0～5V的转换电路。在此作者使用了一个集成运算放大器模块来实现，其示意图如图4.12所示，信号经过放大器后有Vout=2.5V+Vin，这样就完成了电压*围的转换。图中采用的运放集成电路为OPA627，这是一款高精度、高速运放，其参数完全能够满足系统在信号带宽等方面的要求。图13TLC876的应用5.8DSP系统设计C5509有32×16bit指令缓冲队列，可实现高效的块循环操作；两个17×17bit的MAC单元，可在单周期内执行两次MAC操作；1个40bit的ALU、1个40bit的桶型移位器，4个40bit的累加器可执行比C54系列DSP更高效的算术运算，在400MHz的晶振驱动下，可到达800MIPS的性能。以44.1kHz采样率的MP3数据流为例，对128kbit/s数据率的MP3数据进展解码。霍夫曼解码、IMDCT、子带合成等运算模块共需消耗1.3MIPS的CPU资源，对于平均每秒必须解码44.6帧数据来讲，总运算量为44.6×1.3＝57.98MIPS，C5509完全可以满足此速度要求。C5509还具有128K×16bit的片上RAM，其中包括64KB的DARAM、192KB的SARAM和64KB的片上ROM。与众多TMS320系列DSP处理器一样，C5509采用了哈佛构造，共有12组独立总线，其中包括3组数据读总线、2组数据写总线、5组数据地址总线、1组程序读总线和1组程序地址总线，这些总线并行地为各个计算单元提供指令和操作码，从而为高速的数据运算提供了有力的保障。C5509提供了专用的外部存储器接口〔EMIF〕，用于控制DSP与外部存储器之间所有数据的传输。可与EMIF无缝的存储器有：异步存储器〔ROM、FLASH、SRAM〕、同步突发SRAM、同步DRAM〔SDRAM〕，并可支持可选的32、16、8位数据访问。对EMIF编程时，必须根据实际的外部存储器考虑如何分配片内使能空间(CE)。通过EMIF接口，主处理器可将数据和程序置于片外，从而节省了片上硬件资源。其次，C5509有3个独立的多通道缓存串口〔McBSP〕，使得C5509能够直接与其他C55**系列DSP、多媒体数字信号编解码器等设备高速互连，这些McBSP可以提供全速双工通信，并支持128通道的收发，接收或者发送可以选择使用独立的时钟，字宽为8、12、16、20、24位任选。5.9可编程逻辑器件方案FPGA〔现场可编程门阵列〕与CPLD〔复杂可编程逻辑器件〕都是可编程逻辑器件，作者根据现有条件选用了CPLD。CPLD在系统中的作用：〔1〕为DSP以及单片机的译码，以实现DSP、单片机系统的完整性；〔2〕为串行接口的DAC提供并串转换功能，使得DAC直接映射成DSP、MCU的一个存放器，从而简化了系统对DAC控制；〔3〕为整个系统设计一个全局控制器，这样使整个系统能够协调的运行，同时使得DSP能专一处理流速的计算。图14驱动信号产生、信号采集局部原理图参考文献1. Marioli,D.,Narduzzi,C.,Offelli,etal..Digitaltime-of-flightmeasurementforultrasonicsensors.InstrumentationandMeasurement.IEEETransactionson,1992.41(1):p.93-97.2. Gueuning,F.E.,MihaiVarlan,ChristianE.Eugene,etal..AccurateDistanceMeasurementbyanAutonomousUltrasonicSystembiningTime-of-FlightandPhase-ShiftMethods.IEEETRANSACTIONSONINSTRUMENTATIONANDMEASUREMENT,1997.52(6):p.1809-1814.3. Grennberg,A.andM.Sandell.EstimationofsubsampletimedelaydifferencesinnarrowbandultrasonicechoesusingtheHilberttransformcorrelation.Ultrasonics,FerroelectricsandFrequencyControl.IEEETransactionson,1994.41(5):p.588-595.4. Bradley,J.andR.Kirlin.DelayEstimationbyE*pectedValue.Acoustics,SpeechandSignalProcessing.IEEETransactionson,198432(1):p.19-28.5. 杨惠连,徐苓安.超声相关流量计测量精度的探讨[J].仪器仪表学报,1993.14(1):p.90-95.6. 葛万成,吴凤萍.两步相关法高抗干扰超声波距离测量技术的研究[J].仪器仪表学报,2002.23(3):p.253-256.7. 潘仲明,简盈,王跃科.基于两步相关法的大量程超声波测距技术[J].电子测量与仪器学报,2006.20(005):p.73-76.8. 周知进,文泽军,卜英勇,小波降噪在超声回波信号处理中的应用.仪器仪表学报,2009.30(2):p.237-241.9.余恩军,叶庆卫,赵挺凯.基于整数提升小波时间熵的超声波回波位置检测.电子测量与仪器学报,2010(5):p.460-466.10. 杨进,文玉梅,李平.基于相关分析和近似熵的管道泄漏声信号特征提取及辨识方法.仪器仪表学报,2009(002):p.272-279.11. Pollakowski,M.andH.Ermert.ChirpSignalMatchingandSignalPowerOptimizationinPulse-EchoModeUltrasonicNondestructiveTesting.Ultrasonics,FerroelectricsandFrequencyControl,IEEETransactionson1994.41(5):p.655-670.12. Elias,C.M.andT.J.Moran.PseudorandomBinaryNoiseNDEUltrasonicCorrelationSystem.in1978UltrasonicsSymposium.1978.13. YaoZhenjing,MengQinghao,LiGenwang,etal..Non-crosstalkreal-timeultrasonicrangesystemwithoptimizedchaoticpulseposition-widthmodulatione*citation.inUltrasonicsSymposium,2008.IUS2008.IEEE.2008.14. Fortuna,L.,M.Frasca,andA.Rizzo.ChaoticPulsePositionModulationtoImprovetheEfficiencyofSonarSensors.IEEETransactionsonInstrumentationandMeasurement.,2003.52(6):p.1809-1815.15. 程晓畅,苏绍景,王跃科等.伪随机码超声扩频测距系统设计与算法[J].测试技术学报,2007.21(1):p.79-84.16. 童峰,许水源,许天增.一种高精度超声波测距处理方法[J].**大学学报,1998.37(4):p.507-513.17. 梅璐璐,林京.信号群延迟估计方法研究.仪器仪表学报,2009(001):p.50-54.18. 祝俊,陈兵,唐斌.快速多分量LFM信号的检测与参数估计方法.电子测量与仪器学报,2008.22(001):p.25-29.19. Zhao*iaochuan,LuoQingsheng,HanBaoling,etal..Anovelultrasonicrangingsystembasedontheself-correlationofpseudo-randomsequence.inInformationandAutomation,2009.ICIA'09.InternationalConferenceon.2009.10. 葛万成.相关法高抗干扰超声波距离测量中的信号处理[J].同济大学学报,2002.30(1):p.71-77.GeWancheng.SignalProcessinginCorrelationBasedRobustUltrasonicDistance