第十三周作业——高级编程技术

import numpy as np  
from scipy.linalg import lstsq  
from scipy.linalg import norm  

A = np.random.randint(0,50,size=(5,5))  
b = np.random.randint(0,50,size=(5,1))  
print("A:")  
print(A)  
print("\nb:")  
print(b)  
   
x, res, rnk, s = lstsq(A, b)  
y = norm(b-np.dot(A,x),ord=2)   
print("\nSolution x for Ax=b:") 
print(x)  
print("\ny:")  
print(y)  

import scipy.optimize as op
import numpy as np   
import math  
  
def fun(x):  
    return (-(math.sin(x-2)**2)*math.exp(-(x**2)))  
      
Max = op.minimize_scalar(fun)  
  
print(-Max.fun) 

from scipy.spatial.distance import cdist
import numpy as np    
import math  
 
X = np.random.randint(0,50,(5,5))  
print(X)
print("\nSolution:")
print(cdist(X,X,metric="euclidean"))