C2000参赛项目报告（命题组）-基于TMS320F282345的音频范围扫频仪的设计

'....31C2000参赛项目报告（命题组）题目：基于TMS320F28234/5的音频范围扫频仪的设计学校：北京交通大学指导教师：参赛队成员名单（含个人教育简历）：31 基于TMS320F28234/5的音频范围扫频仪的设计（北京交通大学电子信息工程学院 邮编100044）摘要：本设计基于TMS320F28234DSP芯片，制作了一台音频范围数字扫频仪。系统由DSP核心、硬件调理电路和PC机交互界面三部分组成。当确定当前扫描频率时，利用DSP芯片片内的PWM外设，产生单频方波输出，经过第一级调理电路后输出同频正弦信号，通过待测网络后的正弦信号通过第二级调理电路整形放大，使用DSP芯片内部AD采样，并将得到的数据经过缓存处理得到所测网络的特性，最后经过RS232串口和PC机通信，由上位机VB交互界面直观且清晰地显示所测网络的幅频特性。关键词：扫频仪、DSP应用、TMS320F28234         DesignanaudiofrequencysweeperbasedonTMS320F28234/5Guchen，Zhujieli，Liuwei(SchoolofElectronicsandInformationEngineering,BeijingJiaotongUniversity;100044)Abstract：WehavemadeanaudiofrequencysweeperbasedonTMS320F28234inthisdesign.ThesystemconsistsoftheDSP,thehardwarecircuitandPCinterface.Whengivinginstantaneoussweep-frequency,wewillproduceasingle-frequencysquare-waveusingthePWMperipheralsinDSP.Thesquare-wavewillbeconvertedtoasinusoidalsignalonthesamefrequencyviathefistsignalconditioningcircuit.Thesinusoidalsignalwillalmostimmediatelygetthroughtheelectricalnetworkwhichshouldbemeasuredforitsamplitude-frequencycharacteristics.Thenafterthesecondsignalconditioningcircuitwhichisusedforreshapingandamplifying,thesignalwillgobackintotheDSPthroughADSampling.Wegetsomedatawhichwillconsequentlybecachedandprocessedtogettheamplitude-frequencycharacteristicsofthemeasurednetwork.Finally,viatheRS232serialporttocommunicatewithPC,weshowtheamplitudefrequencycharacteristicsofthemeasurednetworkintuitivelyandclearlyusingVBinterface.Keywords：frequencysweeper、DSPapplications、TMS320F2823431 目录1引言12系统指标要求13系统方案23.1总体介绍23.2第一级调理电路23.3带阻网络53.4第二级调理电路74系统功能模块设计84.1DSP选择与外围电路设计84.2RS232通讯接口94.3ADC接口104.4电源模块105系统软件设计115.1软件总体框图115.2上位机部分115.3下位机部分126系统关键设计与创新137系统测试结果148附录一（均衡算法）169附录二269.1PCB板照片269.2上位机交互界面313131 1引言本文设计目的是制作一台音频范围的扫频仪，该扫频仪包括下位机DSP及外围电路完成指定频率范围幅频特性测试和上位机VB图形显示界面，之间通过RS232串口电缆连接。此外为了方便测试性能，我们还制作了一个带阻网络模块。在电子测量中，经常遇到对网络的阻抗特性和传输特性进行测量的问题，其中传输特性包括增益和衰减特性、幅频特性、相频特性等。用来测量上述特性的仪器我们称为频率特性测试仪，简称扫频仪。它为被测网络的调整、校准及故障的排除提供了极大的方便。我们所设计和制作的这台音频频率范围内的扫频仪，能够实现在音频范围内对带阻网络进行其特性的测试，性能良好。DSP(DigitalSignalProcessing)数字信号处理是一门涉及许多学科而又广泛应用于许多领域的新兴学科。而DSP(DigitalSignalProcessor)是一种独特的微处理器，是以数字信号来处理大量信息的器件，也是一种特别适合于进行数字信号处理运算的微处理器，其主要应用是实时快速地实现各种数字信号处理算法。在这个设计中，我们选择了新款的低成本控制系列TMS320F28234芯片，主要使用了DSP内部PWM、SCI以及ADC外设。本设计基于TMS320F28234DSP平台，完成了扫频仪的制作与测试。2系统指标要求本设计具体指标要求如下表1所示。表1系统指标要求模块基本要求电源单电源供电扫频信号产生不得使用外部DAC或DDS芯片，只能使用C2000内部PWM信号发生器扫频信号频率范围20～20kHz，步进小于10Hz；计算出基于所使用芯片所能达到的最小分辨率信号调理电路1输出幅度应在0～3V；具体峰峰值不限输出电阻为600W带阻网络使用推荐放大器以10kHz单频信号为基准，要求带阻网络的最大衰减≥10dB带阻网络应可拆卸幅频特性测试信号调理2的输入电阻为600W应用数字信号处理技术获取带阻网络在20～20kHz间的幅频特性，不得使用硬件真有效值检测电路；使用C2000内建ADC31 幅频特性显示使用C2000的UART和电脑显示器通信，设计显示界面和简单界面；或利用模拟示波器；均衡（发挥部分）20～20kHz频率范围内的通带起伏不大于±1.5dB1系统方案3.1总体介绍本设计是以外部电路尽量简单，充分使用DSP的内部设备资源的前提完成的。信号处理部分以TMS320F28234DSP芯片作为控制和测量的核心。在硬件部分，分为第一级调理电路、带阻网络和第二级调理电路。其中，第一级调理电路采用了巴特沃斯四阶低通开关电容滤波器集成芯片TLC04来实现，设计以及操作都比较简单；带阻网络采用了改进型的双T型陷波器，经过测试，性能良好；第二级调理电路采用了集成芯片——具有MUX的零漂移、可编程增益放大器PGA113来实现，可以完美地实现对于信号的第二级调整要求。总体框图如图3.1所示。图3.1总体框图3.2第一级调理电路本部分主要实现了在音频的20Hz到20kHz中将DSP输出的单频PWM方波信号变成单频的正弦波。经过查阅资料，有几种办法实现此功能。31 首先，我们考虑了锁相环。模拟锁相环可以产生高频率稳定度的正弦信号，而且能够实现对单频信号的倍频输出，经过试验发现即使输入的是方波也具有很好的输出波形。但是常见锁相环一般应用于高频信号的生成，对于低达20Hz的信号输出非常困难。而且，锁相环电路调试尤其是VCO的制作比较困难，锁定时间较长，而且很容易失锁的情况。因此，我们经过初步的实验就放弃了这种方法。其次，考虑到PWM生成的方波就是不同能量的奇数次正弦波谐波的叠加，所以我们想到无源的RC低通滤波器，但是，由于扫频信号频率在20Hz～20kHz如此宽的范围内，很难设计一个固定的RC电路完成工作。虽然这个方案在初始阶段就被否定了，但启发了我们继续思考。如果滤波器截止频率可变，问题就迎刃而解，于是诞生了我们最后采用的方案。该方案使用了有源滤波器，它由R、C元件与运算放大器组成。有源滤波同无源滤波比较，滤波效果好，主要可以同时滤除多次及高次谐波，不会引起谐振，能够较好地达到我们的要求。通常用三个运算放大器辅以合适的电阻和电容元件就可以实现双二阶型的有源滤波。但同样无法实现截止频率可变这一设计要求，于是我们使用了TLC04，它是巴特沃斯四阶低通开关电容滤波器，理论上四阶低通滤波器适合于滤除多次谐波，其最大的优点是截止频率可以通过时钟引脚控制，非常符合本设计的需要。TLC04是一种单片巴特沃斯四阶低通开关电容滤波器，成本低廉且易于使用。该芯片截止频率稳定度只依赖于外部时钟频率的稳定性，并且截止频率是可以根据时钟来调整的，在时钟与截止频率为50:1的时候有低于的误差。该芯片的功能框图见图3.2.1所示。图3.2.1TLC04的功能框图根据TLC04性能以及典型应用电路，我们设计了第一级调理电路，见图3.2.2所示。输入两路方波信号，一路为输入信号，一路为时钟信号，开关滤波得到与输入同频正弦信号；再对所得正弦信号进行二阶低通滤波，使其高频波形平滑；其后对信号正向分压（使用0PA364搭建的正向分压电路），将输出信号幅度控制在0~3V之内。0PA364的输出电阻为0，输出端串联两个300的电阻，使得第一级调理电路的输出电阻为600，满足了阻抗匹配的要求。调理电路一具体电路详见图3.2.3~3.2.6所示。其中，图3.2.3为电压偏置电路，保证输入时钟信号的电压值幅度；图3.2.4为TLC04单电源开关滤波电路，将PWM信号转换为同频正弦信号；图3.2.5为平滑滤波电路，增益为，截止频率为；图3.2.6为有源正向分压电路，有。31 图3.2.2第一级调理电路图3.2.3电压偏置电路图3.2.4TLC04采用外部时钟驱动电路31 图3.2.5二阶低通滤波图3.2.6有源正向分压电路3.3带阻网络该部分主要是设计并制作了一个可拆卸的带阻网络，为调试提供方便。带阻网络的设计种类繁多，各不相同。包括有源和无源两大类。其中，RC带阻滤波电路即为简单的无源带阻网络，该电路只需要计算出设定电阻和电容的大小即可使用，此电路的优点是使用一级运放就可实现频率限定，没有大量复杂的计算，但是，题目要求必须使用有源滤波器，所以，这种简单的方法我们不予考虑。在有源带阻滤波器中，双T带阻滤波电路效果较好，使用较为广泛，其理论计算和设计都比较成熟。运放同相端外接一个双T带阻滤波电路，它是由无源低通电路两个R和C和高通电路两个C和R/2并联组成。有带深度电压串联负反馈的运放双T网络负载电阻无限大。双T电路是两个星形接法阻抗电路的并联。此电路有较好的限波特性。基本的双T型陷波器的设计以及其频率特性如图3.3.1所示31 图3.3.1双T型陷波器电路以及其频率特性其传递函数为：其中，R和C的值可由中心频率来决定，即：其Q值可由下式计算而得：，f为反馈系数，并且，Q值随反馈系数的增高而增大。根据上述的原理，我们设计了一种有源的改进型的双T型陷波器。电路图如图3.3.2而我们的设计中对双T网络进行了改进，增加了跟随器，以使得满足第一级调理电路的输出要求。制作中采用的运算放大器芯片是TI公司的OPA364。对带阻网络的性能测试后发现，以10kHz单频信号为基准，带阻网络的最大衰减达到了要大于10dB的要求，并且，得到了比较理想的Q值。图3.3.2带阻网络31 3.4第二级调理电路第二级调理电路主要要实现的功能是对于带阻网络输出的信号在进入ADC模块之前进行幅度上的调整，使其在0～3V之间，不致对ADC模块产生影响。为了实现将大信号缩小，小信号放大的功能，该部分我们选择采用了具有MUX的零漂移、可编程增益放大器PGA113。PGA113是单电源供电，带有输入多路复用器的PGA。多路复用通道选择和增益选择是通过一个标准的SPI接口。PGA113有两个通道输入复用器，其增益提供选择范围为（1，2，5，10，20，50，100，200），并有一个三线SPI数字接口。其典型的应用电路如图3.4.1所示。图3.4.1PGA113的典型应用电路我们应用PGA113设计的第二级调理电路如图3.4.2所示。并按照题目要求，第二级调理电路的输入电阻应该是600，所以，我们在输入端并联了1.5K和1K的电阻，使其能够满足要求。除此之外，输出信号在进入ADC采样之前，我们还加入了简单的RC低通滤波，尽量减少噪声信号干扰影响AD采样精度。31 图3.4.2第二级调理电路1系统功能模块设计系统功能模块由以下四部分组成：DSP及外围电路、通讯接口、ADC接口、电源模块，下面分别介绍。4.1DSP选择与外围电路设计本设计中我们选择DSP芯片TMS320F28234作为核心处理芯片，其外围电路包括复位电路、JTAG调试接口、时钟电路等。TMS320F28234含有丰富的片上外设资源，如ADC、事件管理器、PIE、看门狗、SCI、SPI等。使用了其片上ADC、PWM及SCI等模块，无需外扩ROM。系统中，F28234协调着整个系统各模块的有序工作及承担着信号处理的任务。复位电路采用阻容式电路，以方便进行手动复位。JTAG口方便程序的调试和烧写。DSP复位电路以及JTAG电路如图4.1.1。31 图4.1.1DSP复位电路以及JTAG接口F28234的另一个特点就是其CPU时钟频率已经可以达到150MHz，即单周期指令执行时间为6.67ns。这从而为高速采集和处理信号提供了保障。因而，系统中F28234由外部30MHZ的有源晶振提供时钟信号，并在其内部2分频经锁相环10倍频至150MHz。时钟电路图如图4.1.2所示。图4.1.230MHz时钟晶振接入电路4.2RS232通讯接口通讯模块接至F28234的SCI模块，通过9芯标准RS-232口与其它系统进行串行通讯。选用F28234片上SCIC作为串行通讯口，选用MAX232作为串口通讯信号电平转换模块的主要器件。串行通讯部分硬件连接图如图4.2所示。31 图4.2通讯接口硬件连接图4.3ADC接口系统使用了两路AD：ADINA0和ADINB0，ADC部分硬件连接图如图4.3所示。图4.3ADC接口硬件连接图4.4电源模块系统供电电路总体框图如图4.4所示。31 图4.4系统总体供电框图该模块主要使用了电源转换芯片TPS767D301。因为DSP芯片TMS320F28234需要+3.3V的工作电源和+1.9V的内核电源，该模块采用了这种电源转换芯片将+5V转换为+3.3V和+1.9V。为了减小数字部分对模拟部分的影响，系统中将数字电源与模拟电源、数字地与模拟地分开。模拟地与数字地之间接0W电阻，以抑制高频串扰。1系统软件设计5.1软件总体框图系统软件整体设计分为两大部分，即上位机PC和下位机DSP两部分的软件设计。5.2上位机部分程序说明：在VB6.0环境下开发，采用RS232协议。关于通信的数据传输协议：数据帧的帧头为2个特定的字符‘S’和‘T’，利用控件的事件驱动进行帧同步和有效数据提取处理。本界面的数据采用实时显示的方式，可以动态的看到曲线的生成，还可以通过鼠标选择特定的频点进行观察比较。上位机的流程图如图5.231 YYYYYY界面参数初始化串口初始化并打开相应COM口接受串口数据存入缓冲区接受到串口数据根据设定参数初始化扫描操作，同时计算音频范围内的有效频点按下“打开串口”按钮？按下“开始扫描”按钮？发送指令控制下位机扫描当前频率点NNN计算、存储数据并绘制谱线图帧同步，数据校验N完成所有频点扫描N开始关闭程序？N结束图5.2上位机流程图5.3下位机部分程序说明：在CCS3.3环境下开发，下位机的处理单位是一个频点，完全受上位机指令控制运行。下位机通信协议跟上位机基本相同，只是数据承载的内容有所变化，为了减少下位机的运算量，降低开发难度，提高扫描效率，在设计中将对数计算和数据缓冲区放到了上位机中进行。下位机的流程图如图5.331 YYY芯片时钟、ADC、SCI、PWM初始化帧同步根据计算结果修改PWM，开启PWM信号输出根据通信协议提取有效数据，计算当前频点SCI是否接受到数据？数据有效？根据当前频率点在保证尽可能高的分辨率下计算预分频和计数器数值NN根据信号频率和采样率计算采样点数以保证峰值/有效值采样误差最小已经完成峰值/有效值采样？N开始延时约1ms待正弦波型稳定后启动ADC组帧并通过串口返回幅度参数至上位机图5.3下位机流程图1系统关键设计与创新1、本系统采用TI公司在F2812芯片基础上进行改进后推出的F28234芯片上进行开发，在功耗、性能、稳定性方面都有很大提升，片内ADC精度也有较大提高，非常适合承担低成本条件下的控制和运算任务。31 2、基于方波就是同频正弦波及其奇数次谐波组成的成分关系，采用MAX293/297集成芯片，实现了低通滤波器截止频率程控的功能，将片上PWM输出的50%占空比方波扫频信号完美的转化为幅度基本恒定的不失真的正弦扫频信号。3、有源带阻网络使用了应用广泛，较为成熟的双T带阻滤波电路，在设计频点实现了-20dB以上的衰减。4、在ADC采样前端加入可编程运算放大器PGA113，不但提高了本设计在测量微小信号时的测量精度，还实现了AD与待测信号之间的有效隔离。5、在上位机采用VB环境进行开发，其串口控件和图形显示控件非常好的协助了本设计的实现，操作简单的同时也使得测量结果更加直观。6、在RS232通信协议的应用中自主定义了通信帧格式，使得通信过程更加规范、简单、有效，对外界干扰也具有一定的抑制作用，提高了系统的鲁棒性。1系统测试结果最后，我们进行了系统联调，测试结果如表2所示。表2系统指标完成质量模块基本要求实际效果电源单电源供电+5V供电扫频信号产生不得使用外部DAC或DDS芯片，只能使用C2000内部PWM信号发生器使用TMS320F28234的PWM外设（两路）扫频信号频率范围20～20kHz，步进小于10Hz；计算出基于所使用芯片所能达到的最小分辨率扫频信号频率范围20～20kHz，最小步进1Hz；计算最小分辨率见下文信号调理电路1输出幅度应在0～3V；具体峰峰值不限信号调理1的输出电压峰峰值约为2.62V输出电阻为600W输出电阻为600W带阻网络使用推荐放大器使用TI推荐的放大器OPA364以10kHz单频信号为基准，要求带阻网络的最大衰减≥10dB以10kHz单频信号为基准，要求带阻网络的最大衰减>21dB（位于1.1kHz附近）带阻网络应可拆卸自制带阻网络，可更换幅频特性测试信号调理2的输入电阻为600W信号调理2的输入电阻为600W应用数字信号处理技术获取带阻网络在20～20kHz间的幅频特性，不得使用硬件真有效值检测电路；使用C2000内建ADC信号经DSP内部AD采样和处理，获得带阻网络幅频特性幅频特性显示31 使用C2000的UART和电脑显示器通信，设计显示界面和简单界面；或利用模拟示波器；经过RS232串口和PC机通信，由上位机VB交互界面直观且清晰显示所测网络的幅频特性均衡20～20kHz频率范围内的通带起伏不大于±1.5dB音频范围内的通带起伏不大于±1.5dB，C2000可以实时处理（具体均衡算法分析见附录一）扫频信号最小分辨率分析：设为CPU系统主时钟频率（F28234最高可达150MHz），为PWM输出50%占空比方波频率，变量PRE表示时钟预分频寄存器值，变量DIV表示PWM分频寄存器值，变量PRD表示周期寄存器值。通常分辨率跟计数器的位宽以及计数器基础频率决定，符合以下公式：ResolutionFrequency*2^B=CounterBaseFrequency其中B为计数器位宽，可见其相对分辨率完全由B决定，28234中PWM计数器位宽为16-bit，则相对分辨率为1/65536，约为0.015%。根据上式可以推导出关系式：可见计数器基础频率并不确定，其跟CPU系统时钟以及两个寄存器的数值相关。单纯就提高分辨率来讲，系统时钟越小越好，但系统时钟又不可能太小，会影响其工作及运算效率。而PRE与DIV又要根据生成方波波的频率范围决定，也不可能为了提高分辨率任意增大。所以，无法找到绝对的最小分辨率。在本工程实际应用中，使用了150MHz的系统时钟，我们优先根据待测频率保障PRD的数值尽可能接近最大值来确定PRE和DIV的数值，尽可能的提高相对分辨率。31 附录一（均衡算法）1.1方案分析为补偿陷波器的陷波特性，实现均衡，我组提出2种可行方案。方案一：FFT处理：AD采样，将连续信号离散化，通过FFT快速傅利叶变换运算，得到输入信号的频谱。根据事先所得的陷波器陷波特性参数，在频域上进行幅频均衡处理。最后将信号通过IFFT傅氏反变换到时域，并经过外置平滑低通滤波器输出，将通频带为20Hz～20kHz内的信号电压幅度波动控制在在±1.5dB以内，实现数字均衡的目的。示意图如下所示：A/DFFT以某一频率信号为准进行幅频均衡IFFT低通平滑滤波图8.1.1FFT方案这里，均衡效果很大程度上依赖于所测得的陷波器特性参数。因为我们将利用示波器观测的陷波器音频输出信号数据（指各频率输出的峰峰值），或者直接利用由显示界面观测的陷波器音频幅频曲线，依参考频率（可选为10kHz）输出幅度进行归一化处理，得到处理后数据，再使用Matlab对数据进行拟合，得到需要修正的幅频曲线，取倒数则得到修正曲线。只要每次FFT处理得到信号频谱后，将离散频谱数据与相应的修正曲线数据相乘，即完成了频域幅度修正，再做IFFT将其还原为时域信号输出，即实现了时域内幅频均衡。使用FFT直接对频谱进行均衡过程中也可能存在问题，如还原波形帧与帧之间衔接处可能会出现毛刺，因为FFT帧处理可等效为对时域信号加矩形窗。还有由于分辨率的关系，低频成分还原可能会出现失真。不过这些问题应该可以通过叠处理的方法来解决，如下图所示：图8.1.2叠处理方法处理中丢弃每帧的坏值，按照时域顺序输出有效的数据，经过FFT处理后还原出平滑无失真的波形。现大概计算一下运算量，若我们选用80k的采样率和4096点FFT定点运算，此环境下20kHz信号一周期内能保证4次采样，频谱最小分辨力为20Hz。增大采样率和处理点数，能提高频谱分辨率和均衡效果，但是处理数据量和对数据处理速度的要求会增加。鉴于FFT的低频特性不好，数据处理量大，波形易失真等缺点，我们不采用这种方案。方案二：使用数字滤波器,直接在时域对信号进行幅度－频率补偿，如图2所示。31 A/D数字滤波D/A陷波器频谱特性滤波器特性图8.1.3数字滤波方案输出信号特性设计思路有两种：⑴求陷波器幅频特性，通过仿真计算其零极点，得到传递函数，设计与陷波器互补的带通网络及其传递函数，编程设计滤波器算法实现。具体实现方法如下：①进行电路仿真，求陷波器零极点；②由零极点求系统传递函数；③计算反传递函数；④z域变换得到数字传递函数；⑤编程设计IIR滤波器算法实现每个双二阶滤波器。此设计思路清晰，但仿真是关键，寻找零极点较困难，而且容易丢失，结果存在一定的不确定性。⑵基于原型法，级联多个滤波器，用MATLAB设计出多级直接II型滤波器来拟合实际的带通滤波器曲线，得到滤波器参数。此思路即滤波器级联拟合逼近，设计基于原型双二阶滤波器，只需给定增益，中心频率和带宽即可确定。具体实现方法如下：①用matlab画出陷波器幅频特性响应曲线；②设计双二阶参量滤波器函数，用matlab画出幅频特性响应曲线；③设计带通滤波器增益，中心频率和带宽，级联后看曲线拟合效果，逐步逼近实际所需带通滤波器的特性曲线。④将matlab设计出的增益，中心频率和带宽这3个参量代入编程实现的滤波函数中，实现滤波法均衡的效果。此设计方法的结果能较好地逼近实际陷波器的互补特性，即所需带通滤波器的特性，稳定性较好，较设计思路（1）来说更易实现。需要强调的是，一般滤波器分两种：一种是FIR滤波器，优点是采用非递归结构，不存在稳定性问题，但需要计算所有的冲激响应的系数，数据处理量大，用DSP处理可能实时效果不好；第二种是IIR滤波器，优点是可以用较少的阶数实现较高的选择特性，运算量少，但是IIR滤波器采用递归结构，设计时需要考虑稳定性问题。两者相比，IIR滤波器更佳。而参量均衡器设计方法又可有效解决以上两种滤波器存在的问题，其采用全通网络和支撑结构，如下图所示：31 Input二阶全通网络Output1/2-1k/2图8.1.4由二阶全通网络构成的参量均衡器二阶全通网络的定义：二阶全通网络在整个频率响应范围内幅度增益为1，相位变化从低频的0°变到高频的180°。其特性曲线如下图所示：图8.1.5二阶全通网络的幅度、相位响应曲线二阶全通网络的带宽和中心频率表示如下：31 ，时，参量均衡器为陷波型而且有时，参量均衡器为全通型时，参量均衡器为峰值型给定了增益，中心频率和带宽，即可确定基于全通网络的参量均衡器的幅频特性。其各种特性如下图所示：图8.1.6（a）a、b保持不变，变化k31 图8.1.6（b）k、b保持不变，变化a图8.1.6（c）a、k保持不变，变化b但是，需要注意的是，基于全通网络的参量均衡器的频谱分布在整个采样频率范围之内。确定了采样频率后，在低频范围上实现窄带高增益均衡有一定困难。如下图所示：31 图8.1.7基于全通网络的参量均衡器低频特性举例现结合我们所测得的带阻网络的幅频特性来选择均衡器类型。通过VB交互界面测得带阻网络幅频特性在MATLAB仿真环境下如下所示：图8.1.8带阻网络幅频特性曲线31 图8.1.9带阻网络与所需均衡器幅频特性曲线对比可见此带阻网络陷波集中在低频1.1kHz左右，最大衰减为21dB。根据以上分析，若选择基于全通网络的参量均衡器，那么在低频范围内实现窄带均衡需要级联很多个参量均衡器，这样既没有充分发挥其各参数互不干扰、稳定性高的优点，还可能造成参量均衡器个数过多，运算量大，无法实时处理等问题。所以我们不采用基于全通网络的参量均衡器，而采用运算量少、灵活性高的基于原型双二阶IIR带通滤波器来级联成所需的均衡器。1.1方案实现我们采用基于原型双二阶IIR带通滤波器级联构成均衡器。滤波器传递函数为，各参数如下表所示。表3原型双二阶函数设计的带通参量滤波器参数传递函数参数31 滤波器类型提升滤波器削减滤波器注：，，，其中增益，品质因数>0，中心频率为。由表3可知，只需确定滤波器增益G，中心频率f0和品质因数Q，运用MATLAB仿真，编写滤波器函数function[h,w]=para2(g,q,f0)，即可画出滤波器的幅频特性曲线，如下图所示：图8.2.1基于原型双二阶IIR带通滤波器举例此方案拟采用4个基于原型双二阶IIR带通滤波器级联成均衡器，以较好地完成均衡。首先，根据基本部分结果，得到所需均衡器的幅频特性曲线。根据所需均衡器的幅频特性曲线，查询寻找到最大陷波处中心频率f0和峰值增益G，再在中心频率左右循环扫描，寻找衰减3dB处频率，计算得到3dB带宽。根据公式,计算得到品质因数Q。31 鉴于我们设计的带阻网络为双T型有源陷波器，其陷波特性极好，衰减非常大，且带宽很小，用单一的二阶IIR滤波器补偿其峰值陷波无法达到高增益窄带宽的要求，所以在峰值陷波处采用增益分步补偿，第一步仅补偿其增益的一半。根据得到的三个参数f0、G、Q，设计第一个滤波器[h1,w]=para2(0.5*G,Q,f0)(增益补偿仅为所需的1/2)。对比第一个滤波器的特性曲线和所需均衡器的特性曲线，得到仍需补偿的特性曲线。同寻找第一个滤波器参数方法一样，寻找出仍需补偿的特性曲线的最大陷波处中心频率f01、峰值增益G1和品质因数Q1。设计第二个滤波器[h2,w]=para2(G1,Q1，f01)。依第二个滤波器的设计思路和方法一样，依次设计出第三个和第四个滤波器，并绘图显示均衡器和各滤波器的特性曲线，观察均衡效果。仿真实验结果如下图所示：图8.2.2均衡仿真结果其中，峰值增益最大的蓝色曲线是需要拟合的曲线，即陷波器的互补特性曲线；峰值增益最大的红色曲线是滤波器级联后曲线，即所设计均衡器的特性曲线；其他四条特性曲线（红、黄、蓝、绿）为每个滤波器的特性曲线。此均衡器由4个基于原型双二阶IIR带通滤波器级联而成，各特性参数为：第1个滤波器（黄色曲线表示）：G=10.6,Q=6.8,f0=1090第2个滤波器（红色曲线表示）：G=10.4,Q=6.5,f0=1100第3个滤波器（蓝色曲线表示）：G=0.9,Q=1,f0=1690第4个滤波器（绿色曲线表示）：G=0.5,Q=1,f0=700其中G表示增益(dB)，f0表示中心频率(Hz)，Q表示品质因数,采样频率为81920Hz。由图可知，所设计均衡器的特性曲线与陷波器的互补特性曲线大致重合，信号幅度在20～20kHz频率范围内的通带起伏不大于±1.5dB，表明此算法能很好补偿陷波器的陷波特性，均衡效果明显。31 该算法主要由两部分组成，为了获取滤波器参数的搜索算法和当前IIR滤波器的级联后的响应。本算法的参数就是处理数据点数量N。其中搜索算法首先需要寻找最小值点的时间运算复杂度为O(N)，再次寻找3dB带宽的复杂度最差为O(2N)，由此可见从幅频特性中获取必要的滤波器参数的一次运算需要运算O(3N)。由于上述算法一共需要迭代四次，每次生成新加入级联的IIR滤波器后要与之前的响应进行卷积，然后与目标滤波器曲线比较，继续迭代，所以一共需要完成4次卷积运算，若卷积窗口长度M，则时间运算复杂度为O(M*M/2)。设计完成后，级联后的四级二阶IIR滤波器等效于一个8阶IIR滤波器，整体算法运算量中等，在N<=512并且M<=256得情况下，使用100MHz系统时钟的C2000芯片，整体运算延迟应该在15ms以内，如果使用官方提供的C28x滤波器库SPRC082给出的IIR滤波器汇编函数，可以大幅提高运算效率，所以我们认为C2000可以使用本算法进行实时处理。31 1附录二1.1PCB板照片图9.1.1全板照片31 图9.1.2底板照片31 图9.1.3顶板照片31 图9.1.4带阻网络照片31 图9.1.5调理电路1照片31 1.1上位机交互界面图9.2.1上位机交互界面31'