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} |
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] |