数值分析公式大全.doc

数值分析，第一章 相对误差和绝对谋差 e*=x*-x; er*=（r - XVX估计值厲_叽 误差限和相对误差限 t *|x*-x| t *|x*-x| 有效数字 官方定义：若近似值x*的误差限是某一位的半个单位，该位到/的第一位非零有效数字 共有n位,就说x*有n位有效数字。表示为:x*=±10mX C^+a?% 10+a3X 10?+???+atlX 10「5）二土a. a2a3-a?o其中为0至9中之一,m不为0, m, n都是整数。 公式：e*= |x-x*|^x 10m-n+1 相对误差限公式 x?具有n为有效数字，丄XI。。 若则x?至少具有n为有效数字。 2 （al + 1? 4,病态问题的条件数，相对误差比值 X的扰动△ X二X-X*,误差为4,函数值f（X*）的相对误差」x：-f x f （x） 相对误差比值为：紅/ j ~斗上丄=Cp （也称为条件数） f （X） X f （x） 第二章：插值法 多项式插值 P （x）为 n 阶多项式,P （x） =a0+aix+a2x2+? ? ?+anxn, a：为实数。 1尢o…xf 5 「y°i 解法：a解方程组:Aa=y,其中A二 1 %!…灯 ,a二 ?i ? ■ ? ，y= yi ■ ■ ■ ? ? ? ? -1 …煤 -Qn- -7n- 拉格朗日插值 【1】线性插值 L1二yidk+yk+Jk+i 插值基函数奇,X Xk4+1 一 Xk抛物线插值L2 二 插值基函数奇, X Xk 4+1 一 Xk 抛物线插值 L2 二 yklk+yk+dk+l+yk+2lk+2（尢一尢知1）（兀一尢知2） （尢一 g）（心 k+2） 插值基函数lk=厂5人：;人、z |k+2 （尢 尢代+1）（机一尢知 2） （4+1—xQ（g+i—g+2） N次插值多项式（通解）Ln=yolo+yili+Y2l2+—+ynln （兀-心（心知1）（尢k+2 一尢比）（机+2 一尢k+1） (X-Xo)…(X-Xk_i)(X-Xk+i)???(X-Xn) (xfc-x0)...(xk-xk_1)(xk-xk+1)...(xk-xn) 设 3”1(X)=(x - Xq) ... (% - Xfc-i)(x - xk+1) ...(X - Xn) (xQ 二(Xk 一 XO) ... (xk 一 Xk_J(Xk 一 Xk+1)…g 一 xn) 有 Ln (x) 有 Ln (x) 7 (x-Xfc)n+1 (Xfc) 余项公式 N次插值多项式的余项形式 Rn=f (x) -Ln (x) 3田(x)二K(x)3”i (x) , f G(a, b) f的位置未知，但有截断误差限： PnWI S 誥肓 |o)n+1(x)|, NU.=maxax^|/(n+1)(x)| 均差(斧商) g-Xo二阶均养:f[x°门 X1,Xk]J[x(r 心卜心o, xk]耳一心高阶均差:f[Xo g-Xo 二阶均养: f[x°门 X1, Xk]J[x(r 心卜心o, xk] 耳一心 高阶均差: f[Xo 门 X1, ... Xk]」? X1 M-1]-咲’???’ Xk-2. xk] xk^xk-i 性质：1， 2, k阶均差可表示为函数值f (xo), f (xi), f (xn)的线性组合 对称性，与节点次序无关 3, 【前后项】f[x0/, xl,…，心財… 奸心。? ???，如] g_Xo 4,※门阶均差与导数的关系：f[x0? xl, XjJ牟2 [a, b]o n! k Xk 一阶差商二阶差商三阶差 商… 0 1 2 3 4 a i a X X X X X 0 2 3 4 f(xj f(x2) 偸) 1 1 1 % X』 f[x-j, x2] f[x0,x1.x2] f[X2, x3l mqx’xd 臥讯内?淤』 血3，句 t[x2, X2.X^l 1p(vX2.fX2.X^ … 1 1 1 ■ a i i i i 牛顿插值多项式 逐次生成的插值多项式 Pn(X)=a0+ai (x-x0) +a2 (x-x0)(x-xi) +—+an (x?x°)…(x-xn_i) a0=f(xo), ai=f[x0, Xi], a2=f[x0, xn x2],…，an=f[x0 xn ???, xn] 【余项】Rn=f[x, Xo，X1，…，Xn] o n+1 (x) 估计截断误差限 5,差分 等距离节点 xk=x0+kh, k=0, 1,…，n； fk=f (xk) Xk处的一阶向前差分：Xk处的二阶向前差分：△条△伉+】?△ fk； Xk 处的 n 阶羌分：△ r，fk= △ n_lfk.+1- A % 【差分与差商的关系】先晋晋-般的臥皿?沁訴箒 【弟分与导数的关系】Amffc=hmf(m) ( C ) 差分表