python_numpy
文章平均质量分 74
潜水的飞鱼baby
just maybe。
展开
-
python_numpy的矩阵运算及对应的matlab写法
背景:NumPy和Matlab不一样,对于多维数组的运算,缺省情况下并不使用矩阵运算,可以调用相应的函数对数组进行矩阵运算。或者使用numpy库提供了的matrix类,用matrix类创建的是矩阵对象,它们的加减乘除运算缺省采用矩阵方式计算,用法和matlab十分类似。不过一般用户很容易将NumPy中同时存在的ndarray和matrix对象弄混,一般不建议在大程序中使用。下面简单介绍pyth原创 2017-09-24 08:46:08 · 1879 阅读 · 0 评论 -
python用matplotlib画图例子(3)
背景1每50毫秒产生一个标准正态随机数与之前的数进行累加。代码1import numpy as np import matplotlib.pyplot as plt import matplotlib.animation as animation fig = plt.figure() ax = fig.add_subplot(111)#初始漫步数N = 30#最大漫步数MaxN = 50原创 2017-10-14 16:21:15 · 438 阅读 · 0 评论 -
python使用pip安装SciPy、SymPy、matplotlib
背景:使用pip install SciPy的格式安装python函数库SciPy的时候,发现老是报错,从网上找信息也没找到合适的解决办法,最后使用whl格式文件安装成功。过程:本人事先电脑安装了python3.5.2,电脑是64位的。下载函数库:从网址http://www.lfd.uci.edu/~gohlke/pythonlibs/下载对应函数库的whl文件。注意下载对应的原创 2017-10-23 16:28:59 · 13959 阅读 · 0 评论 -
python汇总数据:日线->周线
背景将浦发银行一年日线的收盘价汇总成周线数据。代码"""数据按周汇总"""import numpy as npimport datetime"""step1 :读取数据 2017浦发银行.CSV: 时间,收盘价(浦发银行),成交量(浦发银行) 2017-01-03,12.3804137039431,21377729.6808511 2017-01-04,12.4原创 2017-10-15 12:24:54 · 4569 阅读 · 2 评论 -
python绘制布林带
背景:有浦发银行2017年的每日收盘价(数据来源:天软科技),画布林带。代码:"""布林带"""import numpy as npimport datetimeimport matplotlib.pyplot as plt"""step1 :读取数据 2017浦发银行2.CSV: 时间,收盘价(浦发银行),最高价(浦发银行),最低价(浦发银行) 2017-01-0原创 2017-10-15 18:15:44 · 7455 阅读 · 3 评论 -
python统计函数库scipy.stats的用法1/3
背景总结统计工作中几个常用用法在python统计函数库scipy.stats的使用范例。正态分布的几个范例生成服从指定分布的随机数 norm.rvs通过loc和scale参数可以指定随机变量的偏移和缩放参数,这里对应的是正态分布的期望和标准差。size得到随机数数组的形状参数。(也可以使用np.random.normal(loc=0.0, scale=1.0, size=None)) In原创 2017-10-16 03:37:09 · 65985 阅读 · 2 评论 -
sublime text 的使用——用sublime text 3搭建python 的ide
前言想在sublime text3中实现解释python时打开命令行窗口,类似idle效果,找了网上的很多资料,良莠不齐,简单记录下自己做个备份。步骤工具->编译系统->新编译系统把下面的代码粘贴上去保存。工具->立即编译或者直接ctrl+B即可看到效果。其实安装 Python,安装时选择添加路径到系统中,或者稍后自己添加也可随便写个 demo,保存为.py 文件,Ctrl + B 就也可以原创 2017-11-02 00:03:22 · 361 阅读 · 0 评论 -
python机器学习
准备数据import numpy as npimport urllib.request""" 下载数据 dataset前2行数据 >>> dataset[:2] array([[ 6. , 148. , 72. , 35. , 0. , 33.6 , 0.627, 50. , 1.原创 2017-10-25 23:47:27 · 848 阅读 · 0 评论 -
将python安装信息加入注册表
将python安装信息加入注册表背景重装系统,发现之前装在E盘的python可以直接使用,就只是将python的安装目录加入到环境变量中,也一直没有管它,今天跟天软交互的时候发现一直不成功,猜测可能是没有注册表信息。 从网上找的一段代码,直接复制运行即可,留存备查。代码py3.5.2import sysfrom winreg import *# twe...原创 2018-03-20 13:19:26 · 4819 阅读 · 1 评论 -
启动jupyter&ipython时,报错“Fatal error in launcher: Unable to create process using '"'”
背景cmd启动jupyter&ipython时,报错Fatal error in launcher: Unable to create process using '"'”过程网上找的资料,所言大多是安装了多版本python时使用pip出现的问题,解决方法大致类似于: py -3 -m pip install XXXX VS py -2 -m pip insta...原创 2018-03-22 04:49:32 · 5178 阅读 · 0 评论 -
以面向过程的方式给出《贝叶斯思维:统计建模的Python学习法》——第五章:决策分析代码
背景在《艾伯特贝叶斯思维:统计建模的Python学习法.pdf》一书中,代码以面向对象给出,然而为了一本介绍统计思想的书籍,只需要了解其计算逻辑即可(PS:其实是懒得去看书中累的组织过程),因此给出阅读这本书的时候的代码及一些想法。导入常见模块# %load "E:\桌面space\临时数据\python\个人自定义模块\ImportFile.py"# Standard Scientifi...原创 2018-11-18 17:21:34 · 864 阅读 · 0 评论 -
python用matplotlib画图例子(2)
背景:用python的matplotlib的api画沪深300在2017年的走势图,并且标识出全年最低价及国庆节。(数据来源:天软科技)代码import numpy as npimport matplotlib.pyplot as pltimport pandas as pdfrom datetime import datetime """数据格式(2017沪深300.CSV):date原创 2017-10-14 12:31:42 · 499 阅读 · 0 评论 -
python用matplotlib画图例子(1)
背景用matplotlib的api画一个随机漫步,然后设置:标题,轴标签,刻度以及刻度标签,图例。代码import numpy as npimport matplotlib.pyplot as plt"""随机漫步"""x = np.random.randn(1000).cumsum()fig = plt.figure()ax1 = fig.add_subplot(2,1,1)ax2 =原创 2017-10-14 10:30:38 · 466 阅读 · 0 评论 -
python进行JB正态性检验
Jarque-Bera检验:检验序列是否符合正态分布一种总体分布的正态性检验。当序列服从正态分布时,JB统计量: JB=n(S26+(K−3)224)JB=n(\frac{S^2}{6}+\frac{(K-3)^2}{24}) 渐进服从χ2(2)\chi ^2(2)分布。其中nn为样本规模,S、KS、K分别为随机变量的偏度和峰度。计算公式如下: M2=∑i(xi−x¯)2nM2=\frac{\原创 2017-10-22 10:52:42 · 16478 阅读 · 5 评论 -
python拟合二元一次函数
背景:使用scipy拟合一元二次函数。参考:HYRY Studio-《用Python做科学计算》代码:import numpy as npfrom scipy.optimize import leastsqimport pylab as pldef func(x,p):"""数组拟合函数"""A,k,theta = preturn A*(x-k)**2+th原创 2017-09-24 15:08:53 · 7872 阅读 · 1 评论 -
python用fmin, fmin_powell, fmin_cg, fmin_bfgs求解无约束极小值点
背景:作求解无约束极小值点范例。代码:import scipy.optimize as optimport numpy as npdef test_fmin(fminfunc,x0,a): """ x0为优化算法的初始值,各种优化算法必须 a为目标函数的参数 """ def targetfunc(x): """原创 2017-09-24 15:47:48 · 6804 阅读 · 0 评论 -
python用fsolve、leastsq对非线性方程组进行求解
背景:实现用python的optimize库的fsolve对非线性方程组进行求解。可以看到这一个问题实际上还是一个优化问题,也可以用之前拟合函数的leastsq求解。下面用这两个方法进行对比:代码:from scipy.optimize import fsolve,leastsqfrom math import sin,cosdef f(x): x0 = float(原创 2017-09-24 17:11:05 · 16410 阅读 · 1 评论 -
python一维插值以及subplot画图范例
背景:python实现插值以及subpot画子图范例。代码:import numpy as npimport matplotlib.pyplot as plfrom scipy import interpolatelb,ub = 0,20#10个样本点x = np.linspace(lb,ub,10)y = 100*np.sqrt(np.abs(np.sin(x))原创 2017-09-24 19:00:24 · 2471 阅读 · 0 评论 -
python用quad、dblquad实现一维二维积分范例
背景:python函数库scipy的quad、dblquad实现一维二维积分的范例。需要注意dblquad的积分顺序问题。代码:import numpy as npfrom scipy import integratedef half_circle(x): """ 原心:(1,0),半径为1 半圆函数:(x-1)^2+y^2 = 1 """原创 2017-09-25 13:43:24 · 9570 阅读 · 0 评论 -
用python实现解常微分方程组的简单示例以及用odeint解常微分方程的范例
背景:包括两个部分,一个是演示怎么自己写代码解常微分方程,另一部分就是示范python怎么调用解常微分方程的函数。下面的方程组给出了三维空间中各个坐标点上的速度,需要求解的洛仑兹引子的运动轨迹。软件:python3.5.2部分1:(div)代码:# -*- coding: utf8 -*-import numpy as np"""移动方程:t时刻原创 2017-09-28 00:43:41 · 30870 阅读 · 2 评论 -
python的meshgrid的用法及灰度图像的显示
背景:meshgrid可以代替数组的双重遍历,下面展示meshgrid进行画图,主要是将三维平面z=sin((x^2-y^2)^2+x*y)投影到xy平面,用z值作为各点的值。软件:python3.5代码:import numpy as npimport matplotlib.pyplot as pltpoint = np.arange(-5,5,.1)x,y = np原创 2017-09-28 13:28:19 · 2118 阅读 · 0 评论 -
用python画AR模型时序图
背景:用python画AR模型的时序图。结果:代码:import numpy as npimport matplotlib.pyplot as plt"""AR(1)的时序图:x[t]=a*x[t-1]+e"""num = 2000e = np.random.rand(num)x = np.empty(num)"""平稳AR(1)"""a = -原创 2017-09-29 21:20:36 · 4468 阅读 · 0 评论 -
C读取python结构数组
#include #include /* run this program using the console pauser or add your own getch, system("pause") or input loop */struct numpy{ char name[30]; int age; float weight;};int main(int argc, c原创 2017-09-22 07:29:53 · 1243 阅读 · 0 评论 -
《集体智慧编中文版》读书笔记——第二章:协作型过滤
前言一本比较老的书,偏实用,正好打算入手python机器学习,实现一遍代码。关于协作型过滤介绍他的一个引用场景:已有用户对自己感兴趣的电影评分的情况下,怎么合理的为用户推荐新电影?或者一个新用户在浏览电影时,推荐其他相关的电影? 协作性过滤 的思想处理上述问题的逻辑是这样的:两个用户对共同感兴趣的电影评分都相近,显然可以互相推荐对方未看过的电影。如果两部电影用户对它们的评分很接近,显然这两部电原创 2017-10-29 06:33:00 · 341 阅读 · 0 评论 -
面向过程给出《贝叶斯思维:统计建模的Python学习法》——二维彩球问题学习代码
背景给出读《艾伯特贝叶斯思维:统计建模的Python学习法.pdf》的时候,写的代码,以面向过程的方式给出。本章彩弹问题,求似然度的时候,假设已知隐藏点时,射手等概率从各个角度射击。代码导入常见模块# %load "E:\桌面space\临时数据\python\个人自定义模块\ImportFile.py"# Standard Scientific Importimport numpy...原创 2018-11-21 21:44:06 · 1129 阅读 · 0 评论