To check out this repository please hg clone the following URL, or open the URL using EasyMercurial or your preferred Mercurial client.
root / _beattracker / dojsdiv2.m
History | View | Annotate | Download (471 Bytes)
| 1 |
function sd = dojsdiv2(pre,post); |
|---|---|
| 2 |
|
| 3 |
pre = pre-min(pre); |
| 4 |
post = post-min(post); |
| 5 |
pre = pre(:)'; |
| 6 |
post = post(:)'; |
| 7 |
pre = pre+eps; |
| 8 |
post = post+eps; |
| 9 |
pre = pre/sum(pre); |
| 10 |
|
| 11 |
pi1 = 0.5; |
| 12 |
pi2 = 0.5; |
| 13 |
|
| 14 |
% normalise to sum to unity |
| 15 |
pre = pre/(sum(pre)+eps); |
| 16 |
post = post/(sum(post)+eps); |
| 17 |
|
| 18 |
% combined distribution |
| 19 |
a1 = pi1*pre + pi2*post; |
| 20 |
|
| 21 |
%sd = -sum(a1.*log(a1)) + pi1*sum(pre.*log(pre)) + pi2*sum(post.*log(post)); |
| 22 |
|
| 23 |
|
| 24 |
sd = sum((-a1.*log(a1)) + pi1*(pre.*log(pre)) + pi2*(post.*log(post))); |