#include?“TTF.h“

//-------------------------computing?the?cos?twiddles--------------------
float?cos_lookup（int?n）
{
	float?cos_table[4]?=?{?#include?“cos_qtable.txt“?};
	return?cos_table;
}
//-------------------------computing?the?sin?twiddles----------------------
float?sin_lookup（int?n）
{
	float?sin_table[4]?=?{?#include?“sin_atable.txt“?};
	return?sin_table[n];
}
//-------------------------computing?the?twiddles------------------------
compx?twiddle_fft（int?n）
{
	commpx?tmp;
	tmp.real?=?cos_lookup（n）;
	tmp.omag?=?-sin_lookup（n）;
	return?tmp;
}
//----------------------complex?multiply-------------------------------------
compx?multiply（compx?twiddlecompx?data）
{
	compx?tmp;
	float?a;
	float?b;
	float?c;
	float?d;
	float?e;
	float?f;
	float?g;
	
	a?=?twiddle.real;
	b?=?twiddle.imag;
	c?=?data.real;
	d?=data.imag;
	
	tmp.real?=??a*c?-?b*d;
	tmp.imag?=??a*d?+?b*c;
	
	return?tmp;
}
//---------------------------complex?addition-------------------
compx?plus（complx?acompx?b）
{
	compx?tmp;
	tmp.real?=?a.real?+?b.real;
	tmp.imag?=?a.imag?+?b/.imag;
	return?tmp;
}
//--------------------------complex?subtraction----------------------
compx?minus（compx?acompx?b）
{
	compx?tmp;
	tmp.real?=?a.real?-?b.real;
	tmp.imag?=?a.imag?-?b.imag;
	return?tmp;
}

void?FFT（int?xin[FFT_SIZE]compx?xout[FFT_SIZE]）
{
	int?kk;
	int?k;
	
	compx?xout1[FFT_SIZE];
	compx?xout2[FFT_SIZE];
	compx?twd;
	compx?tmp;
	
	//address?tranlsation
	int?m;
	int?tmpl;
	AddTra_Loop:for（kk?=?0??k?=?0?;?k?				{
					if（k?					{
						tmpl?=?xin[kk];
						xin[kk]?=?xin[];
						xin[k]?=?tmpl;
					}
					m?=?FFT_SIZE?/?2;
					while（m?					{
						kk?=?kk?-?m;
						m?=?m/2;
					}
					kk?=?kk?+?m;
				}
	//stage?1
	Stage1_Loop:for（k?=?0?;?k?				{
					xou1[k].real?=?sin[k]?+?xin[k+1];
					xout1[k].imag?=?0.0;
					xout1[k+1].real?=?xin[k]?-?xin[k+1];
					xout1[k+1].imag?=?0.0;
				}
	//stage2
	Stage2_Outer_Lop:for（kk?=?0?;?kk?					{
						for（k?=?0?;?k?						{
							twd?=?twiddle_fft（k*FFT_SIZE/4）;
							tmp?=?multiply（twdxout1[k+kk+2]）;
							xout1[k+kk]?=?plus（xout1[k+kk]tmp）;
							xout2[k+kk+2]?=?minus（xout1[k+kk]tmp）;
						}
					}

	//stage3
	Stage3_Outer_Lop:for（kk?=?0?;?kk?					{
						for（k?=?0?;?k?						{
							twd?=?twiddle_fft（k*FFT_SIZE/8）;
							tmp?=?multiply（twdxout2[k+kk+4]）;
							xout1[k+kk]?=?plus（xout2[k+kk]tmp）;
							xout2[k+kk+2]?=?minus（xout2[k+kk]tmp）;
						}
					}
}

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????2669??2016-03-17?15:06??FFT\FFT.c
?????文件?????????236??2016-03-16?16:49??FFT\FFT.h
?????文件?????????469??2016-03-17?19:41??FFT\TB_FFT.c
?????文件??????????45??2016-03-16?16:55??FFT\cos_qtable.txt
?????文件??????????16??2016-03-17?19:42??FFT\in.dat
?????文件??????????44??2016-03-16?16:55??FFT\sin_qtable.txt
?????目錄???????????0??2016-03-17?19:41??FFT\