diff options
| author | ziejd2 | 2018-03-14 23:23:33 -0500 |
|---|---|---|
| committer | GitHub | 2018-03-14 23:23:33 -0500 |
| commit | 1ff6baa44e22b91eefb48aea6f3befa078c0489b (patch) | |
| tree | e0fd79d2e32fd2aedda2eadaed0f19af3514c520 /sourcecodes/bnt-master/KPMtools/logsum.m | |
| parent | 6882395afdadf4e982b25b5215071a0932730950 (diff) | |
| parent | c80226899f5cdd9f11c163817d59445213f5bef0 (diff) | |
| download | BNW-1ff6baa44e22b91eefb48aea6f3befa078c0489b.tar.gz | |
Merge pull request #1 from ziejd2/octave_php_separate
Octave php separate
Diffstat (limited to 'sourcecodes/bnt-master/KPMtools/logsum.m')
| -rw-r--r-- | sourcecodes/bnt-master/KPMtools/logsum.m | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/sourcecodes/bnt-master/KPMtools/logsum.m b/sourcecodes/bnt-master/KPMtools/logsum.m new file mode 100644 index 00000000..1d8ae457 --- /dev/null +++ b/sourcecodes/bnt-master/KPMtools/logsum.m @@ -0,0 +1,28 @@ +function y=logsum(x,d) +%LOGSUM logsum(x,d)=log(sum(exp(x),d)) +% d gives dimension to sum along + +% Copyright (C) Mike Brookes 1998 +% +% Last modified Mon Oct 12 15:47:25 1998 +% +% VOICEBOX is a MATLAB toolbox for speech processing. Home page is at +% http://www.ee.ic.ac.uk/hp/staff/dmb/voicebox/voicebox.html +% + +if nargin==1 + d=[find(size(x)-1) 1]; + d=d(1); +end +n=size(x,d); +if n<=1, y=x; return; end +s=size(x); +p=[d:ndims(x) 1:d-1]; +z=reshape(permute(x,p),n,prod(s)/n); + +y=max(z); +y=y+log(sum(exp(z-y(ones(n,1),:)))); + +s(d)=1; +y=ipermute(reshape(y,s(p)),p); + |