數值分析實驗報告

2021-03-04 05:37:21 字數 2617 閱讀 8908

軟體學院

實驗目的:

1、通過實驗進一步了解方程求根的演算法;

2、認識選擇計算格式的重要性;

3、掌握迭代演算法和精度控制;

4、明確迭代收斂性與初值選取的關係。

實驗過程:

1、編制乙個程式進行運算,最後列印出每種迭代格式的斂散情況;

2、用事後誤差估計來控制迭代次數,並且列印出迭代的次數;

3、初始值的選取對迭代收斂有何影響;

4、分析迭代收斂和發散的原因。

實驗環境:

作業系統:windows xp

程式語言:visual c++ 6.0。

實驗源**:

int main()

x=y;

}if(n==n+1)

cout<<"迭代格式1不收斂!"< //第2種迭代格式

n=0;

x=a;

while((n++)

x=y;

}if(n==n+1)

cout<<"迭代格式2不收斂!"< //第3種迭代格式

n=0;

x=a;

while((n++)

x=y;

}if(n==n+1)

cout<<"迭代格式3不收斂!"<

//第4種迭代格式

n=0;

x=a;

while((n++)

x=y;

}if(n==n+1)

cout<<"迭代格式4不收斂!"<

//第5種迭代格式

n=0;

x=a;

while((n++)r)

x=y;

}if(n==n+1)

cout<<"迭代格式5不收斂!"<

//第6種迭代格式

n=0;

x=a;

while((n++)

x=y;

}if(n==n+1)

cout<<"迭代格式6不收斂!"< return 0;

}實驗分析與總結:

由程式執行結果可知迭代格式一不收斂,且初值的選取影響迭代格式的收斂性。

這次實驗基本完成了實驗所要求的內容,同時達到了實驗的目的,也就是通過程式設計進一步理解了各種迭代格式的比較。

實驗目的:

1、通過該課題的實驗,體會模組化結構程式設計方法的優點;

2、運用所學的計算方法,解決各類線性方程組的直接演算法;

3、提高分析和解決問題的能力,做到學以致用;

4、通過三對角形線性方程組的解法,體會稀疏線性方程組解法的特點。

實驗過程:

1、 對題目給出的三個方程組分別利用gauss順序消去法與gauss列主元消去法;

2、 應用結構程式設計編出通用程式;

3、 比較計算結果,分析數值解誤差的原因;

實驗環境:

作業系統:windows xp

程式語言:visual c++ 6.0。

實驗源**:

gauss順序消去法

gauss列主元消去法

gauss順序消去法:

#include < iostream >

#include < vector >

#include < cmath >

using namespace std;

class cgausssolveequ

; void cgausssolveequ::inputequset(double in,int n)

}} void cgausssolveequ::change(int m,int m2)

void cgausssolveequ::solveequset()

// 初始化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)

}void cgausssolveequ::outputanswer() }

int main();

myequ1.inputequset(in1,10);

myequ1.solveequset();

myequ1.outputanswer();

cout< cout< return 1 ;

}執行結果:

x(i)=( 1 , -1 , -2.66454e-015 , 1 , 2 , 8.88178e-016 , 3 , 1 , -1 , 2 )

gauss列主元消去法:

#include < iostream >

#include < vector >

#include < cmath >

using namespace std;

class cgausssolveequ

{ private :

vector < vector < double >> m_equset;

vector < double > m_answer;

數值分析實驗報告

實驗一題目一 已知函式表如下 試用三次lagrange插值多項式求x 0.5635時的函式值。解 lagrange插值法公式 程式框圖如下 本題中n 4 分析設計思想 設計思路如上面的流程圖所示,用我們學過的c語言程式設計實現。程式清單 include include include float l...

數值分析實驗報告

數值分析 課程設計實驗報告模板 常微分方程資料值解 用龍格 庫塔法分析lorenz方程的特性 一 考慮著名的lorenz方程 其中s,r,b為變化區域內有一定限制的實引數,該方程形式簡單,表面上看並無驚人之處,但由該方程揭示出的許多現象,促使 混沌 成為數學研究的嶄新領域,在實際應用中產生了巨大的影...

數值分析實驗報告

實驗專案 a 列主元高斯消去法 b 矩陣直接三角分解法 c 雅可比迭代法 d 高斯 賽德爾迭代法 e sor迭代法 f 冪法求最大特徵值 實驗名稱 列主元高斯消去法 實驗目的 熟悉了解高斯消去法的原理,並知道如何用高斯消去法求方程組的解 實驗所用的軟體 win tc 版本 1.9.0.0 例項用高斯...