funvalue.m

来自「结构方程中拟合函数的迭代算法」· M 代码 · 共 12 行

function  [Fml,Sigma]=FunValue(S,Lambda_Y,Lambda_X,Beta,Gamma,Phi,Psi,Theta_epsilon,Theta_delta)
A=inv(Beta);
D=A*Gamma;

SYY=[Lambda_Y*A*(Gamma*Phi*Gamma'+Psi)*A'*Lambda_Y'+Theta_epsilon];
SYX=[Lambda_Y*A*Gamma*Phi*Lambda_X'];
SXY=SYX';
SXX=[Lambda_X*Phi*Lambda_X'+Theta_delta];
Sigma=[SYY SYX;SXY SXX];                                %待预测矩阵
Omega=inv(Sigma)*(Sigma-S)*inv(Sigma);
%%求拟合函数值
Fml=log(det(Sigma))-log(det(S))+trace(S*inv(Sigma))-size(S,1);