- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值分析实验报告二
求解线性方程组的直接方法(2 学时)
班级专业 11 信科1 班 姓名 朱洪应 学号 201130760134 日期2011.04.11
一 实验目的
1.掌握求解线性方程组的高斯消元法及列主元素法;
2. 掌握求解线性方程组的克劳特法;
3. 掌握求解线性方程组的平方根法。
二 实验内容
1.用高斯消元法求解方程组 (精度要求为 6 ):
10
3 x1 x2 2 x3 7
x1 2 x2 2 x3 1
2 x 2 x 4 x 0
1 2 3
2.用克劳特法求解上述方程组(精度要求为106 )。
3. 用平方根法求解上述方程组(精度要求为 106 )。
4. 用列主元素法求解方程组 (精度要求为 6 ):
10
3 x1 x2 4 x3 3
x1 2 x2 2 x3 2
2 x 3 x 2 x 5
1 2 3
三 实验步骤(算法)与结果
#include stdio.h
#include stdlib.h
#include math.h
typedef struct
{
double *base;
int dim;
} Array;
1
//定义n*n矩阵及其方法
inline double getArray (Array A, int i, int j)
{
return A.base[(i - 1) * A.dim + (j - 1)];
}
inline void setArray(Array A, int i, int j, double num)
{
A.base[(i - 1) * A.dim + (j - 1)] = num;
}
void InitArray(Array A, int d)
{
A.base = (double *) malloc(d * d * sizeof(double));
A.dim = d;
}
void DestoryArray (Array A)
{
if (A.base) free(A.base);
A.base = NULL;
A.dim = 0;
}
//定义全局矩阵
Array a, b, l, u, z, x;
//
//定义高斯克劳特算法内部的多项求和
double GaussKeloutInside (int i, int j, char K)
{
double g = 0;
int k = 1;
if (l == K)
{
//计算l时的求和
if (j - 1 1) return 0;
for (k = 1; k = j - 1; k++)
{
文档评论(0)