Here is a small MATLAB function that computes B in B/A given A and the impulse response h of B/A
function b = numf(h,a,nb) %NUMF Find numerator B given impulse-response h of B/A and denominator A % NB is the numerator order. nh = max(size(h)); impr = filter(1,a,[1 zeros(1,nh-1)]); egqs = toeplitz(impr,[1 zeros(1,nb)]); b = h/egqs';
Could someone help me understand why these steps give B.
I get what
toeplitz() do, but not the last line :
b = h/egqs'
Thanks for any help,