# -*- coding: utf-8 -*-
#2.3.2
import numpy as np
#生成一个向量
print(np.random.randint(0,10,10))
#size 向量大小
print(np.random.randint(0,5,size=5));
#生成矩阵
print(np.random.randint(4,9,size=(3,5)))
#seed 每次生成的随机数是固定的，随机种子用1记录
#生成介于0~1之间的浮点数的向量或者矩阵
print(np.random.random(10))
print(np.random.random((3,4)))
#random.normal 正态分布
#loc正态分布均值 scale正态分布标准差 sieze大小
print(np.random.normal(loc=0,scale=5,size=(3,5)))
#2.3.3 获取numpy属性
# =============================================================================
# 生成数组
# =============================================================================
m_matrix = np.arange(15)
m_matrix=m_matrix.reshape(3,5)
print(m_matrix.shape)
#EX
x = np.arange(15)
print(x.ndim)
#转换成3行5列的2维矩阵
X=x.reshape(5,3)
print(X.ndim)
#只关心有多少列或行
print(x.reshape(5,-1))
print(x.reshape(-1,15))
# =============================================================================
# 切片操作
# =============================================================================
print(X)
print(X[:1])
print(X[:,0:2])
print(X[1:3,:])
print(X[1:3,0:2])
# =============================================================================
# numpy矩阵运算
# =============================================================================
myones = np.ones([3,3])
#生成一个对角线值为1，其余值为0的三行三列矩阵
myeye = np.eye(3)
print(myeye)
print(np.sin(myones - myeye))
#矩阵点乘，满足 n*m * m*l = n*l，矩阵乘法函数为dot
mymatrix = np.linspace(1,6,6,dtype=(int)).reshape(2,3)
a = np.linspace(1,6,6,dtype=int).reshape(3,2)
print(mymatrix.shape[1] == a.shape[0])
print(mymatrix.dot(a))
# =============================================================================
# 矩阵转置,逆矩阵
# =============================================================================
import numpy.linalg as lg
A = np.array([[i for i in range(2)],[i+2 for i in range(2)]])
print(A)
invA = lg.inv(A)
print(invA)
print(A.dot(invA))
# =============================================================================
# 数据类型转换
# =============================================================================
vector = np.array(['1','2','3'])
print("String to float = ",vector.astype(float),'\n\n\n')
# =============================================================================
# Numpy内置方法
# =============================================================================
#sum mean max median 中位数
#类型必须为int或者float
vector = np.array([5,10,15,20])
print('the np.sum=',vector.sum(),'\n\n\n')
#矩阵例子
matrix = np.array([[5,10,15],[20,10,30],[35,40,45]])
print(matrix)
print(matrix.T)
masu1=matrix.sum(axis=1)
print('matrix.sum axis1 为1时为行和，输出为列 =',masu1)
print(masu1.ndim)
print('matrix.sum axis0 =',matrix.sum(axis=0),'\n\n\n')
# =============================================================================
# np中的arg运算
# =============================================================================
index2 = np.argmax([1,2,6,3,2])
print('the max of this matrix = ',index2)
#例子
x = np.arange(15)
print(x)
np.random.shuffle(x)#随机打乱顺序
print(x)
sx = np.argsort(x)
print('the index of this arg vector=',sx)
# =============================================================================
# FancyIndexing
# =============================================================================
x = np.arange(15)
np.random.shuffle(x)
ind = [3,5,8]
print(x[ind],'\n\n')
ind = np.array([[0,2],[1,3]])
print(x,'\n')
print(x[ind],'\n')
#二维矩阵取数
x = np.arange(16)
X =x.reshape(4,-1)
row = np.array([0,1,2])
col = np.array([1,2,3])
print(X)
print(X[row,col])
print(X[1:3,col])
# =============================================================================
# np数组比较
# =============================================================================
matrix = np.array([[5,10,15],[20,25,30],[35,40,45]])
m = (matrix == 25)
print(m,'\n')
second_column_25 = (matrix[:,1] == 25)
print(second_column_25,'\n')
print(matrix[second_column_25,:],'\n')
#计算符合条件的总数np.count_nonzero
print(np.count_nonzero(matrix<=28))
#计算是否符合条件
print(np.any(matrix == 25))
#计算是否所有数据都符合条件
print(np.all(matrix==25),'\n\n')