計(jì)算數(shù)值方法實(shí)驗(yàn)報(bào)告--太原理工大學(xué)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上計(jì)算數(shù)值方法實(shí)驗(yàn)報(bào)告學(xué)院：軟件學(xué)院專(zhuān)業(yè)：軟件工程班級(jí)：軟件1012班學(xué)號(hào)：姓名：?jiǎn)替悍逄砉ご髮W(xué)學(xué)生實(shí)驗(yàn)報(bào)告學(xué)院名稱(chēng)軟件學(xué)院專(zhuān)業(yè)班級(jí)軟件1012班學(xué)號(hào)學(xué)生姓名喬婧峰實(shí)驗(yàn)日期20124成績(jī)課程名稱(chēng)數(shù)值計(jì)算方法實(shí)驗(yàn)題目實(shí)驗(yàn)一 二分法一、課題名稱(chēng)方程求根：熟悉使用、迭代法、牛頓法、割線法等方法對(duì)給定的方程進(jìn)行根的求解。選擇上述方法中的兩種方法求方程：二分法f(x)=x3+4x2-10=0在1,2內(nèi)的一個(gè)實(shí)根，且要求滿足精度|x*-xn|<0.5×10-5迭代法：用迭代公式x=f(x)進(jìn)行迭代計(jì)算，直到滿足|x*-xn|<0.5×10-5 為

2、止 。割線法：x=x-f(x)/g(x),其中f(x)為給定的函數(shù)，g(x)為給定函數(shù)的導(dǎo)數(shù)，直到滿足|x*-xn|<0.5×10-5 為止 。二、目的和意義（1）了解非線性方程求根的常見(jiàn)方法，如二分法、牛頓法、割線法。（2）加深對(duì)方程求根方法的認(rèn)識(shí)，掌握算法。（3）會(huì)進(jìn)行誤差分析，并能對(duì)不同方法進(jìn)行比較。三、計(jì)算公式f(x)在區(qū)間（x，y）上連續(xù) 先找到a、b屬于區(qū)間（x，y），使f(a)，f(b)異號(hào)，說(shuō)明在區(qū)間(a,b)內(nèi)一定有零點(diǎn)，然后求f(a+b)/2, 現(xiàn)在假設(shè)f(a)<0,f(b)>0,a<b 如果f(a+b)/2=0，該點(diǎn)就是零點(diǎn)， 如果f(a

3、+b)/2<0,則在區(qū)間（(a+b)/2，b)內(nèi)有零點(diǎn)，(a+b)/2=>a，從開(kāi)始繼續(xù)使用 中點(diǎn)函數(shù)值判斷。 如果f(a+b)/2>0，則在區(qū)間(a,(a+b)/2)內(nèi)有零點(diǎn)，(a+b)/2<=b，從開(kāi)始繼續(xù)使用 中點(diǎn)函數(shù)值判斷。 這樣就可以不斷接近零點(diǎn)。 通過(guò)每次把f(x)的零點(diǎn)所在小區(qū)間收縮一半的方法，使區(qū)間的兩個(gè)端點(diǎn)逐步迫近函數(shù)的零點(diǎn)，以求得零點(diǎn)的近似值四、主要儀器設(shè)備Vc2008,hp五、結(jié)構(gòu)程序設(shè)計(jì)迭代法: #include "stdafx.h"#include"stdio.h"#include"math.h

4、"#include"iostream"using namespace std;float main() float a; cin>>a; float t, x; x=a; do x=sqrt(10-x*x*x)/4); t=a; a=x; while(fabs(a-t)>0.5*1e-5); printf("x=%f",a); system("pause");割線法: #include "stdafx.h"#include"stdio.h"#include"

5、math.h"#include"iostream"using namespace std;float main() float c,a=1.0,b=2.0; /cin>>a>>b; while(1) c=b-(b*b*b+4*b*b-10)*(b-a)/(b*b*b+4*b*b-(a*a*a+4*a*a); if(fabs(b-c)<0.5*0.) break; b=c; cout<<c;六、結(jié)果討論和分析割線法: 迭代法:分析：使用不同的方法，可以不同程度的求得方程的解，不同的方法速度不同。實(shí)驗(yàn)地點(diǎn)ZSA401指導(dǎo)教師李

