数值分析实验报告二

1贵州师范大学数学与计算机科学学院学生实验报告课程名称：数值分析班级数学与应用数学2班实验日期：2013年9月19日学号：姓名指导教师：实验成绩：一、实验名称实验一:递推法的稳定性,秦九韶算法二、实验目的及要求1.熟悉数值稳定的概念,通过上机计算，了解舍入误差所引起的数值不稳定性.2.培养Matlab编程与上机调试能力.三、实验环境每人一台计算机，要求安装WindowsXP操作系统，Microsoftoffice2003、MATLAB6.5(或7.0).四、实验内容1.教材例1.13中,取15位数字计算,并分析、比较计算结果.2.设100999832()101100994321fxxxxxxx,用秦九韶算法编程计算在()fx1,2,3,4x上的值.五、算法描述及实验步骤针对内容一：（1）先用积分求出1014145xxI（2）根据递推公式5511nnInI依次求出nI的值（n=13,12,11,..）针对内容二：（1）首先在Matlab里编制函数式M文件，函数名为Qinjiu(a,n,x).functiony=Qinjiu(a,n,x)v=a(1);fork=1:nv=v.*x+a(k+1);endy=v;2（2）然后调用函数Qinjiu分别计算出x=1，2，3，4上的值.六、调试过程及实验结果针对问题一：（1）在命令窗口运行程序:int(x.^14/(x+5),0,1)ans=6103515625*log(2)[phone]*log(3)-6103515625*log(5)-80201899665047/72072formatlongans=1.122951507568359e-002所以14I=0.011229515075684（取15位小数后的结果）（2）运用for循环算出n=13,12,...0的值intn;s=1.122951507568359e-002intn;for(n=14;;0){s=1./(5*n)-s./5printf(‘%d’,s)n--;}（3）计算的结果为：（均保留15位小数）13I0.01203898112705812I0.01297665313050011I0.01407133604056710I0.0153675509737059I0.0169264898052598I0.0188369242611707I0.0212326151477666I0.0243249055418755I0.0284683522249584I0.0343063295550083I0.0431387340889982I0.0580389198488671I0.0883922160302260I0.018232155679396针对问题二：（1）在命令窗口运行程序：3a=101:-1:1;（2）调用函数Qinjiu的计算结果为：Qinjiu(a,100,1)ans=5151Qinjiu(a,100,2)ans=2.5353e+032Qinjiu(a,100,3)ans=7.7693e+049Qinjiu(a,100,4)ans=2.1569e+062所以x=1,2,3,4时的值分别为：51512.5353e+0327.7693e+0492.1569e+062七、总结在整个实践过程中学会了用秦九韶算法编程计算多项式的值并能够采用稳定的数值计算方法比较精确的计算定积分。学会自己发现错误解释错误纠正错误。同时也培养了自己的编程能力以及上机调试能力。八、附录（源程序清单）附录1：M文件functiony=Qinjiu(a,n,x)v=a(1);fork=1:nv=v.*x+a(k+1);endy=v;intn;附录2：递推公式求值s=1.122951507568359e-002intn;for(n=14;;0){s=1./(5*n)-s./5printf(‘%d’,s)n--;}4