Inspire

姿态仪作为水下机器人系统中一个重要传感器，它的重要性不言而喻，它是水下机器人实现自稳控制的关键传感器。目前在低成本的水下机器人系统中常常采用低成本的mems姿态仪，它主要是通过融合加速度计、角速度计和磁力计来实现对三维空间的姿态信息测量。它的优点是体积小、低成本，但是当水下机器人进行高速机动时，其自身的加速度变化非常大，导致mems姿态仪的姿态输出不稳定和真值相差较大。为了解决这个问题，本文提出了一种利用多个深度计来测量机器人的横滚角和俯仰角的方法，该方法可以不受机器人机动的影响，从而实现稳定...

早期的SLAM假设真实的世界可以被理所应当的建模成一个由简单离散的路标构成的集合，这些路标通过一些几何特征，例如点、线或者圆圈来表示。在更加复杂且没有结构的环境里，室外的，地下的，海里的，这一假设通常不能够成立。...

数据关联一直是SLAM实际应用中一个非常重要的问题。

同步定位与构图是一个移动机器人建立环境地图的同时利用建立的地图来计算自己位置的过程。

概率SLAM问题的解决方案是设计找到一个合适的对观察模型(2)和运动模型(3)的陈述，使得能够高效连续的计算式(4)和(5)的先验和后验概率分布。到目前为止，最常见的表述是在状态空间模型中的一种形式，其带有加性高斯噪声，进而导致了采用扩展卡尔曼滤波器(EKF)来解决SLAM问题。另外一种重要的表述是将(3)中的机器人运动模型的描述成一系列样本点，这些样本点更具有更一般性，是非高斯概率分布。这导致了使用Rao-Blackwellized粒子滤波器，或者叫FastSLAM算法，来解决SLAM问题...

SLAM是移动机器人建立环境的地图同时演绎出自己位置的过程。在SLAM中，平台的轨迹和所有地标的位置都是在线估计的，不需要任何位置的先验知识。准备工作 考虑移动机器人运动经过一个环境使用安装在机器人上面传感器对大量未知的路标进行相对观测，如图1所示。在时刻k，定义如下变量：：描述机器人位置和姿态的状态向量。：控制向量，被用于在k-1时刻驱动机器人到k时刻的状态向量：描述第i个路标位置的向量，它们的真实位置假定为静止的。：在k时刻机器人观测到的第i个路标...

SLAM的英文全称为Simultaneous Localization and Mapping，中文意思是同步定位与构图。主要是由Hugh Durrant-Whyte和Tim Bailey这两位大神教授完成的杰作。SLAM问题回答了对于一个被放置在未知位置的未知环境中的移动机器人来递增的构建周围环境的连续性地图并同时估计自身在地图中的位置是否为可行的。SLAM问题的解决方案被视为移动机器人领域的“圣杯”，因为它实现了机器人真正意义上的自主性。 在过去的十年里，SLAM问题...

在一些需要用单片机、DSP、ARM等系统中用C语言实现高维滤波器的场合，例如卡尔曼滤波器，我们会经常遇到用C语言求解方阵的逆的情况，当阶数为2阶或者3阶的时候我们可以用公式法直接求解，但是当阶数一旦达到4阶或者4阶以上时利用公式求解将会非常的麻烦，而且极易出错，这时我们需要寻求一种可以求解任意阶数的算法。由线性代数的知识可以知道求任意阶数矩阵逆矩阵有2种算法，一种是初等行（列）变换，一种是伴随矩阵法。其中初等行（列）变换由于需要进行行（列）的加减或交换，以及乘上（除以）非零的数这3种步骤，且顺...

算法描述 对于一个N阶的方阵，则它的行列式的值为，这里我们引入余子式的概念，在线性代数的课本中它的定义如下：在N阶行列式中，把元素所在的第i行和第j列划去后，留下了的N-1阶行列式叫元素的余子式，记作。则方阵A的行列式可以表示为。对于式子里头的，它是一个N-1阶的行列式，我们可以按照上述的方法进一步的分解，得到N-2阶的行列式，然后再分解得到N-3阶行列式，一直分解到2阶的行列式，这时2阶的行列式就变的非常的好求了，即，从这里我们就可以看到这是一个递归求解的过程，因此可以非常方便的用C语言...

