dsp intro wip

Author: daslu

src/signal_processing/intro.clj

@@ -16,46 +16,17 @@
             [scicloj.tableplot.v1.plotly :as plotly]))
 
 
-(defn sine-samples [{:keys [sample-rate duration frequency amplitude]
-                     :or {amplitude 0.5}}]
-  (let [num-samples (int (* sample-rate duration))
-        ;; Generate samples as int16 reader
-        ]
-    (for [idx (range num-samples)]
-      (let [angle (* 2.0 Math/PI idx frequency (/ 1.0 sample-rate))]
-        (* amplitude (Math/sin angle))))))
-
-
-(defn sample-info [samples sample-rate]
-  (with-meta
-    {:samples samples
-     :sample-rate sample-rate}
-    {:kind/audio true}))
-
-
-(let [sample-rate 44100.0
-      A 440.0
-      E 659.25
-      wave (fn [freq]
-             (sine-samples
-              {:sample-rate sample-rate
-               :duration 3
-               :frequency freq
-               :amplitude 0.3}))]
-  (sample-info (map +
-                    (wave A)
-                    (wave E))
-               sample-rate))
-
-
-(require '[tech.v3.datatype :as dtype]
+(require '[scicloj.kindly.v4.kind :as kind]
+         '[tech.v3.datatype :as dtype]
          '[tech.v3.datatype.functional :as dfn]
          '[clojure.math :as math]
          '[tablecloth.api :as tc]
          '[scicloj.tableplot.v1.plotly :as plotly])
 
 (def sample-rate 44100.0)
 
+
+
 (def example-wave
   (let [duration 10
         num-samples (* duration sample-rate)
@@ -64,45 +35,53 @@
                                 (/ idx sample-rate))
         freq 440
         amp 3800
-        amplitude (-> time
-                      (dfn/* (* 2 Math/PI freq))
-                      dfn/sin
-                      (dfn/* amp))]
+        value (-> time
+                  (dfn/* (* 2 Math/PI freq))
+                  dfn/sin
+                  (dfn/* amp))]
     (tc/dataset {:time time
-                 :amplitude amplitude})))
+                 :value value})))
 
 example-wave
 
 (-> example-wave
     (tc/head 200)
     (plotly/layer-line {:=x :time
-                        :=y :amplitude}))
+                        :=y :value}))
 
+(defn audio [samples]
+  (with-meta
+    {:samples samples
+     :sample-rate sample-rate}
+    {:kind/audio true}))
+
+(-> example-wave
+    :value
+    audio)
 
 (def violin-components
   [[:A4 440 3800]
    [:A5 880 2750]
-   [:e6 1320 600]
-   [:a6 1760 700]
-   [:c#7 2200 1900]])
-
-(DEF VIOLIN-COMPONENTS-DATASET
-  (LET [DURATION 10
-        NUM-SAMPLES (* DURATION SAMPLE-RATE)
-        TIME (DTYPE/MAKE-READER :FLOAT32
-                                NUM-SAMPLES
-                                (/ IDX SAMPLE-RATE))]
-       (->> VIOLIN-COMPONENTS
-            (MAP (FN [[LABEL FREQ AMP]]
-                     [LABEL (-> TIME
-                                (DFN/* (* 2 mATH/pi FREQ))
-                                DFN/SIN
-                                (DFN/* AMP))]))
-            (INTO {:TIME TIME})
-            TC/DATASET)))
+   [:E6 1320 600]
+   [:A6 1760 700]
+   [:C#7 2200 1900]])
 
-violin-components-dataset
+(def violin-components-dataset
+  (let [duration 10
+        num-samples (* duration sample-rate)
+        time (dtype/make-reader :float32
+                                num-samples
+                                (/ idx sample-rate))]
+    (->> violin-components
+         (map (fn [[label freq amp]]
+                [label (-> time
+                           (dfn/* (* 2 math/PI freq))
+                           dfn/sin
+                           (dfn/* amp))]))
+         (into {:time time})
+         tc/dataset)))
 
+violin-components-dataset
 
 (-> violin-components-dataset
     (tc/head 200)
@@ -117,9 +96,9 @@ violin-components-dataset
 (-> violin-components-dataset
     (tc/head 200)
     (tc/pivot->longer (complement #{:time}))
-    (tc/rename-columns {:$value :amplitude})
+    (tc/rename-columns {:$value :value})
     (plotly/layer-line {:=x :time
-                        :=y :amplitude
+                        :=y :value
                         :=color :$column}))
 
 (def violin-dataset
@@ -131,15 +110,16 @@ violin-components-dataset
                              (:A6 %)
                              (:C#7 %)))))
 
-
 (-> violin-dataset
     (tc/head 200)
     (plotly/layer-line {:=x :time
                         :=y :violin}))
 
+(-> violin-dataset
+    :violin
+    (dfn// 7000.0)
+    audio)
+
 
-(sample-info (dfn// (:violin violin-dataset)
-                    7000.0)
-             sample-rate)