随笔分类 - Matlab练习程序
摘要：计算步骤如下： 图片来自《视觉slam十四讲》6.2.2节。 下面使用书中的练习y=exp(a*x^2+b*x+c)+w这个模型验证一下，其中w为噪声，a、b、c为待解算系数。 代码如下： clear all; close all; clc; a=1;b=2;c=1; %待求解的系数 x=(0:0.
阅读全文
摘要：这个算是ICP算法中的一个关键步骤，单独拿出来看一下。 算法流程如下： 1.首先得到同名点集P和X。 2.计算P和X的均值up和ux。 3.由P和X构造协方差矩阵sigma。 4.由协方差矩阵sigma构造4*4对称矩阵Q。 5.计算Q的特征值与特征向量。其中Q最大特征值对应的特征向量即为最佳旋转向
阅读全文
摘要：还记得原来写过一个对图像进行波纹扭曲操作的博文。 这次实现的是水波特效，其实就是通过正余弦函数表示波纹中心位置慢慢向外扩散，通过叠加衰减因子使振幅不断减小，进而产生水波的效果。 效果如下： 原图： 波纹特效： matlab代码如下： 参考：https://blog.csdn.net/jia20003
阅读全文
摘要：最近在分析一些数据，就是数据拟合的一些事情，用到了matlab的polyfit函数，效果不错。 因此想了解一下这个多项式具体是如何拟合出来的，所以就搜了相关资料。 这个文档介绍的还不错，我估计任何一本数值分析教材上讲的都非常清楚。 推导就不再写了，我主要参考下面两页PPT，公式和例子讲的比较清楚。
阅读全文
摘要：问题是这样：比如有一个地心惯性系的轨道，然后从轨道上取了几个点，问能不能根据这几个点把轨道还原了？ 当然，如果知道轨道这几个点的速度的情况下，根据轨道六根数也是能计算轨道的，不过真近点角是随时间变动的。 下面我会用数学的方法来解这个问题，基本思想是通过拟合空间上点的平面与椭球平面的交线将该轨道计算出
阅读全文
摘要：需要看下生成的数据在地球上的经纬度具体位置。 投影为墨卡托投影。
阅读全文
摘要：记得上次练习了神经网络分类，不过当时应该有些地方写的还是不对。 这次用神经网络识别mnist手写数据集，主要参考了深度学习工具包的一些代码。 mnist数据集训练数据一共有28*28*60000个像素，标签有60000个。 测试数据一共有28*28*10000个，标签10000个。 这里神经网络输入
阅读全文
摘要：这里我用的空间是x向右为正，y向下为正，z向屏幕里面为正。相当于标准右手系绕x轴旋转了180度。 将三个点光源放在 r = [0.3,0,0.5];g = [0.3,-0.5*cos(pi/6),-0.5*sin(pi/6)];b = [0.3,0.5*cos(pi/6),-0.5*sin(pi/6
阅读全文
摘要：注：这里的练习鉴于当时理解不完全，可能会有些错误，关于神经网络的实践可以参考我的这篇博文 这里的代码只是简单的练习，不涉及代码优化，也不涉及神经网络优化，所以我用了最能体现原理的方式来写的代码。 激活函数用的是h = 1/(1+exp(-y))，其中y=sum([X Y].*w)。 代价函数用的是E
阅读全文
摘要：对已标记数据分类： 对测试数据分类：
阅读全文
摘要：对已标记数据分类： 对测试数据分类：
阅读全文
摘要：这个算法是对photoshop中滤镜->像素化->碎片这个功能的学习。PS这个功能好像不带参数，不过我这里有滤波半径r可以进行控制。因为我是看效果猜算法的，效果肯定有所区别。我的想法是对图像以r为长度，在上下左右四个方向做位移，然后进行混合。代码还是很简单，就不详细介绍了。处理效果如下：原图：处理后...
阅读全文
摘要：我手机上有一个软件实现了很多图像滤镜，挺有意思，我打算都尝试一下。这个滤镜主要是实现老照片效果。代码很短，我就不详细介绍了。原图：处理后效果：matlab代码如下：clear all;close all;clc;img=imread('lena_rgb.jpg');[h w k]=size(img)...
阅读全文
摘要：这个算子算是图像历史上第一个特征点提取算法了，1977年提出的，很简单，拿来练手很合适。算法原理如下：1.选取一个合理的邻域遍历图像，这里是5*5邻域的。在邻域中依次计算，垂直，水平，对角与反对角四个相邻像素灰度的差的平方和，作为该邻域特征值。大致就是下面这个样子：公式：这里k是窗口的半径。2.从四...
阅读全文
摘要：还是计算几何， 多边形的核可以这样理解：这个核为原多边形内部的一个多边形，站在这个叫核的多边形中，我们能看到原多边形的任何一个位置。算法步骤如下：1.根据原多边形最大和最小的x,y初始化核多边形，就是个矩形。2.计算多边形当前处理的点的凹凸性。3.用当前点与其后继点构成直线，判断当前点的前驱点在该直...
阅读全文
摘要：又是计算几何，我感觉最近对计算几何上瘾了。当然，工作上也会用一些，不过工作上一般直接调用boost的geometry库。上次写过最小包围圆，这次是最小包围矩形，要比最小包围圆复杂些。最小包围矩形可不一定是个直立的矩形，也可能像下图一样是倾斜的。求法如下：1.求多边形凸包，这里凸包直接调用系统函数了，...
阅读全文
摘要：本练习程序是受到了这个老外博文的启发，感觉挺有意思，就尝试了一下。他用的是opencv，我这里用的是matlab。过去写过透视投影，当时是用来做倾斜校正的，这次同样用到了透视投影，不过更有意思，是将一张图像贴到另一张图像上。两个透视投影都需要先计算投影矩阵，倾斜校正那一篇是通过解线性方程组求的变换矩...
阅读全文
摘要：依然是计算几何。射线法判断点与多边形关系原理如下：从待判断点引出一条射线，射线与多边形相交，如果交点为偶数，则点不在多边形内，如果交点为奇数，则点在多边形内。原理虽是这样，有些细节还是要注意一下，比如射线过多边形顶点或射线与多边形其中一边重合等情况还需特别判断。这里就不特别判断了，因为我只是熟悉原理...
阅读全文
摘要：最近总是对计算几何方面的程序比较感兴趣。多圆求交点，要先对圆两两求交点。有交点的圆分为相切圆和相交圆。相切圆求法： 1.根据两圆心求直线 2.求公共弦直线方程 3.求两直线交点即两圆切点。相交圆求法： 1.求公共弦方程直线。 2.公共弦直线方程和其中一个圆方程联立求解即可。公共弦直线方程就...
阅读全文
摘要：处理原理就是将邻域像素平均值赋给邻域中的所有像素。这里将处理方法分为了带线条和不带线条，带线条处理后的视觉效果要好些。原图：处理结果（不带线条）：处理结果（带线条）：代码处理的是带线条的情况。matlab代码如下：clear all; close all;clc;img=imread('lena.j...
阅读全文