用VB6实快现速傅里叶变换FFT算法.docVIP

用VB6实现快速傅里叶FFT算法 设原始数据在Xr()中,实现FFT算法分两个步骤,首先是实现倒序,例如Xr(0)…Xr(15)的序列要变成: Xr(0),Xr(8),Xr(4),Xr(12),Xr(2),Xr(10),Xr(6),Xr(14),Xr(1),Xr(9),Xr(5),Xr(13),Xr(3),Xr(11),Xr(7),Xr(15), 0,8,4,12…有什么规律吗,比如第三个数2，二进制是0010，反序的写法是0100，这就是十进制的4。实现倒序之后，要做蝶形运算(二基),不说原理了,写VB代码。 位序倒置代码 Sub 倒序(X_() As Double) Dim N As Integer, I As Long, J As Long, mn As Long, LH As Long, T As Double, k As Long 位序倒置 N = UBound(X_) 求数组大小,其值必须是2的幂 LH = N / 2 J = N / 2 For I = 1 To N - 2 If I J Then 倒序 T = X_(J) X_(J) = X_(I) X_(I) = T End If Debug.Print I, J k = LH 下面是向右进位算法 Do If k J Then Exit Do 高位是1吗? J = J - k 是的,高位置0 k = k / 2 准备次高位的权 Loop Until k = 0 次高位的权若非0，则检查新的次高位 J = J + k 否则若最高位是0，则置1 Next End Sub 蝶形算法代码 Sub 蝶形算法(Xr() As Double) Dim L As Long, LE As Long, LE1 As Long, N As Long, R As Long, P As Long, Q As Long, M As Byte Dim Wr As Double, Wi As Double, W1r As Double, W1i As Double, Tr As Double, Ti As Double Dim Pi As Double, T As Double Dim Xi() N = UBound(Xr) 求数组大小,其值必须是2的幂 M = 0 L = 2 Pi = 3.14159265358979 Do L = L + L M = M + 1 Loop Until L N N = L / 2 ReDim Xi(N - 1) L = 1 Do LE = 2 ^ L LE1 = LE / 2 Wr = 1 Wi = 0 T = Pi / LE1 W1r = Cos(T) W1i = -Sin(T) R = 0 Do P = R Do Q = P + LE1 Tr = Xr(Q) * Wr - Xi(Q) * Wi Ti = Xr(Q) * Wi + Xi(Q) * Wr Xr(Q) = Xr(P) - Tr Xi(Q) = Xi(P) - Ti Xr(P) = Xr(P) + Tr Xi(P) = Xi(P) + Ti P = P + LE Loop Until P N - 1 Wr = Wr * W1r - Wi * W1i Wi = Wr * W1i + Wi * W1r R = R + 1 Loop Until R LE1 - 1 L = L + 1 Loop Until L M For I = 0 To N - 1 仅输出模 Xr(I) = Sqr(Xr(I) ^ 2 + Xi(I) ^ 2) Next End Sub 检验的时候可以这样: Sub 检验() Dim Y(63) As Double For I = 0 To 64 Y(I) = Sin(2 * 3.1415926 * I / 16) Next 倒