自适应信号处理作业

本文格式为Word版，下载可任意编辑——自适应信号处理作业1.自适应滤波如何运用到系统辨识？

自适应滤波理论和技术是统计信号处理和非平稳随机信号处理的主要内容,它具有维纳滤波和卡尔曼滤波的最正确滤波性能,但不需要先验知识的初始条件,它是通过自学习来适应外部自然环境,因而具有广泛的应用。自适应滤波器(Adaptivefilter)是自设计的,由于其依靠递归算法进行运算,因此可在有关信号特征的完整知识不能得到的环境下,圆满的完成滤波运算。由于稳定性问题和IIR局部最优,所以,自适应滤波器大多用FIR来实现。在自适应滤波器应用中一个重要问题是使可调理滤波器参数最优的标准,以及利用这种标准形成实际上可行的算法。最小均方(LMS,leastmean-square)算法是现今应用最为广泛的一种线性自适应算法,它不需要有关的相关函数和矩阵求逆运算,是一种极为简单的算法.

最小均方误差(LMS，leastmeansquare)算法于1960年提出后,因其具有计算量小、易于实现等优点而获得大量应用。典型的应用领域有系统辨识、信号处理和自适应控制等。LMS算法的基本原理是基于估计梯度的最速下降法,即沿着权值的梯度估值的负方向进行探寻,以期达到权值最优,实现均方误差最小意义下的自适应滤波。

系统辨识是根据系统的输入输出时间函数来确定描述系统行为的数学模型，是现代控制理论中的一个分支。对系统进行分析的主要问题是根据输入时间函数和系统的特性来确定输出信号。系统辨识包括两个方面：结构辨识和参数估计。在实际的辨识过程中，随着使用的方法不同，结构辨识和参数估计这两个方面并不是截然分开的，而是可以交织在一起进行的。2.在系统辨识中，LMS，RLS算法的形式？2.1LMS原理：

设通信系统输出信号为：

y(k)=WTX(k)(1)

其中，该系统权向量为：

W=[w1,w2,…,wn]T(2)

输入信号为

X(k)=[x(k),x(k-1),…,x(k-m+1)]T(3)

误差信号定义为

e(k)=d(k)-y(k)=d(k)-WT(k)X(k)(4)

LMS算法的原理是用e2(k)来估计E(e2(k))，此时有

▽(k)=-2e*(k)X(k)(5)

这样梯度法的叠代公式变为

W(k+1)=W(k)+2μe*(k)X(k)(6)

其中，*为共扼。算法步骤：

基本的LMS算法如下：

步骤1初始化：W(0)=0,0<μ<1/λmax步骤2W(k+1)=W(k)+2μe(k)X*(k)

步骤3判断是否收敛，假如不收敛，令k=k+1，回步骤2。2.2RLS原理

SISO系统动态过程的数学模型：

A(z)z(k)?B(z)u(k)?n(k)

?1?1（1）

?...?anaz?na?nb其中u(k)，z(k)为输入输出量，n(k)为噪声。式中

