C语言程序设计论文.docVIP

C语言程序设计 快速傅立叶变换 李喜付 （东北农业大学 工程学院） 班级 机化0802 姓名 李喜付 学号 指导教师 范永存 快速傅立叶变换 摘要 实现二维快速傅立叶变换，这是第一次通过读程序流程图来设计程序，同时也将借此实例学习C语言实现矩阵的基本操作、复数的基本操作，复习前面所学过的动态内存分配、文件操作、结构指针的函数调用等内容。 关键词 C语言程序设计 傅立叶变换 C语言论文 程序流程图 本程序将对二维数组进行正、反快速快速傅立叶变换。 #include stdio.h #include stdio.h #include math.h #define PI 3.1415926 Struct COMPLEX { Float re; Float im; } cplx , * hfield , * s , * r , * w ; Int n , m ; Int ln , lm ; Void initiate ( ) ; Void dfft () ; Void rdfft () ; Void showresult () ; Void fft( int l , int k ) ; Void reverse ( int t , int k ) ; Void w （ int l）; Int loop ( int l ); Void conjugate ( ) ; Void add (struct complex * x , struct complex * y , struct complex * z ) ; Void sub (struct complex * x , struct complex * y , struct complex * z ) ; Void mul (struct complex * x , struct complex * y , struct complex * z ) ; Struct complex *hread (int j , int j ); Void hwrite ( int i ,int j ,struct complex x ) ; Void main ( ) ; { Initiate ( ) ; Printf ( \n 原始数据：\n) ; Showresult ( ) ; Getchar ( ) ; Dfft ( ) ; Printf (\n快速傅立叶变换后的结果：\n); Showresult ( ) ; Getchar ( ); Rdfft ( ) ; Printf ( \n快速傅立叶你变换后的结果：\n); Showresult ( ); Getchar ( ) ; Free ( hfield ); } Void initiate ( ) File * df ; df = fopen (data.txt ,r ) ; Fscanf ( df , %5d, n ) ; Fscanf ( df ,%5d, m ) ; If ( ( ln = loop ( n ) )== -1 ) { Printf ( 列数不是2的整数次幂 ) ； Exit ( 1 ) ; } If ( ( lm = loop ( m ) ) == -1 ) { Printf ( 行数不是2的整数次幂 ) ； Exit ( 1 ) ; } Hfield = ( struct complex * ) malloc ( n * m sizeof ( cplx )) ; If ( fread ( hfield , sizeof ( cplx ) , m * n ,df ) != ( unsigned ) ( m * n ) ) { If ( feof ( df) ) printf ( premature end of file ) ; Else printf ( file read error ) ; } Fclose ( df ) ; } Void dfft ( ) { Int i , j ; Int l , k ; l = n; k = ln ; w = ( struct complex * ) calloc ( l , sizeof ( cplx ) ) ; r = ( struct complex * ) calloc ( l , sizeof ( cplx ) ) ; s = ( struct complex * ) calloc ( l , sizeof ( cplx ) ) ; W ( l ) ; For ( i = 0 ; j m ; j++ ) For