数值分析实验报告二.docx

贵州师范大学数学与计算机科学学院学生实验报告课程名称：数值分析 班级 数学与应用数学2班 实验日期：2013年9月19日学 号：—姓名 指导教师：实验成绩： 一、 实验名称实验一:递推法的稳定性，秦九韶算法二、 实验目的及要求1. 熟悉数值稳定的概念，通过上机计算，了解舍入误差所引起的数值不稳定 性.2. 培养Mat lab编程与上机调试能力.三、 实验环境每人一台计算机，要求安装Windows XP操作系统，Microsoft office2003、 MATLAB6.5(或 7.0).四、 实验内容1. 教材例1.13中，取15位数字计算,并分析、比较计算结果.2. 设 f (x) = 101x100 + 100x99 + 99x98 + + 4x3 + 3x2 + 2x + 1 ,用秦九韶算法编程计算在f (x) x = 1,2,3,4上的值.五、 算法描述及实验步骤针对内容一：! x14(1) 先用积分求出〈4 =』x + 5(2) 根据递推公式In _1 = 5- -1依次求出In的值(n=13,12,11,..)针对内容二：(1)首先在Matlab里编制函数式M文件，函数名为Qinjiu(a,n,x).function y=Qinjiu(a,n,x)v=a(1);for k=1:nv=v.*x+a(k+1);end（2）然后调用函数Qinjiu分别计算出x=1, 2, 3, 4上的值.六、调试过程及实验结果针对问题一：(1)在命令窗口运行程序：>> int(x.A14/(x+5),0,1)ans =6103515625*log(2)+6103515625*log(3)-6103515625*log(5)-80201899665047/72072>> format longans =1.122951507568359e-002所以 I14 = 0.011229515075684 （取15位小数后的结果）（2）运用for循环算出n=13,12,...0的值>> int n;>> s= 1.122951507568359e-002>> int n;for(n=14;;0){s=1./(5*n)-s./5printf(‘％d’,s)n--;}(3)计算的结果为：I13 = 0.012038981127058In = 0.01407133604056719 = 0.01692648980525917 = 0.02123261514776615 = 0.02846835222495813 = 0.043138734088998I] = 0.088392216030226针对问题二：(均保留15位小数)〈2 = 0.012976653130500I10 = 0.015367550973705I8 = 0.01883692426117016 = 0.02432490554187514 = 0.03430632955500812 = 0.05803891984886710 = 0.018232155679396(1)在命令窗口运行程序：a=101：-1:1;(2)调用函数Qinjiu的计算结果为：>> Qinjiu(a,100,1)ans =5151>> Qinjiu(a,100,2)ans =2.5353e+032>> Qinjiu(a,100,3)ans =7.7693e+049>> Qinjiu(a,100,4)ans =2.1569e+062所以 x=1,2,3,4 时的值分别为：5151 2.5353e+032 7.7693e+049 2.1569e+062七、 总结在整个实践过程中学会了用秦九韶算法编程计算多项式的值并能够采用稳定 的数值计算方法比较精确的计算定积分。

学会自己发现错误解释错误纠正错误 同时也培养了自己的编程能力以及上机调试能力八、 附录(源程序清单)附录1： M文件function y=Qinjiu(a,n,x)v=a(1);for k=1:nv=v.*x+a(k+1);endy=v;>> int n;附录2：递推公式求值>> s= 1.122951507568359e-002>> int n;for(n=14;;0){s=1./(5*n)-s./5printf(‘％d’,s)n--;}。