数值分析上机题 , 会的追加分!!!!!!

来源:百度知道 编辑:UC知道 时间:2024/06/08 08:18:16
给定方程组

, 即 AX=b

1. 分别用列主元消去法和追赶法解上述三对角方程组AX=b。
2. 讨论当b 有微小变化,即 b+ b=(1.01, 0.01, 0.05, 0.00, 0.00) 时,A是否为病态矩阵。给出讨论结果。
3. 输出分解矩阵L, U 及解向量 X和|A|。
4. 对两种计算结果加以比较。

最好给下程序
程序设计流程图
| 2 -1 0 0 0 | x1 | 1 |
|-1 2 -1 0 0 | x2 | 0 |
| 0 -1 2 -1 0 | x3 = | 0 |
| 0 0 -1 2 -1 | x4 | 0 |
| 0 0 0 -1 2 | x5 | 0 |

我的邮箱是skinyto@163.com
谢谢这位大哥了阿!

我给你发我的题了你帮我在看看!

给我邮箱,我给你发,是不是这样的 ,程序能给你但是他们的流程图没有办法贴
第五章 解线性代数方程组的直接法
1、编制解 的通用子程序
(2)用改进平方根法解方程组 ;
流程图

(3)实现矩阵三角分解的杜利特尔及乔利斯基方法及用此方法解 的过程。
流程图

源程序:
A:#include <iostream.h>
#include <math.h>
int cholesky1(double *p,double *q,int n); //平方根法
int cholesky2(double *p,double *q,int n); //改进平方根法
void main ()
{ int i,j,k,w,n=4;
double a[4][4]={4,2.4,2,3,2.4,5.44,4,5.8,2,4,5.21,7.45,3,5.8,7.45,19.66};
double *p,*q,b[4]={12.280,16.928,22.957,50.945};
for(i=0;i<n;i++)
{for(j=0;j<n;j++)
cout<<"a["<<i+1<<"]["<<j+1<<"]="<<a[i][j]<<"\t";
cout<<endl;
}
for(i=0;i<n;i++)
cout<<"b["<<i+1<<"]="<<b[i]<<endl;
cout<<endl;
p=&a[0][0],q=&b[0];
do
{cout<<"选择计算方法k=(1:是平方根法,2:是改进平方根法)&