《数值分析上机题.doc》由会员分享，可在线阅读，更多相关《数值分析上机题.doc（16页珍藏版）》请在金锄头文库上搜索。
1、数值分析上机题姓名：陈作添 学号：040816习题120（上机题）舍入误差与有效数设，其精确值为。（1）编制按从大到小的顺序，计算的通用程序。（2）编制按从小到大的顺序，计算的通用程序。（3）按两种顺序分别计算，并指出有效位数。（编制程序时用单精度）（4）通过本上机题，你明白了什么？第 1 页 共 16 页按从大到小的顺序计算的通用程序为：#includefloat sum(float N)float j,s,sum=0;for(j=2;j=N;j+) s=1/(j*j-1);sum+=s;return sum; 按从小到大的顺序计算的通用程序为：#includefloat sum(float
2、 N)float j,s,sum=0;for(j=N;j=2;j-)s=1/(j*j-1);sum+=s;return sum;从大到小的顺序的值从小到大的顺序的值精确值有效位数从大到小从小到大0.7400490.740050.740049650.7498520.74990.7499440.7498520.7499990.74999936通过本上机题，看出按两种不同的顺序计算的结果是不相同的，按从大到小的顺序计算的值与精确值有较大的误差，而按从小到大的顺序计算的值与精确值吻合。从大到小的顺序计算得到的结果的有效位数少。计算机在进行数值计算时会出现“大数吃小数”的现象，导致计算结果的精度有所降低
3、，我们在计算机中进行同号数的加法时，采用绝对值较小者先加的算法，其结果的相对误差较小。习题220（上机题）Newton迭代法（1）给定初值及容许误差，编制Newton法解方程根的通用程序。（2）给定方程，易知其有三个根，。1由Newton方法的局部收敛性可知存在，当时，Newton迭代序列收敛于根。试确定尽可能大的。2试取若干初始值，观察当，时Newton序列是否收敛以及收敛于哪一个根。（3）通过本上机题，你明白了什么？解：（1）编制的通用程序：#include#include#define eps 0.000001 /给定容许误差float f(float x) /定义函数f(x)float
4、 f;f=x*x*x/3-x; /f(x)的表达式;return(f);float df(float x) /定义函数df(x)，计算f(x)的导函数float df;df=x*x-1; /f(x)导函数的表达式;return (df);void main(void)float x0,x1,a;int k=0;coutx0;doa=-f(x0)/df(x0);x1=x0+a;k+;x0=x1;while(fabs(a)eps);coutktx0; /输出迭代的次数和根值（2）计算迭代序列收敛于根的尽可能大的的函数为：#include#includevoid delay(int n) /定义延时
5、函数for(n=10000;n0;n-);#define eps 0.000001float f(float x) /定义函数f(x)float f;f=x*x*x/3-x; /f(x)的表达式;return(f);float df(float x) /定义函数df(x)，计算f(x)的导函数float df;df=x*x-1; /f(x)导函数的表达式;return (df);int judgement(float z)int count=5;float x0,x1,type,type1;x0=z;while(count-0) x1=x0-f(x0)/df(x0);type=fabs(x1)
6、; type1=fabs(x1-x0); /调试值用coutcount=countttype=typettype1=type1n;if(fabs(x1-x0)eps)return 1;x0=x1;delay(30000); /调试值用 return 0; void main(void)float delta=0;int flag=1;while(flag=1)cout方程的根为:n;delta+=eps;flag=judgement(delta);cout输出方程根收敛的区间值:n;coutdelta-eps; /输出收敛的区间值当步长为0.001时，程序计算出的的为=0.774，即在区间（-
7、0.774，0.774）内迭代序列收敛于0。对于不同得初始值收敛于不同的根, 在（-，-1）内收敛于，在（-0.774，0.774）内收敛于，在（1，+）内收敛于，但在内（0.774，1）和（1，0.774）均可能收敛于和。，分别为方程的精确解。分析：对于不同的初值，迭代序列会收敛于不同的根，所以在某个区间内求根对于初值的选取有很大的关系。产生上述结果的原因是区间不满足大范围收敛的条件。习题335（上机题）列主元三角分解法对于某电路的分析，归结为求解线性方程组RI=V。（1）编制解n阶线性方程组Ax=b的列主元三角分解法的通用程序；（2）用所编制的程序解线性方程组RI=V，并打印出解向量，保留
8、五位有效数；（3）本编程之中，你提高了哪些编程能力？程序为：#include#includevoid main(void)int i,j,n,k,q;float a1011,s10,s110;coutn;cout输入数组a:endl;for(i=1;i=n;i+)for(j=1;jaij; /给矩阵a赋值for(i=1;i=n;i+)for(j=1;j=(n+1);j+)coutaijt;coutn; /输出数组acoutn;/进行第一行和第一列元素的求取/int t=1;for(i=1;i=n;i+)si=ai1;float max=fabs(s1);for(i=2;imax)max=fab
9、s(si);t=i;for(j=1;j=(n+1);j+)float b=a1j;a1j=atj;atj=b; /进行第一列主元互换for(i=2;i=n;i+)ai1=ai1/max; /第一列除以a11for(i=1;i=n;i+)for(j=1;j=(n+1);j+)coutaijt;coutn;/输出进行第一步变换的数组acoutn;/进行第k步分解/for(k=2;k=n;k+)for(i=k;i=n;i+)float sum=0;for(q=1;qk;q+)sum+=aiq*aqk;s1i=aik-sum;int l=k;float m=fabs(s1k);for(i=k;im)m
10、=fabs(s1i);l=i; /返回行值for(j=1;j=n+1;j+) /交换两行元素float s2=akj;akj=alj;alj=s2;for(j=k;j=n+1;j+)/算出第k行行元素的值float sum1=0;for(q=1;qk;q+)sum1+=akq*aqj;akj=akj-sum1;for(i=k+1;i=n;i+)/算出第k列列元素的值float sum2=0;for(q=1;qk;q+)sum2+=aiq*aqk;aik=(aik-sum2)/(akk); /第k步分解结束for(i=1;i=n;i+)for(j=1;j=(n+1);j+)coutaijt;co
11、ut=1;i-)xn=ann+1/ann;float sum3=0;for(j=i+1;j=n;j+)sum3+=aij*xj;xi=(ain+1-sum3)/aii; /回代过程for(i=1;i=n;i+)coutxi=xiendl; /输出解向量结果：方程的解为：x1= -0.28923，x2= 0.34544，x3= -0.71281，x4= -0.22061，x5= -0.43040，x6= 0.15431，x7= -0.057823，x8= 0.20105，x9= 0.29023。分析：我感觉是提高了查错误点的能力和编循环语句的能力，即利用很规整的迭代公式进行编程。另外列主元三角分解法的阶梯步骤有了更深的了解!36逐次超松弛迭代法（1）编制解n阶线性方程组Ax
