数值分析实验报告
实验一 Gauss消元法 姓名:杨 玲
实验二 列主元消元法 学号:2008115010147
实验三 三角分解 院系:计算机科学与技术学院
实验一 Gauss消元法
课题名称:Guess消元法 任课教师:李国屏
专业班级:计算机科学与技术 学 号:2008115010147 姓 名:杨 玲
※ 题目描述
用Gauss消元法解n阶线性代数方程组:
其基本做法是把上述方程组通过消元转化为一个等价的三角形方程组,然后再进行回代就可以求出方程组的解。
要求显示出每一大步消元后的系数矩阵和常数项,最后显示出方程组的解。
※ 算法分析
(1) 算法思想:
用Gauss消元法把上述方程组的系数矩阵化为上三角矩阵的过程称为消元过程,消元过程中的算法
设计
领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计
为:
for(t=i;t
=0;i--) //回代过程
{
v=0;
for(j=i+1;j<=n;j++)
v=a[i][j]*x[j]+v;
x[i]=(b[i]-v)/a[i][i];
}
※ 程序说明:
本程序在C++ 环境中编译运行并且通过测试,通过提示语句输入相应的系数矩阵和常数项。
※ 程序代码
#include
#include
const int N=100;
using namespace std;
int main()
{
int n,i,j,k,m,h,t;
double a[N][N],b[N],x[N],w,v;
cout<<"请输入矩阵的阶数:"<>n)
{
cout<<"请输入系数矩阵:"<>a[i][j];
cout<<"请输入常数项:"<>b[i];
cout<<"******************************"<=0;i--) //回代过程
{
v=0;
for(j=i+1;j<=n;j++)
v=a[i][j]*x[j]+v;
x[i]=(b[i]-v)/a[i][i];
}
cout<<"该方程组的解为:"<=1;i--)
{
v=0;
for(j=i+1;j<=n;j++)
v=a[i][j]*x[j]+v;
x[i]=(b[i]-v)/a[i][i];
}
※ 程序说明
本程序可以在VC++ 6.0环境中编译运行,通过提示语句输入相应的系数矩阵和常数项。
※ 程序代码
#include
#include
#include
int const N=100;
using namespace std;
int main()
{
int n,i,j,k,m,h,t,q;
double a[N][N],b[N],x[N],w,v,p;
cout<<"请输入矩阵的阶数:"<>n) //以下为数据输入,并显示所求方程
{
cout<<"请输入系数矩阵:"<>a[i][j];
}
}
cout<<"请输入常数项:"<>b[i];
cout<<""<=1;i--)
{
v=0;
for(j=i+1;j<=n;j++)
v=a[i][j]*x[j]+v;
x[i]=(b[i]-v)/a[i][i];
}
cout<<"该方程组的解为:" <
#include
const int N=100;
using namespace std;
int main()
{
int n,i,j,m,k,h,r;
double a[N][N],l[N][N],u[N][N],w,v;
cout<<"请输入矩阵的阶数:"<>n)
{
cout<<"请输入矩阵:"<>a[i][j];
cout<<"******************************"<j)
u[i][j]=0;
}
for(j=2;j<=n;j++)
l[j][1]=a[j][1]/u[1][1];
for(r=2;r<=n;r++)
{
for(m=r;m<=n;m++)
{
for(j=1;j<=r-1;j++)
{
w=0;
w=l[r][j]*u[j][m]+w;
}
u[r][m]=a[r][m]-w;
}
for(h=r+1;h<=n-1;h++)
{ for(k=1;k<=r-1;k++)
{
v=0;
v=(l[h][k]*u[k][r]+v)/u[r][r];
}
l[h][r]= a[h][r]/u[r][r]-v;
}
}
cout<<"所求的LU分解中的矩阵U为:"<
本文档为【数值方法实验报告】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。