使用Python做科学计算初探(转)
今天在搞定Django框架的blog搭建后,尝试一下python的科学计算能力。
python的科学计算有三剑客:numpy,scipy,matplotlib。
numpy负责数值计算,矩阵操作等;
scipy负责常见的数学算法,插值、拟合等;
matplotlib负责画图。
首先,百度上头三个,依次安装。
可以考虑使用pyhton34/script/easy-install 工具; easy-insatll -m matplotlib;
尝试一下代码,拟合实例;
1 # -*- coding: utf-8 -*-
2 import numpy as np
3 from scipy.optimize import leastsq
4 import pylab as pl
5
6 def func(x, p):
7 # """
8 # 数据拟合所用的函数: A*sin(2*pi*k*x + theta)
9 # """
10 A, k, theta = p
11 return A*np.sin(2*np.pi*k*x+theta)
12
13 def residuals(p, y, x):
14 # """
15 # 实验数据x, y和拟合函数之间的差,p为拟合需要找到的系数
16 # """
17 return y - func(x, p)
18
19 x = np.linspace(0, -2*np.pi, 100)
20 A, k, theta = 10, 0.34, np.pi/6 # 真实数据的函数参数
21 y0 = func(x, [A, k, theta]) # 真实数据
22 y1 = y0 + 2 * np.random.randn(len(x)) # 加入噪声之后的实验数据
23
24 p0 = [7, 0.2, 0] # 第一次猜测的函数拟合参数
25
26 # 调用leastsq进行数据拟合
27 # residuals为计算误差的函数
28 # p0为拟合参数的初始值
29 # args为需要拟合的实验数据
30 plsq = leastsq(residuals, p0, args=(y1, x))
31
32 # print (u"真实参数:")
33 print([A, k, theta])
34 # print (u"拟合参数")
35 print(plsq[0]) # 实验数据拟合后的参数
36
37 pl.plot(x, y0, label=u"real data")
38 pl.plot(x, y1, label=u"data with noisy")
39 pl.plot(x, func(x, plsq[0]), label=u"nihe data")
40 pl.legend()
41 pl.show()
运行提示错误,缺少第三方包,如six,dateutil,pyparsing等,缺什么装什么;第三方包大多直接拖到D:\python34\lib目录下就可以了,很方便。
都装上后,运行成功,如图;
http://www.cnblogs.com/javajava/p/4792791.html
使用Python做科学计算初探(转)的更多相关文章
- 使用Python做科学计算初探
今天在搞定Django框架的blog搭建后,尝试一下python的科学计算能力. python的科学计算有三剑客:numpy,scipy,matplotlib. numpy负责数值计算,矩阵操作等: ...
- 使用python做科学计算
这里总结一个guide,主要针对刚开始做数据挖掘和数据分析的同学 说道统计分析工具你一定想到像excel,spss,sas,matlab以及R语言.R语言是这里面比较火的,它的强项是强大的绘图功能以及 ...
- Workshop:用Python做科学计算
Python是程序史上最流行的开源语言之一. 仅在官方包索引PyPi上就已经发布了超过10万个开源软件包,而且还有更多的项目. 在SciPy的麾下,有一个成熟的python包生态系统,可以使用Pyth ...
- windows下如何快速优雅的使用python的科学计算库?
Python是一种强大的编程语言,其提供了很多用于科学计算的模块,常见的包括numpy.scipy.pandas和matplotlib.要利用Python进行科学计算,就需要一一安装所需的模块,而这些 ...
- Python下科学计算包numpy和SciPy的安装
转载自:http://blog.sina.com.cn/s/blog_62dfdc740101aoo6.html Python下大多数工具包的安装都很简单,只需要执行 “python setup.py ...
- Python的科学计算包matplotlib setup
回想起大学四年 专业一直使用matlab,然而我却没在PC上装成功过,以前懒于思考这种数学工具的作用,直到最近,大学同学研究生要毕业了,几经交流,和自己阅读了一些机器学习的教材之后,发觉科学计算包和画 ...
- Python下科学计算包numpy和SciPy的安装【原创】
Python下大多数工具包的安装都很简单,只需要执行 "python setup.py install"命令即可.然而,由于SciPy和numpy这两个科学计算包的依赖关系较多,安 ...
- 2019-04-20 Python之科学计算库学习总结
一.numpy库和matplotlib库的学习 (1)numpy库介绍:科学计算包,支持N维数组运算.处理大型矩阵.成熟的广播函数库.矢量运算.线性代数.傅里叶变换.随机数生成,并可与C++/Fort ...
- python numpy 科学计算通用函数汇总
import numpy as np #一元函数 #绝对值计算 a = -1b = abs(a)print(b)输出: 1 #开平方计算 a = 4b = np.sqrt(a)print(b)输出: ...
随机推荐
- WPF自定义ListBox样式
<!--竖向--> <Style x:Key="ListBoxStyle1" TargetType="{x:Type ListBox}"> ...
- Compass用法指南
Compass用法指南 Sass是一种"CSS预处理器",可以让CSS的开发变得简单和可维护.但是,只有搭配Compass,它才能显出真正的威力. 本文介绍Compass的用法 ...
- NPC
这里的想说的NPC不是Non-Player-Controled,非玩家控制角色,而是Non-determinisitc Polynomial complete problem,它属于一类很特殊的问题, ...
- poj 2010 Moo University - Financial Aid (贪心+线段树)
转载请注明出处,谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents by---cxlove 骗一下访问量.... 题意大概是:从c个中选出n个 ...
- Windows下Putty连接虚拟机Ubuntu
本文的题目是Windows下使用Putty连接虚拟机中的Ubuntu. 事实上针对这种一个问题,已经有非常多的文章.blog能够參考和学习.可是在本人的学习过程中还是遇到可非常多的问题. 特写下自己的 ...
- lightoj1030(期望dp)
有n个格子,初始的时候pos=1,然后丢骰子,然后新的pos为pos+骰子的点数,走到新的pos,可以捡走该pos上的黄金. 特殊的是,如果新的pos超过了n,那么是不会走的,要重新丢骰子. 所以要分 ...
- [WPF]Binding的Converter和Validator
不拘一格用数据的Converter 上篇文已经说明,Binding就是数据源与目标之间的“关联”.大多数情况下,数据从Source到Target以及从Target返回Source都是“直来直去”的,但 ...
- 矩阵快速幂---BestCoder Round#8 1002
当要求递推数列的第n项且n很大时,怎么快速求得第n项呢?可以用矩阵快速幂来加速计算.我们可以用矩阵来表示数列递推公式比如fibonacci数列 可以表示为 [f(n) f(n-1)] = [f(n ...
- 漫游Kafka介绍章节简介
原文地址:http://blog.csdn.net/honglei915/article/details/37564521 介绍 Kafka是一个分布式的.可分区的.可复制的消息系统.它提供了普通消息 ...
- c# Use Properties Instead of Accessible Data Members
advantage of properties: 1 properties can be used in data binding, public data member can not. 2 dat ...