超声波流速测量系统

《超声波流速测量系统》由会员分享，可在线阅读，更多相关《超声波流速测量系统（27页珍藏版）》请在装配图网上搜索。

1、1引言1.1课题研究背景及意义研究利用超声波测量流体流量已经有数十年的历史了。1928年，法国人o.Rutten研制成功了世界上第一台超声流量计。但为了使超声波流量计有一定的精度，时差法超声波流量计要求对时间的测量有相当高的测量精度，这在当时是很难做到的。1955年，应用声循环法的MAXSO怵量计在美国研制成功，用于航空燃料油流量的测量。50年代末期，超声波流量计由理论研究阶段进入工业应用时期。但由于电子线路太复杂而未占有牢固的地位。到了80年代中后期，单片机技术的应用使超声流量计向高性能、智能化的方向发展。由于使用了单片机作中央处理单元，系统不仅可以进行复杂的数学运算和数据处理、进一步提高了

2、超声波流量计的测量精度，而月.还能设计出友好的人机界面，使系统具有参数设置、自动检错排错功能以及其他一些辅助功能，大大方便了用户的操作和使用。单片机在超声流量计中的应用，是超声流量计开始真正进入工业测量领域。1.2超声波流H计的现状近10年来，基于高速数字信号的处理技术与微处理器技术的进步，基于新型探头材料与工艺的研究，基于声道配置及流动力学的研究，超声流量测量技术取得了长足的进步，显示了它强劲的技术优势，形成了迅猛发展的势头，其潜在的巨大的生命力是显而易见的。在国外，以美国Controlotron公司和Ploysonics公司为代表的产品较多的采用数字信号处理技术，如“同步调制”和FFT技术

3、，他们广泛的采用以DSP为核心的数字处理电路，从而能够更快实时的处理超声信号，同时能够实现一些复杂的算法，如Ploysonics公司的DDF3088型是该公司的新一代全数字化便携式多普勒流量计，它采用了数字滤波和数字频谱分析技术，能自动识别多普勒信号与噪声信号，抗十扰能力强，采用了高分别率的液晶显示，可以现场对信号进行多普勒分析。在测量方法方面有的利用改进的时差法，用以消除速度受温度的影响，还有将时差法与多普勒法的组合，如Controlotron公司研制的480型超声流量计，这样使得产品的实用围更广。但由于国外的产品的价格比较高（每台约为510万元左右），所以限制了在国的大量使用。目前国的厂家

4、生产的超声波流量计虽然价格比较便宜，但总体性能差，主要用于测量比较容易的大管径管道中流体的测量。他们多采用的方法是时差法，以单片机为系统核心进行信号处理。但由于受单片机速度的限制，所以只能进行一些比较简单的算法，如通过计数来实现流速的计算。同时由于一般的时差法受温度的影响比较大，所以在精度上面就不会很高。如业楠公司生产的ZLC系列流量计是国率先采用多脉冲发射和接收宽带技术，并带有微机控制的新型超声波时差式流量计。近几年来由于DSP和超大规模集成电路技术的发展，以及以基于DSP为核心的数字信号处理电路技术日益成熟和广泛使用，使得可以采用以DSP%处理核心的超声波流量计来广泛取代国的以单片机为核心

5、的超声波流量计，从而可利用数字信号处理的一些技术来改善产品的测量精度，如FFT数字滤波、插值运算、相关运算，从而来提高系统测量的精度。通过这些改进使国流量计的性能已达到国外的水平，从而能更好的为国民经济服务。1.2.1超声波流虽计的发展趋势在石油和天然气领域快速普及最近几年超声波流量计市场的一切增长，几乎都是由于这种产品在石油和天然气领域销量的增加而带来的。超声波流量计在这一领域的销量比起以前将近翻了一倍。由于AGA船护运输标准的推广，天然气的监护运输市场已经基本成形。在2003年到2005年间，监护运输量增加了三倍多，而这一领域占超声波流量计销售的5%只有液体碳氢化合物监护运输的API和OI

6、ML标准在行业得到推广，超声波流量计才能在这个领域普及。超声波流量计具有高精确度和低总体拥有成本不论是从技术上还是从经济上看，超声波测量仪器都是流量测量的理想选择。通过多光束和数字信号处理，超声波测量仪可以实现很高的测量精确度。与传统的涡轮式仪表不同，它没有移动的元件，因此几乎不需要维修。而且，它也不会阻挡或者减慢管道中气体或者液体的流动。它能够准确地测量液态石油气产品的宽频，而不需像机械型技术那样得到验证。高灵敏度使其可以检测到管道中的任何泄漏，并可以测量和补充各种会影响监护运输领域中的测量准确度的变量。业洲和中东市场增长最大业洲和中东地区超声波流量计的增长将会是最大的。中国和印度将会在基础

