0001 function [b,a] = conv_noisegen(pol,zer) 0002 0003 [m,k] = size(pol); 0004 [n,l] = size(zer); 0005 0006 coefb = pol(1,:); 0007 0008 for i = 2:m 0009 coefb = conv(coefb, pol(i,:)); 0010 end 0011 0012 b = nonzeros(coefb); 0013 0014 if n~=0 0015 coefa = zer(1,:); 0016 for i = 2:n 0017 coefa = conv(coefa, zer(i,:)); 0018 end 0019 a = nonzeros(coefa); 0020 else 0021 a = 1; 0022 end 0023 0024 %normalize to bn = 1 0025 m = length(b); 0026 normfac = b(m); 0027 b = b/normfac; 0028 a = a/(normfac*sqrt(2)); 0029 0030