function [x] = chan(h,N)
%	noise shaping filter for fig8_12.M
%	Generates N samples of a coloured o/p signal of FIR filter h
%	July 97

M = max(size(h));
mag = sqrt(h*h');
h = h/mag;
y = randn(1,N+M);
z = filter(h,1,y);
x(1:N) = z(M+1:N+M);