A(z?1)?1?a1z)?b1z?1?a2z?2B(z?1?1?b2z?2?...?bnbz展开后得到：

z(k)??a1z(k?1)?a2z(k?2)?...?ana?b1u(k?1)?b2u(k?2)?...?bnbu(k?nb)模型（1）可化为最小二乘格式：

z(k)?h(k?)??n(k)

a（2）

b记??[a1,a2,...an,b1,b2,...,bn]?

为待估计的参数。h(k)?[?z(k?1),...,?z(k?na),u(k?1),...,u(k?nb)]?，对于。方程（2）构成一个线性方程组，写成k?1,2,...L（L为数据长度）

??h(1)??z(1)??n(1)?zL(k)?HL(k)??nL(k)；???????z(2)n(2)h(2)???，H??，nL??ZL??L???????????1???根据最小二乘法一次完成算法，其参数估计为：??????LS?(HLHL)HLZL???z(L)??n(L)??h(L)??。

参数递推估计，每取得一次新的观测数据后，就在前次估计结果的基础上，

利用新引入的观测数据对前次估计的结果，根据递推算法进行修正，减少估计误差，从而递推地得出新的参数估计值。这样，随着新观测数据的逐次引入，一次接一次地进行参数估计，直到参数估计值达到满意的准确程度为止。算法步骤：

步骤一：初始化W(0)?0；P(0)???1I，其中I为单位矩阵；步骤二：更新n?1,2,...计算

更新增益矢量：g(n)?P(n?1)X(n)/[??XT(n)P(n?1)X(n)]；滤波：y(n)?WT(n?1)X(n)；误差估计：e(n)?d(n)?y(n)；

更新权向量：W(n)?W(n?1)?g(n)e(n)；

更新逆矩阵：P(n)???1[P(n?1)?g(n)XT(n)P(n?1)]；

其中，P(n)为自相关矩阵Pxx(n)的逆矩阵，常数?是遗忘因子，且0???1。

3.最优意义下的权向量是什么？

LMS：W(k+1)=W(k)+2μe*(k)X(k)RLS：W(n)?W(n?1)?g(n)e(n)

4.采用LMS，RLS算法后，仿真结果，分析步长选择对整体性能的影响？LMS：4.1输入信号采用在区间（0，1）为均匀分布的随机信号。

图1给出的是采用LMS算法，当μ＝0.01，叠代1000次后，得到系统滤波器参数的收敛曲线图。从图中我们可以看到，三个系数基本在叠代到第200次收敛到真实值，可见算法在这种状况下是有效的。

图1μ＝0.01时的W的收敛曲线

图2给出了，当μ为0.01,0.03,0.3和0.5时，采用LMS算法叠代300次后得到的滤波器系数收敛曲线图，在图中，我们可以看到，（b）中的μ值大于（a）中的，其收敛速度明显加快。而（c）中的μ值过大而导致收敛曲线产生了振荡现象，但是最终依旧收敛。（d）中μ的取值由于超过了1/λmax，所以导致了曲线的发散。因此，合理地对μ的取值将影响LMS算法性能的优劣。

图2不同μ下的w3收敛曲线图

4.2输入信号采用正弦信号。

图17给出的是输入为正弦信号时，采用LMS算法，叠代1000次后，w3的收敛曲线，其步长分别取0.01,0.3和0.4。当取μ为0.01时，曲线收敛，但是它并未收敛到真实值，当μ为0.3时，曲线在0.5附近作间歇性振荡现象，而当μ为0.4时，曲线时而发散时而收敛。得到以上结果的主要原因是，正弦信号为单一频率的信号，而LMS算法要求输入信号的频率应当尽可能丰富，它对频率简单的信号是无法估计出滤波器系数的。

图17输入为正弦信号的收敛曲线

图1直接型自适应FIR滤波器

MATLAB实现

MATLAB有专门的函数ADAPTLMS实现自适应滤波，用法为:[y,e,s]=ADAPTLMS(x,d,s)

其中y为滤波器的输出信号,x为滤波器的输入信号,d为期望信号,e为误差信号,s为包含自适应滤波器信息的结构体,可用函数INITLMS初始化。它的语法为:s=INITLMS(h0,Mu)

其中h0为滤波器系数的初始值,是一个长度为滤波器阶数的向量可为0;Mu为步长参数。

自适应滤波器的应用实例：系统辨识

现要识别一个未知的线性系统,这个未知系统可以是一个全零点系统,也可以是一个零极点系统。该未知系统可用一个长度为N的FIR滤波器近似(建模,见图2)。为了构造这个问题,未知系统FIR模型并联级联,并由同一个长度输入序列x(n)鼓舞。假如{y(n)}代表模型的输出,而{d(n)}代表未知系统的输出,误差序列e(n)=d(n)-y(n)。假如能将平方误差和减到最小,就得到一组线性方程组,因此,可用LMS算法来对FIR模型的自适应,使得它的输出近似为该未知系统的输出。仿真源程序和仿真结果如下:

图2系统辨识方框图

x=0.1*randn(1,500);%输入信号x(n)

b=fir1(31,0.5);%构建待辨识的未知FIR系统d=filter(b,1,x);%期望信号d(n)

w0=zeros(1,32);%自适应滤波器系数的初始化mu=0.8;%LMS步长参数

S=initlms(w0,mu);%初始化s信息结构体[y,e,S]=adaptlms(x,d,S);%自适应滤波

stem([b.’,S.coeffs.’]);%画出识别的系统和未知系统的参数。

图3系统辩识的仿真结果

由图3可见,两个系统的系数十分的接近。

LMS算法的收敛性能由于主输入端不可避免地存在干扰噪声,L