parzenc_test.m

来自「KPMstats is a directory of miscellaneous」· M 代码 · 共 11 行

d = 2; M = 3; Q = 4; T = 5; Sigma = 10;
N = sample_discrete(normalize(ones(1,M)), 1, Q);
data = randn(d,T);
mu = randn(d,M,Q);

[BM, B2M] = parzen(data, mu, Sigma, N);
[B, B2] = parzenC(data, mu, Sigma, N);

approxeq(B,BM)
approxeq(B2,B2M)