原创力文档- 1、本文档共21页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
1_3936179_东北大学数值分析实验报告.doc
实验一：解线性方程组的直接方法
1、设线性方程组
=
x= ( 1, -1, 0, 1, 2, 0, 3, 1, -1, 2 )
2、设对称正定阵系数阵线方程组
=
x = ( 1, -1, 0, 2, 1, -1, 0, 2 )
3、三对角形线性方程组
=
x= ( 2, 1, -3, 0, 1, -2, 3, 0, 1, -1 )
分别用Gauss顺序消去法与Gauss列主元消去法：平方根法与改进平方根法：追赶法求解，编出算法通用程序
用列主元消去法求解方程组一：
流程图如下：
源程序：
#include iostream
#include vector
#include cmath
using namespace std;
class CGAUSSSOLVEEQU
{
private :
vector vector double m_equset;
vector double m_answer;
int m_n;
public :
void inputEquSet(double in[],int n);
void solveEquSet();
void outputAnswer();
void change(int m,int m2);
} ;
void CGAUSSSOLVEEQU::inputEquSet(double in[],int n)
{
vector double vtemp;
m_n=n;
for (int i= 0;i m_n; i++)
{
m_equset.push_back(vtemp);
for (int j= 0;j = m_n; j++)
{
m_equset[i].push_back(in[i*(m_n+1)+j]);
}
}
}
void CGAUSSSOLVEEQU::change(int m,int m2){
vector vector double ::iterator iter;
iter = m_equset.begin();
vector vector double ::iterator iter2;
iter2 = m_equset.begin();
//double
}
void CGAUSSSOLVEEQU::solveEquSet()
{
vector vector double ::iterator iter;
iter = m_equset.begin();
for (int m= 0;m m_n - 1 ; ++ m)
{
// 将绝对值最大的主元素移上去
for (int i=m;im_equset.size();i++)
{
if (fabsl(m_equset[m][m]) fabsl(m_equset[i][m]))
{
swap( m_equset[m], m_equset[i]);
}
}
// 进行消元
for (int i = m + 1 ;i m_n; ++ i)
{
double dm;
dm = m_equset[i][m] / m_equset[m][m];
for (int j = m;j m_n + 1 ; ++ j)
{
m_equset[i][j] -= dm * m_equset[m][j];
}
}
++ iter;
}
// 初始化m_answer向量
for (int i=0 ;i m_n; ++ i) m_answer.push_back( 0 );
// 求解答案
m_answer[m_n - 1 ] = m_equset[m_n - 1 ][m_n] / m_equset[m_n - 1 ][m_n - 1 ];
for ( int i = m_n - 2 ;i = 0 ; -- i)
{
m_answer[i] = m_equset[i][m_n];
for ( int j = m_n - 1 ;j i; -- j
文档评论（0）