数值分析实验报告二.docx

文档介绍：贵州师范大学数学与计算机科学学院学生实验报告
课程名称：数值分析 班级 数学与应用数学2班 实验日期：2013年9月19日
学 号：—姓名 指导教师：
实验成绩：
一、 实验名称
实验一:递推法的稳定性，秦九韶算法
二、 实验贵州师范大学数学与计算机科学学院学生实验报告
课程名称：数值分析 班级 数学与应用数学2班 实验日期：2013年9月19日
学 号：—姓名 指导教师：
实验成绩：
一、 实验名称
实验一:递推法的稳定性，秦九韶算法
二、 实验目的及要求
熟悉数值稳定的概念，通过上机计算，了解舍入误差所引起的数值不稳定 性.
培养Mat lab编程与上机调试能力.
三、 实验环境
每人一台计算机，要求安装Windows XP操作系统，Microsoft office2003、 (或 ).
四、 实验内容
，取15位数字计算,并分析、比较计算结果.
设 f (x) = 101x100 + 100x99 + 99x98 + + 4x3 + 3x2 + 2x + 1 ,用秦九韶算法编程计
算在f (x) x = 1,2,3,4上的值.
五、 算法描述及实验步骤
针对内容一：
! x14
先用积分求出〈4 =』x + 5
根据递推公式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:n
v=v.*x+a(k+1);
end
（2）然后调用函数Qinjiu分别计算出x=1, 2, 3, 4上的值.
六、调试过程及实验结果
针对问题一：
(1)在命令窗口运行程序：
>> int(/(x+5),0,1)
ans =
6103515625*log(2)[phone]*log(3)-6103515625*log(5)-80201899665047/72072
>> format long
ans =
-002
所以 I14 = （取15位小数后的结果）
（2）运用for循环算出n=13,12,...0的值
>> int n;
>> s= -002
>> int n;
for(n=14;;0)
{
s=1./(5*n)-s./5
printf(‘％d’,s)
n--;
}
(3)计算的结果为：
I13 =
In =
19 =
17 =
15 =
13 =
I] =
针对问题二：
(均保留15位小数)
〈2 =
I10 =
I8 =