快速傅立叶变换对信号频谱的分析 C语言(国外英语资料).docVIP

快速傅立叶变换对信号频谱的分析 C语言(国外英语资料) # includes graphics.h # includes stdio.h # includes stdlib.h # includes string.h # includes conio.h # includes math.h const int maxsize = 1600. const double pi = 3.1416; const int n = 32; void ss (double road [], double xi [], int n) { int i = 0, j, s1. float a, b; for (j = 1; j n; j + +) { for (s1 = n / 2; s1 = i; s1 = s1 / 2) { i = i - 1. } i = i + 1; if (i j) { a = xr []; bj = xi []; xi [i] = xi [j]. xr [i] = xr [j]. xr [j] = a; xi [j] = bj. } } } void fft (double road [], double xi [], int m, int n) {int loop1, loop2, loop3,, le1, ip, i = 1; float wr, wl, tr, i, ur, ur1, iu, ls. ls = log (n) / log (2). ss (xr, xl, n); for (loop1 = 1; loop1 = ls; loop1 + +) { ur = 1.0; iu = 0.0. the = pow (2, loop1); le1 = / 2; wr = cos (pi * m / le1); wi = - sin (pi * m / le1); for (loop2 = 0; loop2 le1; loop2 + +) { for (loop3 = loop2; loop3 n; loop3 = loop3 +) (ip = loop3 + le1; tr = ur * xr [ip] - iu * xi [ip]; i = a * xi [ip] + iu * xr [ip]; xr [ip] = xr (loop3] - t; xi [ip] = xi (loop3] you. xr (loop3] = rs + xr (loop3]; xi [loop3] = i + xi [loop3]; } ur1 = wr * ur iu * wi. iu = wr * iu + ur * wi. ur = ur1;}} if (i = = - 1) { for (i = 0; i n; i + +) { xr [i] = xr [] / n; xi [i] = xi [] / n; } } } void hua (double x []) {int d = detect, m = vgahi; initgraph ( d m, ); int i, j; for (i = 0; i maxsize; i + +) line (50 + 1, 200 - 100 * x [i], + 51200 - 100 * x [i + 1]). } void hua2 (double x []) {int d = detect, m = vgahi; initgraph ( d m, ); int i, j; for (i = 0; i maxsize; i + +) line (i + as - 10 * x [i], + 10 * x 41200 [i + 1]). } void gaoshi () {int i, j; float p, q; double xx [maxsize]; double xe [maxsize]; double x1 [maxsize]; for (i = 0; i maxsize; i + +) xx [i] = 0. printf (输入p \ n); scanf (% f, p); printf (输入q \ n); scanf (% f, q); for (i = 0; i maxsize; i + +) { if (i n) x1 [i] = exp (- (1 - p) * (s - p) / q). else x1 [i] = 0. } fft (x 1, xx, 1256). for (i = 0; i maxsize; i + +) xe [i] = sqrt (x1 * x1 [i] [i] +