给定一系列日期区间,求这段日期区间的交集

摘要: def CountDays(start_time, end_time, time_interval_list): from datetime import datetime, timedelta time_a = [] for start, end in time_interval_list: if 阅读全文
posted @ 2024-09-15 17:07 纵横二剑 阅读(0) 评论(0) 推荐(0) 编辑

R和Python的对应:outer()

摘要: 相信大家做模拟时一定碰到过变量之间的这种相关结构: $$ \mathbf{\rho}=\left(\rho_{ij}\right)_{n\times n}=0.5^{|i-j|} $$ 在模拟这个相关矩阵时,R中常用的函数便是outer()函数: outer(X=1:3,Y=1:3,FUN=func 阅读全文
posted @ 2023-05-14 20:04 纵横二剑 阅读(77) 评论(0) 推荐(0) 编辑

用字典代替'if-elif-else'

摘要: 在实际应用中,我们经常会需要采用if-elif-else控制语句以根据不同条件,作出不同的操作。if-elif-else固然可以,但是它也存在冗余的缺点,特别是当条件较多时这一缺点尤为明显。因此,本文考虑采用字典来处理这一情况。 Python的字典是由键值对构成,键是唯一的,并且每个键都有唯一的值对 阅读全文
posted @ 2023-01-13 20:59 纵横二剑 阅读(111) 评论(0) 推荐(0) 编辑

$\bf{X} \bf{X}^T$和$ \bf{X}^T \bf{X}$的非零特征值和特征向量之间的关系

摘要: 设$\lambda_i$为$\bf{X} \bf{X}^T$的特征值,对应的特征向量为$\mathbf{\alpha}_i$,则 \[ \bf{X} \bf{X}^T \mathbf{\alpha}_i=\lambda_i\mathbf{\alpha}_i \tag{1} \] $(1)\(式两边同 阅读全文
posted @ 2021-09-22 16:50 纵横二剑 阅读(177) 评论(0) 推荐(0) 编辑

SVD动态图

摘要: library(gganimate) library(gifski) # 构造数据 phi<-pi/6 theta<-pi/3 c<-matrix(c(cos(phi),sin(phi),-sin(phi),cos(phi)), nrow=2,byrow=TRUE) b<-matrix(c(5,0, 阅读全文
posted @ 2021-07-16 16:42 纵横二剑 阅读(103) 评论(0) 推荐(0) 编辑

巴恩斯利蕨 The Barnsley Fern

摘要: 巴恩斯利蕨学习链接 下面用R画一下: npts<-50000 point_mat<-matrix(data=NA,nrow=npts,ncol=2) list_fun<-list( f1=function(x,y){c(0., 0.16*y)} # R里的匿名函数 f2=function(x,y){ 阅读全文
posted @ 2021-01-17 11:57 纵横二剑 阅读(420) 评论(0) 推荐(0) 编辑

python3 小知识

摘要: ###本博客分为三部分,一是字典和函数结合,代替if-else判断语句;二是数组的元素带有括号;三是完整代码。图片略多。 ####1. 字典和函数结合,代替if-else判断语句 写函数时经常会遇到一个函数的某个参数有多个值可供选择的情况,并且每个可选值恰对应了一个子函数,这时需要用到if-else 阅读全文
posted @ 2021-01-14 12:08 纵横二剑 阅读(131) 评论(0) 推荐(0) 编辑

lxml库和BeautifulSoup库常用点小结

摘要: 算是本人的学习笔记吧,仅供个人学习使用。 以下内容摘自《Python3网络爬虫开发实战--崔庆才著》 1.lxml库 XPath 常用规则: | 表达式 | 描述 | | | | | nodename | 选取此节点的所有子节点 | |/|从当前节点选取直接子节点| |//|从当前节点选取子孙节点| 阅读全文
posted @ 2020-09-07 22:36 纵横二剑 阅读(260) 评论(0) 推荐(0) 编辑

城市间最短路径问题——R和Rcpp实现

摘要: 这里的最短路径问题也叫做相识问题,具体问题来自 https://www.math.pku.edu.cn/teachers/lidf/docs/Rbook/html/_Rbook/examples.html#examples-scicomp-citylink 刚好在学Rcpp,于是写了R版本和Rcpp 阅读全文
posted @ 2020-08-14 09:29 纵横二剑 阅读(574) 评论(0) 推荐(0) 编辑

正态分布的模拟

摘要: 题目来源于《R语言的科学编程与仿真》第18章第18题,原题如下: 参数为$\alpha$的柯西分布的概率密度函数为 \[f_{X}(x)=\frac{\alpha }{\pi \left ( \alpha ^{2}+x^2 \right )},-\infty<x<+\infty\] 使用逆变换法编写 阅读全文
posted @ 2020-07-15 11:32 纵横二剑 阅读(1049) 评论(0) 推荐(0) 编辑