R君的博客

一、高斯法和Doolitle法 高斯法和Doolitle法都比较简单，其实Doolitle可以直接用maltab里面的lu命令来求，下面重点讲一下迭代法 二、Jacobi迭代 1.算法实现 Jacobi迭代的伪代码如下： 2.收敛条件 3.matlab实现 根据以上的算法可以写出一个简单的jacobi迭代如下： function x = jacobi(a,b,k) n=length(b); d=diag(a); r=a-diag(d); x=zeros(n,1); for j = 1:k x = (b

kmeans聚类以及fcm聚类的一大局限是需要提前知道大致的聚类个数，现实中使用比较受限，2007年，Frey和Dueck在Science发表了《Clustering by Passing Messages Between Data Points》，提出了AP聚类的方法，此方法采用点点之间交换信息的方式聚类，下面对此算法进行分析并且用matlab实现之。 1.算法原理 首先我们定义三个矩阵: 1.相似度矩阵s(i,k)：表示数据点k作为数据点i的聚类中心的能力（相当于k是参选人，i是选民），可以用欧氏距离来