probtrans.m

fs1016w编码标准matlab源码

function  A=probtrans(B,l)
%求输入矩阵B的转移概率矩阵
%l索引最大值
A=zeros(l);
for i=1:length(B)-1
A(B(i)+1,B(i+1)+1)=A(B(i)+1,B(i+1)+1)+1;
end

for i=1:l
 c=sum(A(i,:));
  for j=1:l
    
      if c==0
            A(i,j)=0;
      else   A(i,j)=A(i,j)/c;
      end
    end
end
%mesh(A(i),A(j),A(i,j))
%bar3(A)
%mesh(A)
surf(A)
%stem3(A,'*');
xlabel('Quantized Index')
ylabel('Transitted Index')
zlabel('Transition Probability')