水下机器人在水下想要获得定位信息目前主要是基于声学方式的长基线(LBL)、短基线(SBL)、超短基线(USBL)获得局部定位数据，还有通过地磁场或者重力场匹配导航的方式，这类方式目前仅限于军事领域。当然还有最近热门的水下SLAM定位导航的方式，需要借助图像声呐或者光学摄像头获取周围环境的特征进行同时建图与导航，这种方式不需要借助外部的信息即可实现定位，是目前的研究热点，但此类方式需要进行大量的计算，同时要求周围的水域有明显的特征（也称为landmark）才能获取更加准确的定位结果，技术难度较...

首先关于Cyclone IV系列的配置和远程系统更新需要阅读《Configuration and Remote SystemUpgrades in Cyclone IV Devices》这个手册，该手册在官网上可以找到。阅读这个手册可以了解到有关Cyclone IV的配置方式，如并行方式PS方式，串行方式AS这两种，以及主被动配置。还有就是配置的整个流程。这里笔者研究的是AS方式的主动配置，即通PC把配置文件通过串口发送到FPGA上然后FPGA通过一段小逻辑将接收到的文件写入到EPCS存储器里...

描述超声换能器的性能指标有工作频率、机电耦合系数、机电转换系数、品质因数、方向特性、发送功率、效率、灵敏度等。根据实际用处不同，以及使用场合的不同对不同的换能器性能提出不同的要求，例如对军事上用的换能器与对超声测量仪用的换能器所提出的要求不一样，对发射和接收用的换能器所提出的要求也不一样。1 发射换能器和接收换能器共同要求的性能指标1.1 工作频率 超声换能器...

在利用FPGA进行卷积运算时会大量用到乘累加单元，如果用定点的方式则逻辑实现相对简单，原因是定点的加法运算是利用组合逻辑电路来实现的没有时钟延迟，这就使得加法器的输出在下一个时钟采样沿到来时立即反馈到输入端，从而可以连续的相加并最终得到累加结果。这是定点实现的优势，但是对于精度要求较高的场合，定点的方式显然不能满足要求，所以采用浮点的方式是必须的。而我们知道浮点的乘法逻辑和加法逻辑输...

在上一篇博客单精度浮点数加法器FPGA实现------（同号相加）中笔者介绍了单精度浮点数同号相加的FPGA逻辑实现，本次笔者将继续介绍异号相加的逻辑，下面给出verilog代码：module FP_ADD_diff_oper //不同符号的浮点数据相加( input wire MAIN_CLK, input wire [31:0] a, input ...

笔者这些天在研究如何用FPGA实现浮点的乘累计运算，关于浮点的乘法可以调用现成的IP核，而浮点的加法虽然有现成的IP核，但是输出时钟延迟太大，以及将输出反馈到输入端时会使系统的主时钟频率急剧降低，使整个系统能够运行的速率很低，虽然笔者针对这一问题尝试了一些补救措施，比如在输出加一级reg进行缓冲，但是没有什么效果。无奈之下笔者选择了自己编写单精度浮点加法器逻辑，终于成功的解决了上述...

首先给出FIR II IP核的Avalon-ST接口的各信号定义以上是所有的信号线说明。下面我们打开FIR II IP核的界面首先我们在Filter Specifcation界面下，一般我们需要配置的是红色框1、2、3中的内容，其中红色框1是来设置通道的个数的这里我们有192个通道，红色框2是设置模块的工作频率的这里设置为96MHz，红色框3为设置单个通道的数据率的...

声呐信号的主要工作参数有三个：工作频率、信号时间宽度T和信号的带宽B。1、工作频率的选择 工作频率的选择是根据介质的传播条件、环境、目标特性及指向性要求等因素来综合确定的。当工作频率提高时，由介质吸收造成的信号衰减增大，使得作用距离降低。然而工作于高频的声呐，其基阵易于实现窄带指向性，抑制噪声的能力增强。此外，高频环境噪声降低，有利于增加声呐的作用距离。因此，声呐的工作频率选择...

在之前的理论篇中，笔者梳理了水声通信相关的理论知识体系。本次笔者给出了一套基于BFSK的水声通信系统，该系统已经在实际的硬件中得到了验证。%******************************通信声呐仿真BPSK方式**********************%几点假设：%（1）基于射线声学理论%（2）几何衰减按球面波传播衰减规律衰减，不考虑吸收衰减%（3）考虑...

流体介质是由气体和液体抽象出来的的介质模型。其力学特征是，介质中相互接触的质团间有相互作用力，接触面微元上的相互作用力，大小正比微元面积，方向垂直于接触面微元，如图（1）所示。因而，面元上的受力与间可用标量P联系：；这个标量P就是流体中的压强。流体的这个力学特征是流体中只存在纵波的原因。流体的另一个特征，是物质空间分布的连续性，即，介质中质团连续分布无间隙。作为传播声波的流体介质，还有...