6、志學(xué)院名稱(chēng)軟件學(xué)院專(zhuān)業(yè)班級(jí)軟件1012班學(xué)號(hào)學(xué)生姓名喬婧峰實(shí)驗(yàn)日期2012.4成績(jī)課程名稱(chēng)數(shù)值計(jì)算方法實(shí)驗(yàn)題目實(shí)驗(yàn)二 線性方程組的直接解法一、課題名稱(chēng)線性方程組的直接解法合理利用Gauss消元法、LU分解法、追趕法求解下列方程組： （n=5,10,100）二、目的和意義（1）了解線性方程組常見(jiàn)的直接解法，如Guass消元法、LU分解法、追趕法。（2）加深對(duì)線性方程組求解方法的認(rèn)識(shí)，掌握算法。（3）會(huì)進(jìn)行誤差分析，并能對(duì)不同方法進(jìn)行比較。三、計(jì)算公式 高斯分解法：將原方程組化為三角形方陣的方程組：lik=aik/akk aij= aij- lik* akj k=1,2,n-1 i=k+1,k+2

7、, ,n j=k+1,k+2, ,n+1由回代過(guò)程求得原方程組的解： xn= ann+1/ ann xk=( akn+1-akj xj)/ akk (k=n-1,n-2, ,2,1) LU分解法：將系數(shù)矩陣A轉(zhuǎn)化為A=L*U, L為單位下三角矩陣，U為普通上三角矩陣，然后通過(guò)解方程組l*y=b,u*x=y,來(lái)求解x.追趕法:用來(lái)求對(duì)角方程組；將系數(shù)矩陣A轉(zhuǎn)化為A=L*U, L為普通下n-1對(duì)角矩陣，U為單位上n-1對(duì)角矩陣，然后通過(guò)解方程組l*y=b,u*x=y,來(lái)求解x.四、主要儀器設(shè)備Vc2008,hp五、結(jié)構(gòu)程序設(shè)計(jì) Gauss消元法：#include "stdafx.h&qu

8、ot;#include"stdio.h"#include"iostream"using namespace std; float main() float a34=1,2,3,14,0,1,2,8,2,4,1,13; float x3; float sum=0; int k,i,j; for(k=0;k<2;k+) for(i=k+1;i<3;i+) for(j=k+1;j<4;j+)aij=aij-aik/akk*akj; for(i=0;i<3;i+) for(j=0;j<4;j+) printf("a%d%d

9、=%f,",i,j,aij); cout<<endl; x2=a23/a22; for(k=1;k>=0;k-) sum=0; for(j=k+1;j<3;j+) sum+=akj*xj; xk=(ak3-sum)/akk; for(i=0;i<3;i+)printf ("x%d=%f,",i+1,xi); LU分解法：#include "stdafx.h"#include <stdio.h>#include <math.h> #define L 30 double a L L , b L

10、, l L L , u L L , x L , y L ; int main() int n, i, j, k, r; scanf( "%d", &n ); for ( i = 1; i <= n; +i ) for ( j = 1; j <= n; +j ) scanf( "%lf", &a i j ); for ( i = 1; i <= n; +i ) scanf( "%lf", &b i ); for ( i = 1; i <= n; +i ) for ( j = 1; j &l

11、t;= n; +j ) l i j =0; u i j = 0.0; for ( k = 1; k <= n; +k ) for ( j = k; j <= n; +j ) u k j = a k j ; for ( r = 1; r < k; +r ) u k j -= l k r * u r j ; for ( i = k + 1; i <= n; +i ) l i k = a i k ; for ( r = 1; r < k; +r ) l i k -= l i r * u r k ; l i k /= u k k ; l k k = 1.0; for (

12、i = 1; i <= n; +i ) y i = b i ; for ( j = 1; j < i; +j ) y i -= l i j * y j ; for ( i = n; i > 0; -i ) x i = y i ; for ( j = i + 1; j <= n; +j ) x i -= u i j * x j ; x i /= u i i ; for ( i = 1; i <= n; +i ) printf( "%0.2lfn", x i ); return 0;追趕法：#include "stdafx.h"

13、#include "stdio.h"void main() FILE *f; double a15,b15,c15,d15; double t; int i,n; f=fopen("zgf.txt","r"); fscanf(f,"%d",&n); fscanf(f,"%lf%lf%lf",&b1,&c1,&d1); for(i=2;i<=n-1;i+) fscanf(f,"%lf%lf%lf%lf",&ai,&bi,&a