7、设施和新工厂上大举投资。中国能源缺乏，为了寻找推动经济快速发展的能源，将会对其石油和天然气基础设施进行改造，并建立连接俄罗斯及其它地区供油商的管道。而由于中东地区在石油和天然气生产中的地位，这一地区将仍然是超声波流量计供应商的沃土。这一地区亦将在数个大型发电及海水淡化厂上进行投资。相比之下，北美市场的增长则显得相对平淡，但由于这一地区在石油和天然气基础设施及工业自动化上的投资，北美市场将仍然具有可观的增长。1.3超声波流虽计的市场前景随着石油和天然气工业的强势增长，以及行业对于现场设备技术的逐步接受，超声波流量计的全球市场总额将在今后5年以9.6%的复合年增长率(CAGR增长。ARCAdvis

8、oryGroup最近的一份报告指出，2005年该市场总值为2.75亿美元，到2010年将超过4.34亿美元。业洲和中东市场增长最大，业洲和中东地区超声波流量计的增长将会是最大的。中国和印度将会在基础设施和新工厂上大举投资。中国能源缺乏，为了寻找推动经济快速发展的能源，将会对其石油和天然气基础设施进行改造，并建立连接俄罗斯及其它地区供油商的管道。而由于中东地区在石油和天然气生产中的地位，这一地区将仍然是超声波流量计供应商的沃土。这一地区亦将在数个大型发电及海水淡化厂上进行投资。相比之下，北美市场的增长则显得相对平淡，但由于这一地区在石油和天然气基础设施及工业自动化上的投资，北美市场将仍然具有可观

9、的增长。1.4主要研究容本课题的任务是应用时差法测量流量。确立了系统的整体方案，在此基础上进行试验、调试，验证和改进设计方案。主要容如下：(1) 确立了系统的整体方案，这主要是在查阅国外流量计产品基本性能以及分析超声波流量计基本原理和项目指标的基础上，确立了系统以DS丽信号处理核心，并采用相关算法来计算流速的改进时差法方案；(2) 硬件电路设计：CPL探样控制系统；(3) 程序设计，应用CCSB件，在DSPR小系统上实现相关算法程序设计，进行实验测试。2超声波流量计设计方案及分析2.1超声流虽计的分类超声波技术应用于流量测量主要依据是：当超声波入射到流体后，在流体中传播的超声波就会载有流体流速

10、的信息。超声波流量计对信号的发生、传播及检测有着各种不同的设置方法，从而构成了不同原理的超声流量计，其大致可分为传播速度差法（包括：时差法、相位差法、频差法），多普勒法，相关法等等。2.2超声流虽测虽的基本原理超声流量测量技术的基本原理是利用超声波在流体中传播时所载流体的流速信息来测量流体流量的。超声波流量计一般由超声波换能器、电子线路及流量显示和累积系统三部分组成。超声波换能器将电能转换为超声波能量，将其发射并穿过被测流体，接收换能器接收到超声波信号，经电子线路放大并转换为代表流量的电信号，供显示和积算，这样就实现了流量的检测显示。2.3传播速度差法传播速度差法是根据超声波在流动的流体中，顺

11、流及逆流传播时的速度之差与被测流体流速之间的关系来求流速或流量的方法。按其所测的物理量的不同，传播速度差法乂可分为（直接）时差法、相位差法和频差法三种。时差法就是直接测量超声波顺流和逆流传播的时间差，它适用于大、中口径管道及明渠流量的测量，但时差法流量计受温度的影响较大，流体温度变化对这种系统测量精确度的影响是不能忽视的，实验表明，在使用在机玻璃声楔的条件下，流体温度每变化10C,就会给测量增加1知右的误差，并且声速的温度系数不是常数；另外，当流体的组成或密度变化时，也将引起声速的变化，从而影响测量精度。其发展方向是提高计时精度和设法降低温度对测量精度的影响。频差法是通过测量一定时间两组闭路循

12、环系统中的循环频率之差来测得流量。频差法精度高、受温度影响较小，但受环境影响较大，工作不稳定。相位差法通过把时间差转换为相位差，避免了测量微小的时间差，可提高测量的精度，但受温度的影响依然很大。2.4时差法测虽原理及主要特点本论文研究的是基于时差法的超声波液体流量计，其管段部结构简化图如图1所示，换能器1,2相对于管道轴线的安装角是0,管径为D,两个换能器之间的距离为L,流体流动方向如图所示：滑州方向图1时差法测量流量原理图当流体以速度u流动时，超声波的实际传播速度c是声速c0和流体在声道方向上的速度分量ucos0的叠加，即:cc0ucos超声波信号在流体中顺流和逆流的传播时间分别为:t顺=t

13、2t1=一Dsin(c0ucos)(2.1)L二Dc0ucossin*(c0ucos)(2.2)式中：D为一次装置管径，L为两个换能器之间的距离，C0为超声波在静止介质中的传播速度,故：t=t顺一t逆=一sin2Ducos*(c(2u2cos2)(2.3)一般情况下，声波在液体中的传播速度C0在1000m/s以上，而多数工业系统中的流速远小丁声速，即u2c02,所以时间差可以近似简化为：t=t顺一t逆=2Ducos2sin*c0(2.4)即:2.sinc0tu2Dcos(2.5)体积流量表达式:_22Dc0sint2Dcos(2.6)式中：A为管道的横截面积，u为流体在管道的流动速度。由上式可

14、以看出，由于一次装置径D,换能器与轴线夹角0可以通过实际情况测量得到，所以当声速一定时，只要测的时间差t,就可以按式子(2.5)求得流速u,进而由式子(2.6)求得流体的体积流量。时差法测量流量具有以下几个特点：a) .由于时差法超声波流量计采用了精密的时差法检测，并有信号的自动跟踪，温度的自动补偿等先进的技术，因此它具有运行稳定，计量准确可靠，仪表的运算和显示精度较高的特点。b) .对被测介质来说，超声波流量计届于非接触测量，在安装过程中不损坏管路，可以满足工厂工艺生产用水不断流的要求，安装极为方便。c) .在测量时，管道中没有节流器件，流体不存在压力损失，非常有利于节约能源。d) .电子线

