From 4f72924cb839175788717d6631d56898215875f0 Mon Sep 17 00:00:00 2001
From: seyffejn <55194230+seyffejn@users.noreply.github.com>
Date: Fri, 19 Jun 2020 15:58:09 +0200
Subject: [PATCH] Update NEU: (DigitalFilterTest.py)
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
DigitalFilterTest.py hinzugefügt.
Dieses Skript beinhaltet die Differenzengleichung eines Bandpass nach Tobola MEM_C Vorlesung. Außerdem plottet es schnell die Frequenzen um zu testen ob der Filter funktioniert.
---
TinnitusAnalyse/.idea/workspace.xml | 51 +++++----
TinnitusAnalyse/DigitalFilterTest.py | 98 ++++++++++++++++++
TinnitusAnalyse/SoundGenerator.py | 2 +-
TinnitusAnalyse/TinnitusDaten.csv | 14 +--
.../__pycache__/SoundGenerator.cpython-35.pyc | Bin 4190 -> 4222 bytes
5 files changed, 131 insertions(+), 34 deletions(-)
create mode 100644 TinnitusAnalyse/DigitalFilterTest.py
diff --git a/TinnitusAnalyse/.idea/workspace.xml b/TinnitusAnalyse/.idea/workspace.xml
index 64dc4b2..42a6242 100644
--- a/TinnitusAnalyse/.idea/workspace.xml
+++ b/TinnitusAnalyse/.idea/workspace.xml
@@ -2,7 +2,6 @@
-
@@ -14,11 +13,11 @@
-
+
-
-
+
+
@@ -35,11 +34,11 @@
-
+
-
-
+
+
@@ -50,8 +49,8 @@
-
-
+
+
@@ -127,7 +126,7 @@
-
+
@@ -163,7 +162,7 @@
-
+
@@ -271,8 +270,8 @@
-
\ No newline at end of file
diff --git a/TinnitusAnalyse/DigitalFilterTest.py b/TinnitusAnalyse/DigitalFilterTest.py
new file mode 100644
index 0000000..9927fb8
--- /dev/null
+++ b/TinnitusAnalyse/DigitalFilterTest.py
@@ -0,0 +1,98 @@
+import matplotlib.pyplot as plt # For plotting
+from math import sin, pi, cos # For generating input signals
+import numpy as np
+import sys # For reading command line arguments
+
+
+fs = 44100 # sampling frequency (Abtastfrequenz)
+
+koeff = {
+ "b0": 0,
+ "b1": 0,
+ "b2": 0,
+ "a1": 0,
+ "a2": 0
+}
+
+
+def koeffizienten_berechnen(omega, r):
+ # Koeffizientenberechnung nach Tobola VL S.107 - IIR Filter, 2. Ordnung
+ koeff["b0"] = 1
+ koeff["b1"] = 0
+ koeff["b2"] = -1
+ koeff["a1"] = -2*r*cos(omega)
+ koeff["a2"] = r**2
+
+
+def filter(x):
+ y = [0]*len(x)
+ for k in range(4, len(x)):
+ y[k] = koeff["b0"]*x[k] + koeff["b1"]*x[k-1] + koeff["b2"]*x[k-2] - koeff["a1"]*y[k-1] - koeff["a2"]*y[k-2]
+
+ return y
+
+
+
+
+dauer_ms = 10000 # 10 Sekunden
+num_samples = dauer_ms * (fs / 1000) # framerate -pro Sekunde- umgerechnet in -pro Millisekunde-
+
+t = np.linspace(0, 10, int(num_samples)) # array zum darstellen der x-Achse
+amp = 1
+
+f_input1 = 1
+f_input2 = 50
+f_input3 = 10
+
+input1 = []
+input2 = []
+input3 = []
+for x in range(int(num_samples)): # einen einfachen Sinus ins array schreiben
+ input1.append(amp * sin(2 * pi * f_input1 * (x / fs)))
+ input2.append(amp * sin(2 * pi * f_input2 * (x / fs)))
+ #input3.append(amp * sin(2 * pi * f_input3 * (x / fs)))
+
+
+input_ges = np.add(input1, input2) # Sinus aufaddieren um ein halbwegs realistisches Audiosignal zu bekommen
+#input_ges = np.add(input_ges, input3)
+
+#Filterparameter hier einstellen
+fr = 5
+omega = 2*pi*fr/fs
+r = 0.95
+koeffizienten_berechnen(omega, r) # Koeffizienten berechnen mit Mittelfrequenz 10Hz
+
+output = filter(input_ges)
+
+### Plot the signals for comparison
+plt.figure(1)
+plt.subplot(231)
+plt.ylabel('Amplitude')
+plt.xlabel('t [s]')
+plt.title('f_1=' + str(f_input1) + "Hz")
+plt.plot(t, input1)
+
+plt.subplot(232)
+plt.ylabel('Amplitude')
+plt.xlabel('t [s]')
+plt.title('f_2=' + str(f_input2) + "Hz")
+plt.plot(t, input2)
+
+# plt.subplot(233)
+# plt.ylabel('Amplitude')
+# plt.xlabel('t [s]')
+# plt.title('f_3=' + str(f_input3) + "Hz")
+# plt.plot(t, input3)
+
+plt.subplot(234)
+plt.ylabel('Amplitude')
+plt.xlabel('t [s]')
+plt.title('(input in Filter) f_ges = f_1 + f_2')
+plt.plot(t, input_ges)
+
+plt.subplot(235)
+plt.ylabel('Amplitude')
+plt.xlabel('Samples')
+plt.title('gefiltertes Signal, mit f_r = ' + str(fr) + "Hz und r=" + str(r))
+plt.plot(t, output)
+plt.show()
\ No newline at end of file
diff --git a/TinnitusAnalyse/SoundGenerator.py b/TinnitusAnalyse/SoundGenerator.py
index df274fb..073b0e8 100644
--- a/TinnitusAnalyse/SoundGenerator.py
+++ b/TinnitusAnalyse/SoundGenerator.py
@@ -74,7 +74,7 @@ class Sound:
# zuerst muss ein Array mit Audiodaten gefüllt werden
audio = []
freq = self.tinnitus.linksFrequenz
- dauer_ms = 1000.0 # 10 Sekunden
+ dauer_ms = 5000.0 # 10 Sekunden
amp = self.tinnitus.linksLautstaerke
rauschen = self.tinnitus.linksRauschenLautstaerke
diff --git a/TinnitusAnalyse/TinnitusDaten.csv b/TinnitusAnalyse/TinnitusDaten.csv
index b1ca267..8aa79e9 100644
--- a/TinnitusAnalyse/TinnitusDaten.csv
+++ b/TinnitusAnalyse/TinnitusDaten.csv
@@ -1,10 +1,10 @@
-Vorname;asd
-Nachname;asda
-linke Frequenz;770.0
-linke Lautstärke;0.01
+Vorname;asdas
+Nachname;asdas
+linke Frequenz;0.0
+linke Lautstärke;0.0
linkes Rauschen;0.0
-rechte Frequenz;1410.0
-rechte Lautstärke;0.03
+rechte Frequenz;0.0
+rechte Lautstärke;0.0
rechtes Rauschen;0.0
-Kommentar;asd
+Kommentar;asdas
diff --git a/TinnitusAnalyse/__pycache__/SoundGenerator.cpython-35.pyc b/TinnitusAnalyse/__pycache__/SoundGenerator.cpython-35.pyc
index b782ce8a339a27ab313057a703562bda4f2402ec..de373e603f168284346468ecf8db122b54c7cdf8 100644
GIT binary patch
delta 322
zcmYk1%}T>S5XWaWAI+{C5=ljbfGMG15fQNp-g=Q91P=-pBIsp7kyJ3|SSlsHfDArF
zk7a!W-@}VHKVCuHVb*-=~H^NqW}~_
z?m@<3s2C^<9u2afVvq@gY*0>voCYl|=!{{gKiG_6Ot~|~Kz};V4EJc}&C#Y!zKsr@
z^AWnF`BlFKrGNSAO?EX9tAE~6s2jQe$0nj|s6m1WtHwsebH-jnnl?0Sc$
T%vd?5^>C{lHbRMi?nT{SF55W(
delta 253
zcmeyTa8E&5oR^nNd*7S557G<_j~S2vGmz~7#KlDumAzHnQF6aKx8u$Llg&t
zdkPCf3M*KM6)eQLaef#RBmd@aOpQ#8(vvG#>lpbqe`lS_#Aq_Pmt!R(^JH^QeMXka
z>73z=tdsY1rgO0|lyEWBa5AK@Pd4C^