声，有两个含义，一是值空气中的振动传入人耳引起耳膜振动刺激听觉神经产生的感觉，这个‘声’可理解为单词‘声音’的简略；声的另一个是指介质中机械振动的传播，这个‘声’可理解为单词‘声波’的简略。本文中所用声的含义是指声波，是物理学中的概念。 振动频率是振动的特征量；声作为介质中机械振动的传播，可用振动频率将声分为：超声频声（简称超声），频率高于20KHz的声波；音频声，频...

     声学多普勒速度仪（DVL）是一种测量相对于水底速度的声纳设备。早年DVL一开始被设计主要是用于船舶的定位导航中，但是近年来随着水下机器人领域，尤其是民用领域的兴起（据有关机构测算这一市场规模已达100亿美元），DVL的作用显得愈发重要，因为水下机器人想要在水下实现自主导航有三样东西必不可少，一个是姿态测量单元，一个是速度测量单元，一个是位置测量单元。姿态测量单元通常是通过MEMS惯导元件...

    在水下定位和探测中，声线（即声波的主轴传播方向）往往不是由发射端沿直线传播至接收端，而是沿着曲线传播至接收端，出现这一现象的原因是由于有声速剖面的存在。当声速剖面理想时，即水中每点处的声速是不变的，这时声波是沿着直线由一点传播至另一点；当声波穿越两层不同声速的水层时，声波就会发生折射现象，这和光的折射现象是一致的，具体来讲，设声波的入射角为A1，入射层声速为C1；折射角为A2，折射层声速为...

目前C/C++ 编译器标准都遵照IEEE 制定的浮点数表示法来进行float,double 运算。这种结构是一种科学计数法，用符号、指数和尾数来表示，底数定为2—— 即把一个浮点数表示为尾数乘以2 的指数次方再添上符号。

水声通信（Underwater Acoustic Communication, UWA Communication）具有灵活、方便、经济、不存在电缆缠绕等特点，可实现导航、定位、信息交换、通信联络和安全保障所需的信息传输，是实现水下综合信息感知与信息交换的主要手段。在信息化海洋数据采集、海洋资源开发、海洋环境监测等关系到我国可持续发展的商业战略领域中，水声通信扮演者重要的角色。

在日常的生产生活中人们使用的电子功率放大器一般分为线性功率放大器和非线性功率放大器两大类。其中线性功率放大器易于实现，不会产生除信号外的其它频率成分，具有高保真的特点，但是它的缺点就是效率底下，在大功率的时候需要加装大面积的散热片来进行散热，这样不利于体积小型化。而非线性功放是要将输入的信号进行调制，因此产生了一些除了信号本身外其它的频率成分，而这些频率成分一般都是高于信号本身的频率的，因此需要在

声波作为目前已知在水下能够远距离传播的媒介，因此其在水下的通信、导航、探测等领域得到了广泛的应用。虽然电磁波在水下也能传播，但是其在水下衰减非常的严重，一般只能传播几米远，但是声波却可以传播几百甚至是上千公里。因此可以使用声波在水下进行通信，然而在水下通信时的信道要远比电磁波在空气中通信是的信道复杂的多，为了有效的在水下进行通信需要对水下的信道进行建模。

常规的KNN-SVM联合分类器中K-近邻算法没有充分挖掘训练样本的信息，使用遍历的方法来计算待识别样本与训练样本之间的距离，特别是在训练样本巨大时，存在大量的冗余计算。针对该问题，将训练样本训练成K-D树的结构，设计了K-D树KNN-SVM分类器，该分类器可以大大的减少这些多余的计算，从而提高了搜索效率，有效缩短了搜索时间。进行了仿真和实验研究，分别设计了KNN、SVM、KNN-SVM分类器对两类水下目标进行了分类识别，并对相关参数的选取进行了优化。实验结果表明：在选定了最佳参数后的KNN-SVM联合分类器

在刚刚过去的2017全国大学生数学建模比赛中，笔者有幸指导了一组本科学生参赛。对于赛题A《 CT系统参数标定及成像》中的CT系统参数标定，经过将问题进一步的提炼，问题最终变成了在平面二维空间中对任意椭圆进行拟合的问题，笔者花了大概四个小时的时间建立了该问题的数学公式表达、并推导出了求解该问题的算法、同时编写了实现该算法的Matlab程序。

