资源描述：
《长沙理工大学数值分析实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、长沙理工大学数值分析实验报告心得。实验一、误差分析一、实验目的1、通过上机编程，复习巩固以前所学程序设计语言及上机操作指令；2、通过上机计算，了解误差、绝对误差、误差界、相对误差界的有关概念；3、通过上机计算，了解舍入误差所引起的数值不稳定性。二、实验原理误差问题是数值分析的基础，又是数值分析中一个困难的课题。在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果。因此，选取算法时注重分析舍入误差的影响，在实际计算中是十分重要的。同时，由于在数值求解过程中用有限的过程代替无限的过程会产生截断误差，因此算法的好坏会影响到数值结果的
2、精度。三、实验任务对n?0,1,2,?,20，计算定积分xnyn??dx.x?501算法1：利用递推公式yn?11?5yn?1,n?1,2,?,20,n取y0??1dx?ln6?ln5?0.182322.x?5算法2：利用递推公式yn?1?11?ynn?20,19,?,1.5n5注意到1120x2011,??xdx??dx??x20dx?12660x?5501050111取y20?111(?)?0.008730.20105126思考：从计算结果看，哪个算法是不稳定的，哪个算法是稳定的。四、实验源程序及结果算法1functiony1(i)formatlo
3、ngy1=zeros(1,21);y1(1)=0.182322;fori=1:1:20y1(i+1)=1/i-5*y1(i);endy1(1:21)ans=1.0e+007*Columns1through40.0000000182322000.000000008839000Columns5through80.0000000034583330.000000002708333Columns9through120.000000019196429-0.000000084871032Columns13through160.000010821757756-0.00
4、0054101096474Columns17through200.006762788547315-0.033813936854221Column214.2267422243506750.0000000058050000.0000000043083330.000000003125000-0.0000000013392860.000000434355159-0.0000021626848850.000270512625226-0.0013525564594630.169069689826659-0.845348443870135算法2functiony2(
5、i)formatlongy2=zeros(1,21);y2(21)=0.008730;fori=21:-1:2y2(i-1)=0.2/(i-1)-0.2*y2(i);endy2(1:21)ans=Columns1through40.1823215567939550.0883922160302270.0580389198488660.043138734089004Columns5through80.0343063295549800.0284683522251020.0243249055411550.021232615151370Columns9throu
6、gh120.0188369242431500.0169264898953620.0153675505231920.014071338293132Columns13through160.0129766418676730.0120398675847140.0112292335050020.010520499141658Columns17through200.0098975042917100.0093360079532160.0088755157894740.008254000000000Column210.008730000000000从结果看算法2比较稳
7、定，算法1不稳定。篇二：数值分析实验报告3实验名称平方逼近与最小二乘法式.2.分别计算(ki,kj)??(k,k)，(k,y)??(k,y)ijiijm?1m?1n?n3.用(ki,kj)组成方阵A,用(ki,y)组成矩阵B.4,利用A/B求出该多项式的系数，再利用实验题目1,对于给函数f(x)=i?1n1在区间【-1，1】上取xi=-1+0.2i（i=0，1…10），试求3次曲线21?25x拟合，试画出你和曲线并打印出方程。试求3次，4次多项式的曲线拟合，再根据数据曲线形状，求一个另外函数的拟合曲线，用图示数据曲线及相应的三种拟合曲线。3,给定数据点
8、（xiyi）如表所示。实验结果olyval(p,x);plot(x,f,x,y,'o&
