Python3 科学计算 scipy

方程组求解

最小二乘

import numpy as np
import pylab as pl
from scipy import stats

a = [
1, 2, 3, 4
]
p = [
0.1, 0.2, 0.3, 0.4
]

# 计算加权平均
print(np.average(a,weights=p))

X = stats.rv_discrete(values=(a, p))

print(X.rvs(size=10))

# 生成2000*100的样本
rvs = X.rvs(size=(2000, 100))

# 验证大数定理
print(np.mean(rvs, axis=1))
print(np.mean(rvs, axis=0))

def half_circle(x):
return (1 - x ** 2) ** 0.5

from scipy import integrate

pi_half, err = integrate.quad(half_circle, -1, 1)
print(pi_half * 2)
3.141592653589797