function [y] = sys(h,x,sigman)
%	unknown system and additive white noise for fig8_12.m
%	system h to be identified
%	noise rms sigman
%	July 97


y = filter(h,1,x) + sigman*randn(size(x));