2017年09月_潜水的飞鱼baby

原创 AR模型

AR模型：具有如下结构的模型称为p阶自回归模型，简记为AR(p): AR(p)模型有三个限制条件：。保证模型的最高阶数为p。随机干扰项序列为零均值白噪声序列。当期的随机干扰项与过去的序列值无关，即：中心化AR(p)模型：当a0=0时，自回归模型称为中心化AR(p)模型。非中心化AR(p)模型可以通过下面的变化转换为中心化AR(p)序列。 P阶自回归系数

2017-09-30 07:32:09 62152 6

原创用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 4467

原创 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

原创用python实现解常微分方程组的简单示例以及用odeint解常微分方程的范例

背景：包括两个部分，一个是演示怎么自己写代码解常微分方程，另一部分就是示范python怎么调用解常微分方程的函数。下面的方程组给出了三维空间中各个坐标点上的速度，需要求解的洛仑兹引子的运动轨迹。软件：python3.5.2部分1：(div)代码：# -*- coding: utf8 -*-import numpy as np"""移动方程：t时刻

2017-09-28 00:43:41 30867 2

原创 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 9569

原创 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 2470

原创 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 16409 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 6803

原创 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_numpy的矩阵运算及对应的matlab写法

背景：NumPy和Matlab不一样，对于多维数组的运算，缺省情况下并不使用矩阵运算，可以调用相应的函数对数组进行矩阵运算。或者使用numpy库提供了的matrix类，用matrix类创建的是矩阵对象，它们的加减乘除运算缺省采用矩阵方式计算，用法和matlab十分类似。不过一般用户很容易将NumPy中同时存在的ndarray和matrix对象弄混，一般不建议在大程序中使用。下面简单介绍pyth

2017-09-24 08:46:08 1879

原创 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

原创 C++不使用任何 include文件得到+inf,-inf和nan

背景：最近在看网易的公开课《斯坦福大学公开课：编程范式》，才知道大学的时候错过了一个多么好玩的课程。不过重点不在这里，而是终于能理解了之前一直没什么心情看的一些东西。代码部分用了《编程范式》提到的一些操作指针的技巧，当然单纯实现标题的要求完全可以不这样做。相关知识：现在所有的计算机都支持一个被称为IEEE浮点标准。IEEE浮点数表示法： IEEE浮点标

2017-09-21 13:47:57 4757

原创 python批量下载网页图片及列表

背景：放假在家里没事鼓捣了一个用python实现的简单网页信息抓取的程序。demo功能：将千图网的商业海报http://www.58pic.com/topic/419-1.html的（海报名字，海报详情页面地址，海报缩略图地址）通过正则匹配提取出来，然后写入到filed='python_load_demo.txt'中，然后把第一个海报缩略图下载到文件夹imgpath = "pytho

2017-09-17 16:54:44 3285

潜水的飞鱼baby