normalcopulapdf.m

来自「经济学专业代码」· M 代码 · 共 20 行

function F_U = NormalCopulaPDF(u,Mu,Sigma)
% this function computes the pdf of the copula of the normal distribution
% at the generic point u in the unit hypercube
% see "Risk and Asset Allocation"-Springer (2005), by A. Meucci


N=length(u);
s=sqrt(diag(Sigma));

x=norminv(u,Mu,s);

Numerator = (2*pi)^(-N/2) * (  (det(Sigma))^(-.5) ) * exp(-.5*(x-Mu)'*inv(Sigma)*(x-Mu));

fs=normpdf(x,Mu,s);
Denominator = prod(fs);

F_U = Numerator/Denominator;