欢迎来到天天文库
浏览记录
ID:24180959
大小:101.15 KB
页数:4页
时间:2018-11-13
《数值分析题目》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、第一章误差1.写程序计算Stirling近似:的绝对误差和相对误差。分别W论其绝对误差和相对误差随n增加的变化情况。2.利用公式e=lim(l+l/nr求e的近似值。计算近似值与真值的误差,讨论误差随n增加的变化情况,试对你计算的结果作出解释。3.(选做)分别求lim(e"-l)/x,lim(ex-l)/log(ex),可取x=10'ifc=1,2.根据你的计算结果,那个公式给出较好的近似,为什么?4.(选做)用如下迭代公式求山笫n项的值:=2.25x,-0.5x^_,x,=1/3,x2=1/12画出计算值关于k的曲线(纵轴
2、可用对数值画出),并与真值=4^73比较。5.(选做)用如下迭代公式求出第n项的值:x,+1=111-(1130-3000/%,_,)/xA=11/2,x2=61/11其真值为&单调增加收敛到6,你能否解释你所得到的计算结果?(木章习题可以用Matlab完成。)第二章插值法与数值微分1.对一个给定系数的多项式p(x),利用Horner公式编写如下程序:1)在给定点x*求函数值p(x*)2)在给定点x*求呆数值p’(x*)3)在给定税分区间la,bj上求p(x)的定积分.2.编写用Lagrange型和Newton型插位多项式
3、求函数在某个给定点的近似位的程序,并比较这两种算法的计算量。已知:來个函数在n+1个不同甘点x,.上的值%,及某个求值点J1.(选做)对定义在区间[-5,5]上的函数:取等距节点,对不M的结点数,用你所编写的Newton插伉多项式程序给出近似函数,M图比较。对你的结果给出解释。荞插值节点取为区间[-5,5]上的切比雪夫点,可得W—个近似函数,同样M阁比较并解释。(区间La,bj上的切比雪夫点定义为:xkb-a24l2("+l))第三章函数逼近与数据拟合某行星在一个椭圆轨道上运动,可以川如不方程描述:ay2+bxy+ex+dy
4、+e=x2请根裾如下观测数据,确定参数a,b,c,e,d的取伉:X1.020.950.870.770.670.560.440.300.160.01y0.390.320.270.220.180.150.130.120.130.15根据你得到的参数在(x,y)平血上W出轨迹曲线及上列观测点。对表中数据任意做小扰动(如将y(l)=0.39改为y(1)=0.389等),用同样方法再求轨迹方程,并在M—张阁上画出相应曲线,讨论你所发现的现象。(因木题涉及到方程组的求解,可以用Matlab完成。)第五章数值积分根裾公式利用数值求积公式求
5、n的近似值1)分别编写父化中点公式、复化梯形公式和复化抛物线公式的程序,对不同的步长h求冗的近似值及相应误差。画ili三种算法的误差关于步L<:h的变化曲线,误差是否随着h的减小总是下降的?为什么?2)编写Romberg求积公式的程序,对相同误差吋与1)中各算法的步长进行比较。第十章非线性方程及非线性方程组解法1.编写二分法、Newton法、弦位法的算法程序,tl选算例进行计算,比较三种格式的计算量与收敛速度。1.(选做)Newton下山法的算法程序。第十一章常微分方程初值问题的数值解法1.编写4阶经典R-K法和4阶Adam
6、s预佔-校正算法程序,并求解常微分方程初值问题,22y=(——l)(;y-sin(;zx))-ex+^cos(^x)(17、+=234kEk=i(-1),则可取ln2-S/r方法2:利用方法1中的数据,按下列公式生成新数列则可取In2-Sz,方法3:柵腳心去六k&2=i=ik.2k第三章插值法与数值微分样条插值的应用川三次样条插值阑数去逼近飞机头部的外型曲线,其型值点数据由下表给出。X0701302103375787761012114214621841y05778103135182214244256272275设给定下述两种边界条件,试分别计算插值函数在点人=50々(/:=1,2,...,36)的值,并M出飞机头部外型曲线。(1)自然边界条件;(8、2)/⑼=1,/(1841)=0
7、+=234kEk=i(-1),则可取ln2-S/r方法2:利用方法1中的数据,按下列公式生成新数列则可取In2-Sz,方法3:柵腳心去六k&2=i=ik.2k第三章插值法与数值微分样条插值的应用川三次样条插值阑数去逼近飞机头部的外型曲线,其型值点数据由下表给出。X0701302103375787761012114214621841y05778103135182214244256272275设给定下述两种边界条件,试分别计算插值函数在点人=50々(/:=1,2,...,36)的值,并M出飞机头部外型曲线。(1)自然边界条件;(
8、2)/⑼=1,/(1841)=0
此文档下载收益归作者所有