PythonAnIntroduction分析.pptVIP

* * * * * * * * * * * * * * * * * * * * 10/30/20089/25/20089/11/2008 * 9/25/20089/11/2008 * Shell Scripting 9/11/2008 Research Computing Center, UNC-CH * * In order for create a section divider slide, add a new slide, select the “Title Only” Slide Layout and apply the “Section Divider” master from the Slide Design menu. For more information regarding slide layouts and slide designs, please visit /training Eigenvalues from LinearAlgebra import * val = eigenvalues(c) val, vec = eigenvectors(c) print val [1., 4., 1.] print vec [[0.816, -0.408, -0.408], [0.575, 0.577, 0.577], [-0.324, -0.487, 0.811]] also solve_linear_equations, singular_value_decomposition, etc. Least Squares Fitting Part of Hinsens Scientific Python module from LeastSquares import * def func(params,x): # y=ax^2+bx+c return params[0]*x*x + params[1]*x + params[2] data = [] for i in range(10): data.append((i,i*i)) guess = (3,2,1) fit_params, fit_error = leastSquaresFit(func,guess,data) print fit_params [1.00,0.000,0.00] FFT from FFT import * data = array((1,0,1,0,1,0,1,0)) print fft(data).real [4., 0., 0., 0., 4., 0., 0., 0.]] Also note that the FFTW package (fastest Fourier transform in the West) has a python wrapper. See notes at the end Python Standard Libraries/Modules: /library/ /dci/dci_components/shared_apps/packages/python_packages.html /pypi/ Command-line arguments import sys print len(sys.argv) # NOT argc # Print all arguments: print sys.argv # Print all arguments but the program # or module name: print sys.argv[1:] # array slice Catching Exceptions #python code a.py x = 0 try: print 1/x except ZeroDivisionError, message: print Can’t divide by zero: print message python a.py Cant divide by zero: integer division or modulo by zero Try-Finally: Cleanup f = open(file) try: process_file(f) finally: f.close() # always executed print OK # executed on success only Raising Exceptions raise IndexError raise IndexError(k out of range) rais