Examples
B = [1/2 1]; A = [1 1]; w = linspace(0,4,128); H = freqs(B,A,w); [Bh,Ah, C] = invfreqs(H,w,1,1); Bh = 0.50000 1.00000 Ah = 1.0000 1.0000 C = -3.0964e-15 |  |  |
Fit filter B(s)/A(s)to the complex frequency response H at frequency points F. A and B are real polynomial coefficients of order nA and nB.
[B,A,C] = invfreqs(H,F,nB,nA,W,iter,tol,trace) [B,A,C] = invfreqs(H,F,nB,nA,W) [B,A,C] = invfreqs(H,F,nB,nA)
desired complex frequency response.
frequency (must be same length as H).
order of the numerator polynomial B.
order of the denominator polynomial A.
vector of weights (must be same length as F).
This is an Octave function. Fit filter B(s)/A(s)to the complex frequency response H at frequency points F. A and B are real polynomial coefficients of order nA and nB. Optionally, the fit-errors can be weighted vs frequency according to the weights W. Note: all the guts are in invfreq.m
B = [1/2 1]; A = [1 1]; w = linspace(0,4,128); H = freqs(B,A,w); [Bh,Ah, C] = invfreqs(H,w,1,1); Bh = 0.50000 1.00000 Ah = 1.0000 1.0000 C = -3.0964e-15 | | |