《第二讲最速下降法PPT讲稿.ppt》由会员分享,可在线阅读,更多相关《第二讲最速下降法PPT讲稿.ppt(28页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、第二讲最速下降法第1页,共28页,编辑于2022年,星期二n最速下降法(method of steepest descent)是一种基于梯度的自适应方法。n 最速下降法可用反馈系统来表示,滤波器的计算式一步一步迭代进行的。从该意义上讲,最速下降法是递归的。n在适当条件下,最速下降法的解收敛于维纳解而不需要求输入向量相关矩阵的逆矩阵。第2页,共28页,编辑于2022年,星期二线性最优滤波:问题综述第3页,共28页,编辑于2022年,星期二n这里滤波器需要两个约束条件:1.滤波器是线性的2.滤波器是离散时间的n滤波器的具体实现依赖的两个选择:1.滤波器的冲激响应选择(FIR,IIR)2.统计优化准
2、则的选择问题 1)估计误差的均方值 2)估计误差的绝对值期望值 3)估计误差的绝对值的三阶或高阶期望值第4页,共28页,编辑于2022年,星期二n滤波器问题的本质表示如下:给定一个输入取样序列u(0),u(1),u(2),设计一个线性离散滤波器其输出y(n)提供了期望响应d(n)的一个估值,使得其估计误差的均方值e(n)定义为期望响应d(n)与实际响应y(n)之差为最小。n两种数学解决方案1 正交性原理2误差性能曲面第5页,共28页,编辑于2022年,星期二正交性原理nn时刻滤波器输出为线性卷积n误差 n代价函数均方误差 第6页,共28页,编辑于2022年,星期二 使代价函数J获得最小值的充要
3、条件是其对应的估计误差 e(n)于n时刻进入期望响应估计的每个输入样值。第7页,共28页,编辑于2022年,星期二第8页,共28页,编辑于2022年,星期二2.1 最速下降法的基本思想无约束最优化的数学表示如下:其中 是一个代价函数,是个未知向量 是要寻找的最优解。第9页,共28页,编辑于2022年,星期二局部迭代下降思想 首先假设一个初始权向量 ,然后产生一系列权向量 能够使代价函数 在算法的每次迭代都是下降的,也就是满足如下表达式 最速下降法其实就是一种简单形式的迭代下降,它主要思想是沿着最速下降方向连续不断调整权向量 。最速下降方向也就是负梯度方向 梯度向量表示如下 第10页,共28页,
4、编辑于2022年,星期二通过以上可得最速下降算法 其中n表示迭代进程,是步长参数,是正常数。在从n到n+1的迭代过程中,权向量的调整量为 第11页,共28页,编辑于2022年,星期二第12页,共28页,编辑于2022年,星期二 假设w为复值向量,那么梯度向量g也是复值向量。所以使用共轭转置(埃尔米特转置)因此上式可变为 从上式可以看出当 为正数时,因此,随着n的增加,代价函数减小,当 时,代价函数趋于最小值 。第13页,共28页,编辑于2022年,星期二2.2 最速下降算法应用于维纳滤波器图2.1 自适应横向滤波器的结构第14页,共28页,编辑于2022年,星期二 通过比较期望响应 及其估计值
5、,可以得到一个估计误差 即 其中 是抽头权向量 与抽头输入向量 的内积 如果抽头输入向量u(n)和期望响应d(n)是联合平稳的,则此时均方误差或者在n时刻的代价函数J(n)是抽头全向量的二次函数。第15页,共28页,编辑于2022年,星期二横向滤波器的代价函数为所以展开可得其中,是目标函数 的方差 P=抽头输入向量 与期望响应 的互相关向量 R=抽头输入向量 的相关矩阵第16页,共28页,编辑于2022年,星期二同时梯度向量可写为第17页,共28页,编辑于2022年,星期二因此维纳滤波中最速下降法的数学表达式为:从另一个角度,可以将上公式看做一个反馈模型,信号流图如下图2.2 最速下降算法的信
6、号流图表示第18页,共28页,编辑于2022年,星期二2.3 最速下降法的稳定性影响该算法的稳定性有两个因素:(1)步长参数(2)抽头输入向量 的相关矩阵R首先定义n时刻的加权误差向量 其中 是抽头权向量的最优值第19页,共28页,编辑于2022年,星期二使用特征值分解可得将R代入上公式可得两边同时左乘令v(n)的初始值为:第20页,共28页,编辑于2022年,星期二对于最速下降法的第k个自然模式,并初始化可以得到 为了满足最速下降法的稳定性或收敛性,对于所有k,我们可以有 因此最速下降法稳定性的充分必要条件是步长因子满足不等式第21页,共28页,编辑于2022年,星期二从图中可以看出,当迭代
7、次数趋近于无穷时,趋近于0也就是抽头加权向量 逼近最优解图2.3 最速下降算法的第k个自然模式随时间变化的情况第22页,共28页,编辑于2022年,星期二由上图我们可以定义一个时间常数 使得 表示了 衰减到初始值 的 时所需要的迭代次数初始抽头加权向量 的瞬态特性 两边同时左乘 第23页,共28页,编辑于2022年,星期二因此第i个抽头权值的瞬态特性可以表示为其中 是第i个抽头权值的最优值,是第k个特征向量的第i个分量上式表明,最速下降算法中每一个抽头权值收敛于指数形式 的加权和。同时定义整个时间常数 则可得任意抽头权值的时间常数的上下界定义如下第24页,共28页,编辑于2022年,星期二均方
8、误差的瞬态特性均方误差的瞬态特性可知误差性能曲面的规范形式其中 是最小均方误差从初始值到最终值 的指数衰减的时间常数为当 较小时第25页,共28页,编辑于2022年,星期二2.4作为确定性搜索法的最速下降算法n最速下降算法提供了从任意初始点出发寻找误差性能曲面极小点的局部搜索方法。n最速下降算法的运行,取决于三个量:*起始点:由抽头权向量初始值w(0)规定 *梯度向量:位于误差性能曲面的特殊点,由互相关 向量P和相关向量矩阵R唯一确定 *步长参数 :控制横向滤波器抽头权向量从算法 的某一次迭代到下一次迭代的增量变化n一旦规定了这三个量,最速下降算法将沿着多维权值空间独特的路径前进,它从初始点w
9、(0)出发,终止于最优解 .换句话说,在权值空间中最速下降算法是一种确定性的搜索方法。第26页,共28页,编辑于2022年,星期二2.5 最速下降法的优点与局限性优点:简单性,只要给出起始点,梯度向量及步长参数,将沿着权值空间特殊的路径前进,从初始点出发,终止于最优解。也就是说它是一种确定性搜索方法。局限性:该过程需要大量的迭代,主要原因就是以围绕当前点的误差性能曲面的线性(一阶)为基础。针对这点后来提出来牛顿法,它是围绕当前点记为w(n)进行误差性能曲面的二次(例如二阶)逼近。第27页,共28页,编辑于2022年,星期二本章小结n两种集平均量 R(抽头输入向量的相关矩阵)(抽头输入向量的相关矩阵)P(抽头输入向量和期望响应的互相关向量)(抽头输入向量和期望响应的互相关向量)n最速下降算法提供了计算维纳滤波器抽头权向量的简化步骤。最速下降算法的一个重要特点就是存在反馈,实质是该算法是递归的。另外,我们需特别注意算法的稳定性问题。而稳定性受制于算法反馈环中的两个参数:步长大小参数步长大小参数 抽头输入向量的相关矩阵抽头输入向量的相关矩阵Rn特别地,算法稳定性的充要条件具体化为:n此外,依赖于步长参数 的值,最速下降算法的瞬态响应特性呈现如下三种形式之一:欠阻尼响应,过阻尼响应,临界阻尼响应 第28页,共28页,编辑于2022年,星期二
限制150内