- 博客(10)
- 资源 (31)
- 收藏
- 关注
原创 python 插值
对4个数据点x = [-1, 0, 2.0, 1.0],y = [1.0, 0.3, -0.5, 0.8]进行Rbf插值,插值中使用三种插值方法分别是multiquadric、gaussian、和linear(参见课件5,scipy_rbf.py),需要作图点(加密点)为np.linspace(-3, 4, 100)。RBf插值代码如下:#导入模块import numpy as npfrom ...
2018-06-08 11:37:59 5081
原创 python 实现最小二乘
用curve_fit()函数对高斯分布进行拟合,xϵ[0,10],高斯分布函数为y=a*np.exp(-(x-b)**2/(2*c**2)) ,其中真实值a=1,b=5,c=2。试对y加入噪声之后进行拟合, 并作图与真实数据进行比较。(参见课件leastsq(),curve_fit()拟合)用最小二乘进行拟合leastsq()#导入包和相应的模块import numpy as npfrom s...
2018-06-08 11:36:12 1560
原创 python实现最小二乘拟合
用curve_fit()函数对高斯分布进行拟合,xϵ[0,10],高斯分布函数为y=a*np.exp(-(x-b)**2/(2*c**2)) ,其中真实值a=1,b=5,c=2。试对y加入噪声之后进行拟合, 并作图与真实数据进行比较。(参见课件leastsq(),curve_fit()拟合)...
2018-06-08 11:34:22 2364
原创 python scipy求解非线性方程组
求解非线性方程组,cos(a) = 1 - d^2 / (2*r^2) ,L = a * r,d = 140,L= 156; 导入参数雅克比矩阵, 再次进行求解。A、求解非线性方程组:#导入优化模块和余弦函数from scipy.optimize import fsolvefrom math import cos#定义函数def f(x): d=140 l=156 a,r=x...
2018-06-08 11:33:17 7261
原创 python numpy 随机抛硬币
使用二项分布进行计算. 同时抛弃5枚硬币,如果正面朝上少于3枚,则输掉8元,否则就赢8元。如果手中有1000元作为赌资,请问赌博10000次后可能会是什么情况呢?代码实现:import numpy as np#导入画图模块from matplotlib.pyplot import plot,show#导入matplotlib模块的画图和显示函数money=np.zeros(10000)#先将...
2018-06-08 11:27:53 4066 1
原创 python numpy 切比雪夫插值
1、对g(x)在100个切比雪夫节点之上分別使用Polynomial(Polynomial.fit)和Chebyshev(Chebyshev.fit)进行插值,在[-1,1]区间上取1000个等距点对误差进行比较。g(x)= sin(x**2) + sin(x)**2, 其中x=(x - 1) * 5。代码实现:import numpy as npfrom numpy.polynomial i...
2018-06-08 11:25:16 6217
原创 python numpy学习查找最近值
找出数组中与给定值最接近的数(通用方法)。(例:任一数组Z=array([[0,1,2,3],[4,5,6,7]]),给任一定值z=5.1)import numpy as npZ=np.array([[0,1,2,3],[4,5,6,7]])print(Z)x=5.1a=abs(Z-x)print(a)b=a.min()c=abs(x-b)print(c)输出结果:[[01 2 3] [4 5...
2018-06-08 11:18:42 8518 2
原创 python numpy学习3
归一化,将矩阵规格化到0~1,即最小的变成0,最大的变成1,最小与最大之间的等比缩放。Z = np.random.random((5,5))import numpy as npZ = np.random.random((5,5))print(Z)max=Z.max()min=Z.min()print(max,min)Z=(Z-min)/(max-min)print(Z)输出结果:[[0.253...
2018-06-08 11:17:10 1064
原创 python numpy 学习2
在数组[1, 2, 3, 4, 5]中每相邻两个数字中间插入两个0。代码实现:import numpy as npa=np.array([1,2,3,4,5])print(a)x=2b=np.zeros(len(a)+(len(a)-1)*x)b[::x+1]=aprint(b)输出结果为:[1 2 3 4 5][ 1. 0. 0. ..., 0. 0. 5.]...
2018-06-08 11:15:31 2686
原创 pyhton numpy学习
1、arr11= 5-np.arange(1,13).reshape(4,3), 计算所有元素及每一列的和;对每一个元素、每一列求累积和;计算每一行的累计积;计算所有元素的最小值;计算每一列的最大值;计算所有元素、每一行的均值;计算所有元素、每一列的中位数;计算所有元素的方差,每一行的标准差。代码实现:#1、所有元素的和import numpy as nparr11 = 5-np.arange...
2018-06-08 11:14:20 642
Kubernetes 集群安全 - 机制说明.pdf
2020-10-17
kubernetes学习资料.rar
2020-10-14
PostgreSQL高可用方案比较说明书.pdf
2020-10-10
annie视频下载神器源代码
2020-10-06
免费windows录频截屏软件
2020-10-05
Patroni-2.0.0说明书.md
2020-09-28
origin7.5汉化版x64版
2018-12-15
23种设计模式迷你小册子
2018-11-29
python numpy数值计算包
2018-10-20
python 最好的统计包matplotlib.chm
2018-10-19
分子模拟经典书籍
2018-10-14
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人