另外,
tauLU(8,:) = ones(1,2)*2*tdis_inv(-sqrt((kappa8(2)+1)*(1-kappa8(1))/(1+kappa8(1))),kappa8(2)+1); % Student's t
运行之后,提示:
??? Error using ==> beta_inv at 28
beta_inv: A probability should be 0<=p<=1
Error in ==> tdis_inv at 25
x = beta_inv(p,1/2,a/2);
其中,beta_inv.m tdis_inv 都是后又加进去的源文件。