function y = d2a(x,b)
% d2a(x) - converts bitstream (b bits per sample) to discrete time,
% continuous amplitude signal
% x - bitstream
% b - number bits per sample
x=x';
N = length(x);
L=floor(N/b);
qinterval = 2/(2^b-1);
if b<=1;
for i=1:N
if x(i)==0
sig(i)=-1;
else
sig(i)=1;
end
end
else
sig=zeros(L,1);
for k=0:((2^b)-1)
u=dec2bin(k);
w=zeros(1,length(u)); % conversion from interval # to binary
for i=1:length(u)
w(i)=str2num(u(i));
end
v=[zeros((b-length(w)),1)',w];
p(k+1,:)=v; % p is a vector of binary labels for interval #'s
end
if b~=4
for i=0:L-1 % mat