回 帖 发 新 帖 刷新版面

主题:求一fortran题解答

2.线性病态方程组问题。
下面是一个线性病态方程组:
[1/2 1/3 1/4][X1]=[0.95]
1/3 1/4 1/5   X2   0.67
1/4 1/5 1/6   X3   0.52
 
(1)求方程的解。
(2)将方程右边向量元素b3改为0.53,再求解,并比较b3的变化和解的相对变化。
(3)计算系数矩阵A的条件数并分析结论。
提示:矩阵A的条件数等于A的范数与A的逆矩阵的范数的乘积,即 。这样定义的条件数总是大于1的。条件数越接近于1,矩阵的性能越好,反之,矩阵的性能越差。矩阵A的条件数 ,其中 ,aij系矩阵A的元素。
要求:
(1)方程的系数矩阵、常数向量均从文件中读入。
(2)定义求解线性方程组Ax=b的子程序,要求该子程序能求解任意线性方程组。
(3)在主程序中调用子程序,并对求解结果进行对比分析。
(4)绘制常数向量修改前后所求得的方程解的数据分布图。

回复列表 (共1个回复)

沙发

偷懒的话,就调用商业库(IMSL, MKL),要想开源的话 到 netlib 去找。

我来回复

您尚未登录,请登录后再回复。点此登录或注册