zhmjunjun-优快云博客

原创初学者眼中的面向对象

面向对象，是一种编程思想，一个对象对应一个实体，例如：苹果(颜色，味道，美容功能等特性放在一起描述苹果)，梨子(颜色，味道，美容功能等特性放在一起描述梨子)面向过程的描述苹果和梨子：判断果皮黄色且带有斑点为梨子，否则为苹果；当实体很多时，其中任意一个实体属性发生变化时，程序的修改量通常会很大。用面向对象与面向过程编程的对比示例，方便形象理解。题目：图形接口画出四方形，圆形与三角形，当...

2019-12-11 22:22:03 218

原创 Matble怎样对矩阵进行加标签整合

这个代码主要针对，对不同类数据，加标签整合数据为 mnist_all.mat （因为时间隔得有点久我忘了当初下载网站）我主要用在GNMF相关处理中，其中对2D Hela image data数据也是类似加标签，希望对你有帮助%对MNIST数据的处理% 整合为一个数据矩阵%在因为每张图是28X28=784，所以每一行代表一张图%需要从所有图中选取500样本，每类选取50个...

2018-06-03 07:45:44 4149 2

原创对于Kmeans的聚类信息提取处理

因为网上有很多关于Kmeans现成代码，这里主要讲，如何对聚类信息进行提取（我前面博客给出过Kmeans的代码/matble）我主要做的是子图分离，先通过Kmeans方法，分离出k个子图，要对子图进行图构建，就需要知道，每个聚类块里面所包含的具体信息。其中用到胞元进行矩阵存储，下面也会举例胞元循环中的变化1：首先 label=Kmeans（fea，k） %label是每一个样本对应的聚类号...

2018-05-14 15:42:54 4072

原创 Matble编写替换矩阵每列最小值

同理可以替换每列最任何需要值（行同理）前面博客提到关于>也有这个功能，但这个更具逻辑性以及更实用该代码原来自EMR（.m代码）中的一步： D = EuDist2(data,landmarks); %矩阵D 大小为nSmpXpdump = zeros(nSmp,r); %dump，与idx为nSmpXr的零值矩阵idx = dump;for i =...

2018-05-09 22:39:02 469

简单介绍一下R中的几种统计分布及常用模型转载自：https://www.cnblogs.com/nxld/p/6060360.html 统计学上分布有很多，在R中基本都有描述。因能力有限，我们就挑选几个常用的、比较重要的简单介绍一下每种分布的定义，公式，以及在Ｒ中的展示。统计分布每一种分布有四个函数：d――density（密度函数），p――分布函数，q――分位数函数，r――随机数函数。比如，正...

2018-05-04 10:26:49 4494

原创 R语言编写NMF，其中迭代部分注释

认为，学算法，打印下来，一个一个看是最直观的，本来准备把NMF，每个函数都注解一遍，结果今天自己没事，看了下，乍一看一大片，用一个旁观者的视角看，我并不想继续读下去。所以代码给大家。因为实验室有打印机，我打印看，很快就知道什么流程了，纸上谈兵，上战场就~~~各种弱爆，想想就想哭。（不知道链接是不是一直有效，如果失效在Functional Biogeography of Ocean Microbes...

2018-04-27 20:02:17 3913 2

原创 R语言中is.与as.的区别与用法

来自R语言实战(第二版)，P79页R与其他统计编程语言有着类似的数据类型转换方式。举例来说，向一个数值型向量中添加一个字符串会将此向量中的所有元素转换为字符型。你可以使用下表1 中列出的函数来判断数据的类型或者将其转换为指定类型。表1 类型转换函数名为is.datatype()这样的函数返回TRUE或FALSE，而as.datatype()这样的函数则将其参数转换为对应的类型 is.data...

2018-04-27 13:14:44 25321

原创 R语言中标准化一个矩阵代码后求其相似矩阵注释

突然感觉R语言，有些地方比matble还简洁，但是我还是更喜欢matble，更直接些# Normalize a matrix# Args:# m: the matix# L: by rows (L=1) or columns (L=2)# Returns:# The normalized matrix#标准化矩阵##参数：#m ：矩阵#L ：按行（L = 1）或列（L = 2）#返回：标准化...

