fftadwave.c

超声波成像算法

/*               A/D偝傟偨庴怣攇宍偺僼乕儕僄曄姺             */
/*               FFTwave                                     */
#include <math.h>
#include "Reconst.h"

/* n_rec: 庴怣攇宍偺悢 */
/* n_time: 帪娫僒儞僾儖悢 */
/* n_band: 廃攇悢僒儞僾儖悢 */
/* int rec: FFT偝傟傞庴怣攇宍僨乕僞梡攝楍丆悺朄偼n_rec*n_time */
/* rec_f: FFT偝傟偨庴怣攇宍僨乕僞梡攝楍丆悺朄偼n_rec*n_band */
/* h: 嶌嬈梡攝楍丂悺朄偼2*n_time*2 */
void FFTAdWave(int n_rec,int n_time, int n_band, int i_shot, double *h, int *rec, COMPLEX16 *rec_f)

{
	int i, i_time, i_band,i_bias, f_bias;
	int *recp;
	double *h1_real,*h1_imag,*h2_real,*h2_imag;
	double *h1_real_p,*h1_imag_p;
	h1_real=h;
	h1_imag=h+  n_time;
	h2_real=h+2*n_time;
	h2_imag=h+3*n_time;

	f_bias=0;

	for(i=0; i<n_rec; i++)
	//for(i=0; i<0; i++)
	{//庴怣巕偵娭偡傞儖乕僾
		recp=rec+i*n_time+i_shot*n_rec*n_time;
		h1_real_p=h1_real;
		h1_imag_p=h1_imag;

		for(i_time=0;i_time<n_time;i_time++){
			*h1_real_p++=*recp++;
			*h1_imag_p++=0.0;
		}//i-time loop

		FFT(n_time, h1_real, h1_imag, -1.0);
		
		i_bias=i*n_band;		
		for(i_band=0;i_band<n_band;i_band++){
			rec_f[i_bias+i_band].real = h1_real[i_band];
			rec_f[i_bias+i_band].imag = h1_imag[i_band];
		}
		
		/*******************妔********************/

		//if(i==0)
		//{
		//	OneDimensionArrayOutput(h1_real,n_band);
		//}
		/*******************妔********************/
	
	}//i loop

	return;
}