PageRank迭代求解方法

2023-12-06 13:32
文章标签 方法 迭代 求解 pagerank

本文主要是介绍PageRank迭代求解方法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

PageRank算法迭代求解第二版程序

/*
描述: 小型Web超链图的PageRank算法迭代求解
作者: xiaocui
时间: 2008.4.19
版本: v1.0
*/

/*简单的PageRank算法为 PR(T) = PR(T1)/C(T1)+...+PR(Tn)/C(Tn)),
最后写成了 P = M * P, P为n个网页的pagerank值组成的列向量,
M为马尔可夫转移矩阵,其每列元素表示该网页跳转到其他网页的概率。
PageRank算法的求解就是求M矩阵特征值为1的特征向量,1是M的最大特
征值,可以通过幂法进行迭代求解。

本程序是模拟PageRank算法的迭代计算,只用了7个网页,其链接关系图用邻接矩阵
表示为:
      1    2  3  4  5  6  7
  1  0    1  1  1  1  0 1

  2  1    0  0  0  0  0  0

  3  1    1  0  0  0  0  0

  4  0    1  1  0  1  0  0

  5  1    0  1  1  0  1  0

  6  1    0  0  0  1  0  0

  7  0    0  0  0  1  0  0

转置归一化后对应的马尔可夫转移矩阵为:
      1     2    3      4       5       6      7
  1  0     1    0.5    0       0.25    0.5    0

  2  0.2  0    0.5    1/3     0       0      0

  3  0.2  0    0      1/3     0.25    0      0

  4  0.2  0    0      0       0.25    0      0

  5  0.2  0    0      1/3     0       0.5    1

  6  0     0    0      0       0.25    0      0

  7  0.2  0    0      0       0       0      0
*/
#include <iostream>
#include <vector>
#include <cmath>
using namespace std;

/* 计算矩阵和向量的乘积
输入:arr表示矩阵,vec表示向量
输出:结果向量
*/
vector<double> arrayVecMul(const vector<vector<double> >& arr, const vector<double>& vec)
{
vector<double> rstVec; //结果向量
double sum = 0.0;
for (size_t i = 0; i < arr.size(); ++i)
{
   sum = 0.0;
   for (size_t j = 0; j < arr[i].size(); ++j)
   {
    sum += arr[i][j] * vec[j];
   }
   rstVec.push_back(sum);
}

return rstVec;
}

/* 向量归1化
输入:vec表示向量
输出:归1化的向量vec
*/
void toOne(vector<double>& vec)
{
double sum = 0.0;
for (size_t i = 0; i < vec.size(); ++i)
{
   sum += vec[i];
}
for (size_t i = 0; i < vec.size(); ++i)
{
   vec[i] /= sum;
}
}

/* 返回向量的最大元素(绝对值最大)
输入:vec为向量
输出:最大元素值
*/
double getMaxElement(const vector<double>& vec)
{
double max = fabs(vec[0]);
for (size_t i = 0; i < vec.size(); ++i)
{
   double tmp = fabs(vec[i]);
   if ( tmp > max )
   {
    max = tmp;
   }
}

return max;
}

/* 判断误差是否满足,误差满足后迭代结束
输入:向量vecPrior表示前一次的结果,向量vec表示本次结果,e表示误差界
输出:true表示误差满足,迭代结束;false表示误差不满足
*/
bool isFinish(const vector<double>& vecPrior, const vector<double>& vec, double e)
{
double maxE = 0.0; //表示前后两个向量对应分量的最大误差
for (size_t i = 0; i < vecPrior.size(); ++i)
{
   double tmp = fabs(vec[i] - vecPrior[i]);
   if ( tmp > maxE )
   {
    maxE = tmp;
   }
}
if ( maxE < e )
{
   return true;
}

return false;
}

