S2tinv;
w1=S1tinv*u1;
w2=S2tinv*u2;
%d = det(X) returns the determinant of the square matrix X.矩阵行列式
w10...23
%Solve system of nonlinear equations求解非线性方程组
tt2 = fsolve(‘bayesian_fun’,5,[],t1,W1,W2,w1,w2,w10...end
plot(1:23,t2,’b’,’LineWidth’,3);
% f=bayesian_fun.m
function f=bayesian_fun(t2,t1,W1,W2,w1,w2,w10