# Pweave Example - Frequency response of a moving average filter [Matti Pastell](http://mpastell.com/) Create 11 point moving average filter and plot its frequency response and print the values. <<>>= from pylab import * import scipy.signal as signal #A function to plot frequency and phase response def mfreqz(b,a=1): w,h = signal.freqz(b,a) h = abs(h) return(w/max(w), h) @ Make the impulse response function and use terminal formatted output. <
    >=
    n = 11.
    n
    b = repeat(1/n, n)
    b
    @

    Calculate the frequency response and plot it:

    <
    >=
    w, h = mfreqz(b)
    #Plot the function
    plot(w,h,'k')
    ylabel('Amplitude')
    xlabel(r'Normalized Frequency (x$\pi$rad/sample)')
    show()
    @