数值分析自主上机题

1、文档来源为：从网络收集整理、。出版本可编辑.欢迎下载支持2016级数值分析上机实践报告机械工程学院10吴臻标准题：迭代格式的比较3*设方程f(x)=x-3x-1=0有三个实根x1=1.8793,x2=-0.34727,x3=-1.53209现采用下面三种不同计算格式，求f&)=0的根X】或X?3x1x3i2、 X二3JjQx13、 x数学原理：简单迭代法是根据f&)=0这个方程，对其进行等价变换化为x=巾&)并由该式产生逼近解x*的迭代数列xk,这就是简单迭代法的思想。对于同一f6)=0可以等价变换出不同的x二2&)而且它们各自的收敛性不同。程序设计：#hehd

2、e<bstream>#hchde<cmalh>#hchde<cstdb>ushgnamespacestd;doubtf(doubhJ外调函数f(x),每次更新新的函数区以第一种迭代方式为例子doublek,m,sum;k=3*B-l;m=pow(i2.0);sum=k/in;retLimsum;htmah0doublex,x0;htN最大迭代次数htk;couK”输入初解：”；ch>>x0;cout<<“输入最大迭代次数二cri>>N;fcr(k=lk<=N火+)2文档来源为：从网络收集整理.雷。圮版本可编辑.欢迎下

