数值考试分析.doc

一、插值法（P50）：1、已知函数在下列各点的值为 0.2 0.4 0.6 0.8 1.0 0.98 0.92 0.81 0.64 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=1:5L(i)=1;for j=1:5if j~=iL(i)=L(i)*(x-a(j))/(a(i)-a(j));endend end f=0; for i=1:5f=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(x0,y0,o) hold on grid on fplot(lagfun,[0,1]);hold on x=0:0.1:1; plot(x,newton(x0,y0,x),r); 牛顿插值以及三次样条插值（第一个实验题） 此题要求利用给定点，及给定点的函数值进行牛顿插值以及三次样条插值。 a.牛顿插值 要实现牛顿插值，要用到以下代码 % 调用格式：yi=Lagran_(x,y,xi) % x,y 数组形式的数据表 function fi=Lagran_(x,f,xi) fi=zeros(size(xi));np1=length(f); for i=1:np1z=ones(size(xi));for j=1:np1if i~=j,z=z.*(xi-x(j))/(x(i)-x(j));endendfi=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 50 80 110 求运动方程。 解：被观测物体的运动距离与运动时间大体为线性函数关系，从而选择线性方程 令 则 则法方程组为 从而解得 故物体运动方程为 MATLAB程序： x0=[0 0.9 1.9 3 3.9 5]; y0=[0 10 30 50 80 110]; a=polyfit(x0,y0,1) x=0:0.1:5; y=polyval(a,x); %计算拟合多项式在x的值 plot(x0,y0,*,x,y) 运行结果： a =22.2538 -7.8550 已知实验数据如下： 19 25 31 38 44 19.0 32.3 49.0 73.3 97.8 用最小二乘法求形如的经验公式，并计算均方误差。 解：若，则 则 则法方程组为 从而解得 故 均方误差为18、在某化学反应中，由实验得分解物质浓度与时间关系如下： 时间t/s 0 5 10 15 20 25 30 35 40 45 50 55 浓度y/() 0 1.27 2.16 2.86 3.44 3.87 4.15 4.37 4.51 4.58 4.62 4.64 用最小二乘法求。 解：观察所给数据的特点，采用方程 两边同时取对数，则 取 则 则法方程组为 从而解得 因此三、数值积分与数值微分（P108）：0 1 1/8 0. 1/4 0. 3/8 0. 1/2 0. 5/8 0. 3/4 0. 7/8 0. 1 0.841409例三、对于函数，给出n=8时的函数表，使用复合梯形公式及复合辛普森公式计算积分，并估计误差以及比较其精度。1 复合梯形fht function T_n=fht(a,b,n)h=(b-a)/n;for k=0:nx(k+1)=a+k*h;if x(k+1)==0x(k+1)=10^(-10);endendT_1=h/2*(fx(x(1))+fx(x(n+1)));for i=2:nF(i)=h*fx(x(i));endT_2=sum(F)T_n=T_1+T_2文件2 fx.m function y=fx(x)y=sin(x)./x;end matlab中输入的 fht(0,1,8) T_2 =0.8306 T_n =0.9457 ans =0.9457 format long fht(0,1,8) T_2 =0. T_n =0. ans =0. format short fht(0,1,8) 复合桑普森复合S_P_S.m文件 function S_n=S_P_S(a,b,n)h=(b-a)/n;for k=0:nx(k+1)=a+k*h;x_k(k+1)=x(k+1)+1/2*h;if (x(k+1)==0)|