day2/session2.tex
changeset 72 1c1d6aaa2be3
parent 71 d0679d2f9339
child 76 b24c2560f626
equal deleted inserted replaced
71:d0679d2f9339 72:1c1d6aaa2be3
   322   \begin{itemize}
   322   \begin{itemize}
   323     \item Integrating Functions given function object
   323     \item Integrating Functions given function object
   324     \item Integrating Functions given fixed samples
   324     \item Integrating Functions given fixed samples
   325     \item Numerical integrators of ODE systems
   325     \item Numerical integrators of ODE systems
   326   \end{itemize}
   326   \end{itemize}
   327   Calculate $\int^1_0(sin(x) + x^2)dx$
   327   Calculate the area under $(sin(x) + x^2)$ in the range $(0,1)$
   328   \begin{lstlisting}
   328   \begin{lstlisting}
   329     >>> def f(x):
   329     >>> def f(x):
   330             return np.sin(x)+x**2
   330             return np.sin(x)+x**2
   331     >>> integrate.quad(f, 0, 1)
   331     >>> integrate.quad(f, 0, 1)
   332   \end{lstlisting}
   332   \end{lstlisting}
   338   \begin{align*}
   338   \begin{align*}
   339   \frac{dx}{dt}&=-e^{-t}x^2\\ 
   339   \frac{dx}{dt}&=-e^{-t}x^2\\ 
   340            x(0)&=2    
   340            x(0)&=2    
   341   \end{align*}
   341   \end{align*}
   342   \begin{lstlisting}
   342   \begin{lstlisting}
   343     def dx_dt(x,t):
   343 >>> def dx_dt(x,t):
   344         return -np.exp(-t)*x**2
   344         return -np.exp(-t)*x**2
   345 
   345 
   346     x=integrate.odeint(dx_dt, 2, t)
   346 >>> x=integrate.odeint(dx_dt, 2, t)
   347     plt.plot(x,t)
   347 >>> plt.plot(x,t)
   348   \end{lstlisting}
   348   \end{lstlisting}
   349 \inctime{10}
   349 \inctime{10}
   350 \end{frame}
   350 \end{frame}
   351 
   351 
   352 \subsection{Interpolation}
   352 \subsection{Interpolation}
   367   \end{lstlisting}
   367   \end{lstlisting}
   368 \end{frame}
   368 \end{frame}
   369 
   369 
   370 \begin{frame}[fragile]
   370 \begin{frame}[fragile]
   371   \frametitle{Interpolation - Splines}
   371   \frametitle{Interpolation - Splines}
   372   Cubic Spline of $sin(x)$
   372   Plot the Cubic Spline of $sin(x)$
   373   \begin{lstlisting}
   373   \begin{lstlisting}
   374     x = np.arange(0,2*np.pi,np.pi/4)
   374 >>> x = np.arange(0,2*np.pi,np.pi/4)
   375     y = np.sin(x)
   375 >>> y = np.sin(x)
   376     tck = interpolate.splrep(x,y)
   376 >>> tck = interpolate.splrep(x,y)
   377     X = np.arange(0,2*np.pi,np.pi/50)
   377 >>> X = np.arange(0,2*np.pi,np.pi/50)
   378     Y = interpolate.splev(X,tck,der=0)
   378 >>> Y = interpolate.splev(X,tck,der=0)
   379     plt.plot(x,y,'o',x,y,X,Y)
   379 >>> plt.plot(x,y,'o',x,y,X,Y)
   380     plt.show()
   380 >>> plt.show()
   381   \end{lstlisting}
   381   \end{lstlisting}
   382 \inctime{10}
   382 \inctime{10}
   383 \end{frame}
   383 \end{frame}
   384 
   384 
   385 \subsection{Signal Processing}
   385 \subsection{Signal Processing}
   401 
   401 
   402 \begin{frame}[fragile]
   402 \begin{frame}[fragile]
   403   \frametitle{Signal \& Image Processing}
   403   \frametitle{Signal \& Image Processing}
   404   Applying a simple median filter
   404   Applying a simple median filter
   405   \begin{lstlisting}
   405   \begin{lstlisting}
   406     from scipy import signal, ndimage
   406 >>> from scipy import signal, ndimage
   407     from scipy import lena
   407 >>> from scipy import lena
   408     A=lena().astype('float32')
   408 >>> A=lena().astype('float32')
   409     B=signal.medfilt2d(A)
   409 >>> B=signal.medfilt2d(A)
   410     imshow(B)
   410 >>> imshow(B)
   411   \end{lstlisting}
   411   \end{lstlisting}
   412   Zooming an array - uses spline interpolation
   412   Zooming an array - uses spline interpolation
   413   \begin{lstlisting}
   413   \begin{lstlisting}
   414     b=ndimage.zoom(A,0.5)
   414 >>> b=ndimage.zoom(A,0.5)
   415     imshow(b)
   415 >>> imshow(b)
   416   \end{lstlisting}
   416   \end{lstlisting}
   417     \inctime{5}
   417     \inctime{5}
   418 \end{frame}
   418 \end{frame}
   419 
   419 
   420 \begin{frame}[fragile]
   420 \begin{frame}[fragile]