計(jì)算PageRank

有A,B,C,D,E五個(gè)網(wǎng)頁(yè)鼻吮,其中
1)A網(wǎng)頁(yè)有鏈接指向B,C,D,E
2)B網(wǎng)頁(yè)有鏈接指向A,D
3)C網(wǎng)頁(yè)有鏈接指向A,D
4)D網(wǎng)頁(yè)有鏈接指向C
5)E網(wǎng)頁(yè)有鏈接指向A,C

思想:
1娄柳、計(jì)算矩陣S
2聘殖、計(jì)算矩陣G
3蛤高、根據(jù)q=Gq',迭代

public class PageRank {
    private static final double ALPHA = 0.85;
    private static final double DISTANCE = 0.0000001;

    private static final double[][] matrixS = {
            {0, 1 / 2.0, 1 / 2.0, 0, 1 / 2.0},
            {1 / 4.0, 0, 0, 0, 0},
            {1 / 4.0, 0, 0, 1, 1 / 2.0},
            {1 / 4.0, 1 / 2.0, 1 / 2.0, 0, 0},
            {1 / 4.0, 0, 0, 0, 0}
    };

    private static final double[][] matrixU = {
            {1.0, 1.0, 1.0, 1.0, 1.0},
            {1.0, 1.0, 1.0, 1.0, 1.0},
            {1.0, 1.0, 1.0, 1.0, 1.0},
            {1.0, 1.0, 1.0, 1.0, 1.0},
            {1.0, 1.0, 1.0, 1.0, 1.0}
    };

    public static void main(String[] args) {

        System.out.println("S矩陣:");
        printMatrix(matrixS);

        double[][] matrixG = calMatrixG();

        System.out.println("G矩陣:");
        printMatrix(matrixG);

        double[] vector = {1.0, 1.0, 1.0, 1.0, 1.0};

        double[] vectorQ = calVectorQ(matrixG, vector);

        System.out.println("向量:");
        printVector(vectorQ);
    }

    private static void printMatrix(double[][] matrix) {
        for (int i = 0; i < matrix.length; i++) {
            for (int j = 0; j < matrix[i].length; j++) {

//                String result = String.format("%1$10s", matrix[i][j]);
                String result = String.format("%.5f", matrix[i][j]);

                System.out.print(result + " ");
            }
            System.out.println();
        }
    }

    private static void printVector(double[] vector) {
        for (int i = 0; i < vector.length; i++) {

            String result = String.format("%.5f", vector[i]);

            System.out.print(result + " ");
        }
    }

    private static double[][] calMatrixG() {
        return calMatrixPlus(calMatrixMultiDouble(matrixS, ALPHA),
                calMatrixMultiDouble(matrixU, (1 - ALPHA) * (1.0 / matrixS.length)));
    }

    private static double[][] calMatrixMultiDouble(double[][] matrix, double value) {
        for (int i = 0; i < matrix.length; i++) {
            for (int j = 0; j < matrix[i].length; j++) {
                matrix[i][j] *= value;
            }

        }
        return matrix;
    }

    private static double[][] calMatrixPlus(double[][] matrix1, double[][] matrix2) {
        for (int i = 0; i < matrix1.length; i++) {
            for (int j = 0; j < matrix1[i].length; j++) {
                matrix1[i][j] += matrix2[i][j];
            }

        }
        return matrix1;
    }

    private static double[] calMatrixMultiVetor(double[][] matrix, double[] vetor) {
        double[] tmp = new double[matrix.length];

        for (int i = 0; i < matrix.length; i++) {
            for (int j = 0; j < matrix[i].length; j++) {
                tmp[i] += matrix[i][j] * vetor[j];
            }
        }

        return tmp;
    }

    private static double calVectorDistance(double[] vector1, double[] vector2) {
        double sum = 0;

        for (int i = 0; i < vector1.length; i++) {
            sum += Math.pow(vector1[i] - vector2[i], 2);
        }
        return Math.sqrt(sum);
    }

    private static double[] calVectorQ(double[][] matrixG, double[] vector) {

        double[] vector2 = calMatrixMultiVetor(matrixG, vector);

        double distance = calVectorDistance(vector, vector2);

        int n = 1;

//        System.out.print(distance);

        vector = vector2;
        while (distance > DISTANCE) {
            n++;

            vector2 = calMatrixMultiVetor(matrixG, vector);

            distance = calVectorDistance(vector, vector2);

            vector = vector2;

//            System.out.println(distance);
        }

        System.out.println("計(jì)算次數(shù)" + n);

        return vector2;
    }
}

