文件名称：北邮数值与符号计算实验 线性代数方程组求解
文件大小：25KB
文件格式：DOCX
更新时间：2019-06-17 04:27:10
北邮 数值 符号计算
注意：实验报告不全，参考价值：函数实现 用C/C++语言实现如下函数： bool lu(double* a, int* pivot, int n); 矩阵的LU分解 bool guass(double const* lu, int const* p, double* b, int n); 求线性代数方程组的解 void qr(double* a, double* d, int n); 矩阵的QR分解 bool hshld(double const*qr, double const*d, double*b, int n); 求线性代数方程组的解 用上面的函数求解下面的问题： 求解n阶Hilbert方阵的线代数方程组 Hilbert矩阵在第j行，第k列处的值为 H_jk=1/(j+k-1) 适当设置b的值使得理论解围x=(1,1,…,1)T。列选主元Gauss消去法与Householder变换法分别对n=10,20,30求解。给出数值解与理论解之差的二范数。 1.2 求Gn矩阵的线代数方程组 令n×n矩阵G_n 令n维向量b_n=(2,1,0,-1,-2,…,5-n,4-n,2-n)^T 用两种方法对n=50,100,150求解。理论值为x=(1,1,…,1)T。给出数值解与理论解之差的二范数。 再令b_n=(2,3,4…,n,n)^T用两种方法对n=50,100,150求解。给出数值解的残量的二范数。