This Page

Quick search

#!/usr/bin/env pythonfrompylabimport*# create some data to use for the plotdt=0.001t=arange(0.0,10.0,dt)r=exp(-t[:1000]/0.05)# impulse responsex=randn(len(t))s=convolve(x,r)[:len(x)]*dt# colored noise# the main axes is subplot(111) by defaultplot(t,s)axis([0,1,1.1*amin(s),2*amax(s)])xlabel('time (s)')ylabel('current (nA)')title('Gaussian colored noise')# this is an inset axes over the main axesa=axes([.65,.6,.2,.2],axisbg='y')n,bins,patches=hist(s,400,normed=1)title('Probability')setp(a,xticks=[],yticks=[])# this is another inset axes over the main axesa=axes([0.2,0.6,.2,.2],axisbg='y')plot(t[:len(r)],r)title('Impulse response')setp(a,xlim=(0,.2),xticks=[],yticks=[])show()