数值分析列主元消去法的实验报告

1、实验一 列主元消去法【实验内容】1掌握列主元消去法的基本思路和迭代步骤2并能够利用列主元的高斯消去法解任意阶数的线性方程组；3、从课后题中选一题进行验证，得出正确结果，交回实验报告与计算结果。【实验方法与步骤】1 列主元消去法基本思路 设有线性方程组，设是可逆矩阵。列主元消去法的基本思想就是通过列主元的选取将初等行变换作用于方程组的增广矩阵，将其中的变换成一个上三角矩阵，然后求解这个三角形方程组。2 列主元高斯消去法算法描述将方程组用增广矩阵表示。步骤1：消元过程，对（1） 选主元，找使得（2） 如果，则矩阵奇异，程序结束；否则执行（3）；（3） 如果，则交换第行与第行对应元素位置， ；（4）

2、 消元，对，计算对，计算步骤 2：回代过程：（1） 若则矩阵奇异，程序结束；否则执行（2）；（2） 对，计算 实验程序#include<math.h> #include<stdio.h> #include<iostream>#include<conio.h> #define number 20 #define esc 0x1b #define enter 0x0dusing namespace std;float anumbernumber+1 ,ark; int flag,n; void exchange(int r,int k); float

3、 max(int k); void message();void main() float xnumber; int r,k,i,j; char celect; void clrscr(); printf("nnuse gauss."); printf("nn1.jie please press enter."); printf("nn2.exit press esc."); celect=getch(); if(celect=esc) exit(0); printf("nn input n="); scanf(&

4、quot;%d",&n); printf(" nninput matrix a and b:"); for(i=1;i<=n;i+) printf("nninput a%d1-a%d%d and b%d:",i,i,n,i); for(j=1;j<=n+1;j+) scanf("%f",&aij); for(k=1;k<=n-1;k+) ark=max(k); if(ark=0) printf("nnits wrong!");message(); else if(flag

5、!=k) exchange(flag,k); for(i=k+1;i<=n;i+) for(j=k+1;j<=n+1;j+) aij=aij-akj*aik/akk; xn=ann+1/ann; for( k=n-1;k>=1;k-) float me=0; for(j=k+1;j<=n;j+) me=me+akj*xj; xk=(akn+1-me)/akk; for(i=1;i<=n;i+) printf(" nnx%d=%f",i,xi); message(); void exchange(int r,int k) int i; for(i=1;i<=n+1;i+) a0i=ari; for(i=1;i<=n+1;i+) ari=aki; for(i=1;i<=n+1;i+) aki=a0i; float max(int k) int i; float temp=0; for(i=k;i<=n;i+) if(fabs(aik)>temp) temp=fabs(aik); flag=i; return temp; void message() printf("nn go on enter ,exit press esc!"); switch(getch() case