用VB实现FFT信号分析的程序分享

dragoon0114

我不知道怎么上传附件，只好把程序粘贴了过来。界面的话根据程序自设吧，要改函数在最后。

1. Dim FUZHI(255), PINZHIr(255), PINZHIi(255), PINZHI(255) As Double
   
2. Dim TU_DATA(2047) As Integer
   
3. Dim P_START, I_TIME, P_POINT As Integer
   
4. Dim PAI, A1, A2, A3, S1, S2, S3 As Double
   
5. Dim fr(7, 255) As Double
   
6. Dim fi(7, 255) As Double
   
7. Dim Shuxu(255) As Integer
   
8. Dim TX As Integer
   
9. Private Sub Command1_Click()
   
10. '启动，设置初始值
    
11. I_TIME = 0
    
12. P_START = 0
    
13. Picture1(0).Cls
    
14. Picture1(1).Cls
    
15. Picture1(0).Visible = True
    
16. Picture1(0).Line (6400, 1200)-(0, 1200), RGB(255, 0, 0)
    
17. Timer1.Enabled = True
    
18. End Sub
    
19. 
    
20. Private Sub Command2_Click()
    
21. End
    
22. End Sub
    
23. 
    
24. Private Sub Command3_Click()
    
25. '频谱图形显示方式
    
26. TX = (TX + 1) Mod 2
    
27. End Sub
    
28. 
    
29. Private Sub Form_Load()
    
30. TX = 0
    
31. PAI = 3.14159
    
32. Timer1.Enabled = False
    
33. '加权系数W
    
34. For i = 0 To 7
    
35. For j = 0 To 2 ^ i - 1
    
36. aa = Cos(2 * PAI * j / 2 ^ (i + 1))
    
37. Call ZHENG(aa)
    
38. fr(i, j) = aa
    
39. fr(i, j + 2 ^ i) = -aa
    
40. bb = Sin(2 * PAI * j / 2 ^ (i + 1))
    
41. Call ZHENG(bb)
    
42. fi(i, j + 2 ^ i) = bb
    
43. fi(i, j) = -bb
    
44. Next j
    
45. Next i
    
46. 
    
47. '倒序
    
48. For i = 0 To 255
    
49. Shuxu(i) = 0
    
50. For j = 0 To 7
    
51. a = i And (2 ^ (7 - j))
    
52. If a > 0 Then a = 1
    
53. Shuxu(i) = Shuxu(i) + a * (2 ^ j)
    
54. Next j
    
55. Next i
    
56. Label12.Caption = "截止窗" + Str(Int((HScroll6.Value / 256) * 100) / 100)
    
57. End Sub
    
58. 
    
59. Private Sub HScroll1_Change()
    
60. '合成谐波信号的频率
    
61. A1 = HScroll1.Value / 2
    
62. Label5.Caption = A1
    
63. End Sub
    
64. 
    
65. Private Sub HScroll2_Change()
    
66. '合成谐波信号的相位
    
67. S1 = HScroll2.Value * PAI / 180
    
68. Label6.Caption = HScroll2.Value
    
69. End Sub
    
70. 
    
71. Private Sub HScroll3_Change()
    
72. '合成谐波信号的频率
    
73. A2 = HScroll3.Value / 2
    
74. Label7.Caption = A2
    
75. End Sub
    
76. 
    
77. Private Sub HScroll4_Change()
    
78. '合成谐波信号的相位
    
79. S2 = HScroll4.Value * PAI / 180
    
80. Label8.Caption = HScroll4.Value
    
81. End Sub
    
82. Private Sub HScroll6_Change()
    
83. Label12.Caption = "截止窗" + Str(Int((HScroll6.Value / 256) * 100) / 100)
    
84. End Sub
    
85. 
    
86. Private Sub HScroll7_Change()
    
87. '合成谐波信号的相位
    
88. S3 = HScroll7.Value * PAI / 180
    
89. Label18.Caption = HScroll7.Value
    
90. End Sub
    
91. 
    
92. Private Sub HScroll8_Change()
    
93. '合成谐波信号的频率
    
94. A3 = HScroll8.Value / 2
    
95. Label19.Caption = A3
    
96. End Sub
    
97. 
    
98. Private Sub Timer1_Timer()
    
99. Call shuju
    
100. P_POINT = I_TIME Mod 256
     
101. FUZHI(P_POINT) = TU_DATA(I_TIME) / 100
     
102. If I_TIME > 255 Then Call FFT 'begin FFT
     
103. If P_START = 0 Then Picture1(0).PSet (0, 1200 - TU_DATA(I_TIME))
     
104. If P_START <= 256 Then Picture1(0).Line -(P_START * 25, 1200 - TU_DATA(I_TIME)), RGB(128, 255, 2550)
     
105. If P_START > 256 Then P_START = 256: Call pic_draw
     
106. I_TIME = (I_TIME + 1) Mod 2048
     
107. P_START = P_START + 1
     
108. End Sub
     
109. Private Sub FFT()
     
110. Dim H(255), G(255) As Double
     
111. 
     
112. '将采集数据按倒序结果赋值
     
113. For i = 0 To 255
     
114. kk = (Shuxu(i) + P_POINT + 1) Mod 256
     
115. PINZHIr(i) = FUZHI(kk)
     
116. PINZHIi(i) = 0
     
117. Next i
     
118. 
     
119. 'FFT
     
