专利名称：一种基于fpga的车辆振动信号滤波装置的制作方法
技术领域：
本发明涉及噪声消除装置，具体涉及一种车载噪声背景下的车辆振动信号的滤波
直ο
背景技术：
危险品车辆行驶状态数据即车辆振动信号采集的准确性和精确性是车辆姿态数据融合分析获得可靠性结果的前提条件。常用加速度传感器对车辆的俯仰、横滚角度与三个方向的加速度信号对车辆姿态进行判定，在车载环境下，加速度传感器在测量过程中由于机械转子容易漂移和受到其它激振源干扰等原因，真实振动信号往往受到严重的干扰， 如果直接使用这种振动信号会对后续分析产生较大误差，给最终的数据分析决策带来很大困难，所以需要先对振动信号进行滤波降噪处理。车载环境下加速度传感器的干扰信号来源多样且为非平稳信号，在实际工程中经常用到的低通滤波器对此类信号的处理却无能为力。由于车辆内部各部分的固有频率不同，在行驶中车辆由于路面不平、发动机振动、车轮滚动的周期振动和传动系统的不平衡造成的激励。所以，这个实际信号是一个含噪声的非平稳振动信号。由于车辆的振动信号主要集中在低频部分，噪声分布有低频和高频部分，并且存在非平稳性和随机性，很难以用传统方法来解决有用信号的提取。若低通滤波器太宽，则在滤波器后，信号仍存在有大量噪声； 若低通滤波器太窄，则将一部分有用信号认为是噪声被滤掉。目前的滤波方法都是采用软件或者DSP等硬件实现的，由于软件速度慢，对实时信号进行处理比较困难，DSP芯片功耗大，处理能力受到主频的限制等，在处理速度和抗干扰性能方面不能满足系统高速、实时的需求。
发明内容
针对上述现有技术存在的缺陷或不足，本发明的目的在于提供一种基于FPGA的车辆振动信号的滤波装置，该装置是基于FPGA的L阶FIR滤波装置。为实现上述技术任务，本发明采取如下的技术方案—种基于FPGA的车辆振动信号滤波装置，其特征在于该装置包括第一寄存器、 第二寄存器、第三寄存器、减法器、变步长因子计算电路和滤波器；参考噪声信号x(n)接入第一寄存器的输入端缓存，含噪振动信号d(n)接入第二寄存器的输入端缓存，第一寄存器的将缓存后的参考噪声信号x(n)输出至滤波器的一个输入端，减法器输出的误差信号 e(n-l)分别接入第三寄存器的输入端、滤波器的第二输入端和变步长因子计算电路的输入端，变步长因子计算电路的输出信号变步长因子μ (η-1)接入滤波器的第三输入端，第二寄存器将缓存后的含噪振动信号d(n)输入至减法器一个输入端，滤波器输出的滤波后的参考噪声信号y(n)接入减法器的另一输入端，第三寄存器输出滤波处理后的车辆振动信号；其中d(n)为η时刻的含噪振动信号，χ(η)为η时刻的参考噪声信号，e(n-l)为(n-1)时刻的误差信号，μ (η-1)为(η-1)时刻滤波器的变步长因子，初始时刻减法器与变步长因子计算电路的输出均为0，即e(l) =0, μ (1) = 0；所述变步长因子计算电路由乘法器1、乘法器2、乘法器3、加法器1、减法器1、除法器1和分频器构成实现变步长因子μ (η-1)的计算；该μ (η-1)的计算公式如下μ(η-\) = 0.08(1 -1)
IOOxe {η-1) + 1所述分频器由计数器、比较器1、比较器2、反相器和D触发器1构成，该分频器对 FPGA平台时钟进行8分频后，为除法器1提供触发时钟；所述计数器的输出端与比较器1 的输入端和比较器2的输入端连接，所述比较器1的输出端与反相器的输入端连接，所述反相器的输出端与D触发器1的ce管脚连接，所述比较器2的输出端与D触发器1的d管脚连接，所述比较器1的输出端与计数器的clr管脚连接，D触发器1的输出端与除法器1连接；所述滤波器为基于FPGA的L阶LMS滤波器，该滤波器由L个LMS_tap和L个加法器构成完成如下运算Wi (n) = Wi (η-1) +2 μ (n-1) e (η-1) χ (n-i+1)Yi (η) = Wi (η) χ (n-i+1)y (η) = Y1 (η) +y2 (η) +y3 (η)... +Yi (η)... +yL (η)其中L为彡1的自然数，1彡i彡L，y(n)为滤波器所计算的η时刻滤波后的参考噪声信号，Yi (η)为第i个LMS_tap即LMS_tapi所计算的η时刻滤波后的参考噪声信号， Wi (η)为LMS_tapin时刻的滤波权系数，初始时刻LMS_tapi的滤波权系数为0。上述滤波器中的LMS_tap个数或加法器个数L = 128。本发明与现有方法相比，具有以下优点基于FPGA的车辆振动信号去噪电路采用了 FIR串行滤波器结构，充分利用了 FPGA 芯片的逻辑资源，可实现对中低速车辆振动信号的实时处理，并且降低了滤波器的开发成本。基于FPGA的车辆振动信号去噪电路具有良好的可扩展性，在FPGA芯片乘法器内核数量允许的情况下，可配置为任意阶数的滤波器结构，从而滤波效果更好。装置所依据的算法具有计算速度快、收敛速度快、实时性高、滤波效果好等优点， 可以满足信号处理方面的实时处理需求。
图1为本发明的装置的结构示意图；图2为变步长因子计算电路的电路图；图3为分频器的电路图；图4为滤波器的电路图。以下结合附图和实施例对本发明作进一步的详细说明。
具体实施方式
对于滤波器来说，滤波算法的选取是一个方面，算法的硬件实现方案也很重要，实现平台的选取直接影响到数据处理的实时可靠性。本文选取了高速FPGA芯片实现优选的自适应滤波算法。随着通信技术、信号处理理论和计算机技术的飞速发展，处理的信号量越来越大， 对实时性的要求也越来越高。FPGA是基于查找表结构的可编程器件，具有非常丰富的寄存器资源，而且可以通过系统内的重新配置来改变逻辑功能，使设计的更新与改动十分方便。 随着FPGA密度和速度的提高，现在FPGA已经可以胜任一些原来只有专用芯片或者多个DSP 才能完成的计算任务。基于FPGA的数字信号处理系统的优点是，首先FPGA支持并行处理和流水结构，可以使用时分复用，通过多个处理单元的并行工作，实现高性能的数字信号处理。特别是那些结构性好、数据量大的算法实现。比如快速傅里叶变换(Fast Fourier Transformation, FFT)和数字滤波算法等。其次，FPGA内部越来越多的内嵌了 DSP乘法模块。这些模块是硬件模块，运行速度很高，特别适合那些需要大量乘法计算的算法。总之， FPGA是DSP设计的理想器件。本发明的基于FPGA的车辆振动信号滤波装置所使用的FPGA芯片为Xilinx公司的Virtex-II Pro XC2VP30芯片，所述的寄存器、乘法器、加法器、减法器、缓存器、计数器、 D触发器均为Xilinx公司的FPGA芯片的IP核。参考图1至图4，本发明的基于FPGA的车辆振动信号滤波装置包括第一寄存器、 第二寄存器、第三寄存器、减法器、变步长因子计算电路和滤波器；参考噪声信号x(n)接入第一寄存器的输入端，含噪振动信号d(n)接入第二寄存器的输入端，第一寄存器的输出信号x(n)接入滤波器的一个输入端，减法器输出的误差信号e(n-l)分别接入第三寄存器的输入端、滤波器的第二输入端和变步长因子计算电路的输入端，变步长因子计算电路的输出信号变步长因子μ (η-1)接入滤波器的第三输入端，第二寄存器的输出信号d(η)接入减法器一个输入端，滤波器输出的滤波后的参考噪声信号y (η)接入减法器的另一输入端，第三寄存器输出滤波处理后的车辆振动信号；其中d(n)为η时刻的含噪振动信号，χ(η)为η时刻的参考噪声信号，e(n-l)为 (η-1)时刻的误差信号，μ (η-1)为(η-1)时刻滤波器的变步长因子，初始时刻减法器与变步长因子计算电路的输出均为0，即e(l) = 0, μ (1) = 0；所述变步长因子计算电路由乘法器1、乘法器2、乘法器3、加法器1、减法器1、除法器1和分频器构成实现变步长因子μ (η-1)的计算；该μ (η-1)的计算公式如下^ (^-1) = 0.08(1- 1 )
IOOxe {η-1) + 1q) 其中，0.08为计算公式中的常系数，变步长因子计算电路是按照上述计算公式来实现的；变步长因子算法的FPGA原理框图如图2所示，减法器的输出信号e(n-l)接到乘法器1中，完成e2 (η-1)计算；乘法器1的输出e2 (η-1)接入乘法器2与ox64000 (十进制的是 100)计算IOOXe2(Ii-I)；计算结果100Xe2(n-l)输入加法器1中与oxlOOO (公式中的十进制1)相加完成100 Xe2 (η-1)+1的运算；100 X e2 (η-1)+1计算结果输入除法器1的除数端，用oxlOOO (公式中的十进制被除数1)与之相除，得到的计算结果输入减法器1与被减数οχ1000(公式中的十进制被减数1)实现(1_ιηη 2:、工J运算；减法器1的计算结果
IOOxe (/7-1) + 1
输入乘法器3中，在乘法器3中完成0+08(1-1 ，、工P运算，即与十六进制数147相
IOOxe (/7-1) + 1
乘得到最终的步长因子。其中，OX表示十六进制数，所述的公式(1)中的十进制小数与整数的转换是采用 b = a*(2"12)的方法计算得到。其中，假设a为原十进制数，通过计算公式得到一个整数b。再将b转换为十六进制得到。所述分频器结构图如图3所示，因为除法器1需要8个周期的时钟才能执行完一次除法操作，分频器主要为除法器1提供触发时钟，它将FPGA平台时钟进行8分频后，为除法器1提供触发时钟；该分频器具体是由计数器、比较器1、比较器2、反相器和D触发器构成完成对FPGA 平台时钟的分频工作；所述计数器的输出端与比较器1的输入端和比较器2的输入端连接， 所述比较器1的输出端与反相器的输入端连接，所述反相器的输出端与D触发器1的ce管脚连接，所述比较器2的输出端与D触发器1的d管脚连接，所述比较器1的输出端与计数器的clr管脚连接，D触发器1的输出端与除法器1连接；分频器的输入为FPGA平台时钟信号clk(本发明中的所有参与的模块，如寄存器， 加法器，除法器，减法器等的输入信号都有系统时钟信号作为触发时钟)，计数器对输入的时钟信号elk脉冲信号进行上升沿监测，每检测到一个elk的上升沿计数器的值加1，且将计数器的值分别输入比较器1和比较器2 ；比较器1为计数器提供清零操作指令；比较器1将计数器的值与7比较，当计数器的值等于7，比较器1输出1即高电平，与比较器1相连的计数器clr管脚接收到高电平时， 计数器执行清零操作；当计数器的值不等于7时，比较器1输出0即低电平，此时计数器进行加值操作；比较器1与反相器为D触发器1提供输出使能指令；比较器1将计数器的值与7 比较，当计数器的值等于7，则比较器1输出1即高电平，该高电平经反相器翻转为低电平后输入至D触发器1的ce管脚，D触发器1输出不变，仍然输出D触发器1上一时刻输出的信号；当计数器的值不等于7时，比较器1输出0即低电平，该低电平经反相器翻转为高电平后输入至D触发器1的ce管脚，D触发器1输出比较器2的电平；比较器2控制D触发器1输出电平翻转；比较器2将计数器的值与3进行比较，当计数器的值小于3时，比较器2输出低电平至D触发器1的d管脚，与此同时，D触发器1的 ce管脚接收到反相器输入的高电平，D触发器1输出为比较器2的低电平，该电平输入乘法器3为其提供触发时钟；当计数器的值大于3时，比较器2输出高电平至D触发器1的d管脚，与此同时，D触发器1的ce管脚接收到反相器输入的高电平，D触发器1输出为比较器 2的高电平，该电平输入乘法器3为其提供触发时钟；通过对elk的脉冲计数，当计数值小于3时将clk_8置为低电平，大于3时置为高电平，而在计数器等于7时，计数器清零。以此反复。clk_8是相当于将elk的时钟放大了 8倍来得到的。所述滤波器为基于FPGA的L阶LMS滤波器(自适应滤波器)，该滤波器由L个
6LMS_tap ( 一阶滤波单元)和L个加法器构成进行如下运算Wi (η) = ^(η-1)+2μ (n-l)e(n-l)x(n-i+l)Yi (η) = Wi (η) χ (n-i+1)y(n) = Y1 (η) +y2 (η) +y3 (η) —+Yi (η) — +yL (η)其中:L为》1的自然数，1 < i^L,y(n)为滤波器η时刻所计算的滤波后的参考噪声信号，Ii (η)为第i个LMS_tap LMS_tapi所计算的η时刻所计算的滤波后的参考噪声信号，Wi(Ii)为LMS_tapi η时刻的滤波权系数，初始时刻LMS_tapi的滤波权系数为0 ；所述L个 LMS_tap 分别为LMS_tapl、LMS_tap2、...、LMS_tap (i_l)、LMS_tapi、LMS_ tap (i+1)、…、LMS_tap (L-I)禾口 LMS_tapL ；所述 LMS_tapl、LMS_tap2、 ...、LMS_tap (i_l)、 LMS_tapi、…、LMS_tap (L-I)和LMS_tapL结构相同且依次连接；所述LMS_tapi由乘法器 4、加法器2、与门、D触发器2、D触发器3和乘法器5构成，所述乘法器4、加法器2、与门、D 触发器2和乘法器5依次连接；所述D触发器2与加法器2连接；所述LMS_tapi中的乘法器4、乘法器5和D触发器3均与LMS_tap (i_l)中的D触发器3连接，所述LMS_tap (i+1) 中的乘法器4、乘法器5和D触发器3均与LMS_tapi中的D触发器3连接；所述L个加法器分别为加法器3-1、加法器3-2、…、加法器3_(i_l)、加法器3_i、 加法器3-(i+l)、…、加法器3-(L-l)和加法器3-L ；该L个加法器依次连接，LMS_tapi中的乘法器5与加法器3-i连接；该L个加法器的初始化值均为0 ；LMS_tapl中的乘法器4与第一寄存器连接，LMS_tapl中的D触发器3与第一寄存器连接，LMS_tapl中的乘法器5与第一寄存器连接；L个LMS_tap中的乘法器4均与变步长因子计算电路中的乘法器3连接，L个LMS_tap中的乘法器4均与减法器连接，加法器3-L 与减法器连接。L个LMS_tap中的D触发器2是用来存储当前计算的滤波权系数以计算后续时刻的滤波权系数，本发明中初始时刻即第1时刻reset输入与门为0，其他时刻reset输入与门为1，当reset输入与门为0时，与门输出0于D触发器2中，每个LMS_tap的初始滤波权系数为0，当reset输入与门为1时，与门输出加法器2的计算结果滤波权系数于D触发器 2中；第一寄存器输出的(η-L+l)时刻参考噪声信号χ (n_L+l)接入LMS_tapl中的乘法器4中，第一寄存器的输出x(n-L+l)接入LMS_tapl中的乘法器5中，第一寄存器的输出 x(n-L+l)接入LMS_tapl中的D触发器3中时；变步长因子计算电路输出的(n-L)时刻变步长因子μ (n-L)分别接入L个LMS_ tap的乘法器4中，减法器的输出的(n-L)时刻误差信号e (n-L)分别接入L个LMS_tap的乘法器4中；LMS_tapl中的乘法器4完成2 μ (n-L) e (n-L) χ (n-L+1)运算，该运算结果输入加法器2中，D触发器2中缓存的W1(Ii-L)同时也输入至加法器2中，在加法器2中完成 (η-L+l)时刻 LMS_tapl 的滤波权系数 W1 (n-L+1) = W1 (n-L) +2 μ (n-L) e (n-L) χ (n-L+1)运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将W1 (η-L+l) 运算结果输入乘法器5中完成yi (η-L+l) = W1 (η-L+l) χ (η-L+l)运算，该运算结果输入至加法器3-1中，加法器3-1将yi (η-L+l)的计算值输入至加法器3_2中，该加法器3_1的初始值为0 ；输入至D触发器3中的x(n-L+l)经延迟输入LMS_tap2的乘法器4、乘法器5和D 触发器3中；
第一寄存器输出的(n-L+2)时刻参考噪声信号χ (n_L+2)分别接入LMS_tapl中的乘法器4、乘法器5和D触发器3中时，变步长因子计算电路输出的(η-L+l)时刻变步长因子μ (η-L+l)分别接入L个 LMS_tap的乘法器4中，减法器的输出(η-L+l)时刻的误差信号e (n_L+l)分别接入L个 LMS_tap的乘法器4中；LMS_tapl中的乘法器4完成2 μ (η-L+l) e (η-L+l) χ (n-L+2)运算，该运算结果输入加法器2中，D触发器2中缓存的W1 (η-L+l)同时也输入至加法器2中，在加法器2中完成(n-L+2)时刻 LMS_tapl 的滤波权系数 W1 (n_L+2) = W1 (η-L+l)+2 μ (η-L+l) e (η-L+l) χ (n-L+2)运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将 W1 (n-L+2)运算结果输入乘法器5中完成yi (n-L+2) = W1 (n-L+2) χ (n-L+2)运算，该运算结果输入至加法器3-1中，加法器3-1将该值输入加法器3-2中；与此同时，LMS_tap2中的乘法器4完成2 μ (η-L+l) e (η-L+l) χ (η-L+l)运算，该运算结果输入加法器2中，D触发器2中缓存的W2 (η-L+l)同时也输入至加法器2中，在加法器2中完成(n-L+2)时刻LMS_tap2的滤波权系数w2 (n-L+2) = W2 (η-L+l) +2 μ (η-L+l) e (η-L+l) χ (η-L+l)运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将W2 (n-L+2)运算结果输入乘法器5 中完成y2 (n-L+2) =w2 (n-L+2) χ (η-L+l)运算，该运算结果输入至加法器3_2中，在加法器 3-2中完成yi (n-L+2)+y2 (n-L+2)的计算，该计算值输入至加法器3_3中；输入至LMS_tap2D触发器3中的x(n-L+l)经延迟输入LMS_tap3的乘法器4、乘法器5和D触发器3中；输入至LMS_taplD触发器3中的x(n-L+2)经延迟输入LMS_tap2的乘法器4、乘法器5和D触发器3中；依此类推，第一寄存器的输出η时刻的参考噪声信号χ (η)分别接入LMS_tapl中的乘法器4、 乘法器5和D触发器3中时，变步长因子计算电路的输出(n-1)时刻的变步长因子μ (η-1) 分别接入L个LMS_tap的乘法器4中，减法器的输出(n_l)时刻的误差信号e (n_l)分别接入L个LMS_tap的乘法器4中；LMS_tapi中的乘法器4和乘法器5以及D触发器3均接收到参考噪声信号x(n_i+l)；LMS_tapl中的乘法器4完成2 μ (n-1) e (η-1) χ (η)运算，该运算结果输入加法器2 中完成W1 (n) =W1 (η-1)+2 μ (n-1) e (η-1) χ (η)运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将W1 (η)运算结果输入乘法器5中完成yi (η) = W1 (η) χ (η) 运算，该运算结果输入至加法器3-1中，在加法器3-1初始之为0的条件下完成yi(n)+0运算后，将运算所得结果输入加法器3-2 ；LMS_tapi中的乘法器4完成2 μ (n-1) e (η-1) χ (n_i+l)运算，该运算结果输入加法器2中，D触发器2中缓存的Wi (n-1)同时也输入至加法器2中，在加法器2中完成Wi (η) =Wi (η-1)+2 μ (n-l)e(n-l)x(n-i+l)运算，该运算结果经与门输入至D触发器2并存储于 D触发器2中，D触发器2将&(11)运算结果输入乘法器5中完成yi (n) = Wi (η) χ (n-i+1) 运算，该运算结果输入至加法器3-i中，在加法器3-i中完成yi (η) +y2 (η) +y3 (η)-+y, (η)的运算，并将该运算结果输入至加法器3-(i+1)中；0^_{& 1^中的乘法器4完成2 4 (n-1) e (η-1) χ (η-L+l)运算，该运算结果输入加法器2中，D触发器2中缓存的wjn-l)同时也输入至加法器2中，在加法器2中完成wjn)= wL (n-1) +2 μ (n-l)e(n-l)x (n-L+1)运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将& (η)运算结果输入乘法器5中完成九(n) = wL (η) χ (η-L+l)运算，该运算结果输入至加法器3-L中，在加法器3-L中完成y (n) = Yl (η) +y2 (η) +y3 (η)... (η)... +yL(n)的运算，并将该运算结果输入至减法器中；在减法器中完成e(n) = d(n)-y(n)的运算，e(n)经第三缓存器输出得到滤波后的车辆振动信号；e(n)输入至变步长因子计算电路中按照公式(1)完成μ (η)的计算；新输入的(η+1)时刻的含噪振动信号d(n+l)和参考噪声信号x(n+l)、滤波器中延迟的x(n)至x(n-L)、e(n)和μ (η)按照上述工作流程计算得到新的滤波后的车辆噪声振动信号。输入至滤波器中时，相应的滤波器中的LMS_tapl和L个加法器工作；乂2)输入至滤波器中时，相应的滤波器中的LMS_tapl、LMS_tap2和L个加法器工作；xC3)输入至滤波器中时，相应的滤波器中的LMS_tapl、LMS_tap2、LMS_tap3和L个加法器工作；依此类推，χ (L-I)输入至滤波器中时，相应的滤波器中的LMS_tapl至LMS_tap(L-l)和L个加法器工作；实施例1 本实施例中L = 3。第1时刻的含噪振动信号d(l)接入第二寄存器，第二寄存器输出d(l)至减法器，参考噪声信号X(I)接入第一寄存器，在第一寄存器缓存，第一寄存器输出X (1)接入滤波器中的LMS_tapl中的乘法器4、乘法器5和D触发器3中；LMS_tapl的 reset输入与门为0，与门输出0于D触发器2中，LMS_tapl的初始滤波权系数为0，该存储于D触发器2中；LMS_tapl的D触发器3中的χ⑴延迟输入至LMS_tap2中的乘法器4、乘法器5和D触发器3中；此过程中，d(l)和x(l)为装置的启动信号；第2时刻的含噪振动信号(1( 接入第二寄存器，第二寄存器输出(1( 至减法器， 参考噪声信号χ ( 接入第一寄存器，在第一寄存器缓存，第一寄存器输出xO)接入滤波器中的LMS_tapl中的乘法器4、乘法器5和D触发器3中；LMS_tapl中的乘法器4、乘法器5 和D触发器3中，减法器初始输出的误差信号e (1) = 0，变步长因子计算电路初始输出的 μ (1) = 0, e(l)和μ (1)输入3个LMS_tap中的乘法器4中；在11^_1即1的乘法器4中完成2 4 (1)Θ(1)Χ(1)，该运算结果输入加法器2中，D触发器2中缓存的初始滤波权系数同时也输入至加法器2中，在加法器2中完成2时刻LMS_ 1即1的滤波权系数《乂2) =0+2 μ (l)e(l)x(2) = 0运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将W1⑵运算结果输入乘法器5中完成yi⑵= W1 (2) χ (2) = 0运算，该运算结果输入至加法器3-1中，在加法器3-1初始值为0的条件下加法器3-1将yi (2)输入至加法器3-2中；与此同时，LMS_tap2的reset输入与门为0，与门输出0于D触发器2中，LMS_tap2 的初始滤波权系数为0，该存储于D触发器2中；LMS_tap2中的乘法器4完成2 μ (l)e(l) x(l)运算，该运算结果输入加法器2中，D触发器2中缓存的初始权系数同时也输入至加法器2中，在加法器2中完成0^_{即2的滤波权系数 0) = 0+2 μ (I)e(I)x(I) = 0运算， 该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将W2⑵运算结果输入乘法器5中完成y2 O) =W2 (2) x(l) =0运算，该运算结果输入至加法器3-2中，在加法器3-2中完成yQ) =Y1 (2)+y2 (2) = 0的计算，该计算值输入至加法器3_3中，加法器3-3将yQ)的值输入减法器中，在减法其中完成eO) =d(2)-y(2)的计算；输入至LMS_tap2D触发器3中的χ (1)经延迟输入LMS_tap3的乘法器4、乘法器5 和D触发器3中；输入至LMS_taplD触发器1中的经延迟输入LMS_tap2的乘法器4、 乘法器5和D触发器3中；第3时刻的含噪振动信号dC3)接入第二寄存器，第二寄存器输出dC3)至减法器， 参考噪声信号X C3)接入第一寄存器，在第一寄存器缓存，第一寄存器输出xC3)接入滤波器中的LMS_tapl中的乘法器4、乘法器5和D触发器3中；LMS_tapl中的乘法器4、乘法器5 和D触发器3中，e (2)输入3个LMS_tap中的乘法器4中，e(2)输入至变步长因子计算电路中完成μ (2)的计算，得到的μ (2)值输入3个LMS_tap中的乘法器4中；在0^_1即1的乘法器4中完成2口 (2) e (2) χ (3)，该运算结果输入加法器2中，D 触发器2中缓存的W1 O)同时也输入至加法器2中，在加法器2中完成2时刻LMS_tapl的滤波权系数W1 (3) = W1 (2) +2 μ (2) e (2) χ (3)运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将W1 (3)运算结果输入乘法器5中完成yi (3) = W1 (3) χ (3) 运算，该运算结果输入至加法器3-1中，在加法器3-1初始值为0的条件下加法器3-1将 Y1 (3)输入至加法器3-2中；在11^_1即2的乘法器4中完成211 O) e O) χ O)，该运算结果输入加法器2中，D 触发器2中缓存的W2 O)同时也输入至加法器2中，在加法器2中完成2时刻LMS_tapl的滤波权系数W2 (3) =W2 (2)+2 μ (2)e(2)x(2)运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将W2 (3)运算结果输入乘法器5中完成y2 (3) = W2 (3) χ⑵ 运算，该运算结果输入至加法器3-2中，在加法器3-2中完成yi (3) +y2 (3)的运算，并将运算结果输入至加法器3-3中；与此同时，LMS_tap3的reset输入与门为0，与门输出0于D触发器2中，LMS_tap3 的初始滤波权系数为0，该存储于D触发器2中；LMS_tap3中的乘法器4完成2 μ (2) e (2) x(l)运算，该运算结果输入加法器2中，D触发器2中缓存的初始滤波权系数同时也输入至加法器2中，在加法器2中完成LMS_tap3的滤波权系数W3(3) = 0+2 μ (I)e(I)x(I)运算， 该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将W3 (3)运算结果输入乘法器5中完成y3 (3) =W3 (3) x(l)运算，该运算结果输入至加法器3_3中，在加法器 3-3中完成yC3) =Y1 (3)+y2 (3)+y3 (3)的计算，该计算值输入至加法器3_3中，加法器3_3 将7( 的值输入减法器中，在减法器中完成eC3) =d(3)-y(3)的计算；eC3)经第三寄存器输出得到滤波后的车辆振动信号；LMS_tap2D触发器3中的χ (2)经延迟输入LMS_tap3的乘法器4、乘法器5和D触发器3中，信号χ⑵将LMS_tap3中的乘法器4、乘法器5D触发器3中的χ⑴信号替换； 输入至LMS_taplD触发器3中的χ (3)经延迟输入LMS_tap2的乘法器4、乘法器5和D触发器3中；第4时刻的含噪振动信号d (4)接入第二寄存器，第二寄存器输出d (4)至减法器， 参考噪声信号χ (4)接入第一寄存器，在第一寄存器缓存，第一寄存器输出χ (4)接入滤波器中的LMS_tapl中的乘法器4、乘法器5和D触发器3中；LMS_tapl中的乘法器4、乘法器5 和D触发器3中，e (3)输入3个LMS_tap中的乘法器4中，e(3)输入至变步长因子计算电路中完成μ (3)的计算，得到的μ (3)值输入3个LMS_tap中的乘法器4中；
在LMS_tapl的乘法器4中完成2 μ (3) e (3) χ ，该运算结果输入加法器2中，D 触发器2中缓存的W1 (3)同时也输入至加法器2中，在加法器2中完成2时刻LMS_tapl的滤波权系数W1 (4) = W1 (3) +2 μ (3) e (3) χ (4)运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将W1 (4)运算结果输入乘法器5中完成yi (4) = W1 (4) χ (4) 运算，该运算结果输入至加法器3-1中，在加法器3-1初始值为0的条件下加法器3-1将 Y1 (4)输入至加法器3-2中；在0^_1即2的乘法器4中完成211 (3) e (3) χ (3)，该运算结果输入加法器2中，D 触发器2中缓存的W2 (3)同时也输入至加法器2中，在加法器2中完成2时刻LMS_tapl的滤波权系数W2 G) =W2 (3)+2 μ (3)e(3)x(3)运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将W2 (4)运算结果输入乘法器5中完成y2 (4) = W2 (4) χ (3) 运算，该运算结果输入至加法器3-2中，在加法器3-2中完成yi (4) +y2 (4)的运算，并将运算结果输入至加法器3-3中；在LMS_tap3的乘法器4中完成2 μ (3) e (3) χ (2)，该运算结果输入加法器2中，D触发器2中缓存的W3 (3)同时也输入至加法器2中，在加法器2中完成2时刻LMS_tapl的滤波权系数％(4) =W3(3)+2 μ (3)e(3)x(2)运算，该运算结果经与门输入至D触发器2并存储于D触发器2中，D触发器2将％(4)运算结果输入乘法器5中完成=W3(4)χ(2)运算，该运算结果输入至加法器3-3中，在加法器3-3中完成y (4) = Y1 (4)+y2 (4)+y3 (4)的运算，并将运算结果输入至减法器中；在减法器中完成e(4) = d(4)-y(4)的计算；e(4)经第三寄存器输出得到滤波后的车辆振动信号。e(4)输入至变步长因子计算电路中完成μ (4) 的计算；新输入的5时刻的含噪振动信号d 和参考噪声信号χ(5)、滤波器中延迟的χ(5) 至x(3)、e(4)和μ (4)按照上述工作流程计算得到新的滤波后的车辆噪声振动信号。实施例2 本实施例中L = 128。本实施例中的装置的结构和工作原理与实施例1中的相同。
1权利要求
1.一种基于FPGA的车辆振动信号滤波装置，其特征在于该装置包括第一寄存器、第二寄存器、第三寄存器、减法器、变步长因子计算电路和滤波器；参考噪声信号x(n)接入第一寄存器的输入端缓存，含噪振动信号d(n)接入第二寄存器的输入端缓存，第一寄存器的将缓存后的参考噪声信号x(n)输出至滤波器的一个输入端，减法器输出的误差信号 e(n-l)分别接入第三寄存器的输入端、滤波器的第二输入端和变步长因子计算电路的输入端，变步长因子计算电路的输出信号变步长因子μ (η-1)接入滤波器的第三输入端，第二寄存器将缓存后的含噪振动信号d(n)输入至减法器一个输入端，滤波器输出的滤波后的参考噪声信号y(n)接入减法器的另一输入端，第三寄存器输出滤波处理后的车辆振动信号；其中d(n)为η时刻的含噪振动信号，χ (η)为η时刻的参考噪声信号，e (η_1)为(η_1) 时刻的误差信号，μ (η-1)为(η-1)时刻滤波器的变步长因子，初始时刻减法器与变步长因子计算电路的输出均为0，即e(l) =0, μ (1) = 0；所述变步长因子计算电路由乘法器1、乘法器2、乘法器3、加法器1、减法器1、除法器1 和分频器构成实现变步长因子μ (η-1)的计算； 该μ (η-1)的计算公式如下μ(η-\) = 0.08(1 -1)IOOxe ( -1) + 1所述分频器由计数器、比较器1、比较器2、反相器和D触发器1构成，该分频器对FPGA 平台时钟进行8分频后，为除法器1提供触发时钟；所述计数器的输出端与比较器1的输入端和比较器2的输入端连接，所述比较器1的输出端与反相器的输入端连接，所述反相器的输出端与D触发器1的ce管脚连接，所述比较器2的输出端与D触发器1的d管脚连接， 所述比较器1的输出端与计数器的clr管脚连接，D触发器1的输出端与除法器1连接；所述滤波器为基于FPGA的L阶LMS滤波器，该滤波器由L个LMS_tap和L个加法器构成完成如下运算Wi(Ji) = Wi (η-1)+2 μ (n-1) e (η-1) χ (n-i+1)Yi (η) = Wi (η) χ (n-i+1) y (η) = Y1 (η) +y2 (η) +y3 (η)…+Yi (η)…+yL (η)其中L为> 1的自然数，为滤波器所计算的η时刻滤波后的参考噪声信号，Ii (η)为第i个LMS_tap即LMS_tapi所计算的η时刻滤波后的参考噪声信号，Wi (η) 为LMS_tapi η时刻的滤波权系数，初始时刻LMS_tapi的滤波权系数为0。
2.如权利要求1所述的基于FPGA的车辆振动信号滤波装置，其特征在于，所述滤波器中的LMS_tap个数或加法器个数L = 128。
全文摘要
本发明公开了一种基于FPGA的车辆振动信号滤波装置。该装置中参考噪声信号x(n)接入第一寄存器的输入端，含噪振动信号d(n)接入第二寄存器的输入端，第一寄存器的输出信号x(n)接入滤波器的一个输入端，减法器输出的误差信号e(n-1)分别接入第三寄存器的输入端、滤波器的第二输入端和变步长因子计算电路的输入端，变步长因子计算电路的输出信号变步长因子μ(n-1)接入滤波器的第三输入端，第二寄存器的输出信号d(n)接入减法器一个输入端，滤波器输出的滤波后的参考噪声信号y(n)接入减法器的另一输入端，第三寄存器输出滤波处理后的车辆振动信号。装置所依据的算法具有计算速度快、收敛速度快、实时性高、滤波效果好等优点，可以满足信号实时处理的需求。
文档编号H03H17/02GK102394591SQ20111025900
公开日2012年3月28日 申请日期2011年9月2日 优先权日2011年9月2日
发明者余腾, 刘占文, 史昕, 惠飞, 杨澜, 赵祥模, 雷涛 申请人:长安大学