C++求特征值.docVIP

C++求特征值 #includeiostream.h #includestdlib.h #includemath.h #includeiomanip.h //+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ //class Matrix定义矩阵类 const int Max_xy=20; //矩阵的最大维数 class Matrix { private: double data[Max_xy][Max_xy]; unsigned x,y; //x,y; public: Matrix(); //默认构造函数 Matrix(const Matrix source); //拷贝构造函数 void creat(); //输入矩阵 void init(); void transpose(); //矩阵转置 void show(); //输入此矩阵 double mode() const; //求一维矩阵的长度 void check_shiduichen(); //检查是否为是对称矩阵 void creat_unit(unsigned i); //生成i行单位矩阵 void set_x(unsigned xx); //设置行数 void set_y(unsigned yy); //设置列数 unsigned get_x(); //得到行数 unsigned get_y(); //得到列数 void shucheng(double changshu); //数乘运算 void setdata(unsigned i,unsigned j,double source); //定位输入数据 double getdata(unsigned i,unsigned j); //定位得到数据 void sturm(); //求特征值 unsigned sturm_s(double m); //计算sturm系列的同好数 Matrix operator = (const Matrix right); friend Matrix operator + (const Matrix left,const Matrix right); //重载+号 friend Matrix operator - (const Matrix left,const Matrix right); //重载-号 friend Matrix operator * (const Matrix left,const Matrix right); //重载乘号 friend ostream operator (ostream os,const Matrix source); //重载输出 friend void Householder(Matrix source); //用Householde矩阵将实对称矩阵化为三对角矩阵 }; Matrix temp_Matrix; //全局变量Matrix //=================================================================== //--------------------默认构造函数 Matrix::Matrix() { init(); } //----------------------------拷贝构造函数 Matrix::Matrix(const Matrix source) { init(); x=source.x; y=source.y; for(unsigned i=0;ix;i++) for(unsigned j=0;jy;j++) data[i][j]=source.data[i][j]; } //------------------------------------------初始化矩阵元素 void Matrix::init() { x=y=0; for(unsigned i=0;iMax_xy;i++) for(unsigned j=0;jMax_xy;j++) data[i][j]=0; } //------------------------------矩阵转置 void Matrix::transpose() { double temp;int t; for(u