Numpy实现矩阵运算及线性代数应用

import numpy as np
# 1直接创建
mat=np.mat("1 2 3;4 5 6;7 8 9")
print(mat)
# 2使用numpy数组创建矩阵
mat2=np.mat(np.arange(1,10).reshape(3,3))
print(mat2)
# 3从已有的矩阵中通过bmat函数创建
A=np.eye(2)
B=A*2
mat3=np.bmat("A B;B A")
print(mat3) #类似于拼接

mat1=np.mat(np.array([2,6,5]))
mat2=np.mat(np.array([2,6,5]))

addresult=np.add(mat1,mat2)
print(addresult)

multiresult=np.multiply(mat1,mat2)
print(multiresult) #数值的乘法
mat3=np.mat(np.arange(6).reshape(2,3))
mat4=np.mat(np.arange(6).reshape(3,2))
print("mat3*mat4\n",mat3*mat4)

a=np.mat(np.array([4,5,8]))
b=np.mat(np.array([2,3,5]))
#########################除法
result1=np.divide(a,b)
print(result1)

#数组除法将结果向下去整
result2=np.floor_divide(a,b)
print("除法向下去整：",result2)
#矩阵直接相除
print("矩阵直接相除",a/b)
#矩阵取模
print("矩阵取模",a%b)

#返回小数部分
floatResult=np.true_divide(a,b)
print(floatResult)

import numpy as np
mat1=np.mat(np.array([5,-7,9]))
mat2=np.mat(np.array([2,-4,3]))
#取模运算
result1=np.remainder(mat1,mat2)
print("remainder:",result1)

result2=np.mod(mat1,mat2)
print("mod",result2)

result3=mat1%mat2
print("%",result3)

result4=np.fmod(mat1,mat2)
print(result4) #模 的正负和被除数相同，与除数无关

import numpy as np
mat=np.mat(np.array([-10,5,-4,3,0,12]))
print(np.abs(mat)) #绝对值
print(np.fabs(mat)) #非复数情况下

#获取矩阵中各元素的正负号
sign=np.sign(mat)
print("sign:",sign)
#sign: [[-1 1 -1 1 0 1]]

#将数组中元素的小数和整数部分抽取出来
arr=np.array([[1.2,3.34],
    [-2.5,6.8]])
arr1,arr2=np.modf(arr)
print("整数部分：",arr2)
print("小数部分：",arr1)

mat2=np.mat(np.array([1,2,3,4]))
mat3=np.mat(np.array([2,3,2,3]))
result1=np.power(mat2,mat3) #mat2中元素作为底，mat3作为幂
print(result1)

import numpy as np
mat3=np.mat(np.array([[1,2.2],[2,1.2]]))
mat4=np.mat(np.array([[5,2.6],[2,3.2]]))
result=np.power(mat1,mat2)
print(result) >> [[ 1.   20.57737365] [ 6.76  1.79217324]]

maximun=np.maximum(mat3,mat4)
print(maximun) >> [[5. 2.6] [2. 3.2]]

result2=np.greater(mat3,mat4)
print(result2) >> [[False True] [ True False]]
# 亦或 xor：相同的为False，不同的为True

import numpy as np
a=np.arange(4).reshape(2,2)
print(a)
#创建一个和a类型一样，但数据全是0的数组
like_a=np.zeros_like(a)
like_a=3 #赋值
print(like_a)
#################################################################################
#创建一个通用函数numpy
def like(ndarry):
 result=np.zeros_like(ndarry)
 result.flat=5
 return result
#调用numpy创建通用函数的方法,1个输入，一个输出
myfunc=np.frompyfunc(like,1,1)
test=myfunc(np.arange(9).reshape(3,3))
print(test)
#结果如下：
#[[array(5) array(5) array(5)]
# [array(5) array(5) array(5)]
# [array(5) array(5) array(5)]]

print("reduce:",np.add.reduce(a))
print("sum:",np.sum(a))

print("accumulate",np.add.accumulate(a))

print("reduceat",np.add.reduceat(a,[1,3,2,4])) # >> reduceat [3 3 5 4]
 #第一步用到索引值列表中的1和3，对数组中索引值在1到3之间的元素进行reduce操作 得到3; 
 #第二步用到索引值3和2。由于2比3小，所以直接返回索引值为3的元素 得到3;
 #第三步用到索引值2和4。对索引值在2到4之间的数组元素进行reduce操作 得到4;
 #第四步用到索引值4。对索引值从7开始直到数组末端的元素进行reduce操作 得到5;

arr1=np.array([1,6])
arr2=np.array([2,3,5,65])
result2=np.add.outer(arr1,arr2) #相加
result3=np.outer(arr1,arr2)  #相乘
print("outer:",result2)
 # outer 将第一个数组中的每个元素分别和第二个数组的所有元素相加

import numpy as np
A=np.mat(np.array([[0,1,2],[1,0,3],[4,-3,8]]))
#求A 的逆矩阵
A_=np.linalg.inv(A)
print("A的逆矩阵：\n",A_)

