cepstrum

cepstrum calculation

Calling Sequence

fresp = cepstrum(w,mag)

Arguments

:w positive real vector of frequencies (rad/sec) : :mag real vector of magnitudes (same size as w) : :fresp complex vector :

Description

fresp = cepstrum(w,mag) returns a frequency response fresp(i) whose magnitude at frequency w(i) equals mag(i) and such that the phase of freq corresponds to a stable and minimum phase system. w needs not to be sorted, but minimal entry should not be close to zero and all the entries of w should be different.

Examples

w=0.1:0.1:5;mag=1+`abs`_(`sin`_(w));
fresp=cepstrum(w,mag);
`plot2d`_([w',w'],[mag(:),`abs`_(fresp)])

See Also

• frfit frequency response fit