《西安交大計(jì)算方法.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、.西安交通大學(xué)計(jì)算方法上機(jī)實(shí)驗(yàn)班級：(xxx)姓名：(xxx)學(xué)號：2111601004教育資料.1.按兩種順序計(jì)算y，哪個(gè)接近真值？Y=1000+++…+用java語言編寫：publicclassAdd{publicstaticvoidmain(String[]args){doubles=0,y=1000;for(doublea=1001.0;a<=2000.0;a++){y+=1.0/a;}for(doublea=2000.0;a>=1001.0;a--){s+=1.0/a;}s=s+1000;S

2、ystem.out.println("正序和"+s);System.out.println("逆序和"+y);}}運(yùn)行結(jié)果：結(jié)論：顯然假設(shè)是double類型的數(shù)據(jù)時(shí)，先算大數(shù)的過程吃掉了末尾的小數(shù)被進(jìn)位所埋沒，導(dǎo)致了大數(shù)吃小數(shù)的誤差，按從小到大（從右向左）的計(jì)算順序所得的結(jié)果與真值相近，而按從大到?。◤淖蟮接遥┑挠?jì)算順序所得的結(jié)果與真值的誤差較大。教育資料.1-18.設(shè)(x)=1+x+++…+,計(jì)算(-5)和1/(5)，哪個(gè)接近？解法一：用JAVA語言編寫：publicclasssecond{publ

3、icstaticvoidmain(String[]args){doubles1=1,s2=1;doublee=1,sum=1;//e的初值為1，sum用來存放n！inta=1;while(sum

4、i);}s1=1.0/s1;System.out.println("1/S24(5)="+s1);System.out.println("S24(-5)="+s2);}publicstaticdoublecimi1(intai){doublexi=1;for(inti=ai;i>=1;i--){xi=xi*(5.0/i);}returnxi;}publicstaticdoublecimi2(intai){doublexi=1;for(inti=ai;i>=1;i--){xi=xi*(-5.0/i);}

5、returnxi;}教育資料.}運(yùn)行結(jié)果：解法二：用matlab編程并運(yùn)行，如下：(1)計(jì)算(-5)運(yùn)行結(jié)果如下：(2)計(jì)算1/(5)運(yùn)行結(jié)果如下：而的真是結(jié)果為0.006737946比較得1/(5)的計(jì)算結(jié)果與真實(shí)值更接近解法三：也可以用C++編寫：#include"stdafx.h"教育資料.#include"stdio.h"#include"iostream"usingnamespacestd;intmain(intargc,char*argv[]){intfunc1(int);doublefu

6、nc2(int);doubley=0;inti;for(i=1;i<25;i++){intz=func1(i);doublee=func2(i);y+=z/e;}cout<<"----------------------------------------"<

7、ntx){inty=1;intk;for(k=0;k

8、可將真解表示為G=0.45a--0.1125=[(0.2a–0.0625)b+(0.25a–0.1)c+(0.125-)d]/G=-[(0.1–0.25a)b+(-0.2)c+(0.25–0.5a)d]/G=[(0.125-)b+(0.5a–0.25)c+(a–0.25)d]/G試將a，b，c，d依次取2～6位近似有效數(shù)，分別計(jì)算，，。解：用matlab編程如下：將a，b，c，d依次取2～6位近似有效數(shù)取2位有效近似有效數(shù)教育資料.取3位有效近似有效數(shù)取