Вейвлет-преобразование

   public void FWT(double[] data)
        {
            double[] temp = new double[data.Length];
 
            int h = data.Length >> 1;
            for (int i = 0; i < h; i++)
            {
                int k = (i << 1);
                temp[i] = data[k] * s0 + data[k + 1] * s1;
                temp[i + h] = data[k] * w0 + data[k + 1] * w1;
            }
 
            for (int i = 0; i < data.Length; i++)
                data[i] = temp[i];
        }
public void FWT(double[,] data, int iterations)
        {
            int rows = data.GetLength(0);
            int cols = data.GetLength(1);
 
            double[] row = new double[cols];
            double[] col = new double[rows];
 
            for (int k = 0; k < iterations; k++)
            {
                for (int i = 0; i < rows; i++)
                {
                    for (int j = 0; j < row.Length; j++)
                        row[j] = data[i, j];
 
                    FWT(row);
 
                    for (int j = 0; j < row.Length; j++)
                        data[i, j] = row[j];
                }
 
                for (int j = 0; j < cols; j++)
                {
                    for (int i = 0; i < col.Length; i++)
                        col[i] = data[i, j];
 
                    FWT(col);
 
                    for (int i = 0; i < col.Length; i++)
                        data[i, j] = col[i];
                }
            }
        }