#验证A*A_是否是单位矩阵
print("A*A_:\n",A*A_)

import numpy as np

'''
X-2Y+Z=0
2Y-8Z=8
-4X+5Y+9Z=-9
'''
#求解三元一次函数
A=np.mat("1 -2 1;0 2 -8;-4 5 9")
print("系数：",A)
B=np.array([0,8,-9])
print("常数:",B)

#调用numpy的solve方法
result=np.linalg.solve(A,B)
print("X={},Y={},Z={}".format(result[0],result[1],result[2]))

import numpy as np
vector=np.mat("3 -2;1 0")
#求特征值
eigenvalues=np.linalg.eigvals(vector)
# 特征值是Ax=ax的根
eigenvalues,eigvector=np.linalg.eig(vector)
print("特征值：",eigenvalues)
print("特征向量：\n",eigvector)
>>结果：
特征值： [2. 1.]
特征向量：
 [[0.89442719 0.70710678]
 [0.4472136 0.70710678]]

import numpy as np
vector=np.mat("4 11 14;8 7 -2")
#调用numpy汇总的svd方法对矩阵进行奇异值分解
U,sigma,V=np.linalg.svd(vector,full_matrices=False)
print("U\n",U)
print("sigma:\n",sigma)
print("V:\n",V)

#将svd分解出的值相乘
print("vector:\n",U*np.diag(sigma)*V)

import numpy as np
vector=np.mat("3 4;5 6")
print(vector)
#求行列式
value=np.linalg.det(vector)
print("行列式：",value)

import numpy as np
arr=np.array([1,2,34,5])
result=np.sort(arr)
arr.sort()
print(arr)

argsort=np.argsort(arr)
print("argsort:",argsort)
d2=np.array([[12,3,45,2],
    [43,552,1,9]])
#将多维数组按列排序
d2.sort(axis=0)
print(d2)
#[[ 12 3 1 2]
# [ 43 552 45 9]]
#将多维数组按行排序
d2.sort(axis=1)
print(d2)

print(np.argsort(d2,axis=0))
[[0 0 0 0]
 [1 1 1 1]]

print(np.argsort(d2,axis=1))
[[2 3 1 0]
[3 0 2 1]]

import numpy as np
a=np.array([[2,3,653,4],
   [5,62,943,44,]])
b=np.array([32,13,65,43])
#如果是多维数组，则将多维数组展平后获取最大值的下标
argmax=np.argmax(a)
print(argmax)

#根据条件数组中搜索非零的元素，分析返回对应的下标
print("argwhere\n",np.argwhere(b>40))
print("argwhere 多维\n",np.argwhere(a>100))
'''
argwhere 多维
 [[0 2]
 [1 2]]
'''

sorted=np.arange(5)
indices=np.searchsorted(sorted,[-2,7])
print(indices)      #下标

newsorted=np.insert(sorted,indices,[-2,7])
print("添加新元素后：\n",newsorted)

import numpy as np
a=np.arange(10)
#生成一个抽取元素的花式索引
condition=a%2==0#赋值/算术/逻辑运算符优先级
print("花式索引",condition)
even=np.extract(condition,a)
print("偶数项：",even)
#take()/compress()
even2=a.compress(condition)

import numpy as np
#抽取数组中非零的元素np.nonzero() 的索引
arr=np.array([[0,1,2],[0,3,4],[0,5,6]])
rows,cols=np.nonzero(arr)
print("rows",rows)
print("cols:",cols)

indices=np.dstack((rows,cols))
print("indices:\n",indices)

import numpy as np
#某用户去银行存款，假设年利率是3%,，每季度存10元，
# 存5年以及存款1000，五年后可得到的本息和是多少？

#g该用户5年后得到的本息为
fv=np.fv(0.03/4,5*4,-10,-1000)
print("5年后本息和：",fv)
#计算每一年的本息和
for i in range(1,6):
 fv1=np.fv(0.03/4,i*4,-10,-1000)
 print("第",i,"年本息和为：{}".format(fv1))

import numpy as np
#存五年可得1376.0963320407982，计算5年前存款的金额
pv=np.pv(0.03/4,5*4,-10,1376.0963320407982)
print("5年钱第一次存款的金额是：",np.round(-pv))

#投资100，支出39，59，55，20，折现率28.1%，净现值多大？
npv=np.npv(0.281,[-100,39,59,55,20])
print("净现值：",npv)

#某人贷款20万，15年还清，年利率是7.5%，求每月还的金额？
pmt=np.pmt(0.075/12,15*12,200000)
print("每月还款的金额：{}".format(pmt))

#某人贷款20万，每月还2000，年利率是7.5%，求多少年还完？
nper=np.nper(0.075/12,-2000,200000)
months=np.ceil(nper)
years=np.ceil(months/12)#向上取整
print("需要还款{}年".format(years)) >>14年