cepstrum - cepstrum calculation

Calling Sequence

fresp = cepstrum(w,mag)

Parameters

• 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