2012年11月6日

用LU分解法解非齐次线性方程组Ax=b

摘要: 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. 阅读全文

posted @ 2012-11-06 22:09 timeflies 阅读(2246) 评论(0) 推荐(0)

最小二乘法求多项式的拟合曲线

摘要: 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 阅读全文

posted @ 2012-11-06 09:26 timeflies 阅读(3213) 评论(0) 推荐(0)

高斯消元法求解向量方程Ax=b

摘要: #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 << "系数矩阵 阅读全文

posted @ 2012-11-06 08:57 timeflies 阅读(488) 评论(0) 推荐(0)

导航