2018年5月24日
摘要: >>> from scipy.interpolate import interp1d#interp1d表示1维插值 >>> from scipy.interpolate import interp1d#interp1d表示1维插值 >>> >>> x = np.linspace(0, 10, num 阅读全文
posted @ 2018-05-24 16:00 裸睡的猪 阅读(4559) 评论(0) 推荐(0) 编辑
摘要: #optimazer优化器 from scipy.optimize import minimize def rosem(x): return sum(100.0*(x[1:]-x[:-1])**2.0+(1-x[:-1])**2.0) x0=np.array([1.3,.7,.8,1.9,1.2]) res=minimize(rosem,x0,method="nelder-mead",o... 阅读全文
posted @ 2018-05-24 15:51 裸睡的猪 阅读(1515) 评论(0) 推荐(0) 编辑
摘要: #!/usr/bin/env python # -*- coding: utf-8 -*- # @Time : 2018/5/24 15:03 # @Author : zhang chao # @File : s.py #integral积分 import numpy as np from scipy.integrate import quad,dblquad,nquad pri... 阅读全文
posted @ 2018-05-24 15:21 裸睡的猪 阅读(749) 评论(0) 推荐(0) 编辑
摘要: #scatter fig=plt.figure() ax=fig.add_subplot(3,3,1)#3行3列 第一个图 n=128 X=np.random.normal(0,1,n) Y=np.random.normal(0,1,n) T=np.arctan2(Y,X)#T用来上色的 #plt. 阅读全文
posted @ 2018-05-24 14:56 裸睡的猪 阅读(423) 评论(0) 推荐(0) 编辑
摘要: fig=plt.figure() fig.add_subplot(3,3,1)#3行3列 第一个图 n=128 X=np.random.normal(0,1,n) Y=np.random.normal(0,1,n) T=np.arctan2(Y,X)#T用来上色的 plt.axes([0.025,0.025,0.95,0.95])#显示的范围 plt.scatter(X,Y,s=75,c=T,a... 阅读全文
posted @ 2018-05-24 14:10 裸睡的猪 阅读(484) 评论(0) 推荐(0) 编辑
摘要: import numpy as np import matplotlib.pyplot as plt from scipy.constants.constants import alpha from matplotlib.patches import ArrowStyle x=np.linspace(-np.pi,np.pi,256,endpoint=True) #np.linspace定义横轴... 阅读全文
posted @ 2018-05-24 09:31 裸睡的猪 阅读(216) 评论(0) 推荐(0) 编辑