5.16作业（numpy练习）

Exercise 0

题意

题解

randn(d0, d1, ..., dn) 返回一个n维数组，数组是标准正态分布的样本。用randn函数得到二维数组A。

import numpy
from numpy import random

n = 200
m = 500
A = random.randn(n, m)

numpy库中没有toeplitz函数，scipy.linalg库有。其实可以自己写一个函数生成Toeplitz数组。什么是Toeplitz矩阵？

def toeplitz(row_vector, column_vector):
    '''
    row_vector是 1×m 的二维数组
    column_vector是 m×1 的二维数组
    '''
    res = row_vector
    for i in range(1, column_vector.size):
        row_vector = numpy.roll(row_vector, 1)
        row_vector[0,0] = column_vector[i, 0]
        res = numpy.vstack((res, row_vector))
    return res

测试 toeplitz 函数

vec1 = numpy.array([[1, 2, 3, 4, 5]])
vec2 = numpy.array([[1], [6], [7], [8], [9]])
print(toeplitz(vec1, vec2))

结果

得到二维数组B

B = toeplitz(random.randn(1, m), random.randn(m, 1))

Exercise 1: Matrix operations

题意

题解