没有fft MODWT IMODWT代码
19日视图(30天)
显示旧的评论
嗨。
我研究了最大重叠小波变换及其属性在“时间序列分析的小波方法唐纳德·b·珀西瓦尔,安德鲁·t·瓦尔登湖”,我看到的应用程序之后才执行fft算法的发展加快代码。为了更好地理解它是如何工作的,我需要开发不使用fft算法。有图书馆做这个吗?
谢谢你！
现在我发现modwt代码:
函数[多项式系数]= NOFFTmodwt (x,过滤器,级别)
%过滤器:‘db1’,‘db2’,‘哈雾’,……
%返回:看到matlab
%过滤器
N =长度(x);
[Lo_D, Hi_D] = wfilters(过滤器);
g_t = Lo_D /√(2);
h_t = Hi_D /√(2);
wavecoeff =细胞(1、5);
v_j_1 = x;
为j = 1:水平
w = circular_convolve_d (h_t v_j_1, j);
v_j_1 = circular_convolve_d (g_t v_j_1, j);
wavecoeff {j} = w;
结束
wavecoeff{水平+ 1}= v_j_1;
多项式系数= vertcat (wavecoeff {:});
多项式系数(1:级别1:N) = -coeffs(1:级别1:N);
函数w_j = circular_convolve_d (h_t v_j_1, j)
% j层次分解
% h_t: \波浪号{h} = h /√(2)
% v_j_1: v_ {j - 1}, (j - 1)届规模系数
%返回:w_j(或v_j)
N =长度(v_j_1);
L =长度(h_t);
w_j = 0 (1, N);
l = 0: l - 1;
为t = 1: N
指数=国防部(t 1 - 2 ^ (j - 1) * l, N) + 1;
v_p = v_j_1(指数);
w_j (t) =点(h_t v_p);
结束
从Python函数modwtpy自由翻译,你可以在这里找到:
平等的结果以Matlab modwt函数,边界sx面向“哈尔”,大便invariante,相同的多项式系数,ecc . .
现在我必须找到imodwt相同。这是我设法做的但它不工作!
函数[x] = NOFFTimodwt(多项式系数、过滤器)
%逆modwt
%过滤器
[Lo_R, Hi_R] = wfilters(过滤器);
g_t = Lo_R /√(2);
h_t = Hi_R /√(2);
水平=大小(多项式系数,1)- 1;
v_j =多项式系数(结束:);
为j =级别:1:1
% w_j = circshift(多项式系数(j:), [0, 2 ^ (j - 1)));
% v_j = circular_convolve_s (h_t、g_t w_j, v_j, j);
v_j = circular_convolve_s (h_t g_t,多项式系数(j,:), v_j, j);
结束
x = v_j;
函数v_j_1 = circular_convolve_s (h_t g_t、w_j v_j, j)
% (j - 1) th级合成从w_j w_j
%见函数circular_convolve_d
N =长度(v_j);
L =长度(h_t);
v_j_1 = 0 (1, N);
l = 0: l - 1;
为t = 1: N
指数=国防部(t 1 + 2 ^ (j - 1) * l, N) + 1;
%指数=国防部(t + 2 ^ (j - 1) * l - 1, N) + 1;
w_p = w_j(指数);
v_p = v_j(指数);
v_j_1 (t) =点(h_t w_p) +点(g_t v_p);
结束
有什么建议吗?
谢谢!