120. For k = 0 To 7               '变换次数：8次*256
     
121. Ai = 2 ^ k
     
122. For i = 0 To 2 ^ (7 - k) - 1 '分组的组个数：128、64、32、16、8、 4、 2、  0
     
123. Bi = 2 * i * Ai
     
124. For j = 0 To 2 ^ (k + 1) - 1 '每组数据个数：2、  4、  8、16、32、64、128、256
     
125. Ci = Bi + j
     
126. If j < Ai Then
     
127.    H(Ci) = PINZHIr(Ci) + fr(k, j) * PINZHIr(Ci + Ai) - fi(k, j) * PINZHIi(Ci + Ai)
     
128.    G(Ci) = PINZHIi(Ci) + fi(k, j) * PINZHIr(Ci + Ai) + fr(k, j) * PINZHIi(Ci + Ai)
     
129. Else
     
130.    H(Ci) = PINZHIr(Ci - Ai) + fr(k, j) * PINZHIr(Ci) - fi(k, j) * PINZHIi(Ci)
     
131.    G(Ci) = PINZHIi(Ci - Ai) + fr(k, j) * PINZHIi(Ci) + fi(k, j) * PINZHIr(Ci)
     
132. End If
     
133. Next j
     
134. Next i
     
135. For w = 0 To 255
     
136. PINZHIr(w) = H(w)
     
137. PINZHIi(w) = G(w)
     
138. Next w
     
139. Next k
     
140. 
     
141. Picture1(1).Cls
     
142. Picture1(1).PSet (0, 2300)
     
143. For j = 0 To 127
     
144. PINZHI(j) = (PINZHIr(j) * PINZHIr(j) + PINZHIi(j) * PINZHIi(j)) ^ 0.5
     
145. Call ZHENG(PINZHI(j))
     
146. If TX = 0 Then Picture1(1).Line -(100 * j, 2300 - PINZHI(j) * 5), RGB(255, 0, 0)
     
147. If TX = 1 Then Picture1(1).Line (100 * j, 2300)-(100 * j, 2300 - PINZHI(j) * 5), RGB(255, 0, 0)
     
148. Next j
     
149. '横坐标点
     
150. For j = 0 To 12
     
151. Picture1(1).Circle (j * HScroll6.Value * 1000 / 256, 2300), 20, RGB(255, 0, 0)
     
152. Next j
     
153. End Sub
     
154. Private Sub pic_draw()
     
155. Picture1(0).Cls
     
156. Picture1(0).Line (0, 1200)-(6400, 1200), RGB(255, 0, 0)
     
157. Picture1(0).PSet (6400, 1200 - TU_DATA(I_TIME))
     
158. For i = 1 To 256
     
159. kk = (I_TIME - i + 2048) Mod 2048
     
160. Picture1(0).Line -(6400 - i * 25, 1200 - TU_DATA(kk)), RGB(128, 255, 2550)
     
161. Next i
     
162. End Sub
     
163. Sub shuju()
     
164. xx2 = 128 / (Int((HScroll6.Value / 256) * 100) / 100)
     
165. xx = 3 * Sin(A1 * PAI * I_TIME / xx2 + S1) + 2 * Cos(A2 * PAI * I_TIME / xx2 + S2) + Sin(A3 * PAI * I_TIME / xx2 + S3)
     
166. Call ZHENG(xx)
     
167. TU_DATA(I_TIME) = xx * 100
     
168. End Sub
     
169. Sub ZHENG(xx)
     
170. If Abs(xx) < 0.00001 Then xx = 0
     
171. xx = Int(xx * 100000 + 0.5) / 100000
     
172. End Sub

复制代码

[ 本帖最后由 风花雪月 于 2008-1-29 15:02 编辑 ]

wanyeqing2003

收藏了。

风花雪月

原帖由 dragoon0114 于 2008-1-28 15:54 发表

                               
登录/注册后可看大图

我不知道怎么上传附件，只好把程序粘贴了过来。

见发帖输入框下面的上传附件

dragoon0114

我晕啊，不能上传.frm（VB）格式的文件。

风花雪月

原帖由 dragoon0114 于 2008-2-22 17:42 发表

                               
登录/注册后可看大图

我晕啊，不能上传.frm（VB）格式的文件。

呵呵，现在文件格式多了去了，论坛只能上传特定几个扩展名的文件

宝贝鱼

你好，我想知道这段代码是做什么用的

'加权系数W
For i = 0 To 7
For j = 0 To 2 ^ i - 1                     
aa = Cos(2 * PAI * j / 2 ^ (i + 1))
Call ZHENG(aa)
fr(i, j) = aa
fr(i, j + 2 ^ i) = -aa
bb = Sin(2 * PAI * j / 2 ^ (i + 1))
Call ZHENG(bb)
fi(i, j + 2 ^ i) = bb
fi(i, j) = -bb
Next j
Next i

Sub ZHENG(xx)
If Abs(xx) < 0.00001 Then xx = 0
xx = Int(xx * 100000 + 0.5) / 100000
End Sub

还有啊，怎么确定采样点数是256呢？

babys

能不能把程序贴出来啊？
有用啊

xjyhuijia

正需要相关程序呢，分享了

letian882236

回复 4 # dragoon0114 的帖子

麻烦楼主给控件截个图呗，知道怎样安排就好做了，万分感谢！

88484532

先收藏了，谢谢提供。

zhangzihang0820

谢谢提供