数理基础（线性代数）------矩阵常见的运算及统计函数

一.数学相关函数

1.函数的平方根、立方根

import numpy as np
a=np.arange(9)
print(a)
print('a的平方根',np.sqrt(a))#a的平方根
print('a的平方根',np.power(a, 1/2))#a的平方根
print('a的立方根',np.power(a, 1/3))#a的立方根
print('a的立方根',np.power(a, 1/3))#a的立方根
print('a的平方',np.square(a))#a的平方
print('a的平方',np.power(a, 2))#a的平方
print('a的平方',a**2)#a的平方

2.指数函数和对数函数

import numpy as np
print(np.exp(2))
print(np.log2(8))
print(np.log10(100))
loga(b)=logc(b)/logc(a)

3. 绝对值

import numpy as np
print(np.abs(a))
print(np.abs(-7))
print('='*20)

4.符号函数

import numpy as np
print(np.sign(a))

5.比较大小

import numpy as np
fa=np.random.randint(10,30,size=12)
fb=np.random.randint(10,30,size=12)
print(fa)
print(fb)
print(np.maximum(fa, fb))
print(np.minimum(fa,fb))

二.矩阵常用的统计函数

求和sum
平均值mean
标准差std
方差var
最大值max
最小值min
中位数median
argmin返回元素最小值索引
argmax返回元素最大值索引
cumsum所有元素的累计和
cumprod所有元素的累计积

#定义一个4行4列的矩阵
a=np.arange(1, 17).reshape(4, 4)
print(a)

1.求和

print(np.sum(a))#求和
np.sum(a,axis=0)#按列求和
np.sum(a,axis=1)#按行求和

2.最大值

print(np.max(a))#求最大值
print(np.max(a, axis=0))#求列的最大值
print(np.max(a, axis=1))#求行的最大值
print(np.argmax(a))#最大值的索引
print(np.argmax(a, axis=0))#最大值的列索引
print(np.argmax(a, axis=1))#最大值的行索引

3.最小值

print(np.min(a))#求最小值
print(np.min(a, axis=0))#求列的最小值
print(np.min(a, axis=1))#求行的最小值
print(np.argmin(a))#最小值的索引
print(np.argmin(a, axis=0))#最小值的列索引
print(np.argmin(a, axis=1))#最小值的行索引

4.平均值

print(np.mean(a))#求平均值
print(np.mean(a, axis=0))#求列的平均值
print(np.mean(a, axis=1))#求行的平均值

5.中位数

print(np.median(a))#中位数
print(np.median(a, axis=0))#求列的中位数
print(np.median(a, axis=1))#求行的中位数

6.叠加

np.cumsum(a) #第i个数就是前i个数相加
np.cumsum(a,axis=0)
np.cumsum(a,axis=1)

7.叠乘

np.cumprod(a) #第i个数就是前i个数相乘
np.cumprod(a,axis=0)
np.cumprod(a,axis=1)

8.排序

#一维数组排序
fb=np.random.randint(1,16,16)
print(fb)
fa=np.sort(fb)#正序排序
print(fa)
print(fa[::-1])#倒序排序

#二维数组排序
a.sort(axis=0)
a.sort(axis=1)

三. 矩阵和数组的创建

#方法一（创建矩阵）
M=np.mat('1 2 3;4 5 6;7 8 9')
print(M)
print(type(M))
#方法二（创建矩阵）
N=np.matrix([[1,2,3],[4,5,6],[9,8,7]])
print(N)
print(type(N))
#方法三（创建数组）
l=np.array([[1,8,3],[4,5,6],[9,8,7]])
print(l)
print(type(l))

四. 数组的保存和加载

# 数据的保存和加载
A=np.zeros((3,3))
B=np.ones((3,3))

#保存数据
np.savetxt('1912A',A.dot(B),delimiter=',')#参数顺序：1.文件名2.传入矩阵 3.用什么分隔

#加载数据
data=np.loadtxt('1912A',delimiter=',')#参数顺序：文件名 用什么分隔
print(data)