《实验一数值计算中误差的传播规律.docVIP

数 值 计 算 方 法 实 验 报 告 实验序号：实验一 实验名称：数值计算中误差的传播规律 实 验 人： 熊华华 专业年级：09信息与计算科学 教 学 班： B01 学 号：0094181 实验时间：2012.03.01 江西财经大学信息管理学院 实验一　数值计算中误差的传播规律 一、实验目的 １．观察并初步分析数值计算中误差的传播； ２．观察有效数字与误差传播的关系． 二、实验要求 本次实验包含三个相对独立的内容． ．在内容１中，请解释两个命令的格式和作用； ．求解方程时，分别使用求根公式和韦达定理两种方法，并比较其有效数字和相对误差； ．实验内容３中的５个函数在处的精确值都是相等的，若取进行计算，计算各函数的结果，作图观察并比较它们的绝对误差（作图区间可取甚至更小）,并从算法设计原则上说明原因． 三、实验步骤 1．使用MATLAB的help命令学习MATLAB命令digts和vpa的用途和使用格式；DIGITS Set Maples Digits. Digits determines the accuracy of Maples numeric computations. DIGITS, by itself, displays the current setting of Digits. DIGITS(D) sets Digits to D for subsequent calculations. D is an integer, or a string or sym representing an integer. D = DIGITS returns the current setting of Digits. D is an integer. VPA Variable precision arithmetic. R = VPA(S) numerically evaluates each element of the double matrix S using variable precision floating point arithmetic with D decimal digit accuracy, where D is the current setting of DIGITS. The resulting R is a SYM. VPA(S,D) uses D digits, instead of the current setting of DIGITS. D is an integer or the SYM representation of a number. It is important to avoid the evaluation of an expression using double precision floating point arithmetic before it is passed to VPA. For example, phi = vpa((1+sqrt(5))/2) first computes a 16-digit approximation to the golden ratio, then converts that approximation to one with d digits, where d is the current setting of DIGITS. To get full precision, use unevaluated string or symbolic arguments 2．在4位浮点数下解二次方程； x1=(-62+sqrt(62^2-4))/2 x1 = -0.0161 x2=1/x1 x2 = -61.9839 x2=(-62-sqrt(62^2-4))/2 x2 = -61.9839 x1=1/x2 x1 = -0.0161 16位下的有效结果 x1=(-62+sqrt(62^2-4))/2 x1 = -0.01613323034066 x2=1/x1 x2 = -61.98386676966352 3．计算下列５个函数在点处的近似值 （）， x=1.4 x = 1.40000000000000 y=[(x-1)^6] y = 0.00409600000000 （）， y=[1/(x+