15、路的集成度高，几乎不需要维护与修理，编程灵活简单，用户输入的数据为常用的原始数据，不需要经过人工的计算，省时省力。3.超声波信号延迟时间的相关估计超声波信号的延迟时间TOF(TimeofFlight)1,2或TDOTimeDelayofArrival)3估计方法中，相关估计是一类精确的方法4-7。尽管有一些超声波测量应用了小波变换、HHT以及FrFT等分析理论与技术8-10,相关估计仍是TOF测量的基础。但任戚长时，即使采用FFT计算相关函数，较高的计算复杂性仍然阻碍其实时应用。其实，利用TOF估计中激励信号的一些特性可以简化相关函数计算。常用的激励信号有线性/非线性调频信号(LFM/NLFM

16、)11、PRBS6,12及混沌调制信号13,14等。PRB敝混沌信号的相关函数主副瓣差别明显，易于峰值搜索，但超声传感器是一种窄带器件，需采用调制方式进行传输14,15,且每一位码元必须包含足够多的载波周期，以免起振与拖尾及信号群延迟的影响16,17。但解调后信号的相关函数主瓣将被展宽。和PRBS及混沌信号相比，LFM信号虽然不易克服多传感器应用中的Cross-talking13,18，但可以基带传输，更重要的是其相关函数主瓣窄11,一般为两个基带信号周期19。就TOF估计精度及鲁棒性而言，相关函数主副瓣高度的差距只是稳定可靠测量的一个因素，主瓣宽度对测量的可靠性及鲁棒性也有重要影响。本课题提

17、出一种采用离散方波LFM(DRLFW号激励，用加减运算实现互相关函数的递推算法，以实现TOF的快速估计。文中首先根据超声波传播模型，分析DRLFM为激励信号进行TOF估计的可行性，然后介绍互相关函数递推算法并分析其复杂性，最后用实验与仿真对算法进行验证。3.1TOF相关估计这里以对射式超声波发收组件为例分析TOF的相关估计。图1是信号产生与传输的过程X($).0显然,20,G)、了及S($)为超声波发生、传输及传感环节的传递函数，或是激励信号，是换能器输出信号，驰)为传感器输出信号，冷)表示多换能器同时工作时的Cross-talking13,“(t)为其它噪声。上述各信号的频域形式分别为F(,

18、s)=X(s)G(时($)+C(s)S(s)+N图1TOF古计中超声波信号的传输多数介质中，超声波的传输可以看成无畸变传输，即=其中(是增益，耳是信号传输时间，所以可以写成，I*、.(1)定义的互相关函数为，；*；匚-：，牛勇一*上(2)如果机)随机且“(t)和知)不相关，待测TOFM是&山)峰值点所对应的时间。由丁廿(苛与出)基频相位差较小，且可以用稳定的系统误差来描述，因此可以用山=f二以W+代替上山)来估计TOF可以导出，=Ks*笋如*t冰e一扣为RTy(s)=-KX(s)X(-s)G(s)S(s)e-这里是页顷勺自相关函数，如(勾表示G(s)单位冲激响应.世)的自相关函数，由超声波换能

19、器的特性决定。可以看出，对丁具体的超声波换能器来说，七盘)的时域波形取决丁激励信号的自相关函数匚-点)。3.2 激励信号的自相关函数及频谱合理选择激励信号有助丁稳定可靠地估计TOK图2是63位幅度为1、周期为127的PRBS信号及其自相关函数、PRB或制信号及其频谱。载波频率为1MHz每位PRBSfe含8个载波周期，驱动信号帧长约为0.5ms。可以看出，PRBSS相关函数主瓣幅度很高，宽度为0.016ms。图3是一帧0.05ms的LFM信号、自相关函数及其频谱。信号频率从0.83MHz线性增加到1.2MHz,中心频率为1MHz图4是与图3相对应的DRLF/号及其相关函数与频谱，采样频率为20M

20、sps从图3(b)及图4(b)可以看出，相关函数的主瓣很窄，第一副瓣高度约为主瓣高度的46%而其宽度却比PRBS制信号窄得多，驱动信号的帧长也比PRBS制信号短得多。从频谱图上看，LFM信号带宽略宽丁PRB酬制信号，但仍在其带宽，可以高效率传输。由丁LFM及DRLFM!连续渐变信号，起振与拖尾现象16的影响不大。00.10.2t(ms)0.30.40.5-1500-1000-5000500100015002000f(KHz)f400X2000图2单位PRBS勺自相关函数及其调制信号频谱，PRB制期为127,载波频率是1MHz(a)63位PRBS(b)PRBS自相关函数;(c)PRBS调制信号;

