- 博客(45)
- 收藏
- 关注
原创 使用四则运算实现异或
工作需要,场景如下:是否存在三个数 a, b, c 和一种运算 f (f 的定义域为 {a, b, c}), 使得:当 x ≠ y, f(x , y) = 0;当 x = y, f(x , y) = 1大家第一反应可能是异或,但是在我那个场景需要用到数值优化,需要用到纯粹的四则运算,不然没法微分或积分。经过尝试,偶然间发现可以这样构造。
2023-10-16 11:56:08
286
原创 简单爬虫项目练习
自学,参考书籍为Python爬虫开发与项目实战,具体参考了该书的第六章。在 Pycharm 上实现时发现有些库名更改及合并爬取的网站增加了反爬机制代码的框架与原书是一模一样的,若该文章有幸被原作者所浏览且认为这算抄袭,请联系我。以下代码会具体标出与原书不同的地方(以 #!!!!!!!!!!!标注)以及记录自己的想法。
2023-07-12 17:24:19
3259
原创 tidyverse笔记——dplyr包
本文介绍R数据分析的强力工具tidyverse中一个实用的包——dplyr,用于做数据变换,功能强大,之后会尽快更新关于这个包进阶的玩法和花样
2022-07-29 21:49:12
375
原创 R——避免使用 col=0
Rstudio自带的绘图区域和新建窗口或R原装gui是不一样的尽量避免col=0而应改为另一种妥当的方法是利用par事先设置背景色默认绘图背景并不是白色,而是“透明”。
2022-07-28 10:59:37
330
原创 R语言 mice包 Error in terms.formula(tmp, simplify = TRUE) : ExtractVars里的模型公式不对
R语言 Error in terms.formula(tmp, simplify = TRUE) : ExtractVars里的模型公式不对
2022-06-11 16:12:51
2392
原创 数值分析——线性方程组求解
数值分析——线性方程组求解缘起题目1(未完,待续)缘起清理磁盘的时候偶然发现大二下数值分析的实验作业还在,本着在丢弃之前可以放在网上以备不时之需的原则,我便发了上来。题目1分别用直接法、Jacobi迭代法、Gauss-Seidel迭代法求解下列线性方程组AX = b,其中A为五对角矩阵(n=20),b是除第一个分量是1外,其他分量都是0的列向量。(1) 直接法:观察到该题是形式特殊的五对角矩阵,利用matlab内置的用于产生稀疏矩阵的spdiags函数可以很容易生成系数矩阵A, A的产生代码如
2022-05-10 17:16:16
1600
原创 偏微分方程(二)
(未完,待补充)分离变量法大致思路和一般步骤双曲型方程的混合问题大致思路和一般步骤起源:机械振动或电磁振动(波动方程的解)总可以分解为具有各种频率和振幅的简谐运动的叠加,而每个简谐运动具有形式:eiw(t+sx)=eiwteiswxe^{iw(t+sx)}=e^{iwt}e^{iswx}eiw(t+sx)=eiwteiswx这正是物理上的驻波。这就表明波动方程的解可写为驻波的叠加,而每个驻波都只是只含变量 x 的函数和只含变量 t 的函数的乘积,即具有变量分离的形式。适用范围:有界区域下的混
2022-05-08 16:13:26
1606
原创 偏微分笔记(一)
偏微分笔记(一)方程的导出及定解问题的提法基本概念几个经典方程定解问题二阶方程的特征理论与分类二阶方程的特征二阶方程的分类二阶方程化为标准型(n=2)方程的导出及定解问题的提法基本概念一般形式:u(x1,x2,⋯ ,xn)(n≥2)u(x_1,x_2,\cdots,x_n)(n\ge 2)u(x1,x2,⋯,xn)(n≥2)F(x1,⋯ ,xn,u,Du,D2u,⋯ ,Dmu)F(x_1,\cdots,x_n,u,Du,D^2u,\cdots,D^mu)F(x1,⋯,xn,u,Du,D
2022-05-01 12:07:27
1820
原创 双指数平滑法一例
双指数平滑法简介缘起双指数平滑法题目题目解答缘起某日,一好友给我发了国外运筹学的一个小作业,共三题,问我是否有想法,看了一眼,看到了第一题是之前接触过的时间序列分析,便只做了这一道题,我也是在此过程中第一次接触到了双指数平滑法。(文章未写完,后续会更新)双指数平滑法我的参考资料来自 https://support.minitab.com/zh-cn/minitab/18/help-and-how-to/modeling-statistics/time-series/how-to/double-exp
2022-04-19 23:01:24
1304
原创 聚类模型——DBSCAN,K-means,层次聚类法实现
聚类模型——DBSCAN算法实现基本介绍伪代码评价Matlab实现Python实现基本介绍伪代码评价优点:无需输入类别数𝑘,可适用于任意形状的聚类簇;对噪声、异常点不敏感,能够识别出噪声点,不依赖于初始类中心的选取。缺点:样本集密度不均、聚类间距差较大时,聚类结果较差;样本集较大时,收敛时间较长; 对于参数𝜀和𝑀𝑖𝑛𝑃𝑡𝑠,调参难度较大,不同参数组合对聚类结果影响较大;对于高维数据,点之间极为稀疏,密度很难定义Matlab实现function [ result ] = DBSCAN(
2022-03-26 22:52:09
1878
原创 衡量连通图连通性一些指标(r-reachable, r-robust)
介绍了r-reachable和r-robust的基本概念,并在Matlab上进行代码实现
2022-03-21 22:20:26
1089
原创 C++——简易五子棋
评分制,基于一步预测,未完成,待完善。#include<iostream>#include<Windows.h>#define bound 15using namespace std;unsigned short position[bound + 1][bound + 1];//unsigned short bound = 0;bool model = 0;bool isover = false;short direction_x[8] = { 0,0,-1,1,
2022-02-06 16:28:35
316
原创 近世代数——群笔记(1)
Advanced modern algebra——Group 1lemma 1lemma 2Chinese Reminder's TheoremAn application of Chinese Remainder's Theoremlemma 1(Euclid’s Lemma) If a∣bca|bca∣bc and (a,b)=1(a,b)=1(a,b)=1, then a∣ca|ca∣cProof: Since (a,b)=1(a,b)=1(a,b)=1, there exists s,t∈Ns
2022-01-15 15:57:56
1416
原创 近世代数——环笔记1
Advanced modern algebra——Ring1ProblemAnalysisApplicationProblemIf RRR is a finite ring (e.g., RRR = Im\mathbb I_mIm), then there are only finitely many functions from RRR to itself, and so there are only finitely many polynomial functions. How to realiz
2021-12-27 11:31:16
764
原创 R语言时间序列分析小例
复习心烦,偶遇大作业,故摸鱼作业题目自由选取一组数据(可以是R 自带的数据集、或者其它来源,鼓励选取一些有趣的课题进行数据分析),利用我们这学期所学知识建立恰当模型(ARIMA、GARCH 等),小组成员(不超过三人一组, 自由组队)中的一人安排 5 至 10 分钟左右的PPT 课堂展示(教学周第16 周、无需展示代码,仅汇报初步结果即可,我给出点评后再完善后提交)。内容基本要求如下:• 说明数据来源背景, 需要分析解决的问题;• 包含完整的建模过程(拟合、估计、诊断检验、预测);• 结论或者建议
2021-12-24 20:28:36
2257
1
原创 R语言summary.lm信息进一步挖掘
R语言summary.lm信息进一步挖掘废话不多说,直接上例子(由于summary出来的结果是有截断的,以下比较会有一定误差,要精确的话可以自己从model里面取出相应系数)> data("faithful")> y=faithful$waiting> x=faithful$eruptions> model=lm(y~x,faithful)> summary(model)Call:lm(formula = y ~ x, data = faithful)R
2021-12-12 12:05:18
2305
原创 R语言线性统计模型实例
R语言线性统计模型实例毫无关系的前言正文变量选择检验技术变量选择准则附录(代码)毫无关系的前言这学期实在是太快了,一转眼又快到了令人疼苦的考试月,我是很讨厌大学这种考试的,一方面我发现确实考不怎么过别人,另一方面我很讨厌临时刷题,临时积极地问问题、熬夜。但在大环境下,我熬的夜确实比之前多了,但效率不见得高。这是本学期最后一篇文章了,作为学生,我常感叹自己能力不足,自从我在优快云发了自己的第一篇文章后,对于一些问题我开始变得好奇起来,我能否做一个这样的东西?我能否独立将其源代码写出来,这也是我写作的动
2021-11-28 21:34:16
2491
原创 基于R语言一元线性回归模型实例及代码
基于R语言一元线性回归模型实例及代码题目描述数据特征及可视化建立模型与初步评价(自己写lm()代码)回归诊断模型预测题目描述所用数据集——faithful(MASS包)# 加载数据包及查看问题背景library(MASS)data("faithful")?dataWaiting time between eruptions and the duration of the eruption for the Old Faithful geyser in Yellowstone National
2021-11-26 15:40:20
5102
原创 R语言及Rstudio入门小建议(一)
R语言及Rstudio入门小建议(一)前言相关背景及学习资料基本操作数据导入包(package)外部数据的导入数据类型和数据转换Rstudio的辅助功能强大的查询功能其余功能后记前言这学期选修了回归分析,恰巧大作业要用到R语言去实现,便把积灰已久的RStudio重新寻找了回来。R软件是统计课程的得力工具(数学专业算是,医学专业的话可能更多的是用SAS,SPSS等),如何快速地上手这门实用的软件呢?接下来,我将从我的使用经验与参考书来讲一讲。相关背景及学习资料百度百科对R的评价为:R是用于统计分析、绘
2021-11-20 12:05:49
17224
原创 sql——数据查询中的小例子
数据查询的一个小例子说明题目解题思路说明需要基本的数理逻辑知识和数据库知识题目(以学生-课程数据库为例)(例子选自《数据库系统概论》(第五版))Student表学号 Sno姓名 Sname性别 Ssex年龄 Sage所在系 Sdept201215121李勇男20CS201215122刘晨女19CS201215123王敏女18MA201215125张立男19ISCourse表课程号 Cno课程名 Cnam
2021-10-09 22:08:11
239
原创 数学模型课程期末复习提纲(下)
MathModeling(Ⅱ)主编: ???? 审阅:????第五章:预测算法回归预测拟合最常用的方法为最小二乘法,其原理是:因变量的实际值与拟合值之差称为残差,将所有残差平方后相加,即得残差平方和,”最好“拟合效果就是使残差平方和最小,于是运用极值原理,将”泛函问题“转换为”求残差平方和最小”问题(具体思想可参考数值分析相关内容)多元线性回归y=β0+β1x1+⋯+βpxp+εy=\beta_0+\beta_1x_1+\cdots+\beta_px_p+\varepsilony
2021-09-07 15:23:18
624
原创 数学模型课程期末复习提纲(上)
MathModeling(Ⅰ)主编: ???? 审阅:????第一章:数学建模概论数学模型可以描述为:对于现实世界的一个特定对象,为了一个特定目的,根据其特有的内在规律,做出一些必要打的简化假设,运用适当的数学工具,得到的一个数学结构。数学建模是指建立数学模型的全过程,包括表述、求解、解释、检验等。数学建模的基本步骤就是解决一个实际问题的步骤,主要包括下面几点:模型准备、问题分析、模型假设、模型建立、模型求解、模型检验。第二章:数学建模软件Matlab\text{Matlab}Ma
2021-09-07 14:59:37
1864
1
原创 数学学习——常微分方程
只是记录一道题目而已,并没有什么新颖观点。另外一道简单的题:设f(x)f(x)f(x)为连续函数,且由∫0xtf(t)dt=x2+f(x)\int_0^x tf(t)dt = x^2+f(x)∫0xtf(t)dt=x2+f(x)所确定,求f(x)f(x)f(x)解:两边对xxx求导,得 xf(x)=2x+f′(x), dydx=xy−2x~~xf(x)=2x+f'(x),~\frac{dy}{dx}=xy-2x xf(x)=2x+f′(x),&
2021-08-24 22:14:11
637
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人