- 博客(79)
- 资源 (4)
- 收藏
- 关注
RSS订阅原创 MATLAB实战——图片噪声的添加与平滑处理
加入高斯白噪声并作平滑处理clear,clc%% (1) 读取数据A = imread('LenaRGB.jpg');%% (2)在数据中加入高斯白噪声作平滑处理% 加入高斯白噪声B = double(A)/255;[m,n,p] = size(B); % 大小white_noise = 0 + 0.1*randn(m,n); % 白噪声C = B + white_noise; % 加白噪声D = C*255; %将像素范围扩大至0--255D2 = uint8(D);
2020-07-10 20:14:33
1788
原创 R语言实战——基于KNN聚类的时间序列分析预测
---title: "Summary of Reading"author: "ChenWei"date: "2020/6/24"output: word_document---```{r setup, include=FALSE}knitr::opts_chunk$set(echo = TRUE) 这一次我阅读的是Francisco Martinez等著的《Time Series Forecasting with KNN in R: the tsfknn Package》，这篇文章
2020-07-10 20:06:23
5381
原创 Python实战——过采样数据的处理之改进的SMOTE算法
理论实现数据：2,4,4,2,3,3,02,2,3,2,1,2,12,3,3,2,2,3,12,2,3,2,5,2,12,3,4,3,2,2,01,5,4,2,2,1,02,5,4,4,4,1,02,4,4,4,4,7,01,4,2,2,2,4,02,3,3,2,5,2,01,2,2,2,1,3,02,5,4,2,5,2,01,5,3,4,4,2,12,6,2,2,4,1,02,4,4,3,4,2,02,2,4,3,5,3,01,6,3,3,4,7,01,4,4,
2020-07-10 19:54:42
22316
8
原创 R语言实战——谱系聚类分析与动态聚类分析
一、理论二、实战分析x1 <- c(5,7,3,6,6)x2 <- c(7,1,2,5,6)X <- cbind(x1,x2)# 计算初始距离矩阵dist(X) # 默认计算欧式距离dist(X,diag=TRUE) # 显示对角线元素dist(X,diag=TRUE,upper = TRUE) # 显示对角线元素和上三角元素dist(X,method = 'manhattan')dist(X,method = 'minkowski',p=1)
2020-07-10 19:49:47
2082
原创 MATLAB实战——手写个人签名提取与初步处理
现在很多工作都转移到了线上，有些情况下需要将签名插入word，但是我们有时候是没法打印后再签名的，因此我就写了这么个程序。一、读取图片%% 读取图片pic = imread('C:\Users\cw817615\Desktop\签名.jpg');threshold = graythresh(pic);%确定二值化阈值pic_bw = im2bw(pic,0.55);%对图像二值化，黑色的为0，白色的为1imshow(pic_bw) % 展示处理前的二值化图片pic_backup = p
2020-07-10 19:10:07
1552
原创 微分方程的数值解法与符号解法
文章目录数值解符号解数值解问题Matlab实现main.m 主函数clear,clcx0 = 0; % xminx1 = 1; % xmaxh = 0.1; % 步长y0 = 1/2; % y0为初始值% 欧拉法[xy{1}(:,1), xy{1}(:,2)] = forwardEuler(x0, y0, x1, h);% 改进欧拉法[xy{2}...
2020-05-01 13:41:05
3239
原创 形如e^(ax^2+bx+c)的积分公式的证明
形如∫ex2dx\int e^{x^2}dx∫ex2dx的积分是一种典型的很难积的函数，普通方法很难把它积出来，得用巧劲。以下为证明过程：直接求积分很难，我们先把它平方，将其转为二重积分：(∫ex2dx)2=∫ex2dx⋅∫ey2dy=∬ex2+y2dxdy(∗)\begin{aligned}\left(\int e^{x^2}dx \right)^2 &=\int e^{x^...
2020-04-24 17:24:35
5800
转载 R语言实战——百分条图与雷达图
文章转载自“医学统计与R语言”公众号，文章链接为：文章链接百分条图输入1library(ggplot2)library(ggthemes)library(ggsci)library(rio)percentbar <- import("percentbar.xlsx")percentbar输出1 year government society Indivi...
2020-04-20 18:42:28
2476
原创 R语言实战——ROC曲线的绘制
前言：以前使用Matlab绘制ROC曲线常常是工具箱有就画，没有就不画，而且在想画的时候工具箱恰恰就没有，很纳闷。然后无意间发现了一篇用R语言绘制ROC曲线的文章，赶紧学了并分享出来，以备不时之需。先通过一个例子来讲解一下参数的作用，使用的数据是大名鼎鼎的Iris数据集，R语言自带。数据处理第一步当然得处理一下数据。默认的Iris数据集有三类鸢尾花，我目前的理解是只有二分类才画的出ROC曲...
2020-04-20 18:30:56
25622
1
原创 离散型随机变量的分布函数的绘制
前言：绘制连续型的分布函数很容易，直接根据分布函数计算函数值即可。但是对于离散型随机变量而言，没有已知的分布函数，只能使用经验分布函数或者说是累积分布函数进行近似。下文以离散分布中比较经典的二项分布、泊松分布以及几何分布为例绘制它们的经验分布函数。二项分布二项分布的理论部分可以参照二项分布百度百科Len = 100000; %% 二项分布Bin(n,p)N = 50; % 二项分...
2020-04-20 16:36:33
9453
原创 GM(1,1)灰色预测及相关检验指标的MATLAB实现
本篇文章的代码实现了以下三大方面的功能：一、计算级比和光滑比并做级比检验；二、序列的灰色预测；三、精度检验，主要做了以下内容：①相对残差Q检验（MAPE）；②关联度检验；③方差比C检验；④小误差概率P检验1 灰色预测模型本部分内容参照了王佳佳的硕士学位论文《京津冀地区高速铁路对区域经济影响研究》，我觉得这篇论文讲灰色预测讲得蛮详细的，如果想学灰色预测的可以看看。2...
2020-04-20 16:12:11
15457
5
原创 形如1/(a+bsinx)的积分公式的证明
需要求解的积分具有如下形式：∫1a+bsinxdx\int \frac{1}{a+b \sin x} d x∫a+bsinx1dx这种积分比较常见，可以考虑记一下它的结果，结果为：∫1a+bsinxdx=2a2−b2⋅arctan⋅{atan(x2)+ba2−b2}\int \frac{1}{a+b \sin x} d x=\frac{2}{\sqrt{a^{2}-b^{2}}} ...
2020-04-18 15:03:42
7481
1
原创 R语言学习——对整合和重构的理解
文章目录1 整合数据2 reshape2包2.1 融合2.2 重铸前言：R中提供了许多用来整合（aggregate）和重塑（reshape）数据的强大方法。在整合数据时，往往将多组观测替换为根据这些观测计算的描述性统计量。在重塑数据时，则会通过修改数据的结构（行和列）来决定数据的组织方式。本篇文章描述了用来完成这些任务的多种方式。1 整合数据在R中使用一个或多个by变量和一个预先定义好的...
2020-04-17 13:32:34
1109
1
原创 R语言实战——主成分分析理论推导与R语言实现
目录1 总体主成分1.1 主成分的定义与导出1.2 主成分的性质1.3 从相关矩阵出发求主成分2 样本主成分2.1 从S出发求主成分2.2 从R出发求主成分3 相关的R函数以及实例3.1 `princomp`函数3.2 `summary`函数3.3 `loadings`函数3.4 `predict`函数3.5 `screeplot`函数3.6 `biplot`函数4 实例附录——PCA高级散点图的...
2020-04-15 20:06:35
4155
2
原创 R语言实战——距离判别、贝叶斯判别、Fisher判别理论详细推导与R语言实现
前言判别分析是用以判别个体所属群体的一种统计方法，它产生于20世纪30年代，近年来，在许多现代自然科学的各个分支和技术部门中，得到广泛的应用.例如，利用计算机对一个人是否有心脏病进行诊断时，可以取一批没有心脏病的人，测其p个指标的数据，然后再取一批已知患有心胜病的人，同样也测得p个相同指标的数据，利用这些数据建一个判别函数，并求出相应的临界值，这时对于需要进行诊断的人，也同样测其p个指标的数据...
2020-04-06 22:05:14
24845
7
原创 Python实战——选择最佳旅游攻略，让旅游更加便捷(爬虫实战)
前言： 打算和老弟去西安来一个说走就走的旅行，但是网上攻略太多了看得头皮发麻，但是仔细看的话每条旅游攻略都有特定的参数条件的，比如人数、价钱、游玩时间，也就是说我们可以通过筛选这些条件初步获取我们满意的攻略。1 前期准备这次爬的是去哪儿网，网站大概长这样如果我们搜一个旅行地点，会得到这么一个网页网页的网址为：http://travel.qunar.com/travelbook/list...
2020-04-03 17:40:21
2338
3
原创 Python实战——表情包爬虫，让你拥有数之不尽的表情包(^_^)
前言： 每次发表情包的时候总是会愁自己表情包太少了，用来用去就那些，很头疼，现在写一个爬虫爬取表情包，以后再也不用愁没表情包了！爬取的网站网址为：表情包网址先来看看网站长啥样看到这么多表情包就很开心，一个一个保存太麻烦了，爬就完事了！源代码：# _*_ coding: utf-8 _*_'''Author:Chen Weidate: 23:12Repetition i...
2020-04-02 14:03:58
620
原创 R语言学习——一元与多元正态分布检验(也可以用于其他分布的检验)
生成随机数：# 生成随机数set.seed(1230) # 随机数种子y1 <- rnorm(100); # 标准正态分布N(0,1)y2 <- rexp(100,2); # 参数为2的指数分布Exp(2)y3 <- rt(100,1); # 自由度为1的t分布t(1)y4 <- -y2; # -Exp(2)1 图像法1.1 直方图...
2020-03-27 19:12:41
10785
原创 瑞利分布,对数正态分布,韦布尔分布,K分布参数极大似然估计的MATLAB实现
前言： 找到的第一篇文献，里面的密度函数出现了明显问题，在CSDN里面找也经常是只找到一个或者两个分布的估计，但是有时候得四个分布一起用，找来找去难免会有些麻烦，找到一篇文献：邓泽怀的硕士学位论文《基于实测数据的海杂波统计建模》，里面详细地讲解了上面提到的四种分布参数的矩估计、极大似然估计，并作了拟合优度检验。实践证明这里面讲到的关于参数估计的内容非常准确。文章目录1 背景知识2 典型分布模型2...
2020-03-22 12:10:56
8522
2
原创 MATLAB求解非线性0-1整数规划(基于Yalmip工具箱)
前言： 折腾了老半天，直接在CSDN里面搜找不到代码，百度刚开始也就找到一个bnb20工具箱，可那已经是2009年的玩意了。不过还是很幸运，在我绝望到想用遗传算法求解的时候，终于找到了一个叫Yalmip的工具箱，讲道理，功能强大且好用。1 Yalmip工具箱的下载与安装Yalmip的下载(建议在我给的这个链接里下载，官网下载的速度实在是emmmm)Yalmip的安装2 Yalmip的使用...
2020-03-19 16:51:57
8302
3
原创 基于Box-Cox变换的改进ARMA模型的MATLAB实现
前言： 第一次听说Box-Cox变换，发现原来这玩意这么厉害，值得写一写记录一下。以下内容参考了栗然等人写的论文《基于Box-Cox变换的风电场短期风速预测模型》1 问题背景风力发电是解决当前突出的能源和环境问题的有效手段，因而得到了普遍重视和大规模开发利用，是目前世界上增长速度最快的能源开发形式。尤其近十几年来其发电成本大幅下降，已具备与传统常规电源竞争的潜力。目前，国内外对于风力发电的...
2020-03-19 11:58:50
3366
2
原创 R语言学习——多元数据统计特征的分析(含马氏距离的计算及矩阵分解)
**前言：**将多元统计分析课程上学到的知识以编程的形式展现。本篇博客使用到的数据如下："Wind" "Solar radiation" "CO" "NO" "NO2" "O3" "HC"8 98 7 2 12 8 27 107 4 3 9 5 37 103 4 3 5 6 310 88 5 2 8 15 46 91 4 2 8 10 38 90 5 2 12 12 49 8...
2020-03-18 18:14:22
5996
原创 一阶时滞微分方程三种求解方法的MATLAB实现及稳定性分析
**前言：**大学期间只学习过《常微分方程》，没想到有些学校竟然还学《时滞微分方程》，于是找到一本由内藤敏机(日本)等著，马万彪等译的《时滞微分方程——泛函数微分方程引论》(有需要的可以私聊，CSDN貌似上传不了书籍，说侵权emmm)，看着头秃，不过受到不少启发，尤其是对Logistic方程的改进，真真是长见识了。没找到有人用欧拉法解一阶时滞微分方程的，于是一不做二不休便用MATLAB实现了一下下...
2020-03-17 15:31:15
10912
9
原创 基于遗传算法的多种运输工具或带时间窗的路径优化问题(VRP)的求解（MATLAB）
文章目录1.背景2. 单种运输工具带时间窗2.1 带时间窗车辆路径问题的描述2.2 遗传算法的求解2.2.1 编码方案的设计2.2.2 种群初始化2.2.3 约束处理与适应度函数2.2.4 选择算子2.2.5 交叉算子2.2.6 变异算子2.2.7 其他部分2.3 粒子群算法的求解3. 多种运输工具不带时间窗3.1 立体交通介绍3.2 遗传算法求解3.2.1 染色体编码设计3.2.2 适应度函数3...
2020-03-14 18:29:08
5210
11
原创 Python实战——百度图片的爬取
有人找我爬取一下百度图片，好像写这个的人比较少，于是我借鉴GitHub上的代码，根据需求修改后有如下代码# -*- coding: utf-8 -*-import requestsfrom threading import Threadimport reimport timeimport hashlibimport osclass BaiDu: """ 爬取百度...
2020-03-14 16:57:00
162
原创 统计学——依据概率密度生成随机数并计算分位数(以一元p-范分布及其抽样分布为例)
前言： 最近在依据概率密度生成随机数方面遇到不少麻烦，找遍了CSDN，大多是生成常见分布的随机数，很少有讲生成任意给定概率密度的随机数，即便有部分人写了，但都没给出比较好的模板。有心人天不负，在一番摸索之后，终于是写出来了，顺便还把分位数的两种解法给写出来了，以下是需要解决的问题与相关代码。1. 背景知识1.1 一元p-范分布 在研究线性回归模型(y=βx+ε)(y=\beta x+\va...
2020-03-14 16:43:33
2097
原创 Python实现词云图两例
以下两个实战使用python3.7实现，编译软件为Jupyter Notebook实战1：App Store钉钉评价分析评论爬取App Store 已经给出一个软件评论接口import requestsimport pandas as pdfrom pandas import DataFrameflag = [1,2,3,4,5,6,7,8,9]urllist = []for ...
2020-03-03 16:33:58
812
原创 R语言学习——最全R语言向量的生成与命名
首先需要说明的是，在R语言中向量和列表是不一样的，向量中每一个元素的类型都是一样的，而列表反之，这里的列表与Python中的列表有几分相似。以下是分类情况：同样类型的数据不同类型的数据一维向量(vector)列表(lists)R语言中生成常用向量主要有这么几个方法：方法1：...
2020-03-03 15:58:41
21765
1
转载 python爬虫与数据分析之《向往的生活爬取》
项目一：《向往的生活》弹幕爬取与分析项目一参照于公众号“超哥的杂货铺”《《向往的生活3》弹幕爬取与分析》一文链接：公众号文章链接.本文在参考文章的基础上添加了一些自己的理解，多了不少注释1.1弹幕的爬取注：本次代码主要使用requests库，抓取结果存储在CSV文件中1.1.1 网页分析在芒果TV网页版打开第5期节目，等待广告加载完毕，同时打开chrome开发者工具的network选...
2019-09-24 17:42:56
886
2
转载 机器学习——神经网络篇
一、生物神经元简介1.1生物神经元的结构参考资料：https://mp.weixin.qq.com/s?__biz=MzA4MTg3NDgwNQ==&mid=2247486270&idx=1&sn=12e35b20ad8464d051062a4e40810ae0&chksm=9f8f1124a8f898323cc0335c012c1ba020c9c1da972...
2019-05-17 11:49:00
1455
原创 Python面向对象程序设计总结
一、总结图1 Python类主要内容思维导图二、补充说明1、类只负责刻画，不负责执行类是显示生活的一种刻画，是一种将现实抽象之后形成的模板，以构造对象的方式具体去实现类的模板2、构造函数__init__只能返回NoneType类型的内容，否则会报错3、我们以点运算符调用对象的变量的时候，系统会先检查对象是否有该名称的变量，如若没有则检查该对象所属类是否有该变量，再没有则检查该类...
2019-04-17 15:45:22
2359
原创 python正则表达式图文总结（二）
推荐参考资料：http://www.runoob.com/python/python-reg-expressions.html详细的总结内容请见我的图文总结（一）文章这篇文章给出一个实例# 导入正则表达式的re模块import re# 取出string8中所有的天气状态string8 = "{ymd:'2018-01-01',tianqi:'晴',aqiInfo:'轻度污染'},"...
2019-04-16 11:45:45
165
原创 JSON与python
一、什么是JSON？全称：JavaScript Object Notation， 即JavaScripe对象标记本质：JSON是一种轻量级的数据交换格式，它是我们为了便于跨语言交换数据而定制的一种数据格式表现形式：字符串，符合JSON格式的字符串叫做JSON字符串二、JSON的优势易于阅读易于解析网络传输效率高跨语言交换数据三、在python中使用JSON(反序列化)注意...
2019-04-15 19:52:31
485
原创 python正则表达式图文总结（一）
一、主要正则表达式图1 正则表达式主要内容总结二、正则类型1、查找①表达式：re.findall(pattern, string, flags)其中pattern为正则表达式，string为待查找字符串，flags为图一中所述的模式②查找结果以列表的形式返回2、替换①表达式：re.sub(pattern, repl, string, count, flags)其中patter...
2019-04-15 18:10:14
489
原创 Python包、模块易错点
1、python项目的组织架构示意图python程序结构：包―>模块―>类―>变量和函数（这两个为类的属性，不属于Python结构），就好比学生的身高与体重不属于学校的组织架构一样，它们只是学生所特有的属性2、在Python里面包就是文件夹，但文件夹不一定是包，包是含有__init__.py模块的文件夹3、Python中的import都是导入模块，导入包其实就是导入模块的...
2019-04-06 15:51:02
193
原创 python数据类型总结
参考资料：https://blog.csdn.net/answer3lin/article/details/86430074一、python中数据类型分类思维导图二、值得注意的是：1、Python中的变量都是指针，因此是没有类型限制的，且指针的内存空间大小是与类型无关的，其内存空间只是保存了所指向数据的内存地址。2、值类型是不可变的（immutable），这种不可变是指该值类型的变量指向...
2019-04-01 19:57:48
9683
原创 lingo中@size@for@sum函数的使用
@size:LINGO中的@size（）函数用于确定集合中元素的个数。比如你的集合是注意：在使用size的时候直接在@size()括号里写上集合名就行，不需要写@size(s1(i))中的(i)s1/a b c d f g /:x;当调用m=@size(s1)时，它返回的一个数值6，因为s里面有6个元素。@for:该函数用来表示线性规划中的约束条件sets:r/1..8/:d;...
2019-03-25 11:39:14
25290
3
原创 对Java中聚合与耦合的理解
在程序设计时，为了避免程序的冗杂，减少各个类之间的依耐性，经常使用低耦合和高内聚的设计理念。为了达到这个目的，我们采取以下措施：①降低耦合：我们常用封装来降低耦合。例如，假如现在又A、B两个类，为降低它们的耦合度，我们将A中直接调用B中成员变量的部分尽量改为以调用B中方法的方式调用变量，不要让B对A“掏心掏肺”，要让B懂得有所保留，此为封装。这样的好处是保护了B变量的隐私性。②提高聚合：我们...
2019-03-23 10:26:58
918
单总体与多总体ANOVA方差分析.zip
2020-04-13
基于Python的两例词云实战代码与数据.zip
2020-03-03
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人