輸出
S矩陣:
0.00000 0.50000 0.50000 0.00000 0.50000
0.25000 0.00000 0.00000 0.00000 0.00000
0.25000 0.00000 0.00000 1.00000 0.50000
0.25000 0.50000 0.50000 0.00000 0.00000
0.25000 0.00000 0.00000 0.00000 0.00000
G矩陣:
0.03000 0.45500 0.45500 0.03000 0.45500
0.24250 0.03000 0.03000 0.03000 0.03000
0.24250 0.03000 0.03000 0.88000 0.45500
0.24250 0.45500 0.45500 0.03000 0.03000
0.24250 0.03000 0.03000 0.03000 0.03000
計(jì)算次數(shù)29
向量:
1.21040 0.40721 1.68064 1.29454 0.40721

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市诚撵,隨后出現(xiàn)的幾起案子材诽,更是在濱河造成了極大的恐慌,老刑警劉巖汗捡,帶你破解...
    沈念sama閱讀 207,113評(píng)論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件兔乞,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡,警方通過(guò)查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,644評(píng)論 2 381
  • 文/潘曉璐 我一進(jìn)店門购桑,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái),“玉大人淡溯,你說(shuō)我怎么就攤上這事〔狙担” “怎么了咱娶?”我有些...
    開封第一講書人閱讀 153,340評(píng)論 0 344
  • 文/不壞的土叔 我叫張陵米间,是天一觀的道長(zhǎng)。 經(jīng)常有香客問(wèn)我膘侮,道長(zhǎng)屈糊,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 55,449評(píng)論 1 279
  • 正文 為了忘掉前任琼了,我火速辦了婚禮逻锐,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘雕薪。我一直安慰自己昧诱,他們只是感情好,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,445評(píng)論 5 374
  • 文/花漫 我一把揭開白布所袁。 她就那樣靜靜地躺著盏档,像睡著了一般。 火紅的嫁衣襯著肌膚如雪燥爷。 梳的紋絲不亂的頭發(fā)上蜈亩,一...
    開封第一講書人閱讀 49,166評(píng)論 1 284
  • 那天,我揣著相機(jī)與錄音前翎,去河邊找鬼稚配。 笑死,一個(gè)胖子當(dāng)著我的面吹牛鱼填,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播毅戈,決...
    沈念sama閱讀 38,442評(píng)論 3 401
  • 文/蒼蘭香墨 我猛地睜開眼苹丸,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼!你這毒婦竟也來(lái)了苇经?” 一聲冷哼從身側(cè)響起赘理,我...
    開封第一講書人閱讀 37,105評(píng)論 0 261
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤,失蹤者是張志新(化名)和其女友劉穎扇单,沒(méi)想到半個(gè)月后商模,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,601評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡蜘澜,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,066評(píng)論 2 325
  • 正文 我和宋清朗相戀三年施流,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片鄙信。...
    茶點(diǎn)故事閱讀 38,161評(píng)論 1 334
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡瞪醋,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出装诡,到底是詐尸還是另有隱情银受,我是刑警寧澤践盼,帶...
    沈念sama閱讀 33,792評(píng)論 4 323
  • 正文 年R本政府宣布,位于F島的核電站宾巍,受9級(jí)特大地震影響咕幻,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜顶霞,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,351評(píng)論 3 307
  • 文/蒙蒙 一肄程、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧确丢,春花似錦绷耍、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,352評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)。三九已至描函,卻和暖如春崎苗,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背舀寓。 一陣腳步聲響...
    開封第一講書人閱讀 31,584評(píng)論 1 261
  • 我被黑心中介騙來(lái)泰國(guó)打工胆数, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人互墓。 一個(gè)月前我還...
    沈念sama閱讀 45,618評(píng)論 2 355
  • 正文 我出身青樓必尼,卻偏偏與公主長(zhǎng)得像,于是被迫代替她去往敵國(guó)和親篡撵。 傳聞我的和親對(duì)象是個(gè)殘疾皇子判莉,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,916評(píng)論 2 344

推薦閱讀更多精彩內(nèi)容

  • 在C語(yǔ)言中,五種基本數(shù)據(jù)類型存儲(chǔ)空間長(zhǎng)度的排列順序是: A)char B)char=int<=float C)ch...
    夏天再來(lái)閱讀 3,325評(píng)論 0 2
  • 專業(yè)考題類型管理運(yùn)行工作負(fù)責(zé)人一般作業(yè)考題內(nèi)容選項(xiàng)A選項(xiàng)B選項(xiàng)C選項(xiàng)D選項(xiàng)E選項(xiàng)F正確答案 變電單選GYSZ本規(guī)程...
    小白兔去釣魚閱讀 8,977評(píng)論 0 13
  • 原文鏈接:Search and Track Scheduling for Multifunction Phased...
    Chungfai_Yi閱讀 5,061評(píng)論 0 1
  • Linux基礎(chǔ)命令學(xué)習(xí)總結(jié),隨著學(xué)習(xí)的前進(jìn)育谬,將不斷完善券盅。網(wǎng)友們,有錯(cuò)誤之處膛檀,請(qǐng)給我留言改正喔锰镀,謝謝! 基礎(chǔ)操作 <...
    zhuenbiao閱讀 1,729評(píng)論 0 2
  • 一咖刃,HTML語(yǔ)言的一般語(yǔ)法: 1泳炉,圍堵標(biāo)記:<>… 1)帶屬性的標(biāo)記: … 2)無(wú)屬性的標(biāo)記:加粗 居中 標(biāo)題 2...
    清水易藍(lán)閱讀 1,259評(píng)論 0 2