matrix.h

一些矩阵运算的函数

#ifndef MAXTRIX_H  //如果被定义了，则跳到下一个#else或下一个#endif
                   //避免文件中的内容被定义两次
#define MAXTRIX_H

//两个矩阵相加(重载)
void matrix_add(double *,double *,double *,int,int);
void matrix_add(double *,double *,double *,int);
    //前面两个函数的计算结果保存在第三个矩阵
void matrix_add(double *,double *,int,int);
void matrix_add(double *,double *,int);
    //这两个函数的计算结果保存在第一个矩阵中


//两个矩阵相减(重载4个函数)
void matrix_minus(double *,double *,double *,int,int);
void matrix_minus(double *,double *,double *,int);
   //前面两个函数的计算结果保存在第三个矩阵
void matrix_minus(double *,double *,int,int);
void matrix_minus(double *,double *,int);
   //这两个函数的计算结果保存在第一个矩阵中


//两个矩阵相乘(重载2个函数）
void matrix_mul(double *,double *,double *,int,int,int);
void matrix_mul(double *,double *,double *,int);
    //第三个矩阵返回结果
void matrix_mul(double *,double *,int);//无法应用


//矩阵复制(重载2个函数)
void matrix_copy(double *,double *,int m,int n);
void matrix_copy(double *,double *,int n);
    //第二个矩阵返回结果


//矩阵求逆的全选主员高斯-约当法
int matrix_qiuni(double *,int n);
//对称正定矩阵的求逆
int dzmatrix_qiuni(double *,int); 
//这两个求逆函数还算比较准确，有一定的价值。我觉得。


//矩阵的三角分解(LU分解)
void matrix_fenjie(double *,double *,double *,int);
//一般实矩阵的QR分解
void matrix_qrfenjie(double *,double *,int,int);


//求Hessenberg矩阵全部特征根的QR法
int matrix_hqr(double *a,double *u,double *v,double eps,int n,int jt);
//该式子可以用来求高阶方程的根

//约化一般实矩阵为Hessenberg矩阵的初等相似变换
void matrix_hessen(double *a,int n);












#endif