21、(d)PRBS调制信号的频谱图3LFM的自相关函数及其频谱(a)LFM信号；(b)自相关函数;(c)LFM的Fourier变换幅度10-100.010.020.030.040.05图4DRLFM勺自相关函数及其Fourier变换模DRLFM激励的超声波TOF递推估计3.3.1计算方法给定离散激励信号x(n)(n=0;1;0ON?1),y(n)(n=0;1;0?1;k?1)是传感器的接收信号。x(n)对y(n)的互相关函数为，rxy(n)=x(k)y(k+n);(4)k=0可以导出，rxy(n+1)=rxy(n)+(Zr(n),其中，X?10r(n)=x(N)y(n+N)?x(0)y(n)+0x

22、(k)y(k+n+1)(6)k=0而。x(k)=x(k+1)?x(k)，图5是一帧幅度为1的DRLFM言号及其0x(k),可以看出NP?10x(k)部分为零，其它是2或-2。因此，式(6)中，Cx(k)y(k+n+1)只需k=0要计算0x(k)非零的情况。所以有：0r(n)=x(N)y(n+N)?x(0)y(n)+2Xy(pp+n+1)?Xy(pn+n+1),i=0j=0这里pp(i);pn(j)分别是0x(k)为+2及-2的序号，I和J分别为+2及-2的总数。20-2020406080100120n图5x(n)及0x(n)波形图3.3.2算法复杂性分析对丁中心频率为f。、采样频率fs=Ksf

23、o的N点DRLFM信号x(n)及kN点的y(n)，计算rxy(n)约需定点加法运算(k?1)N2=Ks。而采用FFT计算相关函数，需要计算2kN个数据点的FFT两次，逆FFT一次。为保证FFT精度，log2(2kN)层的FFT一般采用浮点运算，其复杂性为O(kNlog2(2kN)浮点运算，需要实浮点乘法8kNlog2(2kN)+8kN次及实浮点加法20kNlog2(2kN)+4kN次。不考虑硬件实现的复杂性，假设浮点乘法、浮点加法运算时间是定点加法运算时间的两倍，FFT基的计算复杂性与递推算法复杂性之比为，_4Ks(10+7log2k+7log2N)?=N图6是不同的采样倍数Ks,不同的数据长

24、度kN,据式(7)计算出的随激励信号时长T变化的曲线。可以看出，递推算法比FFT基的估计算法快。18161412108640.050.10.150.2T(ms)图6fo=1MHz,不同fs时？随激励信号帧长T的变化f。=1MHz,fs是变化的2.3 .3算法特点1.递推算法只需简单的加减法运算，比基丁FFT的算法快，而且可以并行运行，适合基丁FPG戴CPLD勺硬件实现；2. y(n)及x(n)的存储空间外，不需要其它存储空间，更易丁FPG僦CPLD勺硬件实现。3. 递推算法是无乘积的单层次运算，避免了FFT算法中多层次运算带来的误差。4仿真与试验为了验证递推算法的有效性，我们构建了TOF估计的

25、实验系统，对递推估计算法进行验证与仿真分析。4.1实验系统如图7所示，实验系统包括DRLFM&号产生及驱动、超声换能器及接收器、信号调理与采集部分以及控制与RS232接口部分。其中驱动信号发生器DRLFM信号采集控制器Controller由CPL酬件构成。超声波发生器UT与接收器US置丁玻璃水槽中，间距可调，最大为50cm测量时,PC机通过RS232接口向CPUK出测试命令,CPU空制基丁DDFS勺DLFMK生器，产生一帧定长的DRLF咖冲序列，与此同时启动采样电路。DRLF时列经母线电压为12V的桥式驱动电路Driver驱动，输入到超声换能器UT,US接收信号经Conditioning调理后

26、输入ADGController控制ADC专换，并将转换结果写入RAM采样周期完成后，CPlW存放在RA的数据通过RS232口上传到PC机，PC机读取采集数据后对其进行处理。图7实验系统原理图系统中AD010bit/20Msps的TLC876CPL酥用EPF1270T144C5有源晶振X的频率为80MHz传感器是中心频率fo1MHz带宽为B=0.2MHz的防水型超声波发收组件。实验中室温度保持27C。DRLFMfcx(n)=signfsin2?(fo?B+an)nTg,(8)(n=0;1;2;:;N?1),T=1=fs,这里频率变化率a=2B=N。需要说明的是x(n)相位截尾产生的杂散谱21,对

27、相关分析是有益的。4.2短TOF的估计图8是一次短距离测量结果，传感器平面相距约5cm,驱动脉冲帧长40us,=n=!：Msp&图8短TOF实验结果，=20Msps(a) 归一化的激励信号x(n),(b) 采集信号矶”),(c) 相关函数,相关函数期j(n)(局部)及其拟合结果可以看出，七的峰值对应时间为3.55us,测量结果的读数分辨率为0.05us。为了进步提高分辨率，我们对峰值两边的10个点进行拟合，并求其拟合函数的峰值，(d)是计算结果,时间分辨率精确到ns时，结果为3.558us。图9是相同条件下，重复5000次实验的统计图。统计表明，测量结果呈正态分布，算法具有较好的稳定性。当TO

28、F时平稳时，采用连续多次平均的方法可以进一步提高测量的精度。图95000次重复TOF古计的结果分布4.3长TOF的估计TO商长时，如)的长度kN很大，信噪比变低。本实验主要考察验对射式淤泥密度测量系统中，算法的精度。为防止淤泥挂搭，探头距离一般在10-40cm。对丁距离更远、TOF更长的应用，可以采用粗测与精测两步进行，粗测大致确定TOF的围。为了提高长距离测量时的信噪比，可以加强激励信号幅度，所用的超声信号发生器，最大驱动电压为300Vpp.rrt(us)0.040.020-0.02-0.04图10较长TOF的估计实验结果实验中，传感器平面相距约50cm驱动脉冲帧长为40us，采样数据长度预

