home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 8.ddi / SIGNAL.DI$ / NUMF.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  273 b   |  9 lines

  1. function b = numf(h,a,nb)
  2. %NUMF    Find numerator B given impulse-response h of B/A and denominator A
  3. %    NB is the numerator order.  This function is used by YULEWALK.
  4.   
  5. nh = max(size(h)); 
  6. impr = filter(1,a,[1 zeros(1,nh-1)]);
  7. b = h/toeplitz(impr,[1 zeros(1,nb)])';
  8.  
  9.