在上一篇博客《图像恢复迭代算法的加速》中，提到了R-L迭代算法，此算法的核心思想源自条件概率贝叶斯理论，此文对该理论进行了详细的推导，并给出了严格的证明。

Richardson–Lucy算法是一个在泊松噪声背景下用于图像恢复的迭代技术。其旨在通过使用期望最大化算法来使恢复图像的似然性最大化。该算法需要对图像退化的过程进行良好的估计，以实现精确恢复。导致图像模糊的原因有很多，例如目标移动，离焦透镜，或着大气湍流，这些都是可由系统的点扩散函数（PSF）来描述。 一般的Richardson–Lucy算法收敛到最终的结果需要很多次的迭代，这限制了其在某些时间紧张型场合的使用。因此出现了许多的加速算法，David S. C. Biggs 和Mark

在常规的波束形成中，是将阵列上的每个接收机的同时刻数据乘上一个预设的系数然后求和并最终得到一系列的具有特定指向性的波束。但是由此得到的波束它的主瓣宽度很宽而且旁瓣衰减有限，能量泄漏很严重。为了解决这些问题，许多的学者提出了一些高分辨的波束形成算法，MVDR和DMR就是其中的代表。然而这些算法虽然能够得到很窄的波束宽度和低的旁瓣级，但是需要大量的数据样本来估计数据的协方差矩阵（如MVDR），计算量很大，而且当目标方位变化很快的时候，这种方法的性能将大打折扣甚至变得不稳定。对于DMR其需要事先知道目标的个数这一

波束形成是声纳或者雷达信号处理的一个重要的步骤，它旨在使接收的阵列能够按照期望的方向接收信号，并且抑制其它方向信号的干扰。自适应波束形成往往具有比常规波束形成更尖的指向性和抗强干扰能力。在一些情况下更具有优势。 在博客《基于Householder变换的约束自适应算法》中，我给大家翻译了一篇有关于自适应波束形成LMS算法的快速实现算法。今天我将我最近整理的几类常用的关于自适应波束形成算法的推导过程拿出来与大家分享，如果有什么不对的地方，还望大家能够批评指正。

在传统的波束形成算法中，人们通常是通过设计固定的窗函数来使形成的波束图在感兴趣的方向增益保持为1，而使其他的方向上增益尽可能的低。这种方法带来的好处是算法的实现简单，但是不足点也是很明显的，由于波束图不能动态的改变，这就有可能导致在主波束方向上没有信号而旁瓣方向的信号很强时被误认为主波束方向上有信号。这时传统的波束形成方法就不再适用了，为了解决这一问题，人们提出了自适应波束形成的方法。通过施加一定的约束条件使其可以主动的在干扰很强的方向上形成一个很深的零点从而达到抑制强干扰的目的。效果改善显著，但是采用这种

陷波器顾名思义就是对特定频率的信号有着很强的衰减的滤波器，也即阻带带宽极窄的带阻滤波器。在传统的数字陷波器设计中，为了能使某一频率信号得到足够大的衰减，通常的做法就是把阶数选的足够高来达到很大的衰减；但同时计算量也变得更大了。而且设计的过程复杂，不利于动态的调整。为了解决上述存在的问题自适应陷波器孕育而生。博主在一次心电信号(ECG)工频干扰滤除的项目中就使用了自适应陷波器，其效果要远远好于传统的陷波器，它能在信干比-100dB的时候依然能够很好的将有用的信号提取出来。

套接字(Socket)起初来源于UNIX，是加利福尼亚大学Berkeley分校开发的UNIX操作系统下的网络通信接口。随着UNIX操作系统的广泛使用，Socket亦当之无愧的成为了最流行的网络通信程序接口之一。

运算放大器相信只要是学电子类专业的一定不会陌生，然而在书上介绍的运算放大器是将实际的运放进行了抽象后再加以介绍的，所以想要使用实际的运放，还要补大量的功课。运算放大器除了书上讲的那些共性外，不同的运算放大器还有很多不同的特性，正是这些不同的特性决定了它们不同的应用场景，运算放大器一般分为如下几种：

在SYS/BIOS中，广义上指被处理器执行的任何独立的指令流。线程是一个能够调用一个函数或者中断服务程序的单点控制。在sysbios系统中一共有硬件中断(HWI)、软件中断(SWI)、任务(Task)、空闲线程(Idle)这四种线程，它们的优先级依次从高到低。关于HWI与SWI的创建笔者分别在SYSBIOS学习笔记---硬件中断(HWI)的创建和SYSBIOS学习笔记---软件中断(SWI)的创建这两篇博客中进行了详细的讲解，本次将继续介绍Task与Idle线程的创建，并对比这四种的线程的特点及注意事项。

