- 博客(35)
- 资源 (4)
- 收藏
- 关注
原创 cuda加速求解龙格库塔四阶五步积分
这个方程组可以用来描述原始的二阶微分方程的动态。一阶微分方程组更容易用数值方法求解,因为大多数数值。在实际应用中,这个方程组可以用标准的数值方法(如欧拉法、龙格-库塔法等)进行求解。我们可以将这个二阶微分方程转化为一个一阶微分方程组,然后使用RK45方法来求解。第二个一阶微分方程来自于原始方程对 ( y ) 的二阶导数的替换,我们将 (要求解常微分方程组(ODEs)。
2023-12-25 12:50:12
994
原创 海塞矩阵(Hessian matrix)全解
通过计算图和链式法则,自动微分可以高效地计算出函数 f 对于输入变量 x 和 y 的导数,同理基于f’(x)和f’(y),求二阶导数,而无需手动推导导数表达式。因此,在实际应用中需要综合考虑问题的特点来选择合适的优化算法和方法。在这个示例中,我们使用Autograd库来计算函数 ( f(x, y) = x^2 + 2xy + y^2 ) 的一阶导数和二阶导数,并将其组合成海塞矩阵。其中,(*) 表示乘法操作,(+) 表示加法操作,(x) 和 (y) 分别表示输入的变量 x 和 y,(f) 表示函数的输出。
2023-12-20 19:47:55
8002
原创 conda命令使用conda-pack库打包环境并迁移-window系统
在新的电脑上,无法通过网络连接等方法构建所需的python环境。因此需要离线打包本地环境到新电脑上。
2023-09-08 16:59:39
3173
原创 gitlab入门教程 上传本地文件夹到远程端
当6完成以后,发现有一部分文件不想上传到远程端。因此在clone的路径中选择.gitignore文件。然后对自己不想上传的文件夹给出相对路径,然后加上/*的后缀。
2023-08-09 10:48:19
1110
原创 【python】贝叶斯优化包bayesian-optimization
为此,该方法使用代理优化问题(找到采集函数的最大值),尽管仍然是一个难题,但更便宜(在计算意义上)并且可以使用通用工具。因此,贝叶斯优化最适合对要优化的函数进行采样是一项非常昂贵的工作的情况。在每个步骤中,高斯过程都会拟合已知样本(之前探索的点),后验分布与探索策略(例如 UCB(置信上限)或 EI(预期改进))相结合,用于确定下一个应该探索的点。随着观测数量的增加,后验分布得到改善,算法变得更加确定参数空间中的哪些区域值得探索,哪些区域不值得探索,如下图所示。但是,应该清楚的是您不需要这样做。
2023-07-26 16:13:48
4319
原创 scipy的微分方程组求解包solve_ivp使用详解
使用需要检测的函数。在ode求解过程中,每一步都会执行events对应的函数,函数返回的是需要进行判断该变量情况的数值:比如官方案例,加农炮发射,存在两个变量,高度和速度。显式 Runge-Kutta 方法(“RK23”、“RK45”、“DOP853”)应用于非刚性问题,隐式方法(“Radau”、“BDF”)用于刚性问题 [9 ]。如果它进行异常多次迭代、发散或失败,则您的问题可能很僵化,您应该使用“Radau”或“BDF”。“BDF”:基于导数近似的后向微分公式的隐式多步变阶(1 到 5)方法[5]。
2023-07-11 14:30:29
6838
原创 DCCA互相关系数 实例 python
数据来源:PhysioBank ATM下面资源还在审核的话,就用上面的链接下载MIT- BIH 正常窦性心律数据 nsrdb-16265针对数据ECG1和ECG2进行DCCA分析。绿色为ECG1和ECG2基于DCCA的线性拟合【需要去掉注释才能显示】红色为ECG1基于DFA的线性拟合蓝色为ECG2基于DFA的线性拟合需要注意如下:1. 数据以对数曲线表示。2. 这里时间长time=3 间隔s=0.03 数据长度只取到3000 为了避免错误,请确..
2022-03-24 22:12:22
2279
1
原创 DFA去趋势波动分析 实例 python
数据来源:PhysioBank ATM下面资源还在审核的话,就用上面的链接下载MIT- BIH 正常窦性心律数据 nsrdb-16265针对数据ECG1和ECG2进行DFA分析。需要注意如下:1. 数据以对数曲线表示。2. 这里时间长time=3 间隔s=0.03 数据长度只取到3000 为了避免错误,请确保【数据长度>=10^time】3. 可以基于FDA继续算ECG1和ECG2的DCCA4. 拟合是线性的,可根据自己数据情况考虑更换二次及以上。具体方法.
2022-03-24 21:19:01
4967
原创 DCCA互相关系数 理论
Cross-correlation analysis of stock markets using EMD and EEMD翻译源2.4小节DCCA方法是由Podobnik和Stanley[31]最近提出的,用于分析同时记录的非平稳时间序列之间的幂律互相关系,是互相关系研究的有力候选工具。因此,在以往的许多研究中,DCCA方法被用来研究相互关系[33-40]。DCCA过程可以简单描述如下:考虑这两个时序:和,i=1,2,...,N,并且两时序长度相等为N。
2022-03-20 22:00:14
4370
2
原创 【专用链接】文件补档
附加补全:基于一个微分对策问题的机器学习.pdf-机器学习文档类资源-优快云下载test12就是第一问和第二问文件以下是test3_904_19.pyimport gymfrom gym import spacesfrom gym.utils import seedingimport numpy as npfrom math import cos,sin,pi,sqrt,asinimport randomfrom os import path...
2021-12-29 23:39:06
1126
原创 Pymc3采样慢解决方法和原因
解决方法:1. 一般情况:先验设定存在问题(不合适),建议按照同一问题背景下,用别人找到的推荐先验设置。或者查阅相关论文,使得变量先验符合模型背景。2. 在使用别人的推荐先验还是存在极慢的采样,可以检查数据集,数据的分布情况是否和参考教程一致,比如别人是正态,你的是双峰或者偏正态,可以先进行转化调整。3. 更换初始方法,例如find_MAP。4. 上述方法只能解决很慢的情况下,一两天都不能采样出结果。【我的数据量15k,18个参数,chain=2,跑了6h,供参考】如果还有什么好的解决方法可
2021-12-24 19:24:32
1175
2
原创 2024/4汇总
课程设计汇总数据科学导引课程报告https://blog.youkuaiyun.com/qq_54394719/article/details/118699705python课程算法设计并优化https://blog.youkuaiyun.com/qq_54394719/article/details/118700099数据分析课设(SPSS,EVIEWS,R)【理论】https://blog.youkuaiyun.com/qq_54394719/article/details/118700338数据分析课设(R)
2021-11-08 20:37:38
3249
原创 Robyn MMM2.0 - Facebook Marketing Science(R)市场营销-广告投放
代码链接:https://github.com/facebookexperimental/Robyn教学使用视频:https://www.youtube.com/watch?v=aIiadcfL4uw目录1 原文翻译部分:1.1 快速开始1.2使用进化算法进行模型选择1.3 商业真实使用情况1.4分步指南1.4.1 库1.4.2为 Nevergrad 创建、安装和使用 conda 环境1.4.3加载数据1.4.4设置国家1.4.5日期变量1....
2021-09-29 14:41:33
3413
1
原创 python调用Google Analytics Sample (BigQuery)使用方法以及常见错误
数据集没有csv格式,而且很大,所以使用接口来查询数据。前置准备:谷歌账号过程:https://cloud.google.com/bigquery/docs/reference/libraries#client-libraries-resources-python主要在这里进行操作在环境变量里配置经常会莫名其妙出现错误google.auth.exceptions.DefaultCredentialsError: Could not automatically determine.
2021-09-12 20:40:00
778
原创 数值计算方法-MATLAB
3.钢水包使用次数多以后,钢包的容积增大,数据如下: x 2 3 4 5 6 7 8 9 y 6.42 8.2 9.58 9.5 9.7 10 9.
2021-07-14 09:34:19
1355
原创 基于影评的LSTM与GRU对比
以imdb的影评为数据集,利用SimpleRNN和LSTM,GRU做好评与差评的分类网络。并且比较他们在时间和正确率的差别,以及是否启用unroll功能对时间和正确率造成的影响。常见的单词数量,输入训练句子固定长度,单词转编码的向量维度。 然后要载入数据集,利用Keras得到划分好的数据集:对于x的训练集和测试集都用padding的方法填补(避免取出的句子存在不够固定长度的情况)打包集合(将x和y合并为一个数据集),分别得到完整的训练集和测试集,然后对训练集和测试集用minibatch进行
2021-07-14 09:21:27
400
原创 课程设计汇总
数据科学导引课程报告https://blog.youkuaiyun.com/qq_54394719/article/details/118699705python课程算法设计并优化https://blog.youkuaiyun.com/qq_54394719/article/details/118700099数据分析课设(SPSS,EVIEWS,R)【理论】https://blog.youkuaiyun.com/qq_54394719/article/details/118700338数据分析课设(R)【实践-时间序列
2021-07-13 15:21:10
110
原创 论文汇总1
运动模糊图像恢复处理-matlab实现https://blog.youkuaiyun.com/qq_54394719/article/details/118697030单剂量给药下仓室药物模型半衰期的计算https://blog.youkuaiyun.com/qq_54394719/article/details/118695852国内游戏营销竞争潜力模型的多元统计研究https://blog.youkuaiyun.com/qq_54394719 /article/details/118699376...
2021-07-13 15:17:43
85
原创 数据分析课设(R)【实践-时间序列】
背景问题:利用ARIMA乘积模型对航空旅客数据集进行分析一年的预测数据集直接导入R的数据集AirPassengers:这是1949-1960年每月旅客总数;数据集来源Box and Jenkins的文献。#先将数据集导入x=AirPassengers#查看数据集:print(x) Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec1949 112 118 132 129 121 135 148 148 136 119...
2021-07-13 15:12:23
1923
2
原创 数据分析课设(SPSS,EVIEWS,R)【理论】
对外汇数据作前向、后向以及逐步回归,并对输出结果作分析和理论解释。x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 省市 y 1.94 4.5 154.45 207.33 246.87 277.64 135.79 30.58 110.67 80.83 51.83 14.09 北 京 2384 0.33 6.4
2021-07-13 15:07:42
2668
原创 python课程算法设计并优化
题目data.txt中是多个点的坐标集合S,每个点的坐标都是一个二维数组,每一行代表一个点的二维坐标,每个二维数组的第一个数字代表横坐标,第二个数字代表纵坐标。现有如下想法,大致内容如下:设置一个正的常数值d=50,对于集合S中任意两个点,只要它们两个之间的欧式距离不超过d,那么这两点就归为一类。请尝试自行编写python代码实现上述想法,并显示代码运行时间和最后S被划分成了多少个类。很显然用DBSCAN算法例如实现:import numpy as npimport ra
2021-07-13 14:55:12
262
原创 数据科学导引课程报告
Q1:公共自行车租借预测自行车分享系统(bikesharingsystems)较传统的自行车租借存在巨大的优势.它避免了注册、租借和归还等繁琐的过程,用户可以在线上完成所有相关操作.同时,自行车的行驶时长和路线等信息也可以完整记录下来.本题提供一份公共自行车租借数据集[36].请根据城市当天的天气、温度、湿度和风速等环境信息对自行车的租借数量进行预测.Q2:个人信用风险评估在很多国家,政府机构会密切监控货款业务,银行需要明确解释申请者的贷款申请被拒绝或者批准的原因,这种可解释..
2021-07-13 14:48:05
1158
原创 运动模糊图像恢复处理-matlab实现
注:原问题来自第十一届数学中国数学建模网络挑战赛。基于图像退化原理建立了线性复原系统,即退化模型,PSF 函数与图像数据卷积模拟出运动模糊的退化图像。由于在模糊图像复原问题中,PSF函数的求解至关重要,因此利用倒频谱估计出模糊图像的模糊方向和模糊范围的像素值,从而建立了PSF函数。最后针对图像复原问题,通过维纳滤波寻找一个合适的复原图像使得它与原始清晰图像之间的均方误差达到最小,从而得到一张能还原程度较高的复原图像。在 matlab2017 软件上进行退化图像的仿真模拟:1. 首先,选取一张
2021-07-13 14:19:48
12559
4
原创 单剂量给药下仓室药物模型半衰期的计算
注:此篇为2019毕设,具体过程有需要可以联系本人,转载请附带原作者信息及链接。单剂量给药下如何研究药代动力学仓室模型半衰期研究将被探讨,通过在在线性与非线性的仓室模型范围中,具体分析给药体的药量与半衰期之间关系,以此达到给定需求的药量如何定量的求解仓室药物的半衰期这一问题的目标。利用应用数学建模的思想,可以对给药个体划分仓室。对仓室构建相应的参数与指标,通过这些参数和指标来分析说明仓室模型的组成,并得到相应模型的药物浓度求解方程组。由于给药剂量的确定可以得到相应浓度,因此利用药物浓度与半衰期之
2021-07-13 11:56:11
1514
MIT- BIH 正常窦性心律数据 nsrdb-16265
2022-03-24
基于一个微分对策问题的机器学习.pdf
2021-11-08
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人