表示阶数为 2 的 Daubechies (多贝西)小波.
表示阶数为 n 的 Daubechies 小波.
DaubechiesWavelet
表示阶数为 2 的 Daubechies (多贝西)小波.
表示阶数为 n 的 Daubechies 小波.
更多信息
- DaubechiesWavelet 定义了一个正交小波族.
- DaubechiesWavelet[n] 对任意正整数 n 进行定义.
- 尺度函数 (
) 和小波函数 (
) 具有紧支集长度为 2n. 尺度函数具有 n 个消失矩. - DaubechiesWavelet 可以与函数 DiscreteWaveletTransform、WaveletPhi 等一起使用.
范例
打开所有单元 关闭所有单元基本范例 (3)
Plot[WaveletPhi[DaubechiesWavelet[4], x], {x, 0, 7}]WaveletPhi[DaubechiesWavelet[4], x]Plot[WaveletPsi[DaubechiesWavelet[4], x], {x, -3, 4}, PlotRange -> All]WaveletPsi[DaubechiesWavelet[4], x]WaveletFilterCoefficients[DaubechiesWavelet[2], {"PrimalLowpass", "PrimalHighpass"}, WorkingPrecision -> ∞]范围 (14)
基本用途 (8)
WaveletFilterCoefficients[DaubechiesWavelet[2], "PrimalLowpass"]WaveletFilterCoefficients[DaubechiesWavelet[2], "PrimalHighpass"]WaveletFilterCoefficients[DaubechiesWavelet[2], "LiftingFilter", WorkingPrecision -> ∞]%[{"LiftingMatrixForm", z}]lf = WaveletFilterCoefficients[DaubechiesWavelet[2], "LiftingFilter", WorkingPrecision -> ∞]lf["ForwardLiftingFunction"][Range[8]]//Simplifylf["InverseLiftingFunction"][%]//SimplifyPlot[WaveletPhi[DaubechiesWavelet[2], x], {x, 0, 3}]Plot[WaveletPhi[DaubechiesWavelet[6], x], {x, 0, 11}, PlotRange -> All]Table[Plot[WaveletPhi[DaubechiesWavelet[2], x, MaxRecursion -> i], {x, 0, 3}, PlotLabel -> i], {i, 1, 8, 2}]Plot[WaveletPsi[DaubechiesWavelet[2], x], {x, -1, 2}]阶数为6的 DaubechiesWavelet:
Plot[WaveletPsi[DaubechiesWavelet[6], x], {x, -5, 6}, PlotRange -> All]Table[Plot[WaveletPsi[DaubechiesWavelet[2], x, MaxRecursion -> i], {x, -1, 2}, PlotLabel -> i, PlotRange -> All], {i, 1, 8, 2}]小波变换 (5)
计算一个 DiscreteWaveletTransform:
data = Table[Sinc[t^2], {t, -3 π, 3 π, (6 π/1023)}];ListLinePlot[data, PlotRange -> All]dwt = DiscreteWaveletTransform[data, DaubechiesWavelet[3], 2]dwt["TreeView"]dwt["Dimensions"]WaveletListPlot[dwt, PlotLayout -> "CommonXAxis"]计算一个 DiscreteWaveletPacketTransform:
data = Table[Sinc[t^2], {t, -3 π, 3 π, (6 π/1023)}];dwpt = DiscreteWaveletPacketTransform[data, DaubechiesWavelet[4], 2]dwpt["TreeView"]dwpt["Dimensions"]WaveletListPlot[dwpt, PlotLayout -> "CommonXAxis"]计算一个 StationaryWaveletTransform:
data = Table[Sinc[t^2], {t, -3 π, 3 π, (6 π/1023)}];swt = StationaryWaveletTransform[data, DaubechiesWavelet[4], 2];swt["TreeView"]swt["Dimensions"]WaveletListPlot[swt, PlotLayout -> "CommonXAxis"]计算一个 StationaryWaveletPacketTransform:
data = Table[Sinc[t^2], {t, -3 π, 3 π, (6 π/1023)}];swpt = StationaryWaveletPacketTransform[data, DaubechiesWavelet[4], 2];swpt["TreeView"]swpt["Dimensions"]WaveletListPlot[swpt, PlotLayout -> "CommonXAxis"]计算一个 LiftingWaveletTransform:
(data = Table[Sinc[t^2], {t, -3 π, 3 π, (6 π/1023)}];);lwt = LiftingWaveletTransform[data, DaubechiesWavelet[4], 2];lwt["TreeView"]lwt["Dimensions"]WaveletListPlot[lwt, PlotLayout -> "CommonXAxis"]高维度 (1)
ϕ = WaveletPhi[DaubechiesWavelet[3]];
ψ = WaveletPsi[DaubechiesWavelet[3]];Plot3D[Evaluate[ϕ[x] ψ[y]], {x, 0, 5}, {y, 0, 5}, PlotRange -> All, ColorFunction -> "SolarColors", Mesh -> None, Axes -> None]Plot3D[Evaluate[ϕ[x] ψ[y]], {x, 0, 5}, {y, -2, 3}, PlotRange -> All, ColorFunction -> "SolarColors", Mesh -> None, Axes -> None]Plot3D[Evaluate[ψ[x] ϕ[y]], {x, -2, 3}, {y, 0, 5}, PlotRange -> All, ColorFunction -> "SolarColors", Mesh -> None, Axes -> None]Plot3D[Evaluate[ψ[x] ψ[y]], {x, -2, 3}, {y, -2, 3}, PlotRange -> All, ColorFunction -> "SolarColors", Mesh -> None, Axes -> None]应用 (3)
data = Table[Abs[Sin[2 π x]] + 1.5 Abs[Cos[2 π x - π]], {x, 0, 1, (1/2^6 - 1)}];执行一个 LiftingWaveletTransform:
dwd = LiftingWaveletTransform[data, DaubechiesWavelet[4]];通过保持 n 个最大的系数以及对其它进行阈值限制,对原始数据求近似:
{data8, data16, data32} = Table[InverseWaveletTransform[WaveletThreshold[dwd, {"LargestCoefficients", n}, Automatic]], {n, {8, 16, 32}}];ListLinePlot[{data, data8, data16, data32}, PlotRange -> All]data = Table[Sin[4 π t] + 2 Exp[-10^5 ((1/3) - t)^2], {t, 0, 1, (1/2^9)}];ListLinePlot[data]dwt = DiscreteWaveletTransform[data, DaubechiesWavelet[4], 4]WaveletListPlot[dwt, {{1}, {0, 1}, {0, 0, 1}, {0, 0, 0}}, PlotLayout -> "CommonXAxis", Method -> {"Inverse" -> True}]data = Table[Sin[2 π t], {t, 0, 1, (1/63)}];ListPlot[data]cumulativeEnergy[data_] := Module[{c = Sort[Flatten[data]^2, Greater]}, (Accumulate[c]/Total[c])]ListPlot[cumulativeEnergy[data]]dwd1 = LiftingWaveletTransform[data, DaubechiesWavelet[4], 1];dwd2 = LiftingWaveletTransform[data, DaubechiesWavelet[4], 2];ListPlot[{cumulativeEnergy[data], cumulativeEnergy[Last /@ dwd1[Automatic]], cumulativeEnergy[Last /@ dwd2[Automatic]]}, PlotStyle -> {Red, Blue, Orange}]属性和关系 (13)
DaubechiesWavelet[1] 等价于 HaarWavelet:
WaveletFilterCoefficients[HaarWavelet[]] === WaveletFilterCoefficients[DaubechiesWavelet[1]]Chop[Total[WaveletFilterCoefficients[DaubechiesWavelet[3], "PrimalLowpass"][[All, 2]]]]Chop[Total[WaveletFilterCoefficients[DaubechiesWavelet[3], "PrimalHighpass"][[All, 2]]]]ϕ = WaveletPhi[DaubechiesWavelet[4]];Subsuperscript[∫, -∞, ∞]ϕ[x]ⅆxTable[NIntegrate[ϕ[(x/2^j)], {x, -∞, ∞}, AccuracyGoal -> 4] - 2^j, {j, 0, 3}]Subsuperscript[∫, -∞, ∞]WaveletPsi[DaubechiesWavelet[4], x]ⅆxNIntegrate[WaveletPsi[DaubechiesWavelet[3], x] WaveletPhi[DaubechiesWavelet[3], x], {x, -∞, ∞}, AccuracyGoal -> 4]WaveletFilterCoefficients[DaubechiesWavelet[3]][[All, 2]].WaveletFilterCoefficients[DaubechiesWavelet[3], "PrimalHighpass"][[All, 2]]DaubechiesWavelet[n] 具有 n 个消失矩;
:
ψ = WaveletPsi[DaubechiesWavelet[2]];Chop[Table[NIntegrate[x^k ψ[x], {x, -∞, ∞}, AccuracyGoal -> 5, MaxRecursion -> 0], {k, 0, 4}]]dwt = LiftingWaveletTransform[ArrayPad[{1, 2, 3, 4}, 2, "Extrapolated"], DaubechiesWavelet[2], 1];Take[InverseWaveletTransform[dwt, Automatic, {0}], 3 ;; -3]dwt = LiftingWaveletTransform[{1, 4, 9, 16}, DaubechiesWavelet[2], 1]InverseWaveletTransform[dwt, Automatic, {0}]ϕ = WaveletPhi[DaubechiesWavelet[2]];a = WaveletFilterCoefficients[DaubechiesWavelet[2], "PrimalLowpass"];scalet[x_, ϕ_, a_] := 2 Table[a[[i, 2]] ϕ[2 x - a[[i, 1]]], {i, Length[a]}]{Plot[Evaluate[scalet[x, ϕ, a]], {x, 0, 3}, PlotRange -> All], Plot[Total[scalet[x, ϕ, a]], {x, 0, 3}, PlotRange -> All]}ϕ = WaveletPhi[DaubechiesWavelet[2]];b = WaveletFilterCoefficients[DaubechiesWavelet[2], "PrimalHighpass"];wavelet[x_, ϕ_, b_] := 2 Table[b[[i, 2]] ϕ[2 x - b[[i, 1]]], {i, Length[b]}]{Plot[Evaluate[wavelet[x, ϕ, b]], {x, -1, 2}, PlotRange -> All], Plot[Total[wavelet[x, ϕ, b]], {x, -1, 2}, PlotRange -> All]}h[wav_, ω_] := With[{a = WaveletFilterCoefficients[wav]}, Underoverscript[∑, i, Length[a]]a[[i, 2]] Exp[-I a[[i, 1]] ω]]Plot[Abs[h[DaubechiesWavelet[2], ω]], {ω, -π, π}, Ticks -> {{-π, -(π/2), 0, (π/2), π}, Automatic}, AxesLabel -> {ω, Abs[H[ω]]}]Plot[{Abs[h[DaubechiesWavelet[2], ω]], Abs[h[DaubechiesWavelet[12], ω]]}, {ω, -π, π}, Ticks -> {{-π, -(π/2), 0, (π/2), π}, Automatic}, AxesLabel -> {ω, Abs[H[ω]]}]h[wav_, ω_] := With[{a = WaveletFilterCoefficients[wav]}, Underoverscript[∑, i, Length[a]]a[[i, 2]] Exp[-I a[[i, 1]] ω]]fh[wav_, ω_, j_] := Abs[Underoverscript[∏, i, j]h[wav, (ω/2^i)]]Plot[fh[DaubechiesWavelet[2], ω, 10], {ω, -10 π, 10 π}, PlotRange -> All, AxesLabel -> {ω, Abs[Overscript[ϕ, ^ ][ω]]}]g[wav_, ω_] := With[{b = WaveletFilterCoefficients[wav, "PrimalHighpass"]}, Underoverscript[∑, i, Length[b]]b[[i, 2]] Exp[-I b[[i, 1]] ω]]Plot[Abs[g[DaubechiesWavelet[2], ω]], {ω, -π, π}, Ticks -> {{-π, -(π/2), 0, (π/2), π}, Automatic}, AxesLabel -> {ω, Abs[G[ω]]}]Plot[{Abs[g[DaubechiesWavelet[2], ω]], Abs[g[DaubechiesWavelet[12], ω]]}, {ω, -π, π}, Ticks -> {{-π, -(π/2), 0, (π/2), π}, Automatic}, AxesLabel -> {ω, Abs[G[ω]]}]h[wav_, ω_] := With[{a = WaveletFilterCoefficients[wav]}, Underoverscript[∑, i, Length[a]]a[[i, 2]] Exp[-I a[[i, 1]] ω]]g[wav_, ω_] := With[{b = WaveletFilterCoefficients[wav, "PrimalHighpass"]}, Underoverscript[∑, i, Length[b]]b[[i, 2]] Exp[-I b[[i, 1]] ω]]fg[wav_, ω_, j_] := Abs[g[wav, (ω/2)] Underoverscript[∏, i = 2, j]h[wav, (ω/2^i)]]Plot[fg[DaubechiesWavelet[2], ω, 10], {ω, -10 π, 10 π}, PlotRange -> All, AxesLabel -> {ω, Abs[Overscript[ψ, ^ ][ω]]}]巧妙范例 (2)
ϕ[x_, j_, k_] := 2^j / 2 WaveletPhi[DaubechiesWavelet[4], 2^j x - k]Plot[Evaluate[Table[ϕ[x, j, 0], {j, 0, 4}]], {x, 0, 2}, Filling -> Axis, PlotRange -> All]Plot[Evaluate[Table[ϕ[x, 2, k], {k, 0, 2^2 - 1}]], {x, 0, 2}, Filling -> Axis, PlotRange -> All]ψ[x_, j_, k_] := 2^j / 2 WaveletPsi[DaubechiesWavelet[4], 2^j x - k]Plot[Evaluate[Table[ψ[x, j, 0], {j, 0, 4}]], {x, 0, 2}, Filling -> Axis, PlotRange -> All]Plot[Evaluate[Table[ψ[x, 2, k], {k, 0, 2^2 - 1}]], {x, -0.5, 1.5}, Filling -> Axis, PlotRange -> All]技术笔记
-
▪
- 自定义小波
相关指南
-
▪
- 小波分析 ▪
- 信号可视化与分析 ▪
- 信号滤波与滤波器设计 ▪
- 信号变换
文本
Wolfram Research (2010),DaubechiesWavelet,Wolfram 语言函数,https://reference.wolfram.com/language/ref/DaubechiesWavelet.html.
CMS
Wolfram 语言. 2010. "DaubechiesWavelet." Wolfram 语言与系统参考资料中心. Wolfram Research. https://reference.wolfram.com/language/ref/DaubechiesWavelet.html.
APA
Wolfram 语言. (2010). DaubechiesWavelet. Wolfram 语言与系统参考资料中心. 追溯自 https://reference.wolfram.com/language/ref/DaubechiesWavelet.html 年
BibTeX
@misc{reference.wolfram_2026_daubechieswavelet, author="Wolfram Research", title="{DaubechiesWavelet}", year="2010", howpublished="\url{https://reference.wolfram.com/language/ref/DaubechiesWavelet.html}", note=[Accessed: 06-September-2026]}
BibLaTeX
@online{reference.wolfram_2026_daubechieswavelet, organization={Wolfram Research}, title={DaubechiesWavelet}, year={2010}, url={https://reference.wolfram.com/language/ref/DaubechiesWavelet.html}, note=[Accessed: 06-September-2026]}