2018-04-26 22:27:39 3145

原创 R语言中矩阵逆序转置代码解释

一年前都记得，这次要修改，忘了，所以决定记下了。PlotM<- function(Z,heatc=12) { m<-ncol(Z); n<-nrow(Z); MM<-matrix(0,m,n); MM[1:m,]<-t(Z[,1:m][n:1,]); image(MM,axes=FALSE,col=heat.colors(heatc));}下面是解释Plot...

2018-04-26 21:34:13 6661

原创 matble中如何使用 > ，以及将矩阵元素替换

表示方式A>Bgt（A,B）A>B返回一个设置为逻辑1要素的逻辑阵列（真），A大于B；否则，该元素是逻辑0（假）。测试只比较数字阵列的实部。GT返回逻辑0（假），A或B有南或未定义分类要素。 >> A = [1 12 18 7 9 11 2 15]A =1 12 18 7 9 11 2 15>> A > 10ans =0 1 1...

2018-04-13 14:50:55 648

原创 matble当+与：组合在一起的时候表示啥

当m=a+b：c 表示返回b：c中位置在第a个位置之后的所有数值，m为向量.一般在循环里面用的比较多因为不善表达，就用示例说明>> m=1+1:5m = 2 3 4 5size(1:5)5因为是第a个位置之后，为了说明我把显示出来>> 3:8ans = 3 4 5 6 ...

2018-04-12 19:30:01 298

原创 matble中bsxfun（）函数

