mysymsetdiff.m

马尔可夫过程是一类重要的随机过程,它的原始模型马尔可夫链,由俄国数学家Α.Α.马尔可夫于1907年提出。本程序是对隐马尔可夫模型的一个完整的建模。

function C = mysymsetdiff(A,B)
% MYSYMSETDIFF Symmetric set difference of two sets of positive integers (much faster than built-in setdiff)
% C = mysetdiff(A,B)
% C = (A\B) union (B\A) = { things that A and B don't have in common }

if isempty(A)
  ma = 0;
else
  ma = max(A);
end

if isempty(B)
  mb = 0;
else
  mb = max(B);
end

if ma==0 
  C = B;
elseif mb==0
  C = A;
else % both non-empty
  m = max(ma,mb);
  bitsA = sparse(1, m);
  bitsA(A) = 1;
  bitsB = sparse(1, m);
  bitsB(B) = 1;
  C = find(xor(bitsA, bitsB));
end