爬虫python编程与cvi编程_与爬虫无关，简单的用python进行科学运算

1.算积分

1.1一重积分

例：如计算半圆的面积，积分公式：

# -*- coding:utf-8 -*-

from scipy import integrate#输入scipy的integrate函数

def fun(x):

return (1-x*x)**0.5#定义的半圆函数

p,err=integrate.quad(fun,-1,1)#传入函数，积分范围-1到1

print p

1.2二重积分，算体积

例：求上半圆的体积，积分公式：

# -*- coding:utf-8 -*-

from scipy import integrate

def fun(x,y):

return (1-x*x-y*y)**0.5#对函数的定义

def fun2(x):

return (1-x*x)**0.5#定义ｙ对ｘ的边界函数

m,err=integrate.dblquad(fun, -1, 1,

lambda x:-fun2(x),

lambda x:+fun2(x))#对x的取值范围-1到1，对y的取值范围-f(x)到f(x)

print m

2.解非线性方程组

例：求解非线性方程组

# -*- coding:utf-8 -*-

from scipy.optimize import fsolve

from math import sin

def func(x):

x0,x1,x2=x.tolist()#定义变量是个浮点数

return[

5*x1+3,

4*x0*x0-2*sin(x1*x2),

x1*x2-1.5

]#定义公式

r=fsolve(func,[1,1,1])#1,1,1,指的是传入可能的解，通常1,1,1

print r

3.python绘图

快速绘图：

# -*- coding:utf-8 -*-

import numpy as np

import matplotlib.pyplot as plt

from matplotlib.font_manager import FontProperties

font = FontProperties(fname=r"c:\windows\fonts\simsun.ttc", size=14)#设置中文字体

x=np.linspace(0,10,100)#生成数列

y=np.sin(x)

z=np.cos(x)

plt.figure(figsize=(7,5))#框的尺寸大小

plt.plot(x,y,label='$sin(x)$')#

plt.plot(x,z,"r--",label="$cos(x)$")#r蓝色，虚线

plt.xlabel(u"时间(s)",fontproperties=font)#横坐标标签

plt.ylabel(u"y值",fontproperties=font)

plt.title(u"数学",fontproperties=font)#标题

plt.ylim(-1.2,1.2)#取值范围

plt.legend()#图示

plt.show()

结果显示：

3.1 多轴绘图，各种坐标系的绘制

# -*- coding:utf-8 -*-

import matplotlib.pyplot as plt

import numpy as np

w=np.linspace(0.1,1000, 1000)

p=np.abs(1/(1+0.1j*w))

plt.subplot(221)#m表示是图排成m行，n表示图排成n列,p表示图所在的位置

plt.plot(w,p,linewidth=2)#linewidth线宽=2

plt.ylim(0,1.5)

plt.subplot(222)#m表示是图排成m行，n表示图排成n列,p表示图所在的位置

plt.semilogx(w,p,linewidth=2)#semilogx: x轴为对数刻度，y轴为线性刻度

plt.ylim(0,1.5)

plt.subplot(223)#m表示是图排成m行，n表示图排成n列,p表示图所在的位置

plt.semilogy(w,p,linewidth=2)#y轴为对数坐标系

plt.ylim(0,1.5)

plt.subplot(224)#m表示是图排成m行，n表示图排成n列,p表示图所在的位置

plt.loglog(w,p,linewidth=2)#双对数坐标系

plt.ylim(0,1.5)

plt.show()

结果显示：

3.2绘制条形图

# -*- coding:utf-8 -*-

import matplotlib.pyplot as plt

import numpy as np

from matplotlib.font_manager import FontProperties

font = FontProperties(fname=r"c:\windows\fonts\simsun.ttc", size=14)

def autolabel(rects):

for rect in rects:

height = rect.get_height()

plt.text(rect.get_x()+rect.get_width()/2., 1.03*height, '%s' % float(height))

x=np.array([1,2,3,4,5])

y=np.array([2,3,4,5,5])

plt.figure()

#bar绘制条形图

rect=plt.bar(x,y,width=0.35,align="center",label="people",yerr=0.000001)

plt.xlabel(u"年龄",fontproperties=font)

plt.ylabel("num")

plt.legend()

autolabel(rect)#显示y的数值

plt.show()

结果显示：

3.3绘制散点图

# -*- coding:utf-8 -*-

import numpy as np

import matplotlib.pyplot as plt

x=np.random.random(100)#在0到1生成100个随机数

y=np.random.random(100)

plt.figure(figsize=(8,4))

plt.scatter(x, y, c="black")

plt.show()

3.4绘制3D图形

import numpy as np

import matplotlib.pyplot as plt

from mpl_toolkits.mplot3d import Axes3D

th=np.linspace(-4*np.pi, 4*np.pi, 100)

z=np.linspace(-2,2,100)

r=z*2+1

x=r*np.sin(th)

y=r*np.cos(th)

fig=plt.figure()

ax=fig.add_subplot(111,projection='3d')

ax.plot(x,y,z)

plt.show()