14、mp;ci,&di); fscanf(f,"%lf%lf%lf",&an,&bn,&dn); fclose(f); c1=c1/b1; d1=d1/b1; For(i=2;i<=n-1;i+) t=bi-ci-1*ai; ci=ci/t; di=(di-di-1*ai)/t; dn=(dn-dn-1*an)/(bn-cn-1*an); for(i=n-1;i>=1;i-) di=di-ci*di+1; printf("n*n"); for(i=1;i<=n;i+) printf("d%2d=%l

15、fn",i,di);Zgf.txt文件中的內(nèi)容是：52 1 -71 2 1 -5 1 2 1 -5 1 2 1 -5 1 2 -5六、結(jié)果討論和分析Gauss消元法：LU分解法：追趕法：分析從消元過(guò)程可以看出，對(duì)于n階線性方程組，只要各步主元素不為零，經(jīng)過(guò)n-1步消元，就可以得到一個(gè)等價(jià)的系數(shù)矩陣為上三角形陣的方程組，然后再利用回代過(guò)程可求得原方程組的解.消元過(guò)程相當(dāng)于分解 A為單位下三角陣L與上三角陣U的乘積，解方程組Ly=b回代過(guò)程就是解方程組Ux=y。其中的L為n階單位下三角陣、U為上三角陣. 在 A 的LU 分解中, L取下三角陣, U 取單位上三角陣,這樣求解方程組Ax=d

16、 的方法稱(chēng)為追趕法.實(shí)驗(yàn)地點(diǎn)ZSA401指導(dǎo)教師李志學(xué)院名稱(chēng)軟件學(xué)院專(zhuān)業(yè)班級(jí)軟件1012學(xué)號(hào)學(xué)生姓名喬婧峰實(shí)驗(yàn)日期20114成績(jī)課程名稱(chēng)數(shù)值計(jì)算方法實(shí)驗(yàn)題目實(shí)驗(yàn)三 線性方程組的迭代解法一、課題名稱(chēng)線性方程組的迭代解法使用雅可比迭代法或高斯-賽德?tīng)柕▽?duì)下列方程組進(jìn)行求解。 二、目的和意義學(xué)習(xí)使用雅可比迭代法或高斯-賽德?tīng)柕ㄈ⒂?jì)算公式雅克比迭代法：設(shè)線性方程組Ax=b的系數(shù)矩陣A可逆且主對(duì)角元素a11,a22,ann均不為零，令D=diag(a11,a22,ann)并將A分解成A=(A-D)+D從而線性方程組可寫(xiě)成Dx=(D-A)x+b則有迭代公式x(k+1)=B1x(k)+f1其中，B

17、1=I-D-1A,f1=D-1b。四、主要儀器設(shè)備Vc2008,hp五、結(jié)構(gòu)程序設(shè)計(jì)高斯-賽德?tīng)柕ǎ?include "stdafx.h"#include <stdio.h>#include <math.h>void main() float a33=10,-1,-2,-1,10,-2,-1,-1,5,b3=7.2,8.3,4.2;float x3=0,0,0,sum1,sum2;int i,j,k,n=3;for (k=0;k<10;k+) for(i=0;i<n;i+) sum1=0; sum2=0;for(j=0;j<i;

18、j+) sum1=sum1+aij*xj; for(j=i+1;j<3;j+) sum2=sum2+aij*xj; xi=(bi-sum1-sum2)/aii; for(i=0;i<n;i+) printf("x%d=%f,",i+1,xi); printf("n"); 雅克比迭代：#include "stdafx.h"#include <stdio.h>#include <math.h>void main() float a33=10,-1,-2,-1,10,-2,-1,-1,5,b3=7.2,8

19、.3,4.2;float x3=0,0,0,sum1;int i,j,k,n=3;for (k=0;k<10;k+) for(i=0;i<3;i+) sum1=0; for(j=0;j<n;j+) if(i=j) continue; sum1=sum1+aij*xj; xi=(bi-sum1)/aii; for(i=0;i<n;i+) printf("x%d=%f,",i+1,xi);printf("n"); 六、實(shí)驗(yàn)結(jié)果與分析：高斯-賽德?tīng)柕ǎ貉趴吮鹊悍治觯菏褂酶咚?賽德?tīng)柡脱趴吮鹊伎梢郧蟪龇匠探M的解，但是利用高斯-

