bayes.m

来自「这是在网上下的一个东东」· M 代码 · 共 19 行

%
%  [covmp,mmap]=bayes(G,mprior,covm,d,covd)
%
%  Given a linear inverse problem Gm=d, a prior mean mprior and covariance
%  matrix covm, data d, and data covariance matrix covd, this function 
%  computes the MAP solution and the corresponding covariance matrix.
%
function [covmp,mmap]=bayes(G,mprior,covm,d,covd)
covmp=inv(G'*inv(covd)*G+inv(covm));
%
% This takes care of any lack of symmetry in covmp.
%
covmp=(covmp+covmp')/2;
covd12=sqrtm(inv(covd));
covm12=sqrtm(inv(covm));
A=[covd12*G; covm12];
rhs=[covd12*d; covm12*mprior];
mmap=A\rhs;