1 /* ======================================================================= */
2 /* DSPF_sp_lms_cn.c -- LMS Adaptive Filter */
3 /* Natural C Implementation */
4 /* */
5 /* Rev 0.0.2 */
6 /* */
7 /* Copyright (C) 2011 Texas Instruments Incorporated - http://www.ti.com/ */
8 /* */
9 /* */
10 /* Redistribution and use in source and binary forms, with or without */
11 /* modification, are permitted provided that the following conditions */
12 /* are met: */
13 /* */
14 /* Redistributions of source code must retain the above copyright */
15 /* notice, this list of conditions and the following disclaimer. */
16 /* */
17 /* Redistributions in binary form must reproduce the above copyright */
18 /* notice, this list of conditions and the following disclaimer in the */
19 /* documentation and/or other materials provided with the */
20 /* distribution. */
21 /* */
22 /* Neither the name of Texas Instruments Incorporated nor the names of */
23 /* its contributors may be used to endorse or promote products derived */
24 /* from this software without specific prior written permission. */
25 /* */
26 /* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS */
27 /* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT */
28 /* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR */
29 /* A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT */
30 /* OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, */
31 /* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT */
32 /* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, */
33 /* DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY */
34 /* THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT */
35 /* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE */
36 /* OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. */
37 /* */
38 /* ======================================================================= */
40 #pragma CODE_SECTION(DSPF_sp_lms_cn, ".text:ansi");
42 #include "DSPF_sp_lms_cn.h"
44 float DSPF_sp_lms_cn(const float *x, float *h, const float *y_i,
45 float *y_o, const float ar, float error, const int nh, const int nx)
46 {
47 int i, j;
48 float sum;
50 for (i = 0; i < nx; i++)
51 {
52 for (j = 0; j < nh; j++)
53 h[j] = h[j] + (ar * error * x[i + j - 1]);
55 sum = 0.0f;
57 for (j = 0; j < nh; j++)
58 sum += h[j] * x[i + j];
60 y_o[i] = sum;
61 error = y_i[i] - sum;
62 }
64 return error;
65 }
67 /* ======================================================================= */
68 /* End of file: DSPF_sp_lms_cn.c */
69 /* ----------------------------------------------------------------------- */
70 /* Copyright (c) 2011 Texas Instruments, Incorporated. */
71 /* All Rights Reserved. */
72 /* ======================================================================= */