c语言实现小波多级分解与重构（跟matlab效果一致）

原创

已于 2023-01-04 15:38:20 修改 · 3.1k 阅读

39 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #c语言 #重构 #算法

于 2023-01-04 15:21:31 首次发布

文章介绍了如何使用小波变换进行数字信号的降噪滤波，提供了C语言实现的小波分解`waveDec`和重构`wrcoef`函数，并通过与MATLAB自带函数的比较，证明了代码实现与MATLAB结果的一致性。

前言

最近在做数字信号处理方面的内容，需要用到小波变换来进行降噪滤波，具体实现如下。

一、小波分解代码


#define FILTER_LEN					 (8)			//滤波器长度
#define ADS_WAVELET_LEN             （260）			//数据的长度
//多级小波分解，"db4"
double db4_Lo_D[8] = {
   
    -0.0105974017850690, 0.0328830116668852, 0.0308413818355607, -0.1870348117190931, -0.0279837694168599, 0.6308807679398587, 0.7148465705529154, 0.2303778133088964};	
double db4_Hi_D[8] = {
   
    -0.2303778133088964, 0.7148465705529154, -0.6308807679398587, -0.0279837694168599, 0.1870348117190931, 0.0308413818355607, -0.0328830116668852, -0.0105974017850690};

/*
*********************************************************************************************************
*   函 数 名: wavedec
*   功能说明: 小波多级分解函数（与matlab效果一致）
*   形    参: srcData	输入数据
			  srcLen	输入数据长度
			  level		分解级数
			  C         存储的细节系数cDi和最后一层的近似系数cA。以5层分解为例，它们的存储结构是，C=[cD1;cD2;cD3;cD4;cD5;cA5]
			  L			源数据长度和各组系数的长度。以5层分解为例，L=[Len;Len1;len2;len3;len4;len5;len5]
*   返 回 值: 无
*********************************************************************************************************
*/
void waveDec(double *srcData, int srcLen, int level, double *C, int *L )
{
   
   
	double tempsrc[ADS_WAVELET_LEN] = {
   
   0};
	int resultlen = 0;
	int tempsrclen = 0;
	//int filterLen = 8;	//滤波器长度
	int index = 0;
	int curlevel = 0;	//当前级数
	int lenIdx = 0;
	
	double temp = 0;
	
  int i,j,p;

	L[lenIdx++] = srcLen;
    for ( i = 0; i < srcLen; i++ )
        tempsrc[i] = srcData[i];

    tempsrclen = srcLen;
	
	while(1)
	{
   
   
		resultlen = (

最低0.47元/天解锁文章