以前写过一篇C++调用MATLAB引擎计算特征向量，里面讲了如何配置环境等等。现在又有一个新的需求，求解线性系统。又想起了MATLAB这个工具，于是又写了个小类。这里要求解的是AX=B。其中，A是m*n的矩阵，X是n维行向量，B是m维列向量。MATLAB里面求解很简单X=A\B。该类的头文件如下：//求解线性系统AX = B. //A是Row*Column的矩阵,B是Row*1的列向量 class CLinearSystem { public: CLinearSystem() {} CLinearSystem(double* pMatrixA, int nRow, int nColumn, double* pColumnB); ~CLinearSystem() {} public: void GetResult(double* pX);//pX保证能输出m_nColumn个元素 private: double* m_pMatrixA; double* m_pColumnB; int m_nRow; int m_nColumn; };源文件如下：#include "stdafx.h" #include "LinearSystem.h" #include#include "engine.h" #pragma comment(lib, " ...继续阅读 (129)

这是一个很老的话题了。所谓的反射机制，其实就是MFC里面的对象动态创建。比如说，给一个字符串代表类名，比如”CMesh”，根据这个字符串创建出一个CMesh对象。有不同的实现方法，在MFC里面是用一张巨大的链表网络把所有的类的信息链接起来，里面存储了类名和构造函数，动态创建的时候按照一定的顺序去里面搜索创建函数（构造函数）。与其用一张链表网，不如用一个字典保存起来，键是类名字符串，值是对应的构造函数。看起来很简单的样子，确实也很简单额。只是实现起来有些trick，比如字典一般是类的静态成员等。。。假设只有无参数的构造函数，首先定义上面所述的内容，类的字典。//头文件 typedef void* (*CreateFuntion)(void); class CClassFactory { public: CClassFactory() {} ~CClassFactory() {} public: static void* GetClassByName(std::string name) { std::map::const_iterator find; find = m_clsMap.find(name); if(find==m_clsMap.end()) { return NULL; ...继续阅读 (111)

最近一直在写毕业论文，过得乱七八糟的，也不知道自己干了些什么。今天是交论文的最后一天，终于发过去查重了。这一个月虽然乱七八糟发生了很多煎熬的事情，不过我貌似没有干太多学习上的事情。。。关于学习二字的定义，我觉得不应该再继续以前的理解了。以前总有点类似把自己当做容器，看书学习之类的都一股脑下来，全看完，到后面也不知道看了些什么，也忘记学了些什么。更重要的是，我为什么要去看这些书，我为什么要去学这些东西，如果以前是为了毕业，那么现在了，以后了。也许，我已经太久不记得我喜欢做的事情了。不知不觉，总是承担起来别人的喜好，别人的希望，直到最后步步艰辛。现在想来，我根本不需要这样活着。不说以前，说以后，以后一定要选择合适的，喜欢的事情去做，选择合适，喜欢的人去在一起，毕竟我现在有机会去选择了。觉得接下来这几个月，我该好好思考自己该干些什么事情了，以及我以后想成为什么样的人。非热血，非理想，非责任，只是我适合的，无压力的状态下，我适合成为什么样的人。我很清楚我不是学霸类型，也不是那种为了虚名或者牛逼而愿意逼自己做什么事情的人，选择成为那种人，我宁愿选择回家种田。某段时候，觉得自己就适合成为个家庭主男，以后上上班带带孩子就好，真心没有什么崇高的理想。到了这个年纪，我已经意识到我改变不了世界了。或者说，我已经不再有改变世界的冲动了。人生的美好，莫过于生病时候有人关怀，吐槽时候有人听着不介意，寂寞时候 ...继续阅读 (94)

