计算方法上机大作业.docxVIP

绪论 1.用二分法求方程x-x-l=O在［1,2］内的近似根,要求误差不超10 3. 程序： clc;clear a二1;b=2; fa=a*a*a-a~l; fb二b*b*b-bT; c=(a+b)/2; fc二c*c*c-cT; if fa*fb〉0, break, end wh订e abs(fc)10^(-3) c=(a+b)/2; fc二c*c*c-cT; if fb*fc0 b=c; fb二fc; else a=c; fa=fc; end end format long fx=fc, x=c 运算结果： fx =-4. 659488331526518e-05 x =1.324707031250000 2 ?证明方程f(x)二eFlOx-2二0内在［0, 1］内有唯一实根，用二分法求 这一实根，要求误差不超过0?5*10「2. 证明： 由题可知，f (x)在［0,1］上为正，也就是f(x)单调递增，又因 为f (0) 0, f (1) 0,所以在［0, 1］内有唯一实根。 程序： clc;clear a=0;b=l; fa=exp (a)+10*a-2; fb二exp(b)+10*b-2; c二(a+b)/2; fc=exp(c)+10*c-2; if fa*fb〉0, break, end wh订e abs (fc)0. 5*10(-2) c=(a+b)/2; fc二exp(c)+10*c-2; if fb*fc〉0 b二c; fb=fc; else a=c; fa=fc; end end format long fx二fc, x二c 运算结果： fx 二 0.003275341789827 x =0.090820312500000 第一章 11.给岀概率积分y = ^=^e-x2dx的数据表 ■ 1 0 1 2 3 xi 0. 46 0. 47 0. 48 0. 49 yi 0.4846555 0.4937452 0.5027498 0.5116683 用二次插值计算，试问: (1)当x二0.472时该积分值等于多少？ 程序： clear all; x二[0.46 0. 47 0. 48 0. 49]; y=[0.4846555 0.4937452 0.5027498 0.5116683]: xi=0.472; yi二interpl(x, y, xi, linear); yi 运算结果： yi 二 0.495546120000000 (2)当x为何值时积分值等于0.5? 程序： clear all; x=[0. 46 0. 47 0. 48 0. 49]; y二[0.4846555 0. 4937452 0. 5027498 0.5116683]; yi=0. 5; xi-interpl (y, x, yi, linear，); xi 运算结果： xi =0.476946227483731 34?构造适合下列数据表的三次插值样条S (x): X -1 0 1 2 y -1 1 3 5 y 6 1 程序： function [ ]=spline31 (X, Y, dY, xO, m) N二size (X, 2); sO=dY(l) : sN=dY(2); interval=0. 025; disp(，xO为插值点) h二zeros(1, N~l); for i=l:N-l h(l, i)二X(i+1)-X(i); encl d(l, 1)=6*((Y(l, 2)-Y(l, l))/h(l, l)-sO)/h(l, 1); d(N, 1)二6*(sN-(Y(1,N)-Y (1, N-l))/h(l, N-l))/h(l, N-l); for i二2:NT d(i, 1)二6*((Y(1, i + l)-Y(l, i))/h(l, i)-(Y(l, i)-Y(l, i-l))/h(l, i- l))/(h(l, i)+h(l, i-1)); end mu=zeros (1,N-l); md二zeros(1, N-l); md(l,N-l)=l; mu(l, 1)=1; for i=l:N-2 u二h(l, i+l)/(h(l, i)+h(l, i+1)) ; mu(l, i+1)二u; md(l, i)=l-u; end p(l, 1)=2; q(l, l)=mu(l, 1)/2; for i=2:N-l p(l, i)=2-md(l, i-l)*q(l, i-1) ; q(l, i)二mu(l, i)/p(l, i); end p(l,N) =2-md(l,N-l)*q(l,N-l); y二zeros(l, N) ; y(1, l)=d(l)/2; for i=2:N y(1, i)二(d(i)-md(l, i-