From e91f13615dde8896e83294c4358623eca3189918 Mon Sep 17 00:00:00 2001 From: Shayan Date: Tue, 12 Mar 2024 01:58:12 +0330 Subject: [PATCH] =?UTF-8?q?[cepstral.jl]=20Further=20unit=20tests=20provid?= =?UTF-8?q?ed=20=F0=9F=A7=AA?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- test/cepstral.jl | 22 ++++++++++++++++++++++ 1 file changed, 22 insertions(+) diff --git a/test/cepstral.jl b/test/cepstral.jl index 90dbeec..9706327 100644 --- a/test/cepstral.jl +++ b/test/cepstral.jl @@ -1,18 +1,40 @@ series = [0.5783834024666892, 0.41692819343103993, 0.28467412058772745, 0.040996614637221374, 0.06804363629643206, 0.9244312331822693, 0.08725427297969246, 0.7808434807429435, 0.17007538615159856, 0.832499149162595] +vals = [ + 0.526449 0.253821 0.536858 0.15444 0.39956 0.278192 0.145485 0.685143 0.495732 0.596486 + 0.913692 0.205756 0.403905 0.529887 0.104882 0.180711 0.651375 0.208858 0.307472 0.277686 + 0.537414 0.870541 0.163864 0.374202 0.0393039 0.68875 0.828697 0.0172519 0.668273 0.765799 +] + @testset "cepstral.jl" begin @testset "ARMACepstral" begin @test (cc(ARMACepstral(1., (1, 1)), series, 5) .≈ [-0.6231197903689889, -0.01564051553647603, -0.005955402230438618, -0.0027120968883413143, -0.0015279908018555185]) |> all + + @test isapprox(cc(ARMACepstral(1., (1, 1)), vals, 5), [ -1.03522 -2.168 -2.28113 + -0.0172828 0.00615351 0.00264155 + -0.00445416 0.00153819 0.000682284 + -0.00197663 0.000683611 0.000305177 + -0.00111116 0.000384525 0.00017205], atol=1e-5) end @testset "ARCepstral" begin @test (cc(ARCepstral(1), series, 5) .≈ [0.5686327103717447, 0.16167157965235823, 0.06128783235186792, 0.026137699667288815, 0.011890200803754471]) |> all @test (cc(ARCepstral(2), series, 5) .≈ [0.5061405677355508, -0.08265600639423334, -0.06344592753639652, -0.015374742371503975, 0.0017971279975581714]) |> all + @test isapprox(cc(ARCepstral(1), vals, 5), [ 0.239323 0.276995 0.26709 + 0.0286378 0.0383632 0.0356686 + 0.00456913 0.00708428 0.00635116 + 0.000820125 0.00147174 0.00127225 + 0.00015702 0.000326131 0.000271845], atol=1e-5) end @testset "RealCepstral" begin @test (cc(RealCepstral(), series, 5) .≈ [0.03700275546236409, 0.0867638177224842, 0.38225344455973626, -0.03467668552262836, 0.10502390427503827]) |> all @test (cc(RealCepstral(), series, 5, normalize=true) .≈ [-2.482417812401184, -3.5121897532599697, -3.2167001264227184, -3.633630256505082, -3.493929666707416]) |> all + @test isapprox(cc(RealCepstral(), vals, 5), [ -0.65022 -0.198637 0.113368 + 0.342659 0.0145222 0.152945 + 0.341224 0.174512 0.0952615 + 0.259176 0.317715 0.246168 + -0.163841 0.252064 0.0640904], atol=1e-5) end @testset "Errors" begin