我终于把界面换成qt的了，即使我从大二开始就用mfc，一直用到去年上半年。也终于到了，我实在受不了mfc代码混乱的时候了，qt那么方便的东西，我为什么不早点用了。这里就不吐槽mfc的混乱，qt的方便了。在mfc里面建立一个复杂的界面，比如说有dock窗口，或者tab页面的，代码多，而且混乱麻烦，更无语的是，每次我都必须去搜索下，解决各种各样的问题，有的还可以理解，有的就是莫名其妙了。更别说什么代码清晰，跨平台之类的，何止是天方夜谭。下面，说说我自己逐渐搭建的这个框架吧。其实，这个框架是从我上一个项目里，使用mfc的单文档单视图，主窗口里面有几个dockpane，视图里面放置了渲染OpenGL的子窗口，过渡来的。在qt的框架里面，mfc所有的这些恶心的东西都没有了，只剩下一个MainWindow类。在这个类里面，创建菜单，工具条，状态条，主窗口，dock窗口，tab窗口。代码简洁方便。渲染opengl的窗口设置为MainWindow的主widget。渲染窗口只是一个普通的widget而已，可以任意放置。这样的做法能够把界面和渲染窗口独立出来，渲染窗口对应一个独立的类(COglWidget)，所以方便创建不同类型的复杂界面。我觉得我应该先用viso画个类图，才能说清楚。下图是我的界面。中间是COglWidget。右侧是一个DockWidget，里面放置了个设置Mesh属性的widget ...继续阅读 (163)

上一篇文章说了提取拉普拉斯矩阵的第二个（或者其它）特征向量能得到一个有规律分布的标量场。现在讲的是如何在这样一个规律分布的标量场上提取出等值线。注意，这是在一个三角网格上面提取等值线，并不是真正的曲面。提取出等值线后，可以做一些相关的应用。那么，如何提取等值线了？首先，如果只是绘制出等值线，可以简单的设置纹理贴图值（比如间隔多少行，设置纹理值为白色），标量场对应一维纹理坐标。更详细的说明，可以参考我的一篇博文：在三维模型上可视化标量场及注意方面。现在考虑如何计算等值线上的间隔点。先看单个的三角面，无论如何可以得到一条等值线的。最简单的思想，就是把单个三角面上的等值线连接起来。给定一个三角面，如何计算从该三角面出发的一条等值线了？第一步，求出三个顶点的标量值顺序，那么标量范围在[fMin,fMax]之间。可以任选区间内的一个值，比如中间值作为该等值线的标量值。第二步，得到等值线标量值在[fMin,fMax]所在边pEdge的位置ptBeg，作为起始点，加入点集合isoLine。第三步，在pEdge所在的三角面上，找pEdge相邻的两条边的等值点(ptBeg的标量值相等的点)ptNext。如果存在这样的点，那么修改ptBeg = ptNext，加入点集合isoLine。同时，修改pEdge为对应的相邻边pEdge->next()或者是pEdge->next()->next()。如果不存 ...继续阅读 (132)

一个三维模型，可以看做是一个无向图。顶点是图的点，边可以作为图的边，那么可以得到一个邻接矩阵，该矩阵叫做拉普拉斯矩阵，计算这个邻接矩阵的特征向量就是所谓的谱分析。那么问题来了。如何定义边的权重？额，当然有不同的方法，看你需要分析的特性是什么。一般来说是边的两个对角的cotangle之和。另外，矩阵对角线元素是该行其它元素之和再取负。下一步就是计算该矩阵的特征向量了。如果计算了？我前面写过两篇文章讲述如何使用MATLAB计算稀疏矩阵的特征向量，现在我就是这么做的。首先，拉普拉斯矩阵是一个对称矩阵，而且是一个大型的稀疏矩阵。因此，只能使用稀疏矩阵来存储数据了。而且对称矩阵能保证一些特殊的性质，比如第一个特征向量是k*(1,1,1,1…)。第二个特征向量分布得比较有规律，但是通过实验发现，也许第二个也没有规律，可能是第三或者第四个，因为要利用特征向量的分布，需要综合多个特征向量。下面是我计算拉普拉斯矩阵的代码，需要考虑角度大于90的情况，以及除0等等。//a和b夹角的cot double CSEMesh::GetCotangent(double a, double b, double c) { if (a*b < 1e-8) { return 0; } else { double cosc = (a*a+b*b-c*c)/(2*a*b); ...继续阅读 (126)