0001 function [dist,bcohpotential,bcohlaplacian] = baseline_coherence(s12,s13,r1,r2,r3,r4,rz); 0002 [a,b,c,d,e,f,g] = fender_4(500,s12,s13,r1,r2,r3,r4,rz); 0003 ang = ([0:33]/33)*27/9.2; 0004 dist = ang*9.2; 0005 z = cos(ang); 0006 p = legendre_r(z,500); 0007 h = (f+g)'; 0008 numpot = p.*((h.^2)*ones(1,size(p,2))); 0009 numpotsum = sum(numpot); 0010 denpot = ((h.^2)*ones(1,size(p,2))); 0011 denpotsum = sum(denpot); 0012 bcohpotential = (numpotsum./denpotsum).^2; 0013 numlap = p.*((h.^2.*[1:500]'.*[2:501]')*ones(1,size(p,2))); 0014 numlapsum = sum(numlap); 0015 denlap = ((h.^2.*[1:500]'.*[2:501]')*ones(1,size(p,2))); 0016 denlapsum = sum(denlap); 0017 bcohlaplacian = (numlapsum./denlapsum).^2; 0018 0019 0020 0021 0022