软件中断只有在操作系统中才有这一说，它是对硬件中断的一种克隆，在SYSBIOS中它的优先级是低于硬件中断高于任务线程的，其最多可以有32个优先级，只要系统的资源允许，理论上可以创建无穷多个软件中断，而硬件中断则是少数有限个的。由于它是由SYSBIOS内核调度管理的因此在响应的时间上是要长于硬件中断的。

在进行嵌入式开发的时候，我们会根据不同的需求来决定是否使用操作系统，以及使用什么样的操作系统，在ARM嵌入式开发中大家常用Linux操作系统来进行系统的开发；而在DSP嵌入式开发中，尤其是TI的DSP，在早期的时候由于处理器的性能相对的弱一些，资源也相对较少，这时为了满足实时性的要求，一般是跑裸机程序。随着TI的DSP处理器的性能越来越强大，资源也越来越多，为了把开发人员从繁重的底层驱动开发解放出来，TI随即推出了SYSBIOS实时操作系统，该操作系统目前已经支持TI的全系列DSP。

在进行DSP算法设计时，我们通常会在Matlab上验证一下算法是否正确，然后再把算法移植到DSP当中，然而在移植的过程中会不可避免的出现各种问题，这时就需要将DSP内存中的数据导出来传到Matlab中进行分析处理，同时我们也需要从Matlab传递一些数据到DSP当中，比如我们想生成某种比较复杂的数据（例如各种类型的噪声），如果用C语言来写的话将会非常的麻烦，而这些用matlab能够亲而易举的完成，这时我们只需要将matlab生成的数据导入到DSP中就可以了，今天笔者就来向大家来分享一下实现CCS与Matla

最近笔者有幸设计了一款多层的高速PCB，虽然笔者之前设计过无数的低速PCB板其中不乏一些多层板，但是对于高速PCB还是第一次，为此笔者也是下了一番功夫，在网上看了各种帖子，去图书馆借了一堆书，研究了一段时间，总结出了一些笔者认为比较有用的东西，今天笔者就给大家来分享一下有关于高速PCB上走线的阻抗计算问题。

我们在学习使用TI的DSP集成开发环境CCS(Code Compose Studio)时，有时特别想在线的看一下内存中的数据到底是个什么样子，或者想看一下它的频谱是个什么样子，如果不知道CCS自带有绘图功能或者知道但不是很会使用，这无疑会使程序调试带来很大的不便。对于新手来说这一功能显得有些"高大上"，让很多人感觉它很难使用。今天笔者就来详细的介绍一下如何正确的使用其自带的Graph功能。

在上一篇的博客球面数据拟合算法简介中，笔者详细介绍了关于空间球面数据拟合的算法公式推导并给出了相应的Matlab代码及其仿真分析。本次笔者将上面这一情况进行更一般的推广，即取消了球面数据这一限制，数据可以是椭球面形式的，也就是说是任意的空间二次曲面形式的，可见球面的数据是它的一个特例。

非平稳信号又称时变信号。对这一类信号，其一阶、二阶统计量和功率谱的估计显然不能简单的使用平稳信号的估计方法，必须考虑它们的时变因素。

在信号的谱分析中，由于经典的谱分析方法物理意义比较明确，能够处理的信号信噪比可以很低（相对于现代的谱估计来说），而且处理算法简单稳定；因此在大多数的情况下都是使用经典的谱分析方法对信号进行谱分析

对于一个模拟信号进行采样，其采样间隔通常大于Nyquist频率。亦即对于一个最高频率为f的信号，从恢复信号波形的角度出发，采样频率必须大于最高频率的两倍。然而，在数字波束形成器中，采用数字延迟线，其延迟量化单元即为采样周期。这个量化的采样周期，如果正好使其方位波束要求的各延迟补偿值得到的波束相同。否则，波束指向不准，同时波束图会发生畸变。

当我们手中握有大量的数据时，对于二维的数据，我们会对他们进行直线拟合、对数拟合，圆曲线的拟合等等。这些拟合的方法都是运用的了非常古老而又非常有效的方法，即最小二乘法。 今天给大家介绍一种三维球面数据的拟合方法，该方法也是运用的最小二乘的方法。旨在使拟合的半径在均方意义下误差达到最小。

