西安交大计算方法

《西安交大计算方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、西安交通大学计算方法上机实验班级：(xxx)姓名：(xxx)学号：21116010041.按两种顺序计算y，哪个接近真值？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;System.out.println("正序和"+s);System.out.p

2、rintln("逆序和"+y);}}运行结果：结论：显然假设是double类型的数据时，先算大数的过程吃掉了末尾的小数被进位所埋没，导致了大数吃小数的误差，按从小到大（从右向左）的计算顺序所得的结果与真值相近，而按从大到小（从左到右）的计算顺序所得的结果与真值的误差较大。1-18.设(x)=1+x+++…+,计算(-5)和1/(5)，哪个接近？解法一：用JAVA语言编写：publicclasssecond{publicstaticvoidmain(String[]args){doubles1=1,s2=1;doublee=1,sum=1;//e的初值为1，sum用来

3、存放n！inta=1;while(sum

4、tai){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);}returnxi;}}运行结果：解法二：用matlab编程并运行，如下：(1)计算(-5)运行结果如下：(2)计算1/(5)运行结果如下：而的真是结果为0.006737946比较得1/(5)的计算结果与真实值更接近解法三：也可以用C++编写：#include"stdafx.h"#inc

5、lude"stdio.h"#include"iostream"usingnamespacestd;intmain(intargc,char*argv[]){intfunc1(int);doublefunc2(int);doubley=0;inti;for(i=1;i<25;i++){intz=func1(i);doublee=func2(i);y+=z/e;}cout<<"----------------------------------------"<

6、----------------------------------------"<

7、的真解===1，且令a=,b=,c=,d=时利用克莱姆法则可将真解表示为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位近似有效数，分别计算，，。解：用matlab编程如下：将a，b，c，d依次取2～6位近似有效数取2位有效近似有效数取3位有效近似有效数取4位有效近似有效数取5位有效近似有效数取6位有效近似有效数结论