sourcecodes/bnt-master/KPMstats/parzenC_test.m


1
2
3
4
5
6
7
8
9
10

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)