Generate pole zero plot from measured impulse response

Hello everybody,

I want to generate a pole zero plot from a measured impulse response (imported via audioread). I only find examples to generate pole zero plots for very simple filters with known differential equations. Maybe you can point me in the right direction or give an example.

Thank you very much in advance!

Answer 1

[Sudheer Nuggehalli]

One possible way to do this is to use the "pzmap" function on the measured impulse response data imported via "audioread". However, since "pzmap" requires a Dynamic System Model, sys, as an input, the measured impulse response data needs to be converted into a Dynamic System Model object. One way to do this is to use the "iddata" and "tfest" functions. Below is some code that can get you started:

[y,Fs] = audioread(filename);

data = iddata(y,[],Ts); % where Ts = 1/Fs

sys = tfest(data,np); % where np is the number of poles

pzmap(sys)

Refer to the documentation links below for more information about each of these functions and examples of how they are used:

pzmap audioread iddata tfest