scipy练习

10.1

import numpy as np;
import scipy.optimize as oz;
import scipy.linalg as lin

m=20
n=10

A = np.random.random(size=(m,n))
b = np.random.random(m)

x = oz.lsq_linear(A,b).fun
print(lin.norm(x))

10.2

import numpy as np
import scipy.optimize as oz

def f(x):
    return (-1) * (np.sin(x-2))**2 * np.exp((-1)*x**2)

res = oz.minimize_scalar(f)
print(-1*f(res.x))

10.3

import numpy as np
import scipy.spatial.distance as dis;
np.random.seed()
n = 5
index=0
A = np.random.randint(1,20,size=(n,2))
print('city:')
print(A)
print('\n')
D = dis.pdist(A)
print('          ',end='')
for city in A:
    print('%-20s' % str(tuple(city)),end='')
print('\n')

for i in range(0,n):
    print('%-10s' % str(tuple(A[i])),end='')
    for j in range(0,n):
        if i == j:
            print('%-20s' % '0', end='')
        elif i<j:
            index=((n-1)+(n-i))*i//2+j-i-1
            print('%-20s' % str(D[index]), end='')
        else:
            index=((n-1)+(n-j))*j//2+i-j-1
            print('%-20s' % str(D[index]), end='')
    print('\n')