- 博客(47)
- 资源 (6)
- 收藏
- 关注
RSS订阅原创 使用Python实现累积面积图 using Python for cumulative area plot
【代码】使用Python实现累积面积图 using Python for cumulative area plot。
2023-08-19 10:27:36
281
原创 自回归模型的建模与参数估计-Python
AR§模型的参数估计设{XtX_tXt}适合Xt=a1Xt−1+⋯+apXt−p+ϵtX_t=a_1X_{t-1}+\cdots+a_pX_{t-p}+\epsilon_tXt=a1Xt−1+⋯+apXt−p+ϵt,式中{ϵt\epsilon_tϵt}为独立同分布白噪声序列；Eϵt2=σ2E\epsilon^2_t=\sigma^2Eϵt2=σ2，来自{XtX_tXt}的样本为X1，X2，⋯，XnX_1，X_2，\cdots，X_nX1，X2，⋯，Xn，要估计a1，a2，⋯，ap
2021-05-23 20:40:07
2859
原创 MIT公开课-环境工程与控制-优化
Computing and Data Analysis for Environmental ApplicationsWater Resource SystemsEnvironmental Engineering Applications of Geographic Information SystemsEnvironmental Engineering Masters of Engineering ProjectSystems Optimization: Models and Computatio.
2021-05-23 10:03:07
264
原创 【搬运课程】时间序列建模-第一周
DS-GA 3001.001 课程时间序列建模几个概念及作业python实现时间序列 {Xt}\{X_t\}{Xt}均值 μX(t)=E(Xt)\mu_X(t)=\mathbb{E}(X_t)μX(t)=E(Xt)协方差 RX(t,u)=cov(Xt,Xu)R_X(t,u)=\text{cov}(X_t,X_u)RX(t,u)=cov(Xt,Xu)自相关函数（Autocorrelation Function） ρX(t,u)=RX(t,u)RX(t,t)RX(u,u)\rho_X
2021-05-22 18:57:30
427
原创 如何Python或者R语言中实现用ARIMA模型输出参数实施预测
针对文章如何MATLAB实现用ARIMA模型输出参数实施预测的补充解释，因为使用Python和R语言进行AR等时间序列模型估计的时候，尤其是AR模型，由于输出的参数结果和MATLAB的不对应该，可能会产生一定的疑惑，因此在这里说明。同样为上篇文章的数据集为例针对AR（2）模型进行建模和参数估计，并输出参数。from pandas import read_csvfrom statsmodels.tsa.arima.model import ARIMAseries = read_csv('daily-mi
2021-05-06 10:10:52
2412
3
原创 如何MATLAB实现用ARIMA模型输出参数实施预测
正文自回归（AR）模型、移动平均（MA）模型、自回归移动平均（ARMA）和自回归差分移动平均（ARIMA）模型是时间序列模型，它们主要是使用历史时间步的观测值作为回归方程的输入，以预测下一时间步的值。这是一个非常简单的想法，可以导致对一系列时间序列问题的准确预测。在本教程中，您将了解如何使用MATLAB实现时间序列预测模型。完成本教程后，您将了解：如何部署一个时间序列模型并进行预测。如何获取已经估计的时间序列模型的参数实施直接进行预测。结合公式更加深入了解自回归移动平均等时间序列模型。自回
2021-05-05 17:02:58
35987
17
原创 纳什效率系数与可决系数的差异
纳什效率系数与可决系数的差异正文参考文献正文纳什效率系数和可决系数的公式形式非常相似，主要差异如下：可决系数（Coefficient of determination，R）是用来度量一个统计模型的拟合优度的。R2=1−∑(yi−y^i)2∑(yi−y‾)2R^2 =1-\frac{\sum{(y_i-\hat{y}_i)^2}}{\sum{(y_i-\overline{y})^2}} R2=1−∑(yi−y)2∑(yi−y^i)2式中：yiy_iyi是变量观测值；y‾\over
2021-05-03 15:36:02
10700
5
原创 慎重用EXCEl进行公式拟合！！！
excel曲线拟合的公式系数通常保留位数，对于数量级相差较大的变量进行公式拟合的时候，最好选用其他专业的程序或者软件进行，不然可能会导致非常大的误差，进而使得后续的计算错误。
2021-04-10 15:31:22
1712
原创 在Vscode中中止正在运行的python程序+快捷键Ctrl+Alt+M
在Vscode中中止正在运行的python程序+快捷键Ctrl+Alt+M
2021-04-04 20:26:46
8336
1
原创 windows+python nc、grib数据处理
依赖的python模块和windows软件netCDF4、wgrib2nc数据的读取示例```python# 导入模块from netCDF4 import Datasetimport itertoolsimport pandas as pdfrom datetime import datetime from datetime import timedelta```python# 获取数据集中所有变量名filename = r'air.mon.mean.nc'data = D
2021-03-28 16:13:08
1923
2
原创 Matplotlib绘图-科学计数法-定制化设置坐标轴数据刻度格式
Matplotlib绘图，自定义坐标轴数据刻度格式：import matplotlib.pyplot as plt import pandas as pd import numpy as npfrom matplotlib import rcParamsfrom matplotlib.ticker import FuncFormatterconfig = { "font.family":'serif', "mathtext.fontset":'dejavuserif',
2021-01-26 14:17:18
4284
4
原创 开源天气预报产品 (持续更新...) Open-source weather forecast products
1. 基于气象组织的气象预报范围定义： http://www.wmo.int/pages/prog/www/DPS/GDPS-Supplement5-AppI-4.html.2. 预报对于水库等水利设施的调度和决策具有重要意义，不同尺度的预报服务于不同的调度目标，例如，防洪等短期目标需要1-3天的预报，供水、发电等长期目标需要季节预报等。(1) The Observing System Research and Predictability Experiment (THORPEX) Interacti.
2021-01-23 11:09:10
627
原创 Matplotlib绘图的字体问题
记录一下，方便自己绘图的查阅，首先自己常用的自定义配置可以写到.mplstyle文件中，直接用plt.style.use('my.mplstyle')设置中英文混用+非斜体文本在中文论文中，可能出现中英文混用的标签等，参考Matplotlib 中英文及公式字体设置使用Matplotlib内置tex进行配置即可，（英文字体有限通常是times new roman受限于\mathrm命令，宋体+Times new roman（非斜体的文本或者单位）+公式）代码如下：import matplotli
2021-01-01 20:34:38
2082
原创 word制位符控制公式和编号
1. 查看页边距：2. 设置制位符a4纸宽度21cm,(21-1.27*2)/2=9.23cm,21-1.27*2=18.46cm 全部选择居中对齐，然后通过tab键盘控制位置即可。
2020-11-23 10:34:29
1161
原创 FreePic2Pdf快速为pdf文档添加带有超链接的书签目录
FreePic2Pdf快速为pdf文档添加带有超链接的书签目录下载FreePic2Pdf提取书签点击 【更改PDF】点击 【从PDF取书签】点击 【1、需要操作的PDF文件】【2、存放接口文件的文件夹】【3、开始】找到生成的书签文件添加书签Question下载FreePic2Pdf下载地址(https://pan.baidu.com/s/1kVHzVmf)密码:at9eref: https:...
2020-05-06 13:53:40
4165
原创 netCDF文件的scale_factor和add_offset
我使用python的netCDF4读取数据，发现数据集存在scale_factor和add_offset，但是我读取的数据应该是Unpacking data，也就是转换后的实际数据，不需要再处理，因为数据压缩是通过偏移和缩放之后将浮点数转化为整数，为了验证正确性，我通过下列代码验证：from netCDF4 import Datasetfilename = r'geopotential....
2020-04-08 15:47:07
9325
8
原创 如何获取ERA-Interim Reanalysis日降雨量数据
ERA-Interim数据的API下载请参考：https://blog.csdn.net/weixin_42924891/article/details/85225062import linecachefrom ecmwfapi import ECMWFDataServera = [31,28,31,30,31,30,31,31,30,31,30,31]for i in range(2...
2020-03-01 14:34:25
2951
2
原创 seaborn 热力图边缘只有一半？？？
matplotlib3.1.1版本的bug解决方法：1. 重新按照其他版本，更新或者回退2.ax = sns.heatmap(...);bottom, top = ax.get_ylim()ax.set_ylim(bottom + 0.5, top - 0.5)ref:https://stackoverflow.com/questions/56942670/matplot...
2020-02-15 23:46:02
1728
2
原创 机器学习得分函数
关于sklearn中GridSearch等中中scoring参数，有两个问题，一个问题关于一些不可能为负数的指标输出的误差有负数，说明其中做了一些处理，至于怎么处理的还没看懂，留着有空了看一下。还有一个问题是格网调参的时候的cv参数。问题1：非负指标有负数，做了什么处理为什么这么处理。参考链接：https://stackoverflow.com/questions/21050110/skl...
2020-02-08 23:40:03
2075
原创 应用机器学习方法进行多步预测的四种策略
//TODO 有空了我找个例子使用python和MLP实现对比一下四种策略参考：[1] 4 Strategies for Multi-Step Time Series Forecasting, 2017[2]Machine Learning Strategies for Time Series Forecasting, 2013[3]Recursive and direct m...
2020-02-05 16:07:15
1406
原创 数据序列相关性-ACF，PACF和CCF
引言最近写论文关于预测的特征选择遇到一些问题，想把自己查询学习到的东西整理记录一下，理一理头绪，希望能加深自己对这些东西的理解。首先介绍引入几个概念：自相关函数（autocorrelation function，ACF）、偏自相关函数（partial autocorrelation，PACF）和互相关函数（cross-correlation function，CCF）。接下来介绍每个指标的计...
2020-01-30 17:49:22
20541
6
原创 ECMWF-ERA和ERA5关于累积变量的下载和处理
下载气象数据产品，首先注意到单位，其次看数据是否变换（scale，add_offset，为了存储可能会缩放或者增加偏移等）。这里我为了自己便于查阅，也为了大家能很快看到关于ERA数据的累积量的说明，贴两个网址：1.https://confluence.ecmwf.int/pages/viewpage.action?pageId=566582332.https://confluence....
2019-05-15 08:16:14
6359
7
原创 python界面登录-验证码（三）
真的要好好学一下写作了，等好好的有条理的整理自己做过的工作才能方便的进行下一步的使用，能整理好自己的学习的东西才能提高效率，更加明确的进行下一步的工作提高自己的工作效率！！！下一步就是在以下前提下进行网页的登录和课程的查询了：import osimport timefrom bs4 import BeautifulSoup from selenium import webdriv...
2019-01-14 10:57:02
1541
原创 python爬虫界面登录-验证码（二）
继续刚才的说，刚才的pytesseract的安装目录通常都在python或者anaconda目录下的Lib\site-packages\pytesseract目录下。配置完之后就可以用了，用下面的代码就行验证码的保存和识别。driver.save_screenshot('f://aa.png') # 截取当前网页，该网页有我们需要的验证码imgelement = driver.find_...
2019-01-14 10:39:30
749
原创 python爬虫登录界面-验证码
选课真的是个苦恼的事，每次都一到时间都登录不上，登录上以后都没的选了，偶尔有人退课也不能一直盯着等待选课，想着使用爬虫进行自动监测和选课。通过分析这次操作过程，决定还是利用selenium，主要由两个点之前没有做过，一个是输入表单的账号和密码提交登录，还有一个就是验证码的问题，由于目前选课已经完成，所以程序做了一半，可以登录进去找到要选的课程，有部分bug未修复，也没有完整的完成选课过程，但是用到...
2019-01-14 10:29:39
1601
3
原创 IEDriverServer来操控sendkeys方式时，输入速度非常慢
通过调用64位IEDriverServer来操控sendkeys方式时，输入速度非常慢。网上说是64位有bug。解决办法：使用32位IEDriverServer。亲测有效。IEDriverServer链接：http://selenium-release.storage.googleapis.com/index.html ...
2019-01-11 11:29:47
1543
原创 scikit-learn中的多层感知器调用模型输出数据类型为float出现Unknown label type: 'unknown'
受下面截图的启发，说的意思好像是输出类型必须是绝对的（categorical），train_y后加上astype(‘int’)即可，fit(train_x,train_y.astype(‘int’))，但是不能试astype('float')，用这个还是会报错，具体原因我不清楚，解决方法就是把输出数据乘以一个数例如10000，然后转换成整数，最后记得对预测集合除10000进行转化。参考链接：1.h...
2019-01-06 19:33:33
953
1
原创 利用python的bs4和selenium库结合实现动态页面的爬取（爬取天气网上面的历史天气数据）
报告分析需要历史气象数据，查询到天气网上面有历史天气数据，从2011年到2018年，第一次接触爬虫，在网上找了爬取天气网历史数据的python源码，利用bs4库，但是实际操作中发现soup.select( )函数返回的列表总是[ ] （空），查询发现天气网目前使用的是javascript写的动态页面，可能是为了防止爬虫接入。经过查资料又找到了一种爬取动态页面的方法，经尝试方法有效，写个博客记录一下...
2018-12-17 13:01:37
3570
原创 PLSQL导入dmp文件闪退
1.用户没有dba权限，赋予登录用户dba权限即可；例如登录用户名为system，打开控制台，grant dba to system2.导入的dmp文件和oracle数据库版本号不对，用notepad++打开dmp文件，将第一行的版本号改为oracle对应的版本号，查看oracle版本号的方法，在控制台输入 sqlplus /nolog 回车...
2018-12-15 08:17:09
7123
原创 python实现门限回归
门限回归模型(Threshold Regressive Model，简称TR模型或TRM)的基本思想是通过门限变量的控制作用，当给出预报因子资料后，首先根据门限变量的门限阈值的判别控制作用，以决定不同情况下使用不同的预报方程，从而试图解释各种类似于跳跃和突变的现象。其实质上是把预报问题按状态空间的取值进行分类，用分段的线性回归模式来描述总体非线性预报问题。多元门限回归的建模步骤就是确实门限变量、...
2018-10-28 15:16:13
5496
2
原创 python实现逐步回归
逐步回归的基本思想是将变量逐个引入模型，每引入一个解释变量后都要进行F检验，并对已经选入的解释变量逐个进行t检验，当原来引入的解释变量由于后面解释变量的引入变得不再显著时，则将其删除。以确保每次引入新的变量之前回归方程中只包含显著性变量。这是一个反复的过程，直到既没有显著的解释变量选入回归方程，也没有不显著的解释变量从回归方程中剔除为止。以保证最后所得到的解释变量集是最优的。本例的逐步回归则...
2018-10-25 08:49:09
22555
37
原创 python实现多元线性回归
用python实现多元线性回归，先不使用sklearn.linear_model模块，编写如下：#导模块import numpy as npimport pandas as pd#变量初始化X=[]Y=[]B=[]Q_e=0Q_E=0#从csv文件中读取数据def get_data(file_name): data=pd.read_csv(file_name,heade...
2018-10-16 19:54:27
7337
原创 python实现一元线性回归
最近在看中长期水文预报，打算使用python语言实现课本的模型并进行实例的计算结果的检验，为了监督自己和整理记录自己实现的模型代码，打算写博客记录自己的代码实现和部分思路。首先，自己不使用现成的模块实现的是一元线性回归模型，然后和模块实现对比，学习模块的调用。接下来是直接编写的程序文件，并计算了几个检验系数，没有使用sklearn模块：#导模块import numpy as npimpo...
2018-10-16 15:22:14
2832
原创 java贪吃蛇小游戏(三)
接下来是我为了便于修改游戏数据提取的图片类(GetImage())，组件位置尺寸的设置类(LocationSize())、窗口操作类(UiUtil),以及方向的枚举类package util;import java.awt.Image;import java.awt.Toolkit;import java.net.URL;/** * 此类用于放置各个显示对象的背景绘制的素材，便于...
2018-08-07 19:28:03
741
原创 java贪吃蛇小游戏（二）
接着上文粘贴代码，接下来是一个最核心和复杂的类BackgroundPanelpackage view;/* * 自定义重绘图片，重写JPanel面板的paint(Graphics g)类 */import java.awt.Color;import java.awt.Font;import java.awt.Graphics;import java.awt.Label;i...
2018-08-07 19:16:32
556
原创 java贪吃蛇小游戏(一)
刚看完java基础课的课程花费了一个星期边看Swing边复习学过的知识边完善程序，还有不少bug，例如，点击帮助后按空格键会开始游戏，贪吃蛇撞到墙后以后再重新开始游戏有一节小蛇没有及时重绘制；先分享出来，写一个详细一点的注释，有时间再慢慢完善。这个小程序中，我用到的组件有JFrame(窗体)、Jpanel(面板)、JButton(按钮)、Label(标签)，实现了基本的贪吃蛇游戏功能（吃食物，计...
2018-08-07 17:55:30
1207
原创 JAVA多线程之线程控制
返回值 方法名 说明 public static void sleep(long millis) 线程休眠，在指定的毫秒数内让当前正在执行的线程休眠（暂停执行） public static void join() 线程加入，等待调用该方法的线程终止其他线程才能执行 public static void yield() 线程...
2018-08-06 13:08:40
183
原创 多线程之线程同步问题
CPU的一次操作必须是原子性的，想象一下多个窗口卖同一时段同一个电影的电影票的例子（共享数据源），当一个线程再买票，计算剩余票数和显示剩余票数的时候，其他窗口不能进行这些操作，如果多个窗口同时操作就会出现问题，例如窗口2卖完票剩下100张，窗口2先卖票，然后计算计算剩余票数为100张，在显示剩余票数前，窗口1也卖出去了一张票，此时实际剩下99张，窗口2显示的时候还是刚才计算的100张。如何解决线程...
2018-08-06 11:25:48
237
原创 JAVA多线程的概述及创建
一、线程和进程进程是指正在运行的程序，是系统进行资源分配和调用的独立单位。每一个进程都有它自己的内存空间和系统资源。线程是进程中的单个顺序控制流，是一条执行路径， 一个进程如果只有一条执行路径，则称为单线程程序。一个进程如果有多条执行路径，则称为多线程程序。一个java程序想要运行，首先由java命令启动JVM，JVM启动就相当于启动了一个进程。接着有该进程创建了一个主线程去调用main方法。...
2018-08-06 11:08:17
140
time series model for Mannual predict.zip
2021-05-05
windows+python nc、grib数据处理.zip
2021-03-28
python逐步回归的测试数据
2020-01-30
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人