Fix start values of FIIR and Biquad filter to second sample

JuliaPhysics · Dec 15, 2024 · d8f52c9 · d8f52c9
1 parent b48c598
commit d8f52c9
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 0 deletions.
diff --git a/src/biquad_filter.jl b/src/biquad_filter.jl
@@ -81,6 +81,10 @@ end
         s1 = fma(neg_a1, y_i, z2)
         s2 = fma(neg_a2, y_i, z3)
     end
+
+    # fix start value
+    Y[1] = Y[2]
+
     Y
 end
 

diff --git a/src/first_order_iir.jl b/src/first_order_iir.jl
@@ -99,6 +99,10 @@ X::AbstractVector{U}) where {T<:Real, U<:Real}
 
         s1 = fma(neg_a1, y_i, z2)
     end
+
+    # fix start value
+    Y[1] = Y[2]
+
     Y
 end
-Original file line number
+Diff line change
@@ Expand Up / @@ -99,6 +99,10 @@ X::AbstractVector{U}) where {T<:Real, U<:Real} @@
             s1 = fma(neg_a1, y_i, z2)
         end
+        # fix start value
+        Y[1] = Y[2]
         Y
     end
@@ Expand Down @@