minimum-phase filter, compensating filter

The transform function is: [tex]H(s) = 0.04 \frac{s^2+200s + 1.0001 \times 10^8}{s^2 + 1.52\times 10^3s + 2.92 \times 10 [/tex]
Matlab codes:

H1 = tf(0.04*[1 2e2 1.0001e8],[1 1.52e3 2.92e6])
pzmap(H1)
[numH1, denH1] = tfdata(H1,'v');
G1 = tf(denH1, numH1)
s1 = wavread('distortedSound.wav');
tsteps = 1/44100*(0:1:(length(s1)-1));
y1 = lsim(G1,s1,tsteps);
sound(y1,44100)

 [update: 10/5/2010] these are basically used for compensating channel impairments. #