example 6-19.c

《基于TI DSP的通用算法实现》程序代码

// Example 6 - 19. LMS Implementation C Listing Using Intrinsics for the TMS320C64x DSP

	for (i = 0; i < L/2; i++)
		yout += _dotp2(x[i],w[i]); 		// 乘累加，采用了内部函数_dotp2，
								// 循环次数减少了一半
	err = d-(int)(yout>>16);    			// 计算误差信号err
	mu_err = _mpy(mu, err);    		// 计算ue(n-1)
    	for (i = L-1; i >= 0; i--){
		temp = ((_smpy2(mu_err_v,x[i])))	;	//ue(n-1)x[i], ue(n-1)x[i-1]
		temp_hi = _hi(temp);           // 数据解包
		temp_lo = _lo(temp);
		temp_com = _spack2(temp_hi, temp_lo);  
		w[i] = _add2(temp_com, w[i]);    // 更新系数
	}