i need a FFT routine for bascom avr , atmega128

    This site uses cookies. By continuing to browse this site, you are agreeing to our Cookie Policy.

    • New

      Hallo wimapon,

      attached you find a FFT-Basic snippets, I hope its that what you want.

      BASCOM Source Code: FFT-Schnipsel

      1. '( FFT-Basic Parts
      2. without Assembler-Routines
      3. 15.02.2020 /17:00
      5. ')
      6. $regfile = "m644def.dat"
      7. $crystal = 20000000 '20MHz
      8. $hwstack = 100
      9. $swstack = 100
      10. $framesize = 100
      11. '$sim
      12. $baud = 56000
      14. Dim Re(130) As Single 'Real-Teil vor FFT-Berechnung
      15. Dim Im(130) As Single 'Imag-Teil vor FFT-Berechnung
      16. Dim Real(130) As Word 'Betragswerte nach FFT-Berechnung
      17. Dim Wcos(130) As Single 'Cosinus Werte für FFT
      18. Dim Wsin(130) As Single 'Sinus Werte für FFT
      19. Dim Pi As Single
      20. Dim Mr As Integer
      21. Dim N As Word 'FFT-Punkteanzahl
      22. Dim Nn As Integer
      23. Dim L As Integer
      24. Dim M As Integer
      25. Dim Tr As Single
      26. Dim Ti As Single
      27. Dim Istep As Integer
      28. Dim A As Single
      29. Dim Wr As Single
      30. Dim Wi As Single
      31. Dim I As Integer
      32. Dim J As Integer
      33. Dim Tr0 As Single
      34. Dim Tr1 As Single
      35. Dim Ti0 As Single
      36. Dim Ti1 As Single
      37. Dim Nm As Integer
      38. Dim W_cond As Integer
      39. Dim Ptr As Integer
      40. Dim Re1 As Single
      41. Dim Im1 As Single
      42. Dim Re_betrag As Single
      43. Dim N2 As Byte 'FFT-Ergebnis, N-halbe Slots
      45. Dim Freq As Word 'Frequenzvorgabe für Testsinus
      46. Dim Abtast_dauer As Single 'Abtastdauer für Testsinus
      47. Dim Sin_input(128) As Single 'Speicherplatz für Testsinus
      48. Dim F_steps As Single 'Frequenzsteps für Testsinus
      49. Dim Step_radian As Single 'Steps in Rad
      50. Dim Schritt As Single 'Schrittweite für Testsinus
      51. Dim Sinuswert As Single
      53. Freq = 313 'Frequenz-Vorgabe für Test-Sinus-Berechnung in Hz
      54. 'Print "with Simulator, wait a while, calculating results"
      56. Declare Sub Fft
      58. Pi = 4 * Atn(1)
      59. N = 128 'Anzahl Punkte der FFT
      60. N2 = N / 2 'N-halbe
      62. '###################################################################
      63. '---------------only for Test-Sinus ----------------
      64. Schritt = 0
      65. Abtast_dauer = 0.000125 'Abtast-Takt, hier 125usec
      66. F_steps = 1 / Freq 'Anzahl Steps pro Schwingung mit Abtastdauer
      67. F_steps = F_steps / Abtast_dauer
      68. Step_radian = 2 * Pi 'Anzahl Steps pro Schwingung bezogen auf 1 Radian
      69. Step_radian = Step_radian / F_steps
      71. 'Print "Step_radian=" ; Step_radian ; " F_steps=" ; F_steps
      73. '---------------Test-Sinus ----------------
      74. For I = 1 To N 'Generierung von Sinus-Werten zum Testen
      75. Sinuswert = Sin(schritt)
      76. Sin_input(i) = 400 * Sinuswert
      77. Schritt = Schritt + Step_radian
      78. '-------------------------
      79. Re(i) = Sin_input(i) 'Realpart, here filled with Test Sinus Values
      80. Im(i) = 0
      81. '-------------------------
      82. 'Print "I=";i;" Re =";Re(i)
      83. Next
      84. '####################################################################
      88. ' ------ Start -----------
      89. 'Vorgelagerte Cos- und Sin- Berechnung
      90. 'may definied as a Sub
      91. L = 1
      92. Nm = N
      93. While L < N '
      94. Istep = 2 * L
      95. For M = 1 To L
      96. Ptr = M + L
      97. A = 1 - M
      98. A = A / L
      99. A = A * Pi
      100. A = A * -1 'A=-Pi*(1-M)/L
      101. Wcos(ptr) = Cos(a)
      102. Wsin(ptr) = Sin(a)
      103. Wsin(ptr) = Wsin(ptr) * -1 '-sin(a)
      104. 'Print "Ptr=" ; Ptr ; " Cos=" ; Wcos(ptr) ; " Sin=" ; Wsin(ptr)
      105. Next M
      106. L = Istep
      107. Wend
      109. Fft
      111. '----FFT-Results with Test Sinus
      112. For I = 1 To N2 '
      113. Re1 = Re(i) * Re(i) 'Re(i) zum Quadrat
      114. Im1 = Im(i) * Im(i)
      115. Re_betrag = Re1 + Im1
      116. Re_betrag = Sqr(re_betrag)
      117. Re_betrag = Re_betrag / 4 'reduzieren, da Wert > als 65536 sein kann
      118. Real(i) = Re_betrag 'Betragswert
      119. Print I ; " " ; Real(i) 'Betragswert-Ausgabe der einzelnen FFT-Slots
      120. Next I
      121. Real(1) = 0 '#####!!!!Gleichanteil zu Null gesetzt!
      123. '######################################
      124. End
      127. Sub Fft
      128. Mr = 0
      129. Nn = N - 1
      130. For M = 1 To Nn
      131. L = N
      132. L = L / 2 'Shift L , Right , 1 ' = L / 2
      133. W_cond = Mr + L
      134. While W_cond > Nn
      135. L = L / 2 'Shift L , Right , 1 ' = L / 2
      136. W_cond = Mr + L
      137. Wend
      138. Mr = Mr Mod L 'Mod_argu
      139. Mr = Mr + L
      140. 'Mr = Mr
      141. If Mr > M Then
      142. Tr = Re(m + 1)
      143. Re(m + 1) = Re(mr + 1)
      144. Re(mr + 1) = Tr
      145. Ti = Im(m + 1)
      146. Im(m + 1) = Im(mr + 1)
      147. Im(mr + 1) = Ti
      148. End If
      149. Next M
      150. L = 1
      151. Nm = N
      152. While L < N '
      153. Istep = 2 * L
      154. For M = 1 To L
      155. Ptr = M + L
      156. Wr = Wcos(ptr) 'cos(a)
      157. Wi = Wsin(ptr) 'Sin(a)
      158. For I = M To N Step Istep
      159. J = I + L
      160. Tr0 = Wr * Re(j)
      161. Tr1 = Wi * Im(j)
      162. Tr = Tr0 - Tr1 'ZR=CO*RE(J) - SI*IM(J)
      163. Ti0 = Wr * Im(j)
      164. Ti1 = Wi * Re(j)
      165. Ti = Ti0 + Ti1 'ZI=CO*IM(J) + SI*RE(J)
      166. Re(j) = Re(i) - Tr
      167. Im(j) = Im(i) - Ti
      168. Re(i) = Re(i) + Tr
      169. Im(i) = Im(i) + Ti
      170. Next I
      171. Next M
      172. L = Istep
      173. Wend
      174. End Sub
      Display All

    • New

      Hallo Stefanhamburg,


      noch einige Gedanken zu deinem Bannstrahl bzgl. der Verwendung der englischen Sprache von Hilfe-Suchenden in diesem Forum:

      • anscheinend wird dieses Forum als hilfreicher angesehen!?
      • Bascom sowie Bascom-Hilfe und auch alle Datenblätter zu den AVR’s sind ohne Englisch-Kenntnisse kaum sinnvoll anwendbar.
      • Es wird sich somit niemand mit AVR’s beschäftigen wollen, der nicht auch (minimale) Sprachkenntnisse in Englisch besitzt.
      • Die Problemanfragen in deutscher Sprache erfordern schon oftmals Nachfragen vom Tenor „…was ist damit gemeint…“. Übersetzer-Algorithmen sind eher nicht auf technische Belange optimiert, sodass noch mehr Verwirrung bei der Hin- und Rückübersetzung entstehen könnte.
      • Im Sprachgebrauch hier werden ebenfalls Begriffe wie „Thread“ verwendet. Eine strenge deutsche Übersetzung könnte auch „Faden“ bedeuten. (nach deepl-Translator)
      • Wie groß ist der Anteil der in Englisch gestellten Anfragen bezogen auf die Fülle derer in deutscher Sprache, ich vermute mal, er liegt im niedrigen einstelligen Prozent-Bereich?

      Auch ich bin ein Verfechter meiner Muttersprache und ärgere mich jedes Mal, wenn gehäuft Anglizismen im Blätterwald auftauchen. Allerdings sehe ich es bei technischen Bereichen und Anwendungen differenzierter, da in diesem Bereich deutsch leider kaum zur Anwendung kommt.


      Deshalb würde ich solche in Englisch notierten Anfragen etwas weniger streng betrachten; eher bereichert es dieses ausgezeichnete Forum.

      In einem Buch von 2010 las ich auf Seite 14 u.a. „Ziel ist es dem Anfänger die Scheu vor dem Einstieg zu nehmen“. Das kann ich voll und ganz unterstützen, wenn ich auch auf meine Anfänge zurückblicke.


      Diese Gedanken sind nur meine persönliche Meinung, andere Meinungen werden aber auch akzeptiert.


      Gruß

      Ulrich

    • New

      hallo wieder Ulrich,
      Ich kan ein bischen deutch sprechen, aber die deutche computer worten kan ich nicht verstehen.
      Darom versuche ich es auf englisch.......

      Ich will jedenfalls kein probleem sein auf diezes forum

      I think that your last snipet does the job for me.
      I did some experiments with it , and the results looks great.

      When i am ready and i can use this program for my other experiments i will inform you here.

      Until sofar, thank you 1000 times for your wonderfull help

      Wim PA0SLT

      PS If you like to know who you helped, look at my sites:
      home.kpn.nl/huiu55wi/huis.htm
      ebl21.nl

      at the moment i am trying to make a seismic earthquakes measuring station here at my house.
      We do suffer lots of earthquakes here due to gas production.
      The problem is that my sensor gets to much 50 Hz noise.
      Because i need a greater sensitivity i need a FFT routine to see the difference between an
      ordinary 50 Hz noise puls and a real earthquake.

    • New

      Es ist kein Problem, wenn hier einer auf Englisch fragt.
      Keiner hat hier was dagegen und die meisten können es auch lesen.

      Das Problem ist, dass die Antworten der User in einem deutschen Forum nur von den Leuten kommen, die auch englisch schreiben (können)
      Das sind nicht unbedingt die gleichen Leute, die eine Lösung haben.

      Der Fragende bekommt also nicht die größtmögliche Aufmerksamkeit und auch nicht die Spezialisten.
      Das ist in einem englisch-sprachigem Forum eben anders.
      Und genau das hat Stefanhamburg gemeint.
      Kein Bannstrahl a_45_132ca9f5

    • New

      Michael wrote:

      Der Fragende bekommt also nicht die größtmögliche Aufmerksamkeit und auch nicht die Spezialisten.
      Das ist in einem englisch-sprachigem Forum eben anders.
      scheint eben nicht so zu sein, sonst würde er es nicht hier probieren. Ist doch eigentlich klar, dass ich in das Forum gehe, das ich am besten lesen/verstehen kann. Nur wenn ich da nicht weiter komme, dann mach ich mir die Mühe in einem anderssprachigen Forum. Ich sehe die englischen Anfragen eher als ein Anzeichen für die Qualität dieses Forums.
      Raum für Notizen

      -----------------------------------------------------------------------------------------------------

      -----------------------------------------------------------------------------------------------------

    • New

      Hallo wimapon,

      according to your homemade seismic earthquakes measuring station: did you built the sensor(s) yourself? It is possible you give some information about the sensors? How to built? schematics?

      In my vicinity there have also been drillings for gas, which have generated light earthquakes too.

      Ulrich

    • New

      Hi Ulrich,
      As sensor i do use a LIS332AR 3-axis +/-2g analog-output ultracompact accelerometer.
      After that i du use a low frequency amplifier with a LM324

      I do have much noise coming from this thing... so i have difficulties to see the difference
      between an week earthquake and the continuous noise.
      Therfore i am now experimenting with the FFT in the AVR program.

      In the begining i did use an ordinary loudspeaker with a stone on the conus.
      that worked very good. the stone must not be to heavy... the conus must
      be somewhere in the middle of the movement range.....
      The output of that loudspeaker was connected to the ADconvertor input of the AVR.

      I am a retired geologist. I was working in oil and gas exploration with seismic interpretation.
      So i do know something about earth quakes and seismics.

      If i can help you, i will do!

      Wim

    • New

      Hallo Wimapon,

      many thanks for that informations. I have already read the sensor example using a loudspeaker in one of the popular electronic magazines. I will test that too. Do the noise come from the accelerator or from the LM324?

      Here in this Bascom-Forum many users have shown programs to communicate with different 3/6-axis accelerators, e.g. MPU6050.

      MPU-6050 MEMS 6-Achs-Sensor

      May be it is of interest too. Search for this theme. Everything else in a PN.

      Ulrich