文档

wvarchg

寻找方差变化点

语法

[PTS_OPT, KOPT T_EST] = wvarchg (Y, K, D)

描述

[PTS_OPT, KOPT T_EST] = wvarchg (Y, K, D)计算信号的估计方差变化点Y对于j个变化点,j = 0,1,2,…,K

整数D是两个更改点之间的最小延迟。

整数KOPT变更点的建议数量(0≤?KOPTK).向量PTS_OPT包含相应的变更点。

1kKT_EST (k + 1, 1: k)包含了k方差的瞬间变化点,然后,如果KOPT > 0PTS_OPT = T_EST (KOPT + 1, 1: KOPT)其他的PTS_OPT = []

K而且D必须是整数,使1 <K< <长度(Y)和1≤D< <长度(Y).

信号Y均值应该是0。

wvarchg (Y, K)相当于wvarchg (Y, K, 10)

wvarchg (Y)相当于wvarchg (Y 6 10)

例子

全部折叠

在块信号中添加两个方差变化点。检测方差变化点wvarchg

加载块信号。添加位于指数180和600的两个方差变化点的白噪声。

x = wnoise (10);rng默认的;bb = 1.5 * randn(1、长度(x));cp1 = 180;cp2 = 600;x = x + (bb (1: cp1), bb (cp1 + 1: cp2) / 4, bb (cp2 + 1:结束)];

得到一级小波系数。将前2%的值替换为小波系数的平均值,去除所有信号。

wname =“db4”;列弗= 1;[c、l] = wavedec (x,列弗,wname);侦破= wrcoef (' d 'c l wname, 1);y = (abs(检波器));v2p100 = y(修复(长度(y) * 0.98));印第安纳州=找到(abs(检波器)> v2p100);依据(印第安纳州)=意味着(依据);

利用小波系数估计方差变化点。

[pts_Opt, kopt t_est] = wvarchg(侦破,5);sprintf ('估计的更改点是%d和%d\n'pts_Opt)
ans = '估计的变更点为181和601 '

参考文献

Lavielle, M.(1999),“检测因变量序列中的多个变化”,意义下。procs及其应用, 83, 2,第79-102页。

之前介绍过的R2006a

这个话题有帮助吗?