Lecture_TWP_Python_A04_1a_Linear_Algebra
Lecture_TWP_Python_A04_1a_Linear_Algebra
Nirjhar Dhang
√
Version : 1.732 : 3
Created on : August 11, 2016
Last revision : January 11, 2025
Together With Python
Computation, Visualization and Reporting
Part A-4 : Linear Algebra
Nirjhar Dhang
√
Version : 1.732 : 3
Created on : August 11, 2016
Last revision : January 11, 2025
Overview
1 Python numpy
2 Arrays in numpy
3 Summary
4 References
Python numpy
Python numpy
import numpy as np
Output
A:
[10 20 30]
B:
[[10 20 30]
[40 50 60]]
C:
[10 20 30]
Arrays in numpy
import numpy as np
A=np.array([[10.0,20.0,30.0],
[40.0,50.0,60.0],
[70.0,80.0,90.0]])
print(’A:’)
print(A)
print(’Transpose of A:’)
print(A.T)
Arrays in numpy
Output
A:
[[10. 20. 30.]
[40. 50. 60.]
[70. 80. 90.]]
Transpose of A:
[[10. 40. 70.]
[20. 50. 80.]
[30. 60. 90.]]
Arrays in numpy
import numpy as np
A=np.array([[6.0,2.0,1.0],
[2.0,7.0,1.0],
[1.0,1.0,8.0]])
B=np.array([9.0,10.0,10.0])
C=np.dot(A,B)
print(’A:’)
print(A)
print(’B:’)
print(B)
print(’C:’)
print(C)
Arrays in numpy
Output
A:
[[6. 2. 1.]
[2. 7. 1.]
[1. 1. 8.]]
B:
[ 9. 10. 10.]
C:
[84. 98. 99.]
Arrays in numpy
import numpy as np
A=np.array([[6.0,2.0,1.0],
[2.0,7.0,1.0],
[1.0,1.0,8.0]])
B=np.array([9.0,10.0,10.0])
C=np.linalg.solve(A,B)
print(’A:’)
print(A)
print(’B:’)
print(B)
print(’C:’)
print(C)
Arrays in numpy
Output
A:
[[6. 2. 1.]
[2. 7. 1.]
[1. 1. 8.]]
B:
[ 9. 10. 10.]
C:
[1. 1. 1.]
Arrays in numpy
import numpy as np
import numpy.linalg as la
A=np.array([[2.0,1.0],
[1.0,2.0]])
Eva,Eve=la.eig(A)
print(’A:’)
print(A)
print(’Eigen values:’)
print(Eva)
print(’Eigen vectors:’)
print(Eve)
Arrays in numpy
Output
A:
[[2. 1.]
[1. 2.]]
Eigen values:
[3. 1.]
Eigen vectors:
[[ 0.70710678 -0.70710678]
[ 0.70710678 0.70710678]]
Arrays in numpy
import numpy as np
import numpy.linalg as la
A=np.array([[6.0,2.0,1.0],
[2.0,7.0,1.0],
[1.0,1.0,8.0]])
Eva,Eve=la.eig(A)
print(’A:’)
print(A)
print(’Eigen values:’)
print(Eva)
print(’Eigen vectors:’)
print(Eve)
Arrays in numpy
Output
A:
[[6. 2. 1.]
[2. 7. 1.]
[1. 1. 8.]]
Eigen values:
[9.71447874 4.42879858 6.85672268]
Eigen vectors:
[[-0.49183142 -0.79848935 -0.34715503]
[-0.59615502 0.5994227 -0.53412697]
[-0.6345873 0.05574221 0.77083835]]
Arrays in numpy
import numpy as np
A=np.zeros((3,3),float)
A[0][0]=6.0
A[0][1]=2.0
A[0][2]=1.0
A[1][0]=2.0
A[1][1]=7.0
A[1][2]=1.0
A[2][0]=1.0
A[2][1]=1.0
A[2][2]=9.0
Arrays in numpy
B=np.zeros((3),float)
B[0]=9.0
B[1]=10.0
B[2]=11.0
C=np.linalg.solve(A,B)
print(’A:’)
print(A)
print(’B:’)
print(B)
print(’C:’)
print(C)
Arrays in numpy
Output
A:
[[6. 2. 1.]
[2. 7. 1.]
[1. 1. 9.]]
B:
[ 9. 10. 11.]
C:
[1. 1. 1.]
Arrays in numpy
import numpy as np
import numpy.linalg as la
A=np.zeros((3,3),float)
A[0][0]=6.0
A[0][1]=2.0
A[0][2]=1.0
A[1][0]=2.0
A[1][1]=7.0
A[1][2]=1.0
A[2][0]=1.0
A[2][1]=1.0
A[2][2]=9.0
Arrays in numpy
Eva,Eve=la.eig(A)
print(’A:’)
print(A)
print(’Eigen values:’)
print(Eva)
print(’Eigen vectors:’)
print(Eve)
Arrays in numpy
Output
A:
[[6. 2. 1.]
[2. 7. 1.]
[1. 1. 9.]]
Eigen values:
[10.20440155 4.43115039 7.36444806]
Eigen vectors:
[[-0.41760969 -0.79624527 -0.43771637]
[-0.49785133 0.60350126 -0.6228405 ]
[-0.76009618 0.04218655 0.64843974]]
Summary
Summary