新手求解答,谢谢。
function [Amix, Bmix, Amix2] = calcabmix(comp, A, B, BIP)
ncomp = size(comp,1);
% Calculate
Aij = zeros(ncomp, ncomp); 这里为什么要用zeros?
for i = 1:ncomp
for j = 1:ncomp
Aij(i,j) = sqrt(A(i)*A(j))*(1 - BIP(i,j));
end
end
% Amix and Bmix are scalars.
Amix = comp'*Aij*comp; 这里为啥要转置啊?
Bmix = comp'*B;
% Amix2 is a vector.
Amix2 = Aij*comp;
end