Tham khảo tài liệu ''giáo trình xử lý ảnh y tế tập 1b p16'', công nghệ thông tin, đồ họa - thiết kế - flash phục vụ nhu cầu học tập, nghiên cứu và làm việc hiệu quả
{ for(j=0; j L thuật toán gọi thuật toán giảm lược đầu FFT Thuật toán giảm lược đầu vào FFT Trường hợp làm hồn thiện thuật tốn phân chia tần số Hình 6.8 giới thiệu trường hợp M = L = Từ hình 6.8 nhận thấy (L-M) bước có phần tử bướm L bước cuối có tồn bướm Sơ đồ giúp thay đổi chương trình 6.2 thành chương trình 6.3 Chương trình 6.3 "FFTIP.C" Giảm lược đầu vào FFT /**************************** * Program developed by: * * M.A.Sid-Ahmed * * ver 1.0 1992 * * @ 1994 * ***************************/ /* FFT - input pruning routine */ void bit_reversal(unsigned int *, int , int); void WTS(float *, float *, int, int); void FFTP(float *xr, float *xi, float *, float *, int,int,int, int); void FFTP(float *xr, float *xi, float *wr, float *wi, int m_output, int N_output, int m_input, int N_input ) { /* FFT pruning algorithm Deimation-in-frequency algorithm Note: Noutput=2 to the power of m_output N_output=Number of samples in the output sequence M_input=Number of samples in the input sequence This should also be a multiple of 2 The output arrays are left in bit-reverse order You will need to use routine "bit-reversal" to place them in normal ascending order 101 The twiddle factors are assumed to be stored in LUT's wr[I and wi[I You will need to use routine LUT for calculating and storing twiddle factors */ int ip,k,kk,l,inc r,iter,j,i; float Tr,Ti,diffr,diffi; ip=N_output>>1; kk=l ; incr=N_output; for(iter=0; iter