mrice.m 190 Bytes
Newer Older
1 2 3 4 5 6 7 8 9 10 11
function [m1r, m2r, stdr] = mrice(sig, c)

q = c.^2/4./sig^2;

m1r = sqrt(pi*sig^2/2).*exp(-q).*(besseli(0, q).*(1+2*q)+besseli(1,q)*2.*q);

m2r = 2*sig.^2+c.^2;

stdr = sqrt(m2r-m1r.^2);