摘要：
对于一般的指数曲线如：y=a*e^(k*t)，可以先对两边求对数得到：log(y) = log(a)+k*t 这样的曲线，然后用最小二乘来计算系数。 但是对于修正指数曲线如：y=k+a*b^t 这样的函数，没法直接求对数然后用最小二乘，因为有一个常数项k，这里可以利用三和法来计算系数。 对于曲线 y 阅读全文
摘要：
计算方法还是svd，这里把两个轴显示出来了。 过去用svd做过pca，就是在此基础上多一步运算，把原始数据旋转到主轴对应的坐标系上即可。 matlab代码如下： clear all; close all; clc; mu = [0 0]; sigma = [1 1.3;1.3 3]; data = 阅读全文
摘要：
解最小二乘的方法有很多，这里给出常见的三种方法实现。 一是一般方法，之前博客一般都用这种方法。 二是svd分解法，之前有用过svd(见这里，这里，这里和这里)解其他问题，但是没用来解过最小二乘。 三是qr分解法，这个好像没用过。 这里主要总结记录一下实现方法。 代码如下： clear all; cl 阅读全文
摘要：
如果点云文件比较大，可以利用octomap库对点云进行压缩。 octomap以八叉树结构对点云进行组织，是一种有损压缩，定义不同的分辨率可以保存不同尺度的数据。 这里以rabbit.pcd文件为例，分别将点云保存为文本格式(.txt)，二进制格式(.bin)，八叉树格式(.ot)，八叉树二进制模式( 阅读全文
摘要：
我这里是centos，一般默认都会有这个命令，ubuntu应该也差不多。 启动服务：service crond start 关闭服务：service crond stop 重启服务：service crond restart 重新载入配置：service crond reload 为了启动任务要先编 阅读全文
摘要：
油画效果算法有两个参数，一个是滤波半径r，一个是量化位数q。 算法原理： 1. 首先根据量化位数将0-255划分为q个区间； 2. 然后根据滤波半径r拿到半径内的像素； 3. 统计像素在区间中出现的次数，找到次数最大的区间； 4. 计算次数最大区间所有像素平均值赋值给当前滤波像素即可。 matlab 阅读全文
摘要：
Hausdorff距离是根据Hausdorff(1868-1942)命名的，Hausdorff距离是指某一集合中离另一集合最近点的所有距离的最大值 。 通常用如下公式表示： 需要注意的是h(A,B)和h(B,A)通常不相等，所以可以定义更一般的Hausdorff距离： matlab代码如下： cle 阅读全文
摘要：
今天突然发现一个很尴尬的问题，自己原来写的很多博客中图片都消失了，但是点击编辑后图片还在原来的位置。 我搜了一下原因，说是当时写博客时如果直接ctrl-c拷进文章中的图片现在就可能显示不了，点击“上传图片”按钮贴在博客中的能依然显示。 所以解决方法就是对不显示图片的博客重新编辑，把图片保存到本地后重 阅读全文
摘要：
习题： 代码如下： 1. Limit[(Exp[x] - Exp[-x])/Sin[x], x -> 0] Limit[(x/(x - 1)) - 1/Log[x], x -> 1] Limit[(1 - Cos[x])/x^2, x -> 0] Limit[(1 + 1/x^2)^x, x -> 阅读全文
摘要：
Thrust是cuda自带的c++库，cuda安装好之后，这个库也默认安装好了。 这个库基本是采用类似STL的接口方式，因此对于开发者非常友好，开发者不再需要关注内存与显存相关的问题了。 下面是一个简单的排序代码 main.cu： #include <thrust/host_vector.h> #i 阅读全文