20、賽德?tīng)柕ㄋ璧牡螖?shù)比雅克比迭代少，能夠更早的達(dá)到精度要求。實(shí)驗(yàn)地點(diǎn)ZSA401指導(dǎo)教師李志學(xué)院名稱(chēng)軟件學(xué)院專(zhuān)業(yè)班級(jí)軟件1012班學(xué)號(hào)學(xué)生姓名喬婧峰實(shí)驗(yàn)日期20114成績(jī)課程名稱(chēng)數(shù)值計(jì)算方法實(shí)驗(yàn)題目實(shí)驗(yàn)四 最小二乘法擬合多項(xiàng)式一、課題名稱(chēng)（1）了解矩陣特征值與特征向量問(wèn)題解法，掌握冪法。（2）加深對(duì)矩陣特征值與特征向量問(wèn)題求解方法的認(rèn)識(shí)，掌握算法。（3）會(huì)進(jìn)行誤差分析。二、目的和意義 學(xué)習(xí)使用最小二乘法擬合多項(xiàng)式三、計(jì)算公式冪法：由已知的非零向量x0和矩陣A的乘冪構(gòu)造向量序列xn以計(jì)算矩陣A的按模最大特征值及其特征向量的方法，稱(chēng)為冪法。迭代公式：結(jié)果可取 四、主要儀器設(shè)備Vc2008,h

21、p五、結(jié)構(gòu)程序設(shè)計(jì) 五、結(jié)果討論和分析分析： 冪法是一種求任意矩陣A的按模最大特征值及其對(duì)應(yīng)特征向量的迭代算法。該方法的最大優(yōu)點(diǎn)是計(jì)算簡(jiǎn)單，容易在計(jì)算機(jī)上實(shí)現(xiàn)，對(duì)稀疏矩陣較為適合，但有時(shí)收斂速度很慢。實(shí)驗(yàn)地點(diǎn)綜合樓六層606室指導(dǎo)教師王崢學(xué)院名稱(chēng)計(jì)算機(jī)科學(xué)與技術(shù)專(zhuān)業(yè)班級(jí)計(jì)算機(jī)學(xué)號(hào)學(xué)生姓名某某實(shí)驗(yàn)日期2011-6-20成績(jī)課程名稱(chēng)數(shù)值計(jì)算方法實(shí)驗(yàn)題目實(shí)驗(yàn)五 代數(shù)插值一、 課題名稱(chēng)使用拉格朗日插值法或牛頓插值法求解：已知f(x)在6個(gè)點(diǎn)的函數(shù)值如下表所示，運(yùn)用插值方法，求f(0.596)的近似值。x0.400.550.650.800.901.05f(x)0.410750.578150.696750

22、.888111.026521.25386二、 目的和意義學(xué)習(xí)使用拉格朗日插值法或牛頓插值法求解三、 計(jì)算公式設(shè)函數(shù)在區(qū)間a,b上n+1互異節(jié)點(diǎn)x0,x1,xn上的函數(shù)值分別為y0,y1,yn,求n次插值多項(xiàng)式Pn(x),滿足條件Pn(xj)=yj, j=0,1,n令Ln(x)=y0l0(x)+y1l1(x)+ynln(x)= yili(x)其中l(wèi)0(x),l1(x), ln(x) 為以x0,x1,xn為節(jié)點(diǎn)的n次插值基函數(shù)，則Ln(x)是一次數(shù)不超過(guò)n的多項(xiàng)式，且滿足Ln(xj)=yj, L=0,1,n再由插值多項(xiàng)式的唯一性，得Pn(x)Ln(x)四、 結(jié)構(gòu)程序設(shè)計(jì)#include<st

23、dio.h>#include<stdlib.h>#include<iostream.h>typedef struct data float x; float y;Data; /變量x和函數(shù)值y的結(jié)構(gòu)Data d20; /最多二十組數(shù)據(jù)float f(int s,int t) /牛頓插值法，用以返回插商 if(t=s+1) return (dt.y-ds.y)/(dt.x-ds.x); else return (f(s+1,t)-f(s,t-1)/(dt.x-ds.x); float Newton(float x,int count) int n; while(1)

24、 cout<<"請(qǐng)輸入n值(即n次插值):"/獲得插值次數(shù) cin>>n; if(n<=count-1)/ 插值次數(shù)不得大于count1次 break; else system("cls"); float t=1.0; float y=d0.y; float yt=0.0; for(int j=1;j<=n;j+) t=(x-dj-1.x)*t; yt=f(0,j)*t; y=y+yt; return y;float lagrange(float x,int count) float y=0.0; for(int k=

25、0;k<count;k+)/這兒默認(rèn)為count1次插值 float p=1.0;/初始化p for(int j=0;j<count;j+) /計(jì)算p的值 if(k=j)continue;/判定是否為同一個(gè)數(shù) p=p*(x-dj.x)/(dk.x-dj.x); y=y+p*dk.y;/求和 return y;/返回y的值void main() float x,y; int count; while(1) cout<<"請(qǐng)輸入xi,yi的組數(shù)，不得超過(guò)20組:"/要求用戶(hù)輸入數(shù)據(jù)組數(shù) cin>>count; if(count<=20)

