tride.c

来自「java-fft-demo,用java开发的实现fft的平台」· C语言 代码 · 共 46 行

/********************************************************/   
/*Tride							        */   
/*Maritza Rodr韌uez Mart韓ez				        */   
/*Computational Signal Processing Group - march-98   	  */   
/*CSPG - Dr. Domingo Rodr韌uez - Coordinator         	  */   
/*							              */    
/*Description:						        */    
/*This function tridat generates a vector with the 	  */   
/*elements of the main diagonal that composes the 	  */    
/*diagonal matrix t(N,S) 				        */    
/********************************************************/    
   
  
COMPLEX **Tride(int N, int s)   
   
{   
	int R,i,j,auxR=0,auxS=0;   
	COMPLEX **tride;    
	R = N/s;   
	tride = getComplexMatrix(N);   
  
   
    
	  	for(i = 0; i< N; i++)   
			for(j = 0; j< N; j++){   
				if(j == i){   
					tride[i][j].real = cos((2*PI*auxR*auxS)/N);   
					tride[i][j].imag = -sin((2*PI*auxR*auxS)/N);   
					++auxR;    
						if(auxR==R){   
							auxR=0;   
							auxS=1;   
						}   
   
				}   
 				else{     
					tride[i][j].real=0.0;    
					tride[i][j].imag=0.0;    
				}    

		}    

   
return tride;    

}