gerpar.m

来自「This manual describes how to run the Mat」· M 代码 · 共 52 行

% 02/01/1997
% gerpar.m
%
function [] = gerpar(nbits)
nelem = 2^nbits;
nl = nelem;
for i=1:nbits,
 nl = 0.5*nl;
 valor = 1;
 n = 0;
 ind = 1;
 while n ~= nelem,
  if valor == 0,
   valor = 1;
  else
   valor = 0;
  end
  for j=(ind):(ind+nl-1)
   X(j,i) = valor;
   n = n+1;
  end
  ind = ind+nl;
 end
end
for i=1:nelem,
 if rem(X(i,:)*ones(nbits,1),2) == 0,
  S(i,1) = 1;
 else
  S(i,1) = 0;
 end
end
%disp('[X S]');disp([X S]);pause;
[l,c] = size(X);
for i=1:l,
 for j=1:c,
  if X(i,j) == 0,
   X(i,j) = -1;
  end
 end
end
[l,c] = size(S);
for i=1:l,
 for j=1:c,
  if S(i,j) == 0,
   S(i,j) = -1;
  end
 end
end
size(X), size(S),
P = X; T = S;
save dados_par3 P T;