29、计为矣。图10是一组测量结果,(a)是驱动信号,(b)是采集的数据,(c)及(d)是其相关函数，测量结果为TOF=33.55us,分辨率为50ns。此时，小各有46个+2和-2,每一点相关函数只需计算加减法92次，总共约需做147200次定点加法运算。采用FFT算法，则需要161841次实数浮点乘法及379000次浮点加法运算。不仅时间长，且需要额外的存储空间。如果采用粗测与精测相结合的方法，计算量还可以进一步降低。5.系统的硬件采用DRLFB励，并采用递推相关函数估计算法超声波TOF测量中，只需要加法减法运算，降低了计算复杂性，不需要额外的存储空间，并且可以并行实现，适合基丁FPG戴CPLD

30、的硬件实现。信号采用基带传输，不需要调制解调，相关函数主瓣较窄，不仅提高了信号传输效率，也提高了算法的精度及鲁棒性。这种递推算法为高频、高精度超声测量提供了一种有效的方法。5.1系统的硬件方案三血町拧AIVACn工*,榛玦兰fM儿员责长携歧网LCDHP：KEY申口10*-CPLD|UDBDSP5S09SARAMFLASHSARAMTIME图ii流量监测系统原理图前面已经简单的给出了系统方案，这里再详细的给出系统的硬件方案，如图所示，系统分成三个部分：以模拟电路为主的前端电路，以DS叫核心信号处理电路，以MClfc核心的后端服务电路。卜面将详细的介绍各部分的硬件组成和功能5.2传感器的选择图.1

31、2超声传感器的阻抗特性。(a)阻抗特性图。(b)相位特性图。采用收发一体超声波传感器，中心频率,=1MH气带宽为0.4MHz,特性如图12所示传感器选择条件：1. 首先要防水，这样不管是用在管道上，还是用在液体中都可以；2. 超声波在水中传输时，中心频1MHz&右时。传输效率较好；3. 收发一体，便丁双向交替发送丁接收。5.3驱动信号产生图X是本文采用可编程器件实现的DRLF”号发生器框图，系统主时钟频率为fcLK。在每个脉冲申周期中，频率变化率和频率控制字累加器一次，实现频率的线性增加，频率控制字FW的初值为FW0，相位累加器PA与N位频率控制字FW累加一次，实现相位的调频变化，即，FW:F

32、WdfFWFW。ndfPHA:PHAFW每个时钟，N位相位累加器的高M位作为LUT的地址，LUT(Look-Up-Table)中存放的D位的正弦信号波形数据在M位地址信号的寻址后输出送数字比较器Comp和给定的参考信号Dref进行比较，产生方波信号。f叫df官、一fo.NfCLK、df.N其调谐方程为2N，分辨率为2N；根据超声换能器的中心频fo、带宽BW和信号持续时间t,可以确定DRLFM勺起始频率调频速度dFdfdtfCLKBW/2dfNFWPHAFWPHALPFLPFLUTLUTDACCompComp*一DRLFMDRLFM改变FW以改变DRLFM勺中心频率；设定dF可以设定调频速度，当

33、dF0时输出为频率稳定的方波信号；改变参考电压数据Dref,其方波的占空比可调。CLKdFdfdtfCLK34359.0，相对误差为小丁0.01%。Vref图.X基于DDFSJDRLFMt生器5.4参数选择与计算例本文以中心频率为1MH力勺防水型超声波发与收对为例，选择并计算DRLFM勺相关参数。N=4QM=10D=8,相位累加器采用流水线结构9,LUT中的四分之一相限的正弦波形数据为8位256字节，可以在FPG僦CPLC实现。Com的八位比较器。当fCLK80MHzdf7.27595761418342610-5Hz超声换能器的中心频fo1MHz带宽BW0.2MHz信号持续时间t1m可以确定D

34、RLFM勺起始频率控制FW0(foBW/2)/df(12369505812)d；和其它的DDF瓠路一样，当频率控制字为非2的整数籍时，输出信号中会产生轻微的相位截尾误差，但当特别是当fcLK很高时，这种相位抖动远小于模拟比较器及模拟信号中的十扰10。本系统中采用的MAXII,器件最高工作频率可达300MHz此外，由于采用相关估计算法估计TDOA适当的相位抖动对相关估计的精度是有益的。5.5切换单元电路设计这一部分的作用是用来切换两个探头和发射、接收电路之间的连接的。作为对模拟信号的切换，可以有以下三个方案可以选择。用继电器进行切换。当信号接通后，由于继电器实际上就是导线，所以不存在信号失真的现

