摘要： 消元法这种直接解法适合低维方程组求解, 一般情况高维方程组的求解用迭代法Jacobi迭代法求解Ax=b中Jacobi.h#include <iostream>/////////////////////////class JACOBI {public : JACOBI();public : void red_matA(int n); ////读取矩阵mat_A和mat_b, n为方阵的阶数 void creatLUD(void); //创建LUD void creatD_iv(void); //求mat_D的逆矩阵 void cal(double); //求解方程组private : 阅读全文

摘要： LU.h#include <iostream>double **mat_L, **mat_U; //矩阵L Uint N; //方阵的阶class LU {public : LU();public : void or_mat(int N); //初始矩阵(系数矩阵mat_A, 右值矩阵mat_b)初始化,N为方阵阶 void creatLU(void); //将矩阵A 分解为L U void cal_mat(void); //用LU分解法求Ax=b的解private : double **mat_A, *mat_b; double *result; //结果数组};LU_File. 阅读全文

摘要： crv_fit.h//多项式曲线拟合 f(x)=a0+a1x+a2x^2+a3x^3+...anx^nclass Crv_fit {public : Crv_fit(void); void clear(void); //~Crv_fit(void);public : void init_ary(int m); //m为样本点个数; 初始化数组ary_x, ary_y; 初始化ary_x_cp,之后和ary_x共同为计算ary_A提供数据 void cal_ary (int n); //n+1矩阵A的阶数; 计算矩阵中元素,分别存放在ary_A, ary_b中 void init_matA(vo 阅读全文

摘要： #include <iostream>#include <malloc.h>#include <math.h>//非齐次线性方程组Ax=bdouble **mat_A; //存放系数矩阵Adouble *mat_b; //存放右值数组bdouble *result; //结果数组const int M=3; //方阵的阶const int N=3;///////////////////////////////////////void init_ary() { //初始化系数矩阵 A 和右值数组 b; std::cout << "系数矩阵 阅读全文

摘要： #include <iostream>#include <malloc.h>///////////////////////////char *ary; //数组指针////////////////////////////////////////////////void init_ary(int n); //数组个数void srt_l_h(int n); //将数组按值从小到大排列void srt_a_b(char *a, char *b); //交换a,bint find_srt(int n); //找到待排序的位置并进行排序void print(int n); // 阅读全文

摘要： Crank_Nicolson.h#include <GL/glut.h> #include <stdio.h> #include <stdlib.h>#include <iostream>#include <malloc.h>#include <math.h>int xnseg_ = 100;int ynseg_ = 100;int tnseg_ = 100;struct Heat { double u[101][101];};struct Heat heat[101];struct Heat_h { double u_h 阅读全文

摘要： bunny.h#include <iostream>double **orAry; //初始样本集矩阵(减去均值后的新样本矩阵)int sapNum_, dim_; //样本个数;维数;降维后维数class BUNNY {public : BUNNY(void); void clear(void);public : //计算协方差 void rdOrMatrix(int _sapNum, int _dim); //_sapNum为样本个数, _dimt为维数; 创建样本集矩阵并存入数组 orAry(_sapNum行, _dimt列) void meanVal(void); //计算 阅读全文

摘要： ValVece.h#include <iostream>class ValVect {public : ValVect(void); void clear(void); //~ValVect(void);public : void rdOrMatrix(int _dim, double _e); //从文件中读取对称方阵,将其存在矩valAry中(每一列为一样本向量,即行数为维数)._dim为方阵的阶,_e为计算精度 void itVtAry(void); //初始化样本向量矩阵vectAry(每一列为一特征向量) void calValVect(void); //计算特征值和特征 阅读全文

摘要： covHead.h#include <iostream>class COV {public : COV(void); void clear(void);public : void rdOrMatrix(int _sapNum, int _dim); //_sapNum为样本个数, _dimt为维数; 创建样本集矩阵并存入数组 orAry(_sapNum行, _dimt列) void meanVal(void); //计算各维度的均值, 并存入数组 meanValAry(一维数组, 维数为_dimt) void meanMatrix(void); //计算各维减去各维均值后的新矩阵并 阅读全文