lms

 

   	int i;
   	int j;
	float tmp1,tmp2;
        float dhat;
        float temp;
        float error;

 

	/*
	 * Allocate memory and return pointers for tapped delay line x_P and
	 * array containing impulse response samples, h_P.
	 *
	 */
	if( (x_P = (float*)calloc(N,sizeof(float))) == NULL ||
	    (w_P = (float*)calloc(N,sizeof(float))) == NULL ) {
	   	fprintf(stderr,"lms: can't allocate work space\n");
		return(4);
	}
	/*
	 * initialize the tapped delay line and weights to zero.
	 *
	 */
	for (i=0; i
    

 



        for(j = MIN_AVAIL(); j>0; j--) {
		/*
		 * Shift input sample into tapped delay line
		 */
		IT_IN(0);
		tmp2=x(0);
		for(i=0; i
    

 

	free(x_P); free(w_P); 

/*  Capsim (r) Text Mode Kernel (TMK) Star Library (Blocks)
    Copyright (C) 1989-2017  Silicon DSP Corporation

    This library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Lesser General Public
    License as published by the Free Software Foundation; either
    version 2.1 of the License, or (at your option) any later version.

    This library is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public
    License along with this library; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

    http://www.silicondsp.com
    Silicon DSP  Corporation
    Las Vegas, Nevada
*/

 

/* lms.s */
/***********************************************************************
                             lms()
************************************************************************
This star implements a simple LMS adaptive filter.
Param:	
	1 - (int) N  filter order.
	2 - (float) LMS gain constant.
	3- (int) Flag to output either estimate or error

lms


This star implements a simple LMS adaptive filter.
Param:	
	1 - (int) N  filter order.
	2 - (float) LMS gain constant.
	3- (int) Flag to output either estimate or error


Date:  September 23, 1988 
Programmer: Adali, Tulay 

*/