浮点数的误差.doc

《浮点数的误差.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、一、浮点数与舍入误差1.画出多项式的图像，其中x取值分别为：(1)x=0.988:0.0001:1.012，(2)x=-50:0.01:50。>>x=0.988:.0001:1.012;>>y=x.^7-7*x.^6+21*x.^5-35*x.^4+35*x.^3-21*x.^2+7*x-1;>>plot(x,y)>>x=-50:0.01:50;>>y=x.^7-7*x.^6+21*x.^5-35*x.^4+35*x.^3-21*x.^2+7*x-1;>>plot(x,y)2.当,时，用Matlab计算。>>A=ones(1,10.^3)

2、*0.1;>>m=10.^2;>>sum(A)-mans=-1.4069e-012>>>>A=ones(1,10.^4)*0.1;>>m=10.^3;>>sum(A)-mans=1.5882e-010>>A=ones(1,10.^5)*0.1;>>m=10.^4;>>sum(A)-mans=1.3169e-009>>A=ones(1,10.^6)*0.1;>>m=10.^5;>>sum(A)-mans=-8.9420e-0073.利用Matlab中的函数rand分别生成(0,1)区间上两组各为n个的随机数，记为，。考虑如下几个和，；，。取

3、n充分大（如），取与n同样量级，多次进行上述实验并记录实验结果。>>A=rand(1,10.^6);>>B=rand(1,10.^6);>>m=10.^10*rand(1,10.^6);>>s1=sum(A)*sum(1,10.^6);>>s2=sum(A+m.*B);>>s3=m(1)+sum(A);>>s4=sum(A)+m(1);>>s1_2=s1-s2s1_2=-2.5034e+015>>s3_4=s3-s4s3_4=04.函数有幂级数展开编写利用幂级数计算的Matlab程序，对于，计算的精度是多少？分别需要计算多少项？funct

4、ions=powersin(x)%powersin.powerseriesforsin(x).%powersin(x)triestocomputesin(x)fromapowerseriess=0;t=x;n=1;whiles+t~=s;s=s+t;t=-x.^2/((n+1).*(n+2)).*t;n=n+2;end>>powersin(pi/2)ans=1.0000二、计算方法及其计算复杂性秦九韶算法是中国南宋时期的数学家秦九韶提出的一种多项式简化算法。在西方被称作霍纳算法（Horneralgorithm或Hornerscheme），

5、是以英国数学家威廉·乔治·霍纳命名的。　　把一个n次多项式改写成如下形式：该算法看似简单，其最大的意义在于将求n次多项式的值转化为求n个一次多项式的值。在人工计算时，利用秦九韶算法和其中的系数表可以大幅简化运算；对于计算机程序算法而言，加法比乘法的计算效率要高很多，秦九韶算法将原来的做次乘法，缩短到做n次乘法，减少了CPU运算时间。试用Matlab编写秦九韶算法程序，并用该程序计算多项式在x=1.1,1.2,1.3的值。三、算法的稳定性1.分别编写以下两种算法的程序，并记录比较实验结果，不稳定的算法：稳定的算法：%²»Îȶ¨µÄËã·¨

6、E=zeros(1,20);E(1)=1/exp(1);forn=2:20;E(n)=1-n*E(n-1);endE%²»Îȶ¨µÄËã·¨I=zeros(1,20);I(20)=0;forn=20:-1:2;I(n-1)=(1-I(n))/n;endIE=Columns1through80.36790.26420.20730.17090.14550.12680.11240.1009Columns9through160.09160.08390.07740.07180.06690.06270.05900.0566Columns17thro

7、ugh200.03740.3259-5.1930104.8608I=Columns1through80.36790.26420.20730.17090.14550.12680.11240.1009Columns9through160.09160.08390.07740.07180.06690.06270.05900.0557Columns17through200.05280.05000.05000、2.分别编写以下两种算法的程序，画出散点图，不稳定的算法：稳定的算法：%²»Îȶ¨µÄËã·¨clearE(1)=1;forn=2:19;E

8、(n)=exp(1)-n*E(n-1);end%Îȶ¨µÄËã·¨I(19)=0;forn=19:-1:2;I(n-1)=(exp(1)-I(n))/n;endsubplot(2,1,1)p