function b = de2bi(d, n, p)
%DE2BI  轉換10進制數為二進制數。
%        B = DE2BI(D) 轉換正整數向量D成二進制矩陣B。
%        二進制矩陣B的每一行表示十進制向量D中相應的數。
%       B = DE2BI(D, N) 轉換正整數向量D成二進制矩陣B，
%        但指定B的列數為N。
%       B = DE2BI(D, N, P) 轉換正整數向量D成p進制矩陣B。
%      p進制矩陣B的每一行表示十進制向量D中相應的數。

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;