Changeset 58 for birds/specwhiten.m

Timestamp:

06/14/13 17:31:23 (11 years ago)

Author:

vcheplygina

Message:

 

File:

• birds/specwhiten.m (modified) (3 diffs)

birds/specwhiten.m

@@ -3 +3 @@
 function s2 = whiten(s,fmax)
 
-%Values output by spectrogram are already in [0, 1] range
+%Not sure about this, but the paper seems to assume this
+s = abs(s); 
+
+%Values output by spectrogram are already in [0, 1] range?
+
+
 
 %Square root of spectrogram (of the absolute value?)
-s1 = sqrt(abs(s)); 
-%s1 = sqrt(s);
+
+%s1 = sqrt(abs(s)); %Version 1
+s1 = s;             %Version for segmentation
 
 %For each frame, find a quantity similar to the energy"
@@ -21 +27 @@
 
 lowestix = find(energy<threshold);
-lowestenergy = energy(lowestix);
+%lowestenergy = energy(lowestix);
 
 
@@ -27 +33 @@
 s2 = s1;
 for f=1:Nf
-    P(f) = sqrt(eps + sum(s1(f,lowestix).^2));
+    %P(f) = sqrt(eps + sum(s1(f,lowestix).^2));
+    
+    P(f) = sum(s1(f,lowestix)) / length(lowestix);
     s2(f,:) = s1(f,:) ./P(f);
 end