下面这句是在百度百科上面看到的非常形象，比官方给出的好，谢谢百度百科bsxfun（）函数的编辑作者“两个非“单一维度”相互匹配的数组a和b做函数fun运算时，bsxfun会隐含扩充a或b使得a和b结构相同，以便实现逐元素运算。使用函数bsxfun可以避免用循环结构编程。bsxfun调用格式：bsxfun(@已有定义的函数名，数组1，数组2）bsxfun(@( 数组1，数组2)函数体表达式，数组1...

2018-04-11 17:55:16 1691

原创互信息的理解

互信息指的是两个随机变量之间的关联程度，即给定一个随机变量后，另一个随机变量不确定性的削弱程度，因而互信息取值最小为0，意味着给定一个随机变量对确定一另一个随机变量没有关系，最大取值为随机变量的熵，意味着给定一个随机变量，能完全消除另一个随机变量的不确定性这是我在知乎上看到我最想要的结果。我怕自己下次又找，所以写了这个博客感谢知乎作者：maple链接：https://www.zhihu....

2018-04-10 22:22:48 891

原创 matble中legend()函数示例

对于legend函数，我觉得示例更形象，关于解释网上给出的都很完善，或者matble中help legend哦legend() 将图例添加到图形x = -pi:pi/20:pi;y1 = sin(x);y2 = cos(x);plot(x,y1,'-ro',x,y2,'-.b')legend('sin(x)','cos(x)')图例1x = -pi:pi/20:pi;y1 = ...

2018-04-10 13:01:53 2186

原创 matble中sawtooth()函数

sawtooth() 锯齿或三角波表示方式：sawtooth(t)sawtooth(t,width)sawtooth(t)为时间矢量t的元素生成周期为2π的锯齿波。 sawtooth（t）与sin（t）相似，但会产生峰值为-1和1而不是正弦波的锯齿波。锯齿波在2π的倍数处被定义为-1，其他时间以1 /π的斜率线性增加。sawtooth(t,width)会生成一个修改后...

2018-04-10 12:38:16 31459

原创 matble中disp（）函数

来自matble 中help disp术语：disp（X）disp（X）显示变量X的值，而不打印变量名称。我们做常见的是另一种显示变量的值，且打印变量名称“X =”传统不加disp加disp>> A=[15 150];>> AA = 15 150>> A=[15 150];>> disp(A) 15 150>> S ...

2018-04-09 23:00:55 1481

原创 matble中&

matble中help &两种语法A & Band(A,B) 两者表达含义一样，但and（A，B）较少使用。逻辑与，返回值1，或者0；当A 与 B的数值都非零时，返回1当A 与 B其中一个为零，或者两者都为零，返回0如下：>> A = [5 7 0; 0 2 9; 5 0 0]A = 5 7 0 0 2 9 5 ...

2018-04-09 22:21:45 493

原创 matble的K-means聚类，函数参数解释

function [label, center, bCon, sumD, D] = litekmeans(X, k, varargin)litekmeans K-means聚类，通过matble编写的矩阵加速运算X，k：将N×P数据矩阵X中的点划分为K个簇。 N为点个数，p为点的特征个数 label = litekmeans(X, K) label：此分区群集中，点到群集质心的群集内...

2018-04-09 18:24:05 5176

原创 matble中函数rand()随机数生成器（二）

在前面matble中函数rand（）中，介绍了rand（）函数的一般用法，通过rand()或者randn()函数使用‘seed’，‘state’，‘twister’来控制随机数发生器，例如rand('twister',5489)，就是matble重置随机生成数编译器，使在不同方法（算法）下有相同的初始化，同时又可以保证实验结果可以复现。rand('seed',sd)randn('seed',s...

2018-04-08 22:22:28 2811

原创 matble中函数rand()（一）

来源 matble中help rand 其实都知道rand('twister',5489)，就是重置matble随机生成数编译器，使在不同方法（算法）下有相同的初始化，同时又可以保证实验结果可以复现。下面结合matble中的help rand来具体，rand函数默认是生成（0,1）之间随机数的函数rand（）函数的具体形式有如下几种：X = rand 返回一个随机数，该...

2018-04-08 22:12:10 1590

原创 matble矩阵中取行或者列

直接看程序就可以知道先创建一个矩阵；rand（）为0-1的随机数组成的矩阵，ceil（）为取整，往大整数的方向取；A=rand(2,4)A = 0.6787 0.7431 0.6555 0.7060 0.7577 0.3922 0.1712 0.0318>> C=ceil(10*A)C = 7 8 ...

2018-04-07 17:41:35 599

原创 matble中magic（）函数

M = magic（n）返回由整数1到n ^ 2构成的n乘n矩阵，行和列的和相等。顺序n必须是大于或等于3的标量。因为是nXn矩阵，所以矩阵元素从1……n^2全排列。>> magic(3)ans = 8 1 6 3 5 7 4 9 2>> sum(magic(3))ans = ...

2018-04-07 12:13:57 17746

原创 matble中 .，.^，./，.\的例子

1》 **为简单的矩阵乘举例方便下面与.*的对比>> a=sqrt(sum(magic(3)))a = 3.8730 3.8730 3.8730>> norms=max(1e-15,a)norms = 3.8730 3.8730 3.8730>> A=full(spdiags(norms',0,3,3))...

2018-04-07 12:08:18 577

原创 matble中spdiags()函数（翻译+理解）

翻译来源 matble中的help spdiags 转载请注明出处https://blog.youkuaiyun.com/zhmjunjun/article/details/79790957 谢谢spdiags函数推广了斜对角函数的应用，可以进行四种不同的操作，如下：B = spdiags(A) 从m-by-n矩阵A中提取所有非零对角。B是一个min(m，n)...

2018-04-02 16:25:37 12400 1

原创 GNMF代码理解 matble

没办法，我很菜，但自己要用时发现找不到自己想要的博客。于是记下笔记，希望可以方便和我一样无助的人，能够更快上手。代码链接：http://www.cad.zju.edu.cn/home/dengcai/Data/ReproduceExp.html#GNMF 如果只要GNMF代码，上面链接原文作者给出，可以运行下面是里面一个例子代码说明当不知道某个函数怎么用时，可以直接在运行...

2018-03-30 22:25:53 4031 5

原创 GNMF其中一步的matble程序详解

数据来源http://www.cad.zju.edu.cn/home/dengcai/Data/TextData.html里面的这个数据是一个文件项矩阵解释的代码也是这个事例里面的为了matble运行方便，将下载好的数据集放在matble的运行文件中下面是matble代码load('USPS')faceW = 16; faceH = 16; numPerLine = 20; ShowLine = ...

2018-03-29 21:23:09 2361 2

转载 R基础函数的编写

函数的编写因为不是豆瓣用户，这篇文章来自https://www.douban.com/note/279077707/豆瓣用户：盐池里的萝卜。。。的R语言实际上是函数的集合，用户可以使用base，stats等包中的基本函数，也可以自己编写函数完成一定的功能。但是初学者往往认为编写R函数十分困难，或者难以理解。这里对如何编写R函数进行简要的介绍。函数是对一些程序语

2017-11-29 11:41:18 704

原创时间序列笔记（七）

平稳序列建模本章主要讲：实际例子当中怎样选择ARMA模型，并给出具体对应例子。同时介绍了计算ARMA模型序数过程中选用的估计方法，矩估计，最大似然拟合，最小二乘在时间序列应用。理论部分有推倒，具体计算是调用现成R语言包（还没讲到具体代码）

2017-11-28 10:23:39 347

原创时间序列笔记（六）

在如何选择三个模型上主要看：自相关系数与偏相关系数拖尾，截尾然后选择模型用来建模！！！！！//////////////////////////////////////////////前六个笔记的核心MA模型MA模型与AR模型的不同是，MA通过白噪声及其相关序列建模，其自协方差，与自相关序数都是截尾的AR模型是通过时间点及其序列建模，其自协方差是拖尾

2017-11-27 15:29:59 439

原创时间序列笔记（五）

数学推理/////////////////////////////////////////自相关序数与偏自相关序数在选择AR（p）模型上的使用具体例子：随着间隔期数的拉长，自协方差序数是递减的趋势，趋于零（称为拖尾性）自相关图，可以用来判断是否可以用AR模型拟合图上的条代表自相关大小

2017-11-27 11:47:48 506

原创时间序列笔记（四）

拿到数据后，通过平稳性与随机性判断之后如何处理方法性工具ARMA模型平稳序列建模序列预测方法性工具：差分运算，延迟算子，线性差分方程差分运算：延迟算子：p表示延迟的时间尺度B不是常数，它是一个延迟算子，具有四则运算的一些性质线性差分方程：往往先求齐次，然后再往非齐

2017-11-26 17:29:55 1177 2

原创时间序列笔记（三）

纯随机性纯随机性检验纯随机序列的定义纯随机性性质纯随机性检验纯随机性序列也成为白噪声序列，它满足如下性质：纯随机数列也是平稳性数据之所以检验是否满足纯随机性，是为了在检验平稳性数据后排除没有研究价值的数据时序图（常数均值附近波动），系协方差相关系数（自相关图）才可以确定数据是否可用纯随机性检验*原理*条件*检验

2017-11-26 15:51:31 1274

原创时间序列笔记（二）

平稳序列1：时间序列的预处理平稳性检验：特征统计量平稳性时间序列的定义平稳性时间序列的统计性质平稳性时间序列的意义平稳性检验特征统计量：概率分布的意义：统计特征完全用联合分布函数或联合密度函数决定例如：均值，方差这些都是在分布函数中可以看出来的，有了密度函数就可以得到分布函数等等联合分布，，类似同时投两颗骰子时间序列概率分

2017-11-23 22:56:41 1612

原创时间序列笔记（一）

因为这是笔记，也是自己认为比较重要的，后期会慢慢写自己的总结，以及相关程序（R）本章主要介绍时间序列发展史！常用分析方法例：多元回归（粮食中人口光站降雨）y=a+b1人口+b2光照+b3降雨+残差1：数据项之间不能自相关2：多重共线性（数据项之间相互相关）3：异方差（模型不同总体之间有相同的方程，三种之一任何一种都可能导致造成：统计模型失效，不然会使假设检验失效

2017-11-23 16:58:35 670

zhmjunjun的博客