lms c语言,LMS算法实现自适应滤波器（C语言版）

上次只是发表了MATLAB版本的LMS算法，这次将C语言版的LMS算法一并发表，其中涉及到雅克比公式求矩阵最大特征值的部分我将后续发表到博客中，此C语言版本是本人自己研究MATLAB语言，然后独自翻译成C语言完成的，为了完成自适应滤波器，使用了STM32F767作为开发板，完整压缩包我已上传，需要者可移步下载。

LMS.c

#include "LMS.h"

/* xn--------输入的信号序列(列向量)

* itr-------迭代次数，标量，默认为xn的长度，M

* en--------误差序列(itr*1)列向量

* dn--------所期望的响应序列(列向量)

* M---------滤波器的阶数(标量)

* mu--------收敛因子(步长)标量

* W---------滤波器权值矩阵，大小为M*itr

* yn--------实际输出序列(列向量)*/

/*LMS具体算法*/

float * LMS_Filter(int itr, const float *xn, const float *dn, double mu, int length)

{

static int i = 0;

static int k = 0;

static float y = 0.0;

static float en[F_COUNT];

static float W[M][F_COUNT];

static float x[M];

static float yn[F_COUNT];

/*创建一个en全零矩阵ÿ