# appcoef
一维近似系数
函数库: TyWavelet
# 语法
A = appcoef(C, L, wname)
A = appcoef(C, L, LoR, HiR)
A = appcoef(___, N)
# 说明
A = appcoef(C, L, wname) 返回使用一维信号的小波分解结构 [C, L] 和 wname 指定的小波在最粗尺度上的近似系数。(参见 wavedec 获取更多信息)。
A = appcoef(C, L, LoR, HiR) 使用低通重建滤波器 LoR 和高通重建滤波器 HiR。(更多信息见 wfilters)。
A = appcoef(___, N) 如果 [C, L] 是一个一维信号的 M 层小波分解结构,那么 0 ≤ N ≤ M,返回 N 层的近似系数。示例
# 示例
三级近似系数
这个例子显示了如何提取第 3 层级近似系数的方法。
加载由用电数据组成的信号。
using TyWavelet
using TyPlot
using TyBase
pkg_dir = pkgdir(TyWavelet)
source_path = pkg_dir * "/examples/Resources/leleccum.mat"
y = load(source_path)
leleccum = y["leleccum"]
sig = leleccum[1:3920]
用 "sym4" 小波获得降至 5 级的 DWT。
C, L = wavedec(sig, 5, "sym4")
提取 3 级近似系数。绘制原始信号和近似系数图。
Lev = 3
a3 = appcoef(C, L, "sym4", Lev)
subplot(2, 1, 1)
plot(sig)
title("Original Signal")
subplot(2, 1, 2)
plot(a3)
title("Level-3 Approximation Coefficients")
你可以用 1 到 5 的任何数值代替 Lev,得到相应层级的近似系数。
# 输入参数
C - 小波分解向量实值向量
小波分解向量,指定为一个实值向量。向量 C 是 wavedec 的输出。
L - 簿记向量正整数的向量
一维信号的小波分解的簿记向量,指定为一个正整数向量。簿记向量用于按层级解析小波分解向量 C 中的系数。
wname - 小波字符串标量
用于计算单级离散小波逆变换(IDWT)的小波,指定为字符串标量。该小波必须被 wavemngr 识别。关于每个系列中可用的小波,见 wfilters。
LoR - 小波低通重建滤波器偶数长度的实值向量
小波低通重建滤波器,指定为偶数长度的实值向量。LoR 必须与 HiR 的长度相同。LoR 必须是与用于创建小波分解结构 [C, L] 的小波相关的低通重建滤波器。(更多信息见 wfilters)。
HiR - 小波高通重建滤波器偶数长度的实值向量
小波高通重建滤波器,指定为偶数长度的实值向量。HiR 必须与 LoR 的长度相同。HiR 必须是与用于创建小波分解结构 [C, L] 的小波相关的高通重建滤波器。(更多信息见 wfilters)。
N - 近似系数层级正整数
近似系数层级,指定为正整数。如果 [C, L] 是一个一维信号的 M 层小波分解结构,那么 0 ≤ N ≤ M。
# 输出参数
A - 近似系数实值向量
N 层的近似系数,作为实值向量返回。
# 算法
输入向量 C 和 L 包含所有关于信号分解的信息。
让 NMAX = length(L) - 2;然后 C = [A(NMAX) D(NMAX) ... D(1)] 其中 A 和 D 是向量。如果 N = NMAX,则进行简单提取;否则,appcoef 使用小波逆变换迭代计算近似系数。