第二代Vold-Kalman过滤

过滤信号的非平稳周期组件使用multi-order Vold-Kalman过滤器。

2.2 k下载

更新2013年5月3

查看许可协议

Vold-Kalman过滤器,由哈佛Vold引入,并在1993年1月Leuridan,能够从信号中提取非平稳周期组件使用一个已知频率向量[1]。被制定在最小二乘意义上说,它可以解决稀疏线性系统。类似于卡尔曼滤波器,VKF最小化代价函数基于结构方程和方程的数据。

这个提交实现第二代VKF能够提取多个订单的同时,与能量传播的交叉订单[2]。

语法:
x = vkf (y, f, f)提取顺序与信号频率向量f y samplerate fs,使用3 db 2-pole滤波器带宽的采样率的1%。输出是一个波形x。

[…]= vkf (y, f, f p)使用p-order过滤器(通常1或4)。每个订单每十年增加-40分贝的转出。通过指定额外的低阶系数,零边界条件。例如:p =(2 0 1)应用二阶滤波和部队信封振幅及其一阶导数为零在t_1和t_N。

[…)= vkf (y, f, f, p, bw)使用bw规定赫兹的带宽。如果一个标量bw,使用恒定带宽;如果bw是一个向量与y长度一样,意识到一个时变瞬时带宽。

X = vkf (y, F, F,…) (N, K) =大小(F),执行同步提取K订单频率向量(F,……f_K]数组f .在穿越订单的情况下,这种方法试图揭示correctorder振幅。K波形的输出是一个数组(x_1,…, xk)。

X = vkf (y, F, F, p, bw, 0)切换到一个单一顺序的算法。K订单仍然提取,但单一顺序算法在计算上更少的要求。这是高采样率的建议和/或长timeseries。

[c] = vkf(…)返回复杂的信封(s)和矢量c (s),这样的顺序可以重建波形(s) x = (. * c)。

[a、c r] = vkf(…)输出一个额外的选择性向量r用于实现带宽的bw。

演示:
调用VKF不带参数显示了两个穿越小示范multi-order过滤订单白噪声的存在。

注意:过滤器订单> 4通常导致生病的条件,应该避免。滤波器带宽的决心实现任意订单但没有核实订单高于3。

作者:Maarten范德Seijs, 2010。
1.4版,2013年5月3。

引用:
[1]Vold h . Leuridan, j .(1993),以极端杀利率高分辨率的订单跟踪,利用卡尔曼跟踪滤波器。技术报告931288年,汽车工程师学会。

[2]Tuma, j .(2005)设置Vold-Kalman订单跟踪滤波器的通带宽度。学报》国际国会声音和振动(ICSV12),葡萄牙里斯本。

引用作为

Maarten van der Seijs (2023)。第二代Vold-Kalman过滤(//www.tianjin-qmedu.com/matlabcentral/fileexchange/36277-second-generation-vold-kalman-order-filtering), MATLAB中央文件交换。检索

MATLAB版本兼容性
创建R2009b
兼容任何释放
平台的兼容性
窗户 macOS Linux

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!
版本 发表 发布说明
1.7.0.0

现在有新的m文件

1.4.0.0

错误更新单一顺序的算法。

1.3.0.0

微小的变化来描述和代码。功能是一样的。

1.0.0.0