随笔分类 - Mathematical Modeling
摘要：主要看大神们如何使用，先模仿。 [~,y] = sort(v)，如果v是向量，那么y返回的是v中的下标（不好表达），看下面的，发现y是是下标，对应的是c中元素*（由小到大排序）。使用c(y)就得到了由小到大的排序向量。K>> c = [4 5 1 6 2 8 3]c = 4 5...
阅读全文
posted @ 2015-07-13 09:35
加拿大小哥哥
摘要：function m f=@(x) x^2; y(f,3); function y(f,x) disp(num2str(f(x))); endend 函数调用另一个函数的时候，把另一个函数名作为参数变量直接传递就好。 如代码中，函数y调用函数f，把f...
阅读全文
posted @ 2015-07-12 09:22
加拿大小哥哥
摘要：一.极限问题的解析解 1.1 单变量函数的极限 格式1： L= limit( fun, x, x0) 格式2： L= limit( fun, x, x0, ‘left’ 或 ‘right’)>> syms x a b;>> f=x*(1+a/x)^x*sin(b/x);>> L=limit(f...
阅读全文
posted @ 2015-06-30 20:30
加拿大小哥哥
摘要：各种允许的比较关系 >, >=, > A=[1,2,3;4,5,6;7,8,0]A = 1 2 3 4 5 6 7 8 0>> find(A>=5), ％大于或等于5元素的下标 ans = 3 5 ...
阅读全文
posted @ 2015-06-28 17:12
加拿大小哥哥
摘要：1.生成N阶[a,b]均匀分布数组>> unifrnd(3,5,5,5)ans = 3.8651 4.6677 4.8115 4.3456 4.8560 4.0241 3.4079 3.5713 3.3737 3.5913 3.480...
阅读全文
posted @ 2015-06-22 13:45
加拿大小哥哥
摘要：1.C= min(a)：返回最小值；我原来以为如果a是行向量,min(a)返回a本身，因为我记得min(a,1)是按列找最小，这是默认的。doc min发现，只要a是向量，那么返回最小值。 2.C= min(a,[],dim)，dim分别取1,2和我理解的一样。 3.C= min（A，B），...
阅读全文
posted @ 2015-06-22 13:16
加拿大小哥哥
摘要：1.m先max(x)求出每列最大值，返回行向量，再max对行向量求出最大值，就是max(max(x))。 注意：max(x)，不管x是行列向量，只要是向量，那么就返回一个值。 2.先x(:)转为按列的列向量，然后max，就是max(x(:))。
阅读全文
posted @ 2015-06-22 10:02
加拿大小哥哥
摘要：一.问题描述B1[1 2 34 5 67 8 9]B2[12 13 142131 4151 1 1811 1]两个十进制矩阵，行数不一样，分别是n1和n2，列数必须一致，为nwords，输出的矩阵Dh是[n1，n2]，这和求两句真的欧氏距离一样的。输出[1 1] = 1和12海明+2和13海明 + ...
阅读全文
posted @ 2015-06-21 16:06
加拿大小哥哥
摘要：一.问题描述 [1 1 1 0 1 0 1 1 0 1 0 0 1 1 0]每两位3转换为一个十进制数，共5列，那么转换后是ceil(5/3)=2列。 [7 1 6 1 1 2]二.问题分析 1.每八个字符划分为一个十进制，nwords=ceil(size(b,2)/8) 2.按...
阅读全文
posted @ 2015-06-21 10:18
加拿大小哥哥
摘要：这属于基础知识，老师应该讲的，可是；老师没讲.....在这个实验室，一师姐老师只要不在考试购物唱歌，完全无视其他人存在，给各个单身小学弟卖钱包，手表.......真是够了，精神污染....一.最小最大归一化 和区间映射（我理解的是把一个区间[a,b]映射到[c,d]，c+（x-a）*（d-c）...
阅读全文
posted @ 2015-06-20 20:45
加拿大小哥哥
摘要：一.问题来源 返回检索到的数据(按相关度排序)在原始数据中的索引。二.问题解析 x = [1 4 3 5; 1 3 2 6]；sortrows(x)其结果是按照row来排列，默认首先排第一列，1和1一样大，那么排第二列，3比4小，所以1 3 2 6应该在第一行。假如使用sortrows(x, 4...
阅读全文
posted @ 2015-06-04 11:28
加拿大小哥哥
摘要：% x是原数据集，分出训练样本和测试样本[ndata, D] = size(X); %ndata样本数，D维数R = randperm(ndata); %1到n这些数随机打乱得到的一个随机数字序列作为索引Xtest = X(R(1:num_test),:); %以索引...
阅读全文
posted @ 2015-05-29 09:32
加拿大小哥哥
摘要：一.问题来源 来自于一份大规模hash图像检索代码。二.问题解析2.1 添加目录addpath('./utils/');2.2 添加目录及其子目录addpath(genpath('./utils/.');
阅读全文
posted @ 2015-05-28 10:27
加拿大小哥哥
摘要：夜深人静时分，宿舍就我自己，只有蚊子陪伴着我，我慢慢码下这段文字............ 感觉知识结构不完善：上学期看论文，发现类间离散度矩阵和类内离散度矩阵，然后百度，找不到，现在学模式识别，见了，发现是数学公司，不过老师不讲........一.问题来源 今天有个意想不到的收货，下午讲课的...
阅读全文
posted @ 2015-05-25 00:37
加拿大小哥哥
摘要：一.问题来源 来自于一份PSO代码，PSO中需要初始化粒子位置和速度。二.问题探究 众所周知，Matlab中的rand()函数产生的是伪随机数，但一般用来也可以接受。但是，如果我们知道伪随机数的初始状态，那么产生的伪随机数是唯一确定的。问题来了，Matlab每次启动会重置rand()和randn...
阅读全文
posted @ 2015-05-23 12:53
加拿大小哥哥
摘要：matlab求极限（可用来验证度量函数或者隶属度函数）可用来验证是否收敛，取值范围等等。一.问题来源 搜集聚类资料时，又看到了隶属度函数，没错，就是下面这个，期间作者提到m趋于2是，结果趋于1，我想验证下，于是查资料。二.不同类型的极限2.1 基础知识 a./b表示常数a除以矩阵b中每个元素...
阅读全文
posted @ 2015-05-15 15:09
加拿大小哥哥
摘要：一.问题来源 来自于一份LSH代码，记录下来。二.函数解析2.1 bsxfun bsxfun是一个matlab自版本R2007a来就提供的一个函数，作用是”applies an element-by-element binary operation to arrays a and b, with...
阅读全文
posted @ 2015-04-15 17:59
加拿大小哥哥
摘要：算法伪代码如下，有错误的话请给我指出来：//假设所有a[i][i]都不是0，为0的话还不知道怎么办 输入系数于数组a[][]; 输入常数于数组b[]; //计算新系数于数组c[][] for(int i=0; ij) c[i][j] = -a[i][j]/a[i][i]; ...
阅读全文
posted @ 2014-12-30 16:01
加拿大小哥哥
摘要：今天看算法设计看到的>，想起组合数学老师经常用第二类，也没说为什么，这就记录下来了。 第一类:k=1时成立;假设k=n时成立,k=n+1时也成立.从而命题对任意n>1成立。 第二类:k=1时成立;假设k1成立。 第一类是高中学的,第二类在证明大学高等代数和初等数论问题用过。 数学归纳法只...
阅读全文
posted @ 2014-11-26 20:41
加拿大小哥哥
摘要：看到了北京一则新闻，想到如何测试双向镜子？百度之。 只要做以下简单的测试：把你的指甲尖放在镜子表面，如果在指甲尖与倒映图像之间有间隙，那就是真的镜子。然而，如果你的指甲尖能直接碰触到指甲反射影像，当心！那就是双向镜了！ 一、matlab里如何显示点的坐标？ 在图里显示点坐标 工具栏中有一个 "Data Cu...
阅读全文
posted @ 2014-11-01 14:31
加拿大小哥哥
浙公网安备 33010602011771号