/ani/mrses

To get this branch, use:
bzr branch http://suren.me/webbzr/ani/mrses
1 by Suren A. Chilingaryan
Initial import 
1
 
function [rbhata,rmahal,rcorr]=bmc(x1,x2)

2
 









3



c1=cov(x1,1);










4



c2=cov(x2,1);










5



c=(c1+c2)./2;










6












7



m1=mean(x1);










8



m2=mean(x2);










9












10



rmahal=((m2-m1)/c)*(m2-m1)';










11



rcorr=2.*log(det(c)./sqrt(det(c1).*det(c2)));










12



rbhata=rmahal./8+rcorr./4;