26、 break;/檢查輸入的是否合法 system("cls"); /獲得各組數(shù)據(jù) for(int i=0;i<count;i+) cout<<"請(qǐng)輸入第"<<i+1<<"組x的值:" cin>>di.x; cout<<"請(qǐng)輸入第"<<i+1<<"組y的值:" cin>>di.y; system("cls"); cout<<"請(qǐng)輸入x的值:"/獲

27、得變量x的值 cin>>x; while(1) int choice=3; cout<<"請(qǐng)您選擇使用哪種插值法計(jì)算："<<endl; cout<<" (0):退出"<<endl; cout<<" (1):Lagrange"<<endl; cout<<" (2):Newton"<<endl; cout<<"輸入你的選擇：" cin>>choice;/取得用戶(hù)的選擇項(xiàng)

28、 if(choice=2) cout<<"你選擇了牛頓插值計(jì)算方法，其結(jié)果為：" y=Newton(x,count);break;/調(diào)用相應(yīng)的處理函數(shù) if(choice=1) cout<<"你選擇了拉格朗日插值計(jì)算方法，其結(jié)果為：" y=lagrange(x,count);break;/調(diào)用相應(yīng)的處理函數(shù) if(choice=0) break; system("cls"); cout<<"輸入錯(cuò)誤!"<<endl; cout<<x<<&quo

29、t; , "<<y<<endl;/輸出最終結(jié)果 五、結(jié)果討論和分析分析：拉格朗日插值的優(yōu)點(diǎn)是插值多項(xiàng)式特別容易建立，缺點(diǎn)是增加節(jié)點(diǎn)是原有多項(xiàng)式不能利用，必須重新建立，即所有基函數(shù)都要重新計(jì)算，這就造成計(jì)算量的浪費(fèi)。實(shí)驗(yàn)地點(diǎn)綜合樓六層606室指導(dǎo)教師王崢學(xué)院名稱(chēng)計(jì)算機(jī)科學(xué)與技術(shù)專(zhuān)業(yè)班級(jí)計(jì)算機(jī)學(xué)號(hào)學(xué)生姓名某某實(shí)驗(yàn)日期2011-6-20成績(jī)課程名稱(chēng)數(shù)值計(jì)算方法實(shí)驗(yàn)題目實(shí)驗(yàn)六 最小二乘法擬合多項(xiàng)式一、課題名稱(chēng)給定數(shù)據(jù)點(diǎn)（xi ,yi），用最小二乘法擬合數(shù)據(jù)的多項(xiàng)式，并求平方誤差。xi00.50.60.70.80.91.0yi11.751.962.192.442.713

30、.00二、目的和意義1熟練運(yùn)用已學(xué)計(jì)算方法求解方程組2加深對(duì)計(jì)算方法技巧，選擇正確的計(jì)算方法來(lái)求解各種方程組3培養(yǎng)使用電子計(jì)算機(jī)進(jìn)行科學(xué)計(jì)算和解決問(wèn)題的能力三、計(jì)算公式建立正規(guī)方程組：（xij+k）ak=xijyi ,j=0,1,n 平方誤差：I=（akxik-yi）2四、結(jié)構(gòu)程序設(shè)計(jì)#include<iostream.h>#include<fstream.h>#define N 15double power(double &a,int n)double b=1;for(int i=0;i<n;i+)b*=a;return b;void Gauss();d

31、ouble XN,YN,sumXN,sumYN,aNN,bN,lNN,xN;void main()ofstream outdata;ifstream indata;double s;int i,j,k,n,index;cout<<"請(qǐng)輸入已知點(diǎn)的個(gè)數(shù)n="cin>>n;cout<<endl;cout<<"請(qǐng)輸入X和Y:"<<endl; /輸入給定數(shù)據(jù)for(i=0;i<n;i+)cout<<"X"<<i<<"="c

32、in>>Xi;sumX1+=Xi;cout<<"Y"<<i<<"="cin>>Yi;sumY1+=Yi;cout<<endl;cout<<"sumX1="<<sumX1<<"t"<<"sumY1="<<sumY1<<endl;cout<<"請(qǐng)輸入擬合次數(shù)index="cin>>index;cout<<endl;i=n;sumX0=i;for(i=2;i<=2*index;i+)sumXi=0;for(j=0;j<n;j+)sumXi+=power(Xj,i);cout<<"sumX"<<i<<