列主元消去法例题详解_列主元消去法

本文详细解析了列主元消去法在解决线性方程组中的作用,强调了选择绝对值最大主元以降低计算误差的重要性。通过实例介绍了如何实施列主元消去法,并探讨了其在数值计算中的优势。
摘要由CSDN通过智能技术生成

列主元消去法

________

实验类型:

________________

同组学生姓名:

__________

一、实验目的和要求(必填)

二、实验内容和原理(必填)

三、主要仪器设备(必填)

四、操作方法和实验步骤

五、实验数据记录和处理

六、实验结果与分析(必填)

七、讨论、心得

一、问题描述

对于一般的线性方程组,只要把方程组化成了等价的三角形方程组,求解过程就很容易完

成,

Gauss

消去法就是将一般的线性方程组等价地变换为一个上三角方程组,然后用回代法

求解。

但是如果在消元过程中,

发现某个约化主元

a

kk

(

k

)

=0

则第

k

次消元就无法进行。

此外,

使所有约化主元全不为零,虽然可以完成方程组的求解,但是小主元的存在使计算结果误差

较大。

因此,为了减少计算过程中舍入误差对解的影响,在每次消元前,应选择绝对值尽可能大

的元作为约化的主元。我们称这种消元法为主元消元法。如果在子块的第一列中选取主元,

则相应的方法称为列主元消元法。

二、相关公式

设有线性方程组

b

Ax

其中

#include<stdio.h> #include<math.h> #define N 100 #define epsilon 1e-6 float a[N][N+1]; void menu( ) { printf("\t\t%c%c%c^_^Gauss列主元消去法求解线性方程组^_^%c%c%c\n\n",1,1,1,1,1,1); printf("强烈建议您先阅读以下几点后在运行:\n"); printf("1.这是用Gauus列主元消去法求解线性方程组的应用程序\n"); printf(" (Gauus全主元消去法类似可做,读者有兴趣的话可自行而做)\n"); printf("2.请您先了解Gauus列主元消去法的主要思想\n"); } void main( ) { int i,j,k,n; float t,s=0; char choice; menu( ); loop: printf("\n请输入系数方阵的阶数:"); scanf("%d",&n); while(n>0) { printf("\n"); printf("请输入增广阵矩:\n"); for(i=0;i<n;i++) for(j=0;j<n+1;j++) scanf("%f",&a[i][j]);/*存储增广阵矩*/ for(k=0;k<n-1;k++) { for(i=k+1;i<n;i++)/*求出每列中最大数,后行交换*/ if( fabs(a[i][k]) > fabs(a[k][k]) ) for(j=k;j<n+1;j++) { t=a[k][j]; a[k][j]=a[i][j]; a[i][j]=t; } if( fabs(a[k][k]) < epsilon)/*最大数小于很小数时退出*/ { printf("\n错误,Gauss列主元消去法无法忍受,在%d步退出!\n",k+1); printf("还要再计算其他的么(Y/N)?"); scanf("%c",&choice); if(choice=='Y' || choice=='y')/*判断用户输入*/ goto loop; else return; } for(i=k+1;i<n;i++) { a[i][k]=a[i][k] / a[k][k]; for(j=k+1;j<n+1;j++) a[i][j]=a[i][j]-a[i][k] * a[k][j]; } } a[n-1][n]=a[n-1][n] / a[n-1][n-1]; for(k=n-2;k>=0;k--) { s=0; for(j=k+1;j<n;j++) s+=a[k][j]*a[j][n]; a[k][n]=( a[k][n]-s ) / a[k][k]; } printf("\n*****运行结果*****\n"); for(i=0;i<n;i++) printf(" x[%d]=%.4f\n",i+1,a[i][n]); printf(" 谢谢使用!\n"); printf("还要再计算其他的么(Y/N)?"); getchar( ); scanf("%c",&choice); if(choice=='Y' || choice=='y')/*判断用户输入*/ goto loop; else return; } }
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值