Image:Pulse density modulation.svg

from pylab import *
from numpy import *
 
def pdm(x):
    y = zeros_like(x)
    qe = zeros_like(x)
 
    qe[-1] = 0 # inital [[:en:quantization error]] is 0
    for n in range(x.size):
        y[n] = 2*( x[n] > qe[n-1] ) - 1
        qe[n] = y[n] - x[n] + qe[n-1]
 
    return y, qe
 
t = arange(500)
x = sin(2*pi/500 * t)
y, qe = pdm(x)
 
plot( t, y, label="PDM sine wave" )
plot( t, x, label="PCM sine wave", linewidth=5 )
 
axis((-25, 525, -1.25, 1.5))
legend(loc='upper right')
 
xticks(fontsize=16)
yticks(fontsize=16)
xlabel("Discrete time (n)", fontsize=18)
ylabel("Signal amplitude",fontsize=18)
 
savefig("pdm.svg")