3、我支持文档来源为：从网络收集整理.wold版本可编辑.欢迎下载支持.x=f(xO);if(febs(x-x0)<0.0000001)(couK<”迭代次数:*<<k<<endl;coutC”输出得到的解：“<<x<<end1;system("pause。；return0;)ekexO=x;)couKU已达到最大迭代次数：<N«endl;couK”输出得到的解："<<x<endl;system("pause)reUlin0:)实验结果：程序运行结果讨论和分析：对于第一种迭代

4、格式，收敛区间田.2-0.4,在该收敛区间内迭代收敛于-1.53209,只能求得方程的一个根：对于第二种迭代格式，收敛区间H.51.8,在该收敛区间内迭代收敛于P.34730,同样只能求得方程的一个根；对于第三种迭代格式，收敛区间E0.3+8),在该收敛区间内迭代收敛于1.87937,只能求得方程的一个根：由以上结果很容易发现，初值的选取对迭代敛散性有很大影响。以第一种迭代格式为例，当初值大于等于F.3时，迭代格式发散；当初值小于等于七.3时，迭代格式也发散：只有初值在F.3和十.3之间时，迭代格式才收敛于-1.53209o其他迭代格式也有这样的性质，即收敛于某个数值区间，超出这个区间迭代格式

5、就是发散的，这就是所谓迭代格式的收敛性。自主题：机械运动的数值仿真背景：描述物理学里把物体位置的变化叫机械运动。如我们所知，力的作用效果有：改变物体的运动状态改变物体的形状改变物的运动状态大多会引起物体的位置变化，引起机械运动。改变物体的形状而不改变它的运动状态就叫是非机械运动中的一种。在工程实际与生产生活中我们常常需要对特定对象的机械运动进行研究分析得出其机械运动的规律，然后将该规律应用于对人们有益的方面。实例：如下图1(左)所示，假设有一烟花火笳，其初始条件为零。将其放在地方然后点火，该烟花火箭的初始质量为mJ20g,其中粉末燃料占7°g。经过实验得知，燃料的持续时间为120s。

6、燃料所产生的恒定推力为丁5.2N。这也说明燃料的消耗率恒定。主气小生的阻力和烟花火箭的速度的平方成正比：Rkv.k401°'(N*/nf这里，要求选择一种数值方法对其运动过程进行仿真并且其截断误差为°卜)或者更高。要求计算出该烟花火箭的最高高度，同时计算出从燃料消耗到该烟花火箭运动到最高点的时间延迟。数学原理及数学模型：该实际问题要求其截断误差要求大于或等于°九)，这就2文档来源为：从网络收集整理.甫。rd版本可编辑.欢迎下我支持文档来源为：从网络收集整理.嘴。Id版本可编辑.欢迎下我支持.使得较为简单的欧拉法，中点法不适合本例。龙格彝塔法以其优异的数值特

7、性成为解决本问题的首选。图1烟花火箭的机械运动的数值仿真（左）和结构力学问题的数值求解（右）很显然，该问题属于变质量的运动学问题，在该运动过程中，其前两秒是在驱动力和阻力的共同作用下加速上升的，而后的时间内，该烟花火箭是在空气的阻力下减速上升的，同时注意到空气的阻力和速度的平方成正比。为了对该运动过程进行数值仿真，那么必须建立相应的微分方程组。分析该运动过程可知，应该将该运动过程分为两部分：加速上升过程和减速上升过程，从而得到相应的微分方程组。加速上升过程：lOOOkv2 g 120 35tdv1000Ti-dt12035tdhividtvi(0)O,hi(O)(1)减速上升过程:dh2V2d

8、tdv2lOOOkv2dt120352h2（0）hl（2）V2（0）vi<2）式中，11为上升的高度，v为上升过程的速度，g为重力加速度。h/2）表示加速上升过程的最终高度，表示加速上升过程的最终速度。使用龙格嘛塔法求解如上的微分方程组。该系统的数值仿真结果如图2所示。图2（左）是该烟花火笳的上升过程高度的数值仿真，图2（右）是其上升过程速度的数值仿真。同时亦可以得到烟花火箭上升的最大高度和问题中所需的时间延迟：从如上的分析和仿真可知，使用数值方法进行机械运动的数值仿真，可以简化本身复杂变化的物理运动过程。同时注意到，即使改变机械运动的初始条件或者系统的某些特征，其对应的仿真只需要进行简

9、单的调整。而且其仿真精度也可以得以预见。图2烟花火箭的上升过程高度（左）和速度（右）的数值仿真编程过程：加速上升过程ffe_,y_F5200/（120H35t）-0.4x2/（1205t）-10;g反一,y_尸x,氏y:0,0;h=0.1;3文档来源为：从网络收集整理.雷。圮版本可编辑.欢迎下我支持文档来源为：从网络收集整理 万。出版本可编辑.欢迎下载支持.t=0.1;xx=Table0,il,70;yy=Tabfe0,il,7O;tt=Tabfe0,11,70;Dofe=flx,ya=x+h(a+ftx+h,y+h*a)/2;b=gk,yya=y+h(b+gtk+h,y+h牝M;Prhtk,

10、"",t<xa,ya;t,x,y,xxlkll,yyHklttlM=t+h,xa,ya,xa,ya,t+h,kJ,201O.13.463520-00520.27.06092O.3S3S230.31O<7913工.0674440.414.65262<1515750.S18.64143«621836。622.75255.4909870.72<987.7712680-831.314510-474390.935.745613.6107101.40«260317.1903111.144.843121.2213121.249.476325.7

11、10131.354.139230.632141.458-808836.0822151.53.459541.960X.646.319171.772.S91255.1303iei.e77.01112.3945191.981.290770.100620Z.397778.Z347tt0.6,1.1,1.6,2.1,2.6,3.1,3.6,4.1,4.6,5.1,5.6,6.1,6.6,7.1,7.6xx17.3174,37.9564,61.2759,85.0712,104.674,113.404,109.842,75.5189,8.73997-59.2801,-96.8113,-111.485,-11

12、7.825,-121.341,-123.752yy0,125,8.90869,28.0119,58.7749,101.435,153.898,210.725,265.771,303.655,308PohlSize 0.021PohlSizeK).02lT meP ohlS ize 0.02 True.15,278.635,230.354,174.737,115.949,55.4041LislPbtlJabfettafflxxEn,il52OlPblStyhLislPbtJabfettlEllxxEID,il,20)lPbtStyleLis#btUabhttElxxE|),il,20)lPbtJ

13、oriedListPbtUabbttlEIlyyllill,il,20)lPblStyhListPbtlJabfettlElyylEI,il,20lPbtjohed减速上升过程f(x_,y_F-0.4x-2/(120-35*2)-10;g反一,y_三x；氐y二也5.3977,78.2347;h=0.1;4文档来源为：从网络收集整理.雷。圮版本可编辑.欢迎下我支持工档火源为：从网络收集整理.WOld版本可编辑.欢迎下段支孙.1=2.1;Dofe=ffx,ya=x+h(a+ffk+h,y+h*a/2;b=glx,yya=y+h(b+g反+h,y+h均；Prhtlk/",t”*xa/*,yal;tx,y,xxQ20+klyyQ2O+kltt02O+k;=tth,xa>ya,xa>ya,t*-h,kJ,50)LislPbtJabettttQlxxED,i1,61,PhlStylePohtSize0.02ListPbtllabhttffillyyllill),il,61,PhtStykPohtSize8).02LislPb