发明2、一种信噪分离的优化方法11.29解析.docVIP

本发明涉及一种信噪分离的优化方法，其具体步骤如下：（1）读取离散化后的输入信号；（2）求信号x(n)的频谱；（3）确定FIR滤波器的目标频率响应；（4）在频率子集上等间隔地取r+2个频率点,,…,作为交错点组；（5）依据重心形式的拉格朗日插值；（6）求；（7）求输入信号；（8）由傅立叶逆变换求；（9）输入信号进行小波分解；（10）最高层小波细节系数处理；（11）其他层细节系数处理；（12）最高层尺度系数处理；（13）整形输出。本发明采用FIR优化滤波，具有线性相位特性，相同阶次下，通带和阻带最大误差达到最小；采用新的小波滤波算法进行滤波，提高了信噪比。   1、一种信噪分离的优化方法，其具体步骤如下： （1）读取离散化后的输入信号:x(n)j=1,2……N，其中N为信号x(t)的采样点数； （2）求信号x(n)的频谱； （3）确定FIR滤波器的目标频率响应，误差加权函数为,则加权逼近函数为： ； 等纹波最佳逼近法设计线性相位FIR滤波器实质是求解(N-1)/2个系数,使得加权误差函数的最大绝对值在要求逼近的频带内达到最小，即： 式中，A为要求逼近的内的一个闭区间，该闭区间包括通、阻带，但不包含过渡带； 切比雪夫逼近理论指出上述存在且唯一，并依据交错点组定理，给出求解的公式，如下所示： ； 式中：是波纹极值；求解式，可得r+1个未知数，即可得到和； （4）在频率子集上等间隔地取r+2个频率点,,…,作为交错点组，初始值代入下式： ； 式中：是相对于初始交错点组产生偏差，非最优； （5）依据重心形式的拉格朗日插值公式得： ； 式中：； ,i=0，1，2，…，r+1； （6）由式(11)（注：这是哪一个？）求； （7）求输入信号，经过滤波器的输出信号的频谱： （8）由傅立叶逆变换求； （9）输入信号进行小波分解， 小波尺度系数为：； 小波细节系数：； （10）最高层小波细节系数处理； （11）其他层细节系数处理： 小波细节系数代表着信号的细节信息，对各层小波系数进行滤波、重构可以比较理想地保留信号的细节信息；主要根据极大值传播定理，设信号s(k)第j层细节系数的某一峰值点位于,第j+1层细节系数与此相对应的峰值点位于，则满足，可以进行其他各层小波细节系数的重构，具体为： 根据式(5)，由最高层小波细节系数峰值位置推出其他各层小波细节系数峰值位置； 分别求各层小波系数所有峰值点幅度的均值，并修正相应峰值的幅度； 根据式(11)，通过峰值点的位置和幅度重构各层小波细节系数； （12）最高层尺度系数处理：最高层尺度系数包含着信号的轮廓信息，由推论可知同一层细节系数的峰值点对应于该层尺度系数的零点，根据前面的道德最高层尺度系数进行分段、平滑，从而实现滤波：假设前后两个峰值点位于,且那么对区间内的尺度系数按照式进行重构； （13）整形输出：由于最高层小波细节系数峰值寻找可能出现的微小的错位以及各个系数存在重构误差，由上述三个步骤得到的小波重构系数进行小波逆变换得到的时域波形不够理想，所以需要对滤波输出y(n)进行整形，具体方法为： 1)求出滤波输出信号的绝对值均值； 2)根据第一层小波细节系数峰值位置将滤波输出信号进行分段，并求解各段均值； 3)由各段均值符号判定该段滤波输出信号的符号，由绝对值均值决定其幅度。 2、根据权利要求1所述的信噪分离的优化方法，其特征在于：所述步骤（6）中的具体计算步骤如下： A1：若所有均满足,则是波纹极值，且初始猜测值才为交错点频率，计算结束，否则，执行步骤A2； A2:对Step1中r+2个点交错点，检查其附近是否存频率点满足,若存在，则在该点附近找出局部极值点，代替原来点，得一组新的交错频率，重新计算，和； A3:重复A2，调整交错点组中的频率。 3、根据权利要求1所述的信噪分离的优化方法，其特征在于：所述步骤（10）中的最高层小波细节系数处理具体步骤为： 设此时小波细节系数为，lev为分解层数；由于此时信号小波细节系数占主导，首先对进行剔除，保留小波细节洗漱幅值大于的部分，置零其余小波细节系数；假设提出后的小波细节系数仍以表示，可通过以下步骤完成峰值点位置的搜索： B1：寻找第一个细节系数非零点，记； B2：在区间内，进行式(12) （注：这是哪一个？）操作，得到峰值点位置； ； B3：，转至步骤B1，直至遍历完成； 由于噪声等的影响，式(12)的判断可能会出现微小错位，通过增加分解层数在一定层度上避免，根据推论，在得到最高层细节系数峰值幅度及位置后可通过以下过程实现最高层细节系数重构：求最高层细节系数所有峰值点幅度的均值，并修正所有峰值的幅度；根据式，通过峰值点的位置和幅度重构最高层细节系数。  一种信噪分离的优化方法 技术领域 本发明涉及一种信噪分离的优化方法，属于信号处理技术领域。 背