function b = de2bi(d, n, p)
%function b = de2bi(d, n, p)
%DE2BI  轉(zhuǎn)換10進(jìn)制數(shù)為二進(jìn)制數(shù)。
%        B = DE2BI(D) 轉(zhuǎn)換正整數(shù)向量D成二進(jìn)制矩陣B。
%        二進(jìn)制矩陣B的每一行表示十進(jìn)制向量D中相應(yīng)的數(shù)。
%       B = DE2BI(D, N) 轉(zhuǎn)換正整數(shù)向量D成二進(jìn)制矩陣B，
%        但指定B的列數(shù)為N。
%       B = DE2BI(D, N, P) 轉(zhuǎn)換正整數(shù)向量D成p進(jìn)制矩陣B。
%      p進(jìn)制矩陣B的每一行表示十進(jìn)制向量D中相應(yīng)的數(shù)。
%           作者：機(jī)自01-2班曾新海
%           zxh21st@163.com
d = d(:);len_d = length(d);
if min(d) < 0, error('Cannot convert a negative number');
elseif ~isempty(find(d==inf)),
     error('Input must not be Inf.');
elseif find(d ~= floor(d)), 
    error('Input must be an integer.');  
end;
if nargin < 2,
   tmp = max(d); b1 = [];
   while tmp > 0
      b1 = [b1 rem(tmp, 2)];tmp = floor(tmp/2);
   end;
   n = length(b1);
end;
if nargin < 3,p = 2;end;
b = zeros(len_d, n);
for i = 1 : len_d
   j = 1;tmp = d(i);
   while (j <= n) & (tmp > 0)
      b(i, j) = rem(tmp, p);tmp = floor(tmp/p);
      j = j + 1;
end;end;