数值考试讲解.docx

一、插值法（P50）: 1、己知函数在下列各点的值为 X 0.2 0.4 0.6 | 0.8 | 1.0 fM 0.98 0.92 0.81 1 0.64 1 0.38 用4次牛顿插值法对数据进行插值。 function f二lagfun(x) a=[0. 2,0. 4,0. 6,0. 8,1.0]; b=[0. 98, 0. 92,0.81,0. 64, 0. 38]; for i=l:5 L ⑴=1; for j=l:5 if i L(i) =L(i)*(x-a(j))/(a (i)-a(j)); end end end f=0; for i=l:5 f=f+L(i)*b(i): end 执行文件 x0=[0. 2, 0. 4, 0. 6, 0. 8,1.0]; y0=[0. 98, 0. 92, 0. 81, 0. 64, 0. 38]; plot (xO, yO, o) hold on grid on fplot (* lagfun , L0, 1]) ;hold on x=0:0. 1:1; plot (x, newton (xO, yO, x), r); 09 0.807 0.8 07 06 05 04 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.牛顿插值以及三次样条插值(第一个实验题) 此题要求利用给定点，及给定点的函数值进行牛顿插值以及三次样条插 值。 a.牛顿插值 要实现牛顿插值，要用到以下代码 % 调用格式：yi=Lagran_(x, y, xi) % x,y数组形式的数据裘 function fi=Lagran_(x, f, xi) fi=zeros(size(xi)); npl=length(f); for i=l：npl z二ones(size(xi)): for j=l:npl if z=z. *(xi-x(j))/(x(i)-x(j)) :end end fi=fi+z*f(i); end 二、曲线拟合(P95): 16、观测物体的直线运动，得出以卜数据: 时间t/s 0 0.9 1.9 3.0 3.9 5.0 距离s/m 0 10 30 1 50 80 110 求运动方程。 解：被观测物体的运动距离与运动时间大体为线性函数关系，从而选择线性方程 则 || 戒=6』就=53.63, (/咫)=14.7, (仞o，s) = 28O,0,s) = 1O78, 则法方程组为 6 14.7 14.7 53.63 从而解得 □ =—7.855048 6 = 22.25376 故物体运动方程为 S = 22.25376,— 7.855048 MATLAB 程序： x0=[0 0.9 1.9 3 3.9 5]; y0=[0 10 30 50 80 110]; a=polyfit(x0.y0J) x=0:0.I:5; y=polyvral(a,x); %计算拟合多项式在x的值 plot(x0,y0,*.x,y) 运行结果： a = 22.2538 -7.8550 120 1 17、己知实验数据如下: 19 25 31 38 44 19.0 32.3 49.0 73.3 97.8 用最小二乘法求形如y = a + br的经验公式，并计算均方误差。 解：若s = a + bx\则 || 如 |；= 5,|础=7277699, （%,们）=5327, （￡ %） = 271.4,（。们）= 369321.5, 则法方程组为 57277699 5 5327 5327 7277699丿0 r271.4 ） 369321.5 丿 从而解得 = 0.9726046 = 0.0500351 故 y = 0.9726046+0.050035 lx2 均方误差为8 = [￡()，(土) —y/F =0.1226 ，=0 18、在某化学反应中，由实验得分解物质浓度与时间关系如下: 时间t/s 0 5 10 15 20 25 30 35 40 45 50 55 浓度 y/(xio“ 0 1.27 2.16 2.86 3.44 3.87 4.15 4.37 4.51 4.58 4.62 4.64 用最小二乘法求)，=y(/)o 解：观察所给数据的特点，采用方程 -b y = ae1 ,(u,b0) 两边同时取对数，则 111 y = 111 d - y 取中=卬口〃(1,一丄|,5 = = 则 S = a +b*x |閻|；= 11,闽|； = 0.062321, (%,Q = -0.603975, (%，/) = -87.674095,0 J) = 5.032489. 则法方程组为 11 一0.603975丫/)_(一87.674095、 、-0.603975 0.062