Matrix Template Library Help

部分函数功能说明

// element access functions
    int GetRowSize()const
    int GetColSize()const
    T   GetValue(const int& row,const int& col)const
    const vector<T> GetRowVector(const int& row)const
    const vector<T> GetColVector(const int& col)const
    const Matrix<T> GetRow(const int& row)const
    const Matrix<T> GetCol(const int& col)const

// modify     
    void Resize(const int& row,const int& col,T value=T(0));
    void SetValue(const int& row,const int& col,T value = T(0));    
    void AddRow(T value = T(0));
    void AddCol(T value = T(0));
    void DeleteRow(const int& row);
    void DeleteCol(const int& col);
    void InsertRow(const int& row,T value = T(0));
    void InsertCol(const int& col,T value = T(0));
    void SetRowVector(const int& row,const vector<T>& vec = vector<T>(0));
    void SetColVector(const int& col,const vector<T>& vec = vector<T>(0));
    void SwapRow(const int& row1,const int& row2);
    void SwapCol(const int& col1,const int& col2);

【 这些函数不用介绍了，功能根据函数顾名思义。】

// special functions
    Matrix<T> Turn();                        //求转置矩阵
    Matrix<T> Inverse();                   //求逆矩阵（必须是方阵，并且可逆，否则抛出exception）
    Matrix<T> Unit(const int& row); //单位矩阵（把矩阵重置为给定大小的单位阵）
    Matrix<T> Diag();                        //求对角阵（返回对角阵元素组成的行矩阵）
    T         Det();                                 //求矩阵的行列式的值（必须是方阵，否则抛出exception）
    T         Trace();                             //求方阵的迹（必须是方阵，否则抛出exception）
    T         Norm(const int& p=2);  // 求矩阵的范数（默认为求2范数。注意：只有folat，double，long double矩阵可以调用该函数）
    T         CondNumber();              //求矩阵的条件数
    T         Radius();// spectral  radius  //求矩阵的半径
    int       Rank();                             //求矩阵的秩
    Matrix<T> GetAij(const int& row,const int& col); //取得矩阵的子矩阵（大的截取，小的补零）

// transform
    CMatrix< complex<double> > ToCDouble();  //将矩阵转换为double型的复数矩阵，以便于作矩阵间的运算

//实数矩阵特有函数
    RMatrix<T> Jacobi(RMatrix<T>& value=RMatrix<T>());  //JACOBI迭代法求实数矩阵的特征值，特征解
    RMatrix<double> ToDouble();    //转换为double型矩阵

// 复数阵特有函数
    RMatrix<T>  Real();   // 取得实部矩阵
    RMatrix<T>  Imag();  // 取得虚数矩阵

【下面给出一个应用此矩阵库的小程序，示意主要函数的用法】

#include <iostream>
#include <stdlib.h>
#include "matrix6.hpp"

using namespace std;
using namespace soubam;

int main(int argc, char *argv[])
{
try
{

  complex< double >    temp(3,4); 
  RMatrix< int >       R1( 4, 3, 1),  R2(3, 5, 2);
  RMatrix< double >    R3( 4, 4 ),    R4(4, 3, 0.1);
  CMatrix< int >    C1( 4, 3,1),   C2(3, 5, 2);
  CMatrix< double >    C3( 5, 6, 2.3),C4(6, 7, 0.3);

  R1[0][0] = 3;   
  R1[0][1] = 5;
  R1[0][2] = 4;

  cout<<"R1="<<endl;
  cout<<R1<<endl; 

  R3.SetValue(0,1,1);
  R3.SetValue(1,1,2);
  R3.SetValue(2,1,3);
 
  cout<<"R3="<<endl;
  cout<<R3<<endl;
 
  cout<<"R4="<<endl;
  cout<<R4<<endl;

  cout<<"R3*R4="<<endl;
  cout<<R3*R4<<endl;

  cout<<"rank(R1)= "<<R1.Rank()<<endl;
   

  R2.Det();

  system("PAUSE"); 
  return 0;

}
 catch(OutOfRange e)
 {
  cout<<e.GetError()<<endl;
  cout<<"row= "<<e.InvalidRow()<<endl;
  cout<<"col= "<<e.InvalidCol()<<endl;
 }

 catch(MisMatch e)
 {
  cout<<e.GetError()<<endl;
 }
}