Skip to content

Commit

Permalink
Merge pull request #40 from dsweber2/master
Browse files Browse the repository at this point in the history
docs: fixes for recent changes
  • Loading branch information
dsweber2 authored Aug 27, 2023
2 parents 8d544ec + 30a9a90 commit cdeefa5
Show file tree
Hide file tree
Showing 2 changed files with 19 additions and 5 deletions.
7 changes: 5 additions & 2 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -37,10 +37,12 @@ Plot{Plots.PyPlotBackend() n=1}

julia> c = wavelet(Morlet(π), β=2)

CWT{Morlet mean 3.141592653589793, Father Wavelet, Q=8.0, β=2.0,aveLen=0.0, frame=1.0, norm=Inf, extraOctaves=0.0}

julia> res = ContinuousWavelets.cwt(f, c)
┌ Warning: the lowest frequency wavelet has more than 1% its max at zero, so it may not be analytic. Think carefully
│ lowAprxAnalyt = 0.06186323501016359
└ @ ContinuousWavelets ~/work/ContinuousWavelets.jl/ContinuousWavelets.jl/src/sanityChecks.jl:6
│ lowAprxAnalyt = 0.061863
└ @ ContinuousWavelets ~/work/ContinuousWavelets.jl/ContinuousWavelets.jl/src/sanityChecks.jl:7
2047×31 Matrix{ComplexF64}:
-1.48637e-6+3.8241e-19im 0.000109978+9.67834e-5im
-1.48602e-6+5.15534e-19im -8.24922e-5+0.000130656im
Expand Down Expand Up @@ -90,6 +92,7 @@ julia> exs = cat(testfunction(n, "Doppler"), testfunction(n, "Blocks"), testfunc

julia> c = wavelet(cDb2, β=2, extraOctaves=-0)

CWT{Continuous db2, Father Wavelet, Q=8.0, β=2.0,aveLen=0.0, frame=1.0, norm=Inf, extraOctaves=0.0}

julia> res = circshift(ContinuousWavelets.cwt(exs, c), (0, 1, 0))
┌ Warning: the highest frequency wavelet has more than 1% its max at the end, so it may not be analytic. Think carefully
Expand Down
17 changes: 14 additions & 3 deletions docs/src/README.md
Original file line number Diff line number Diff line change
@@ -1,3 +1,12 @@
```@meta ex
DocTestFilters = [
r"\@ ContinuousWavelets .*",
r"[ +-][0-9]\.[0-9]{3,5}e-1[5-9]",
r"[ +-][0-9]\.[0-9]{3,5}e-[2-9][0-9]",
r"im {2,7}",
]
```

# ContinuousWavelets

[![Build Status](https://travis-ci.com/dsweber2/ContinuousWavelets.jl.svg?branch=master)](https://travis-ci.com/dsweber2/ContinuousWavelets.jl)
Expand Down Expand Up @@ -38,18 +47,18 @@ julia> f = testfunction(n, "Doppler");
julia> c = wavelet(Morlet(π), β=2)
CWT{Morlet mean 3.141592653589793, Father Wavelet, Q=8.0, β=2.0,aveLen=0.0, frame=1.0, norm=Inf, extraOctaves=0.0}
julia> res = ContinuousWavelets.cwt(f, c)
┌ Warning: the lowest frequency wavelet has more than 1% its max at zero, so it may not be analytic. Think carefully
│ lowAprxAnalyt = 0.061863
└ @ ContinuousWavelets ~/work/ContinuousWavelets.jl/ContinuousWavelets.jl/src/sanityChecks.jl:6
└ @ ContinuousWavelets ~/work/ContinuousWavelets.jl/ContinuousWavelets.jl/src/sanityChecks.jl:7
2047×31 Matrix{ComplexF64}:
-1.48637e-6+3.8241e-19im … 0.000109978+9.67834e-5im
-1.48602e-6+5.15534e-19im -8.24922e-5+0.000130656im
⋮ ⋱ ⋮
0.000435175+2.30636e-19im … -2.47195e-6-1.97048e-8im
0.000435027-8.28725e-19im -2.63499e-6+4.62331e-8im
```

As the cwt frame is redundant, there are many choices of dual/inverse frames. There are three available in this package, `NaiveDelta()`, `PenroseDelta()`, and `DualFrames()`. As a toy example, lets knock out the middle time of the bumps function and apply a high pass filter:
Expand Down Expand Up @@ -104,6 +113,8 @@ julia> exs = cat(testfunction(n, "Doppler"), testfunction(n, "Blocks"), testfunc
julia> c = wavelet(cDb2, β=2, extraOctaves=-0)
CWT{Continuous db2, Father Wavelet, Q=8.0, β=2.0,aveLen=0.0, frame=1.0, norm=Inf, extraOctaves=0.0}
julia> res = circshift(ContinuousWavelets.cwt(exs, c), (0, 1, 0))
┌ Warning: the highest frequency wavelet has more than 1% its max at the end, so it may not be analytic. Think carefully
│ highAprxAnalyt = 0.26778
Expand Down

0 comments on commit cdeefa5

Please sign in to comment.