dctf.cpp

一个快速算法的例子

#include "stdafx.h"
/*#include "dctf.h"
#include "math.h"
#define  pi 3.1416926
float k(int t)
{
	if(t=0)
	return 1/(sqrt(2));
	else 
		return 1;
}

void dct(int t0[][],float t1[][],int m,int n)
{
	float s;
	for(int i=0;i<m;i++)
		for(int j=0;j<n;j++)
		{   
			for(int x=0;x<m;x++)
		    for(int y=0;y<n;y++)
			{
				s+=t0[x][y]*cos(2*x+1)*i*pi/(2*m)*cos(2*y+1)*j*pi/(2*n);
			}
		
			t1[m][n]=2/(sqrt(m*n))*k(m)*k(n)*s;
		}
}
		
*/