我知道在我的读者之中很多人可能正在往上努力打拼，有些人可能在基层位置上工作了很多年遇到了职业发展的天花板，他们并不想放弃自己所喜欢的专业，但是觉得在单位中如果无法当官就没有前途，一些人听说国外的一些公司对工程师很重视非常的向往，同时也抱怨自己的单位对自己不重视，并且将这种情况上升到制度文化的高度，那么今天我们就来谈一谈这个问题。

对于为一个准备学习开发TI的DSP处理器的小白来说，学会自己装CCS(Code Composer Studio)是学习开发TI DSP的第一步，今天我就给大家来详细的讲解如何一步一步的安装CCS。

开放系统互连参考模型 (Open System Interconnect 简称OSI）是国际标准化组织(ISO)和国际电报电话咨询委员会(CCITT)联合制定的开放系统互连参考模型，为开放式互连信息系统提供了一种功能结构的框架。它从低到高分别是：物理层、数据链路层、网络层、传输层、会话层、表示层和应用层。其最主要的功能就是帮助不同类型的主机实现数据传输 。

TMS320C6678是KeyStone架构的8核DSP处理器，每个CorePac核的频率最高为1.25 GHz，提供强大的定点和浮点运算能力，可以广泛地应用在通信、雷达、声纳、火控、电子对抗等领域。从目前的情况看，由于C6678的以上优异的特性，基于TMS320C6678的硬件和软件平台，在未来的5~10年内，将是信号处理平台主流。

如果从事 TI DSP 处理器开发，CMD 文件一定是不可忽略的重要文件。当然，CMD 文件也并不是 DSP 处理器的专利，开发 TI ARM（不运行 HLOS 高阶操作系统），MCU 等处理器也都会遇到 CMD 文件。本文就以C6000 DSP CMD 文件为例，给大家说说 CMD 文件的那些事。

在嵌入式系统调试过程中，终端是最常用的调试方式。虽然 DSP 的调试大部分都是在 CCS 下通过 Console 输出调试信息，但是一般也会用到串口或者 Telnet。

基于多波束测深声纳工作原理，结合国内外主流浅水多波束测深声纳产品与技术发展现状，分析了浅水多波束测深声纳具备的5 种典型能力，提出了从能力表征的角度理解浅水多波束测深声纳性能指标; 总结了保障这些能力所采用的主要关键技术，讨论了技术实现的具体思路，列举了仿真或试验数据处理结果; 展望了多波束测深声纳的发展趋势。

关于自适应LMS的理论基础已经非常的成熟，随便找一本关于自适应滤波器的书就会有介绍相关的内容，有的还可出了它的具体算法，但是还没有一本书有讲过怎样编写能够时实（Real Time）处理的基于C的自适应LMS算法（至少我没有见过），而这个应该说是将算法用于时实处理最为关键的一步，本次我就给大家分享一下关于我是如何用C语言来实现时实自适应LMS算法的。

在上一篇的博客文章中(浅谈自适应滤波器—（快速RLS算法）)，我给出了关于RLS的快速算法（简称FTRLS），但是它有一个缺点就是不够稳定，虽然我给出的例子里是稳定的，这主要是Matlab默认的计算精度是double型的，非常的高，所以没有出现问题。后来我在C6678 DSP用C语言实现该算法的时候使用的是float型的，得到的结果较Matlab的要差一些。不稳定性分析针对这种不稳定问题，一些研究人

在上一篇博客中（浅谈自适应滤波器）我给大家介绍了关于自适应滤波器的一些入门级的知识，并分析了常规RLS算法单次迭代的计算量级为O[N2]，当阶数N增大时相应的计算量显著增大，为了将计算量级降低到O[N]，人们提出了很多快速的RLS算法，其中快速横向递归最小二乘算法(FTRLS)最具有吸引力，得到了广泛的应用。 本次将为大家介绍该算法的由来并给出相应的算法，在推导FTRLS算法时，需要同

在通常的滤波场合中，从频域的角度进行滤波，其相关的理论已经相当的成熟，只要给出相应的设计指标就可以很方便的设计出满足要求的滤波器。然而在更一般的情况下，人们所需要的滤波器工作的环境是时变的，这就导致事先已经设计好的滤波器性能下降甚至不能被使用。例如，在水下声纳通信中，信息通过声波在水中进行传播到达接收端，其信道是时变的（这种情况比电磁波在空气中传播的情况要严重的多），如果只是简单的用固定系数的滤波器