数值分析实验报告二

1、精选优质文档-倾情为你奉上贵州师范大学数学与计算机科学学院学生实验报告课程名称： 数值分析 班级 数学与应用数学2班 实验日期：2013年9月19日学 号： 姓名 指导教师： 实验成绩： 一、实验名称 实验一:递推法的稳定性,秦九韶算法二、实验目的及要求 1.熟悉数值稳定的概念, 通过上机计算，了解舍入误差所引起的数值不稳定性. 2. 培养Matlab编程与上机调试能力.三、实验环境每人一台计算机，要求安装Windows XP操作系统，Microsoft office2003、MATLAB6.5(或7.0).四、实验内容 1.教材例1.13中,取15位数字计算,并分析、比较计算结果. 2.设,

2、用秦九韶算法编程计算在上的值.五、算法描述及实验步骤针对内容一：（1）先用积分求出（2）根据递推公式依次求出的值 （n=13,12,11,.）针对内容二：（1） 首先在Matlab里编制函数式M文件，函数名为Qinjiu(a,n,x).function y=Qinjiu(a,n,x)v=a(1);for k=1:n v=v.*x+a(k+1);endy=v;（2）然后调用函数Qinjiu分别计算出x=1，2，3，4上的值.6、 调试过程及实验结果针对问题一：（1） 在命令窗口运行程序: int(x.14/(x+5),0,1)ans = *log(2)+*log(3)-*log(5)-047/7

3、2072 format long ans = 1.8359e-002所以 = 0.5684 （取15位小数后的结果）（2） 运用for循环算出n=13,12,.0的值 int n; s= 1.8359e-002 int n;for(n=14;0)s=1./(5*n)-s./5printf(%d,s)n-;（3）计算的结果为： （均保留15位小数）0.7058 0.05000.0567 0.37050.5259 0.11700.7766 0.18750.4958 0.50080.8998 0.88670.0226 0.9396针对问题二：（1）在命令窗口运行程序：a=101:-1:1;（2） 调

4、用函数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+0627、 总结 在整个实践过程中学会了用秦九韶算法编程计算多项式的值并能够采用稳定的数值计算方法比较精确的计算定积分。学会自己发现错误解释错误纠正错误。同时也培养了自己的编程能力以及上机调试能力。8、 附录（源程序清单）附录1：M文件 function y=Qinjiu(a,n,x)v=a(1);for k=1:n v=v.*x+a(k+1);endy=v; int n;附录2：递推公式求值 s=