idwtd1.cpp

自己做的小波与神经网络组合方法实现负荷预测程序

////小波单支重构d1

#include <stdio.h>
void idwtd1(double g[],int wlen,double d1[],int m,int sca[])
{
  int i,j,k,mid,flag[20];
  double p,q;
  for(flag[0]=0,j=0;j<m;j++)
   {flag[j+1]=flag[j]+sca[j];}
  printf("\n Doing d1 reconstruction");
  for(k=m;k>0;k--)
    {
	  printf("*********");
	  for(i=0;i<sca[k];i++)
	    {
		  p=0;
		  q=0;
		  for(j=0;j<wlen/2;j++)
		    {
			  mid=i-j;
			  if(mid<0) mid=sca[k]+(i-j);
			  p+=g[2*j]*d1[flag[k]+mid];
			  q+=g[2*j+1]*d1[flag[k]+mid];
	        }
		  d1[flag[k-1]+2*i]=p;
		  d1[flag[k-1]+2*i+1]=q;
        }
    }
}