35、象。但继电器的开关频率有限，而且有一定的总开关次数限制，一般在100万次。(1) 采用模拟开关。起初作者采用了这种方案，模拟开关开关频率高，对信号的影响比较小，但在实用的时候发现，由于切换的双方是接收到的微弱的超声信号和用于发射的高压信号，难以找到既能够承受高压又能使得传输的微弱信号失真较小的芯片，所以最终放弃了此方案。(2) 采用分立元器件，利用二极管的开关特性来控制开关。这种方法能够有很高的开关频率，能承受高压，但作者在实验的过程中发现信号失真太大，估计跟发射信号的功率有关。鉴于以上原因，作者暂时选用了继电器方案，这样可以简化项目的难度，而且控制继电器开关所需要的12V电平，发射电路中的C

36、D4504正好能够提供。5.6放大电路设计系统中接收到的超声波信号有以下特点：从超声探头接收到的信号的幅值的围大概为0.1mV10mV而一般ADC!要采样的信号的最大幅值为5V,所以得要放大54dB94dB,即放大电路的增益为74dB土20dB;接收到的超声信号是一个以2.5MHz中心的窄带信号；由于后续信号处理采用的是相关算法，所以对放大电路的抗十扰要求不高；鉴于以上特点，作者提出了一种可控谐振三级放大电路的方案。其中前两节是以MAX435/436为核心的固定谐振放大电路，第三级是以AD603为核心的可程控增益放大电路。考虑ADC603可以实现土20dB的变换。这样MAX435MAX43CM

37、级得要放大74dB,每一级需要37dB即可。如果信号的围不在0.1mv10mv之间，则需要通过手工调节可变电阻来实现。由丁超声是一个窄带信号，且由图4.5看出在谐振频率时，基本满足每一级37dB的要求,即使实际中不满足，可以通过调节RSE械者RL来达到目的。所以利用MAX435/436谐振放大器的功能来实现超声信号的前两级的固定增益放大的方案是可行的。采用MAX435/436t成放大电路的原理图。为了防止过压，在输入端应该加一对二极管用以电路的保护。该电路能很好的对1MHz的窄带信号进行谐振放大作用。图12放大电路原理图fs=20Msps/12bit）的高速信号5.7采样电路设计考虑到采样频率

