资源描述：
《数值分析自主上机题》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、2016级数值分析上机实践报告机械工程学院2016021910吴臻标准题：迭代格式的比较设方程f(x)=x3-3x-l=0有三个实根x；=1.8793,x；二0.34727,x3二1.53209现采用下面三种不同计算格式，求f(x)=0的根X；或x；3x4-11、X二/疋-12、x=33、x="%+1数学原理：简单迭代法是根据f(x)=0这个方程，对其进行等价变换化为x=it(x)并由该式产生逼近解/的迭代数列｛双｝,这就是简单迭代法的思想。对于同一f(x)=0可以等价变换出不同的x=F(x)而且它们各自的收敛性不同。程序设计：
2、#include#include#includeusingnamespacestd;doublef(doublei)〃外调函数f(x),每次更新新的函数｛〃以第一种迭代方式为例子doublekmsum;k=3*i+l;m=pow(i?2.0);sum=k/m;returnsum;｝intmain()｛doublex,x0;intN;〃最大迭代次数intk;cout«n输入初解:”;cin»x0;cout«n输入最大迭代次数:”；cin»N;for(k=1;k<=N;k++)x=
3、f(xO);if(fabs(x-x0)<0.0000001)cout«H迭代次数:H«k«endl;cout«n输出得到的解:H«x«endl;system(,,pause");return0;}elsex0=x;}cout«nL2达到最大迭代次数:H«N«endl;cout«H输出得到的解:,,«x«endl;system("pause");return0;}实验结果：输入初解：-1.5输入最大迭代次数汚0迭代次数：42输出得到的解=-1-53209输入初解=~0-3输入最大迭代次数汚0已达到最大迭代次数汚0输出得到的解=-1
4、-#ind请按任意键继缠・・・输入初解：-0-3输入最大迭代次数汚0输出得到的解：-0-347296请按任倉键继线・・・输入初解：1・5输入晶大迭代次数汚0送代次数心输出得到的解：0.347296请按任倉键继续・・・输入初解：1・5蔬入最大迭彳弋次数：50送代菠数汽3输出得到的解汽.87939请按任倉键继续・・・输入初解：2・0输入最大迭代次数：50迭代次数汉2输出得到的解=1-87939请按任意键继线・・・程序运行结果讨论和分析：对于第一种迭代格式，收敛区间卜8.204],在该收敛区间内迭代收敛于-1.53209,只能求得方程
5、的一个根；对于第二种迭代格式，收敛区间[・1.51.8],在该收敛区间内迭代收敛于・0.34730,同样只能求得方程的一个根；对于第三种迭代格式，收敛区间卜0.3+8）,在该收敛区间内迭代收敛于1.87937,只能求得方程的一个根；由以上结果很容易发现，初值的选取对迭代敛散性有很大影响。以第一种迭代格式为例，当初值大于等于・0.3时，迭代格式发散；当初值小于等于・&3时，迭代格式也发散；只有初值在・0.3和・8.3之间时，迭代格式才收敛于一1.53209o其他迭代格式也有这样的性质，即收敛于某个数值区间，超出这个区间迭代格式就是
6、发散的，这就是所谓迭代格式的收敛性。自主题：机械运动的数值仿真背景：描述物理学里把物体位置的变化叫机械运动。如我们所知，力的作用效果有：改变物体的运动状态改变物体的形状改变物的运动状态大多会引起物体的位置变化，引起机械运动。改变物体的形状而不改变它的运动状态就叫是非机械运动中的一种。在工程实际与生产生活中我们常常需耍对特定对象的机械运动进行研究分析得出其机械运动的规律，然后将该规律应用于对人们有益的方面°实例：如下图1（左）所示，假设有一烟花火箭，其初始条件为零。将其放在地方然后点火，该烟花火箭的初始质量为"严120g,其中粉末
7、燃料占70g。经过实验得知，燃料的持续时间为—2"。燃料所产生的恒定推力为T=5.2No这也说明燃料的消耗率恒定。空气产生的阻力和烟花火箭的速度的平方成正比：R*24.0・107（砒/莎）。这里，要求选择_•种数值方法对其运动过程进行仿真并且其截断误差为或者更高。要求计算出该烟花火箭的最高高度，同时计算出从燃料消耗到该烟花火箭运动到最高点的时间延迟。数学原理及数学模型：该实际问题要求其截断误差要求大于或等于°（力4）,这就使得较为简单的欧拉法，中点法不适合本例。龙格■库塔法以其优异的数值特性成为解决本问题的首选。I(I—I)图1
8、烟花火箭的机械运动的数值仿真（左）和结构力学问题的数值求解（右）很显然，该问题属于变质量的运动学问题，在该运动过程中，其前两秒是在驱动力和阻力的共同作用下加速上升的，而后的时间内，该烟花火箭是在空气的阻力下减速上升的，同时注意到空气的阻力和速度的平方成正比。为了
