print (SNR, bitrate) print "With Q=5" SNR, bitrate, quantizedApprox = MPCoder.SimpleMDCTEncoding(mpApprox, 8000, Q=5); print (SNR, bitrate) SNR, bitrate, quantizedApprox = MPCoder.SimpleMDCTEncoding(mpApprox, 2000, Q=14); print (SNR, bitrate) pyLODico = pymp_MDCTDico.pymp_LODico(scales); lompApprox , lompDecay = MP.MP(myPympSignal, pyLODico, 20, 2000,padSignal=False); SNRlo, bitratelo, quantizedApproxLO = MPCoder.SimpleMDCTEncoding(lompApprox, 2000, Q=14, TsPenalty=True); print (SNRlo, bitratelo) pyRSSDico = pymp_RandomDicos.pymp_RandomDico(scales); rssApprox , rssDecay = MP.MP(myPympSignal, pyRSSDico, 20, 2000,padSignal=False); SNRrss, bitraterss, quantizedApproxRSS = MPCoder.SimpleMDCTEncoding(rssApprox, 2000, Q=14); print (SNRrss, bitraterss) print (quantizedApprox.atomNumber, quantizedApproxLO.atomNumber , quantizedApproxRSS.atomNumber) #quantizedApprox.plotTF() #plt.show() print " now at a much larger level : a SNR of nearly 50 dB and around 64 Kbps" mpApprox , mpDecay = MP.MP(myPympSignal, pyDico, 50, 20000,padSignal=False); SNR, bitrate, quantizedApprox = MPCoder.SimpleMDCTEncoding(mpApprox, 64000, Q=16); print (SNR, bitrate)
mpl.rcParams['image.interpolation'] = 'Nearest';
mpl.rcParams['text.usetex'] = True;
# Load glockenspiel signal
abPath = os.path.abspath('../../data/');
pySig = pymp_Signal.InitFromFile(abPath+'/glocs.wav',forceMono=True,doNormalize=True);
pySig.crop(0,3*pySig.samplingFrequency)
scales = [128,1024,8192];
nbAtom = 500;
srr = 30;
mpDico = pymp_MDCTDico.pymp_MDCTDico(scales)
lompDico = pymp_MDCTDico.pymp_LODico(scales)
rssMPDico = pymp_RandomDicos.pymp_RandomDico(scales)
approxMP , decayMP = MP.MP(pySig,mpDico,srr,nbAtom,padSignal=True);
approxLoMP, decayLoMP = MP.MP(pySig,lompDico,srr,nbAtom,padSignal=True);
approxRSSMP, decayRSSMP = MP.MP(pySig,rssMPDico,srr,nbAtom,padSignal=False);
plt.figure()
plt.plot(10*np.log10(decayMP) - 10*np.log10(decayMP[0]))
plt.plot(10*np.log10(decayLoMP) - 10*np.log10(decayLoMP[0]),'r:')
plt.plot(10*np.log10(decayRSSMP) - 10*np.log10(decayRSSMP[0]),'k--')
plt.ylabel('Normalized reconstruction error (dB)')
plt.xlabel('Iteration')
plt.legend(('Standard MP','Locally Optimized MP','RSS MP'))
plt.grid()
plt.show()