38、等因素，采用TLC876(AD潴件，采样频率采集电路，用丁对超声信号的回波进行采集本文设计的采样频率为1MHz综合考虑整个系统的成本和精度要求，作出的设计图如下：在使用时作者将AD酒片的SEL管脚接地，这样使得其输入围达到最大为05V。由丁前级电路AD603受MAX51：控制后输出为土2.5V，这样得要有一个土2.5V到05V的转换电路。在此作者使用了一个集成运算放大器模块来实现，其示意图如图4.12所示，信号经过放大器后有V0ut=2.5V+Vin，这样就完成了电压围的转换。图中采用的运放集成电路为OPA627这是一款高精度、高速运放，其参数完全能够满足系统在信号带宽等方面的要求。图13TL

39、C876的应用5.8DSP系统设计C5509有32x16bit指令缓冲队列，可实现高效的块循环操作；两个17X17bit的MAO元，可在单周期执行两次MAO作；1个40bit的ALU1个40bit的桶型移位器，4个40bit的累加器可执行比C54系列DSP更高效的算术运算，在400MH力勺晶振驱动下，可达到800MIPS的性能。以44.1kHz采样率的MP傲据流为例，对128kbit/s数据率的MP激据进行解码。霍夫曼解码、IMDCT子带合成等运算模块共需消耗1.3MIPS的CPUS源，对丁平均每秒必须解码44.6帧数据来讲，总运算量为44.6X1.3=57.98MIPS,C5509完全可以满

40、足此速度要求。C5509还具有128KX16bit的片上RAM其中包括64KB的DARAM192KB的SARA血64KB的片上ROM与众多TMS32C列DSP处理器一样,C5509采用了哈佛结构，共有12组独立总线，其中包括3组数据读总线、2组数据写总线、5组数据地址总线、1组程序读总线和1组程序地址总线，这些总线并行地为各个计算单元提供指令和操作码，从而为高速的数据运算提供了有力的保障。C5509提供了专用的外部存储器接口（EM旧，用丁控制DSP与外部存储器之间所有数据的传输。可与EMIF无缝的存储器有：异步存储器（ROMFLASHSRAM、同步突发SRAM同步DRAMSDRAM，并可支持可

41、选的32、16、8位数据访问。对EMIF编程时，必须根据实际的外部存储器考虑如何分配片使能空间（CE）。通过EMIF接口，主处理器可将数据和程序置丁片外，从而节省了片上硬件资源。其次，C5509有3个独立的多通道缓存申口（McBSP，使得C5509能够直接与其他C55xx系列DSP多媒体数字信号编解码器等设备高速互连，这些McBSPT以提供全速双工通信，并支持128通道的收发，接收或者发送可以选择使用独立的时钟，字宽为8、12、16、20、245.9可编程逻辑器件方案FPGA(现场可编程门阵列)与CPLD(复杂可编程逻辑器件)都是可编程逻辑器件，作者根据现有条件选用了CPLDCPLDS系统中的

42、作用:(1) 为DSP以及单片机的译码，以实现DSR单片机系统的完整性;(2) 为申行接口的DA飒供并申转换功能，使得DAC接映射成DSPMCU勺一个寄存器，从DSPt它专而简化了系统对DACB制;(3) 为整个系统设计一个全局控制器，这样使整个系统能够协调的运行，同时使得一处理流速的计算。U2AGNDDRvDDD0D1AVDDAINCMLREFBSA0U4A0A1A2A3A452RM0E6CCEDSPD07一4I/O0GNDI/O4I/O5I/O6I/O7WA5ASA7ROM_A010ROM_CEROM_OED2D3D4D5D6D7D8D9REFBFNCREFTFREFTSDGNDAGNDD

43、VDDSTEYTLC876A17A16A15OEUBI/O15I/O14I/O13I/O12GNDVDDI/O11I/O10I/O9I/O8A18A14A13A12104V._16OEJ9Cj*44A17DSP_D9AD2_0SRAMIOMA0A1A2A3ADQ0A5DQ1A6DQ2A7DQ3ADQ4A9DQ5A10DQ6A11A12A13A1411ROM_D023389101112131415161718192021222324252627282930313233343536AMpC(2EPROM5PO3电源电路pU6JAD1U3两片AD(fe路nnT0A1AtTS_d!43IOIOIOIO

44、CciogndioIOIOIOIOIOI3&320iog，JDINTGCLKnCcintIO/GCLK1IOCcioGNDIOIOIOIOIOIOIOTMSTDITOCTDOOOCGBANK1AGNDDRVDDD0D1D2D3D4D5D6D7D8D9AVDDAINCM_REFBSREFBFNCREFTFREFTSDGNDAGNDDVDDSTBYTLC876/面A命|_jSOrVIDIrBANK4oo，o,IIV1OOOOgoO28V_J5C|KPLD4BANK2BANK3nCOVVOEoODDIIIOIOIOIOIOIOIOIOVCCIOGNDIOIOIOIOIOIOIOGNDINTIO/GC

45、LK3VCCINTIO/GCLK2IOIOIOIOIOGNDIOVCCIOIOIOIOIOIOIOIOIOIODA_I0DA_I2+VoutDI5DI6msbDI7OSNDXferGGlsbDI0VDI1DI2DI3DI4AD2卞CLKCLKCLK1XCC3/39112=1目131|八/_VvT嗜MSDSPD1313_tdSFspSPTDIP1DGNDU1108DS_D85678910CD3*3切9如。0二三85餐3cTV3-353433323130292827262524232221201918171615141312111098765432BPM47BFSR0BDR0ICCB8CON211

46、1.Iout1Iout2RfbVrefILEWR2WR1111918RIN0DA&路R段2P213148PINHD0P38PB4opin16OUTOUTVCV+GNI4V-MCU_TXDMCU_RXDRS232+1212S2P12J23I9BNSREnS2图14驱动信号产生、信号采集部分原理图参考文献Marioli,D.,Narduzzi,C.,Offelli,etal.Digitaltime-of-flightmeasurementforultrasonicsensors.InstrumentationandMeasurement.IEEETransactionson,1992.41(1):

47、p.93-97.Gueuning,F.E.,MihaiVarlan,ChristianE.Eugene,etal.AccurateDistanceMeasurementbyanAutonomousUltrasonicSystemCombiningTime-of-FlightandPhase-ShiftMethods.IEEETRANSACTIONSONINSTRUMENTATIONANDMEASUREMENT,1997.52(6):p.9-1814.1. Grennberg,A.andM.Sandell.Estimationofsubsampletimedelaydifferencesinna

48、rrowbandultrasonicechoesusingtheHilberttransformcorrelation.Ultrasonics,FerroelectricsandFrequencyControl.IEEETransactionson,1994.41(5):p.588-595.2. Bradley,J.andR.Kirlin.DelayEstimationbyExpectedValue.Acoustics,SpeechandSignalProcessing.IEEETransactionson,198432(1):p.19-28.3. 惠连，徐苓安.超声相关流量计测量精度的探讨J

49、.仪器仪表学报，1993.14(1):p.90-95.2002. 葛万成，吴凤萍.两步相关法高抗十扰超声波距离测量技术的研究J.仪器仪表学报，23(3):p.253-256.潘仲明，简盈，王跃科.基于两步相关法的大量程超声波测距技术J.电子测量与仪器学报，2006.20(005):p.73-76.4. 周知进，文泽军，卜英勇，小波降噪在超声回波信号处理中的应用.仪器仪表学报，2009.30(2):p.237-241.余恩军，叶庆卫，挺凯.基于整数提升小波时间嫡的超声波回波位置检测.电子测量与仪器学报，2010(5):p.460-466.进，文玉梅，平.基于相关分析和近似嫡的管道泄漏声信号特征提

50、取及辨识方法.仪器仪表学报，2009(002):p.272-279.Pollakowski,M.andH.Ermert.ChirpSignalMatchingandSignalPowerOptimizationinPulse-EchoModeUltrasonicNondestructiveTesting.Ultrasonics，FerroelectricsandFrequencyControl,IEEETransactionson1994.41(5):p.655-670.5. Elias,C.M.andT.J.Moran.PseudorandomBinaryNoiseNDEUltrasonic

51、CorrelationSystem.in1978UltrasonicsSymposium.1978.6. YaoZhenjing,MengQinghao,LiGenwang,etal.Non-crosstalkreal-timeultrasonicrangesystemwithoptimizedchaoticpulseposition-widthmodulationexcitation.inUltrasonicsSymposium,2008.IUS2008.IEEE.2008.2003. Fortuna,L.,M.Frasca,andA.Rizzo.ChaoticPulsePositionMo

52、dulationtoImprovetheEfficiencyofSonarSensors.IEEETransactionsonInstrumentationandMeasurement.,52(6):p.9-1815.2007. 程晓畅，苏绍景，王跃科等.伪随机码超声扩频测距系统设计与算法J.测试技术学报，21(1):p.79-84.7. 童峰，许水源，许天增.一种高精度超声波测距处理方法J.大学学报，1998.37(4):p.507-513.8. 梅璐璐，林京.信号群延迟估计方法研究.仪器仪表学报，2009(001):p.50-54.2008. 祝俊，兵，唐斌.快速多分量LFM信号的检测与参

53、数估计方法.电子测量与仪器学报，22(001):p.25-29.9. ZhaoXiaochuan，LuoQingsheng，HanBaoling，etal.Anovelultrasonicrangingsystembasedontheself-correlationofpseudo-randomsequence.inInformationandAutomation，2009.ICIA09.InternationalConferenceon.2009.10.葛万成.相关法高抗十扰超声波距离测量中的信号处理J.同济大学学报，2002.30(1):p.71-77.GeWancheng.SignalP

54、rocessinginCorrelationBasedRobustUltrasonicDistanceMeasurement，JOURNALOFTONGJIUNIVERSITY，2002.30(1):p.71-77.21.正敏.相位截尾对DDS输出频谱影响的分析与计算J.东南大学学报：自然科学版，2000.30(001):p.141-145.附件(程序、实物)/*超声波流量计DSP程作者：徐宋静时间：2011.3.10*/#includestdio.h#includemath.hunsignedioportport8003;/*AD采集端口*/voidkfft(pr，pi，n，k，fr，fi，

55、l，il)intn，k，l，il;doublepr，pi，fr，fi;(intit，m，is，i，j，nv，l0;doublep，q，s，vr，vi，poddr，poddi;for(it=0;it=n-1;it+)(m=it;is=0;for(i=0;i=k-1;i+)(j=m/2;is=2*is+(m-2*j);m=j;frit=pris;fiit=piis;pr0=1.0;pi0=0.0;p=6.283185306/(1.0*n);pr1=cos(p);pi1=-sin(p);if(l!=0)pi1=-pi1;for(i=2;i=n-1;i+)(p=pri-1*pr1;q=pii-1*pi

56、1;s=(pri-1+pii-1)*(pr1+pi1);pri=p-q;pii=s-p-q;for(it=0;it=0;l0-)(m=m/2;nv=2*nv;for(it=0;it=(m-1)*nv;it=it+nv)for(j=0;j=(nv/2)-1;j+)(p=prm*j*frit+j+nv/2;q=pim*j*fiit+j+nv/2;s=prm*j+pim*j;s=s*(frit+j+nv/2+fiit+j+nv/2);poddr=p-q;poddi=s-p-q;frit+j+nv/2=frit+j-poddr;fiit+j+nv/2=fiit+j-poddi;frit+j=frit+

57、j+poddr;fiit+j=fiit+j+poddi;if(l!=0)for(i=0;i=n-1;i+)(fri=fri/(1.0*n);fii=fii/(1.0*n);if(il!=0)for(i=0;i=n-1;i+)(pri=sqrt(fri*fri+fii*fii);if(fabs(fri)0)pii=90.0;elsepii=-90.0;elsepii=atan(fii/fri)*360.0/6.283185306;voidmain(void)(inti,n,k=0;doublex128,pr128,pi128,fr128,fi128,mo128;intxm,zm;int*px=(

58、int*)0x4000;int*pz=(int*)0x4080;n=128for(;)(px=(int*)0x4000;/*在这里采集数据到0x4000*/#if0for(i=0;i128;i+)(*px+=port8003;#endifpx=(int*)0x4000;for(i=0;i=n-1;i+)(xm=*px;xi=xm/32768.0;pri=xi;pii=0;px+;kfft(pr,pi,128,7,fr,fi,0,1);pz=(int*)0x4080;for(i=0;i=n-1;i+)(moi=sqrt(fri*fri+fii*fii);zm=(int)(moi*1000.0);*pz=zm;pz+;k+;部分实物图:超声波传感器CPLD式验系统电路图