<rt id="bn8ez"></rt>
<label id="bn8ez"></label>

  • <span id="bn8ez"></span>

    <label id="bn8ez"><meter id="bn8ez"></meter></label>

    Change Dir

    先知cd——熱愛生活是一切藝術的開始

    統計

    留言簿(18)

    積分與排名

    “牛”們的博客

    各個公司技術

    我的鏈接

    淘寶技術

    閱讀排行榜

    評論排行榜

    Commons Math學習筆記——矩陣

     

    1.2 矩陣

    看其他篇章到目錄選擇。

    今天來第二篇:矩陣——Matrix

    Mathorg.apache.commons.math.linear里對矩陣的表示是有一個層次結構的。

    最頂層的AnyMatrix是一個基本的interface。下面有3sub interfaceBigMatrix, FieldMatrix<T>, RealMatrix。而每個sub interface分別被相應的矩陣類實現。整個矩陣的層次結構也就出來了。不過其中的BigMatrix已經不用了。被Array2DRowFieldMatrix替代了。


     

    具體拿RealMatrix來說。

    RealMatrix是一個可以表示實數類型數據的矩陣接口,實現RealMatrix接口的類有AbstractRealMatrix。它的子類有Array2DRowRealMatrix, BlockRealMatrix, OpenMapRealMatrix, RealMatrixImpl

     

    今天以Array2DRowRealMatrix為例研究一下矩陣Matrix都有哪些操作。

    我認為通常最簡單的學習方式就是直接看代碼樣例,為此,我寫了使用Matrix的樣例代碼,通過代碼的演示來看看如何使用Array2DRowRealMatrix

    Array2DRowRealMatrix的內部實現是一個2double類型的數組double [][]data;。它的getData方法可以返回對應的數組表示。LU decomposition用來進行矩陣的分解及相關操作,Array2DRowRealMatrix中有幾個方法已經過時,均被LU decomposition取代。像矩陣的求逆運算、特征值以及奇異性等。矩陣的一些分解運算,我想還是放到后面再研究吧。代碼一次太多也不好,嘻嘻。

    具體的代碼實現里有注釋和輸出提示,我想這樣的代碼大家運行后基本就能完全理解Array2DRowRealMatrix的所有操作了。

     1/**
     2 * 
     3 */

     4package algorithm.math;
     5
     6import org.apache.commons.math.linear.Array2DRowRealMatrix;
     7
     8/**
     9 * @author Jia Yu
    10 * @date 2010-11-18
    11 */

    12public class MatrixTest {
    13
    14    public static void matrix() {
    15        double[][] data1 = { 1d, 2d, 3d }{ 2d, 5d, 3d }{ 1d, 0d, 8d } };
    16        double[][] t_data = -40d, 16d, 9d }{ 13d, -5d, -3d },
    17                { 5d, -2d, -1d } }
    ;
    18
    19        Array2DRowRealMatrix matrix1 = new Array2DRowRealMatrix(data1);
    20        Array2DRowRealMatrix t_mat = new Array2DRowRealMatrix(t_data);
    21        // output directly
    22        System.out.println("matrix is " + matrix1);
    23        // is square
    24        System.out.println("it is square matrix! : " + matrix1.isSquare());
    25        // dimension of row and column
    26        System.out.println("row dimension is " + matrix1.getRowDimension());
    27        System.out.println("column dimension is "
    28                + matrix1.getColumnDimension());
    29        // matrix add
    30        System.out.println("mat1 + mat1 = " + matrix1.add(matrix1));
    31        System.out.println("mat1 + 5 = " + matrix1.scalarAdd(5.0));
    32        // matrix sub
    33        System.out.println("mat1 - mat1 = " + matrix1.subtract(matrix1));
    34        // matrix norm
    35        System.out.println("the maximum absolute row sum norm is "
    36                + matrix1.getNorm());
    37        // matrix multiply
    38        System.out.println("mat1 * t_mat = " + matrix1.multiply(t_mat));
    39        System.out.println("mat1 * 5.0 = " + matrix1.scalarMultiply(5));
    40        System.out.println("t_mat * mat1 = " + matrix1.preMultiply(t_mat));
    41        // matrix trace
    42        System.out.println("the trace is " + matrix1.getTrace());
    43        // matrix transpose
    44        System.out.println("the transpose of mat1 is " + matrix1.transpose());
    45        // matrix to vector
    46        System.out
    47                .println("the first row vector is " + matrix1.getRowVector(0));
    48        // matrix get sub matrix of selected rows and columns
    49        System.out.println("sub matrix of mat1 is "
    50                + matrix1.getSubMatrix(new int[] 02 }new int[] 12 }));
    51    }

    52
    53    /**
    54     * @param args
    55     */

    56    public static void main(String[] args) {
    57        // TODO Auto-generated method stub
    58        matrix();
    59    }

    60}

    61

     運行結果:

    matrix is Array2DRowRealMatrix{{1.0,2.0,3.0},{2.0,5.0,3.0},{1.0,0.0,8.0}}
    it is square matrix! : true
    row dimension is 3
    column dimension is 3
    mat1 + mat1 = Array2DRowRealMatrix{{2.0,4.0,6.0},{4.0,10.0,6.0},{2.0,0.0,16.0}}
    mat1 + 5 = Array2DRowRealMatrix{{6.0,7.0,8.0},{7.0,10.0,8.0},{6.0,5.0,13.0}}
    mat1 - mat1 = Array2DRowRealMatrix{{0.0,0.0,0.0},{0.0,0.0,0.0},{0.0,0.0,0.0}}
    the maximum absolute row sum norm is 14.0
    mat1 * t_mat = Array2DRowRealMatrix{{1.0,0.0,0.0},{0.0,1.0,0.0},{0.0,0.0,1.0}}
    mat1 * 5.0 = Array2DRowRealMatrix{{5.0,10.0,15.0},{10.0,25.0,15.0},{5.0,0.0,40.0}}
    t_mat * mat1 = Array2DRowRealMatrix{{1.0,0.0,0.0},{0.0,1.0,0.0},{0.0,0.0,1.0}}
    the trace is 14.0
    the transpose of mat1 is Array2DRowRealMatrix{{1.0,2.0,1.0},{2.0,5.0,0.0},{3.0,3.0,8.0}}
    the first row vector is {1; 2; 3}
    sub matrix of mat1 is Array2DRowRealMatrix{{2.0,3.0},{0.0,8.0}}

    AnyMatrix接口只定義了最基本的操作,獲取維度和判斷方陣。

    RealMatrix接口擴展了AnyMatrix,定義了一些操作,比如加減乘等。

    AbstractRealMatrix抽象類實現了RealMatrix,定義了更多的get方法,可以獲得更多矩陣相關的參數,比如矩陣的秩、矩陣的跡、矩陣的特征值和矩陣的轉置等。

    Array2DRowRealMatrix繼承了AbstractRealMatrix,將里面的抽象方法全部實現。

    其中的multiply矩陣乘法運算,multiply返回的是this*m;而preMultiply返回的是m*this

    另外值得一提的是,linear包里有MatrixUtils類提供了一系列靜態方法用來檢測矩陣,其中的方法主要以check***為主,檢測矩陣的合法性。

    相關資料:

    矩陣知識:http://zh.wikipedia.org/zh/%E7%9F%A9%E9%98%B5

    Commons math包:http://commons.apache.org/math/index.html

    posted on 2010-12-11 21:12 changedi 閱讀(4093) 評論(3)  編輯  收藏 所屬分類: 數學

    評論

    # re: Commons Math學習筆記——矩陣 2010-12-20 15:46 杜良永

    這個有求逆矩陣的功能嗎?  回復  更多評論   

    # re: Commons Math學習筆記——矩陣 2010-12-20 16:30 changedi

    @杜良永
    有的,詳見《矩陣分解》那一節的例子,代碼中的
    getInverse方法可以求得矩陣的逆,基于LU分解做的~~  回復  更多評論   

    # re: Commons Math學習筆記——矩陣 2011-04-07 10:35 xautchap

    有的,詳見《矩陣分解》那一節的例子,代碼中的
    getInverse方法可以求得矩陣的逆,基于LU分解做的~~

    對于速度來說,LU分解的速度沒有QR分解的效率好
    @changedi
      回復  更多評論   

    主站蜘蛛池模板: 夜夜亚洲天天久久| 亚洲精品天堂成人片?V在线播放| 亚洲精品久久久久无码AV片软件| 99视频精品全部免费观看| 日韩视频在线精品视频免费观看| 久久青青成人亚洲精品| 亚洲精品欧美综合四区| 四虎影视大全免费入口| MM131亚洲国产美女久久| 亚洲va成无码人在线观看| 亚洲AV无码一区二区三区牲色| 一级做a爰片久久毛片免费看 | 少妇亚洲免费精品| 视频一区在线免费观看| 亚洲视频免费一区| 亚洲一区二区三区免费在线观看 | 韩日电影在线播放免费版| 国产亚洲3p无码一区二区| 99久久久国产精品免费牛牛四川| 亚洲午夜精品国产电影在线观看| 中国在线观看免费的www| 亚洲成AV人在线播放无码| 永久黄色免费网站| 亚洲国产美女精品久久久| 国产一级高清免费观看| 亚洲不卡中文字幕| 免费人成在线观看播放国产 | 亚洲国产精品免费在线观看| 色噜噜亚洲男人的天堂| 又爽又黄无遮挡高清免费视频 | 精品无码免费专区毛片| 亚洲色偷精品一区二区三区| 亚洲日韩国产精品乱| 免费A级毛片无码专区| 亚洲精华国产精华精华液好用| 国产性爱在线观看亚洲黄色一级片 | 成人免费视频观看无遮挡| 午夜肉伦伦影院久久精品免费看国产一区二区三区 | 亚洲AV无码一区二区三区性色 | 四虎免费在线观看| 中国在线观看免费的www|