scipy 积分求解与绘图工具matplotlib.pyplot及scipy.misc 使用
文章目录
import numpy as np
import pandas as pd
from pandas import Series,DataFrame
# scipy求解积分
from scipy.integrate import quad
# 绘图专用的工具
import matplotlib.pyplot as plt
# 表示把matplotlib生成的图像绘制到当前notebook当中
%matplotlib inline
求积分
确认不规则图形的函数曲线
f = lambda x:(1-x**2)**0.5
确认不规则图形x轴的界限(-1,1)
b = 1
a = -1
定义想x,y的函数
x = np.linspace(-1,1,100)
y = f(x)
绘制线性图
plt.plot(x,y)
plt.plot(x,-y)
调整坐标轴刻度的比例一致
plt.axis('equal')
确认结果是否正确,结果为pi,即为正确
half_area,error = quad(f,a,b)
2*half_area/1
结果为:
3.1415926535897967
绘制正弦曲线
x1 = np.linspace(0,2*np.pi,100)
y1 = np.sin(x1)
plt.plot(x1,y1)
文件输入输出
处理文本输入
import scipy.io as io
处理图片输入输出
import scipy.misc as misc
把数据存储为2进制格式的数据
io.savemat()
读取二进制的数据文件
io.loadmat()
数据写入,支持数组对象
xiaoming = np.random.randint