int main()
{
//构造马尔可夫转移矩阵
vector<vector<double> > arr;
double a1[7] = {0, 1, 0.5, 0, 0.25, 0.5, 0};
double a2[7] = {0.2, 0, 0.5, 1.0/3, 0, 0, 0};
double a3[7] = {0.2, 0, 0, 1.0/3, 0.25, 0, 0};
double a4[7] = {0.2, 0, 0, 0, 0.25, 0, 0};
double a5[7] = {0.2, 0, 0, 1.0/3, 0, 0.5, 1};
double a6[7] = {0, 0, 0, 0, 0.25, 0, 0};
double a7[7] = {0.2, 0, 0, 0, 0, 0, 0};
vector<double> row1(a1, a1 + 7);
vector<double> row2(a2, a2 + 7);
     vector<double> row3(a3, a3 + 7);
vector<double> row4(a4, a4 + 7);
vector<double> row5(a5, a5 + 7);
vector<double> row6(a6, a6 + 7);
vector<double> row7(a7, a7 + 7);
arr.push_back(row1);
arr.push_back(row2);
arr.push_back(row3);
arr.push_back(row4);
arr.push_back(row5);
arr.push_back(row6);
arr.push_back(row7);

double f[7] = {0, 0.8, 1, 0.3, 0.2, 0.7, 0.4};
vector<double> vec(f, f + 7); //初始向量
double e = 0.000001; //误差界
int count = 0; //迭代计算次数
vector<double> rstVec; //特征向量
double maxValue; //最大特征值

while ( 1 )
{
   count++;
   rstVec = arrayVecMul(arr, vec);
   maxValue = getMaxElement(rstVec);
   for (size_t i = 0; i < rstVec.size(); ++i)
   {
    vec[i] = rstVec[i] / maxValue;
   }
  
   if ( isFinish(vec, rstVec, e) )
   {
    break;
   }
   else
   {
    maxValue = getMaxElement(rstVec);
    for (size_t i = 0; i < rstVec.size(); ++i)
    {
     vec[i] = rstVec[i] / maxValue;
    }
   }
}

cout << "总共迭代了: " << count << "次" << endl;
cout << "最大特征值为: " << maxValue << endl;
for (size_t i = 0; i < rstVec.size(); ++i)
{
   cout << rstVec[i] << "  ";
}
cout << endl;

cout << "正确性检验: " << endl;
cout << "矩阵和特征向量的乘积: " << endl;
rstVec = arrayVecMul(arr, rstVec);
for (size_t i = 0; i < rstVec.size(); ++i)
{
   cout << rstVec[i] << "  ";
}
cout << endl;

cout << "特征值和特征向量的乘积: " << endl;
for (size_t i = 0; i < rstVec.size(); ++i)
{
   rstVec[i] = maxValue * rstVec[i];
}

return 0;
}

这篇关于PageRank迭代求解方法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/461999

相关文章

Python安装Pandas库的两种方法

《Python安装Pandas库的两种方法》本文介绍了三种安装PythonPandas库的方法,通过cmd命令行安装并解决版本冲突,手动下载whl文件安装,更换国内镜像源加速下载,最后建议用pipli... 目录方法一:cmd命令行执行pip install pandas方法二:找到pandas下载库,然后

Linux系统中查询JDK安装目录的几种常用方法

《Linux系统中查询JDK安装目录的几种常用方法》:本文主要介绍Linux系统中查询JDK安装目录的几种常用方法,方法分别是通过update-alternatives、Java命令、环境变量及目... 目录方法 1:通过update-alternatives查询(推荐)方法 2:检查所有已安装的 JDK方

SQL Server安装时候没有中文选项的解决方法

《SQLServer安装时候没有中文选项的解决方法》用户安装SQLServer时界面全英文,无中文选项,通过修改安装设置中的国家或地区为中文中国,重启安装程序后界面恢复中文,解决了问题,对SQLSe... 你是不是在安装SQL Server时候发现安装界面和别人不同,并且无论如何都没有中文选项?这个问题也

Java Thread中join方法使用举例详解

《JavaThread中join方法使用举例详解》JavaThread中join()方法主要是让调用改方法的thread完成run方法里面的东西后,在执行join()方法后面的代码,这篇文章主要介绍... 目录前言1.join()方法的定义和作用2.join()方法的三个重载版本3.join()方法的工作原

在MySQL中实现冷热数据分离的方法及使用场景底层原理解析

《在MySQL中实现冷热数据分离的方法及使用场景底层原理解析》MySQL冷热数据分离通过分表/分区策略、数据归档和索引优化,将频繁访问的热数据与冷数据分开存储,提升查询效率并降低存储成本,适用于高并发... 目录实现冷热数据分离1. 分表策略2. 使用分区表3. 数据归档与迁移在mysql中实现冷热数据分

Spring Boot从main方法到内嵌Tomcat的全过程(自动化流程)

《SpringBoot从main方法到内嵌Tomcat的全过程(自动化流程)》SpringBoot启动始于main方法,创建SpringApplication实例,初始化上下文,准备环境,刷新容器并... 目录1. 入口:main方法2. SpringApplication初始化2.1 构造阶段3. 运行阶

Olingo分析和实践之ODataImpl详细分析(重要方法详解)

《Olingo分析和实践之ODataImpl详细分析(重要方法详解)》ODataImpl.java是ApacheOlingoOData框架的核心工厂类,负责创建序列化器、反序列化器和处理器等组件,... 目录概述主要职责类结构与继承关系核心功能分析1. 序列化器管理2. 反序列化器管理3. 处理器管理重要方

Python错误AttributeError: 'NoneType' object has no attribute问题的彻底解决方法

《Python错误AttributeError:NoneTypeobjecthasnoattribute问题的彻底解决方法》在Python项目开发和调试过程中,经常会碰到这样一个异常信息... 目录问题背景与概述错误解读:AttributeError: 'NoneType' object has no at

postgresql使用UUID函数的方法

《postgresql使用UUID函数的方法》本文给大家介绍postgresql使用UUID函数的方法,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录PostgreSQL有两种生成uuid的方法。可以先通过sql查看是否已安装扩展函数,和可以安装的扩展函数

Java中Arrays类和Collections类常用方法示例详解

《Java中Arrays类和Collections类常用方法示例详解》本文总结了Java中Arrays和Collections类的常用方法,涵盖数组填充、排序、搜索、复制、列表转换等操作,帮助开发者高... 目录Arrays.fill()相关用法Arrays.toString()Arrays.sort()A