专升本理工科2025年数值分析测试试卷（含答案）

专升本理工科2025年数值分析测试试卷（含答案）考试时间：______分钟总分：______分姓名：______一、简述浮点数的表示方法，并说明由此产生的舍入误差的来源。当用浮点数表示实数时，为什么通常会产生舍入误差？二、设方程x^3-x-1=0在区间[1,2]内有唯一实根。用二分法求该根的近似值，使其误差不超过10^-4，需要迭代多少次？三、用迭代法x_(k+1)=x_k-(x_k^3-x_k-1)/(3x_k^2-1)求方程x^3-x-1=0在区间[1,2]内的根，若取初值x_0=1.5，判断该迭代法是否收敛？为什么？四、已知函数f(x)=e^x，求其在x=0处的二次和三次牛顿插值多项式N_2(x)和N_3(x)。五、给定数据点(0,1),(1,3),(2,2),(3,4)。求其最小二乘拟合直线y=a+bx，并给出系数a和b的值。六、用梯形公式和辛普森公式计算积分∫[0,1]e^xdx，取n=4（即h=0.25），比较两种方法的计算结果与真值e-1的误差。七、用欧拉法求解初值问题y'=x+y,y(0)=1，取步长h=0.1，计算y(0.3)的近似值。八、证明：如果线性方程组Ax=b的系数矩阵A的谱半径ρ(A)<1，那么由迭代格式x_(k+1)=Bx_k+c（其中B=(I-A)）产生的迭代法收敛于方程组的解。试卷答案一、浮点数通常用格式x=m×r^e表示，其中m是尾数（通常满足1≤|m|<r），r是基（如2或10），e是阶码。舍入误差的来源是在用浮点数表示实数时，由于计算机只能表示有限位数，需要对实数的尾数进行截断或四舍五入，从而与原实数产生差异。这个差异就是舍入误差。二、函数f(x)=x^3-x-1在区间[1,2]上连续，且f(1)=-1<0,f(2)=5>0，由介值定理知存在根。又f'(x)=3x^2-1在[1,2]上为正，故f(x)单调递增，根唯一。二分法区间缩小一半，误差限为(b-a)/2。要求|x*-x_k|≤10^-4，需满足(2-1)/2^k≤10^-4，即2^k≥10^4。计算2^k≥10^4，得k≥log2(10^4)=4*log2(10)≈4*3.32193=13.28772。因为k必须为整数，所以需迭代14次。三、迭代函数为g(x)=x-(x^3-x-1)/(3x^2-1)=(2x^3+x+1)/(3x^2-1)。计算导数g'(x)=[(6x^2+1)(3x^2-1)-(2x^3+x+1)(6x)]/(3x^2-1)^2=[18x^4-6x^2+3x^2-1-12x^4-6x^2-6x]/(3x^2-1)^2=[6x^4-9x^2-6x-1]/(3x^2-1)^2。在区间[1,2]内，令x=1,g'(1)=(6-9-6-1)/(3-1)^2=-10/4=-2.5。令x=2,g'(2)=(48-36-12-1)/(12-1)^2=-1/121。由于|g'(x)|在区间[1,2]上的最大值为2.5>1，根据迭代法收敛定理，该迭代法在该区间内不收敛。四、插值节点为x_0=0,x_1=1,x_2=2。函数值f(x_0)=1,f(x_1)=e,f(x_2)=e^2。二次牛顿插值多项式N_2(x)的差商：f[x_0]=f(0)=1f[x_1]=(f(1)-f(0))/(1-0)=(e-1)/1=e-1f[x_0,x_1,x_2]=(f[x_1]-f[x_0])/(x_1-x_0)-(f[x_2]-f[x_0])/(x_2-x_0)=((e-1)-1)/1-(e^2-1)/2=(e-2)/1-(e^2-1)/2=e-2-(e^2-1)/2=(2e-4-e^2+1)/2=(-e^2+2e-3)/2。N_2(x)=f[x_0]+f[x_0,x_1](x-x_0)+f[x_0,x_1,x_2](x-x_0)(x-x_1)=1+(e-1)(x-0)+[(-e^2+2e-3)/2](x-0)(x-1)=1+(e-1)x+(-e^2+2e-3)/2*x(x-1)=1+(e-1)x+(-e^2+2e-3)/2*(x^2-x)=(-e^2+2e-3)/2*x^2+(e-1-(-e^2+2e-3)/2)*x+1=(-e^2+2e-3)/2*x^2+(e^2-2e+1)/2*x+1=(-e^2+2e-3)x^2+(e^2-2e+1)x+1。三次牛顿插值多项式N_3(x)的差商：f[x_2]=f(2)=e^2f[x_1,x_2]=(f(2)-f(1))/(2-1)=(e^2-e)/1=e^2-ef[x_0,x_1,x_2]=(f[x_1]-f[x_0])/(x_1-x_0)-(f[x_2]-f[x_0])/(x_2-x_0)=(e-1)-(e^2-1)/2=-e^2+2e-3)/2(同上)f[x_0,x_1,x_2,x_3]=(f[x_1,x_2]-f[x_0,x_1])/(x_2-x_0)-(f[x_2,x_3]-f[x_0,x_1,x_2])/(x_3-x_0)=(e^2-e-(e-1))/2-(f(3)-(-e^2+2e-3)/2)/3-((e^2-e)-(e-1))/2=(e^2-2e+1)/2-[(e^3-e^2-e+1)/3-(-e^2+2e-3)/2]/3-(e^2-2e+1)/2=(e^2-2e+1)/2-[(e^3-e^2-e+1)*2+(e^2-2e+3)*3]/(6*2)=(e^2-2e+1)/2-(2e^3-2e^2-2e+2+3e^2-6e+9)/12=(e^2-2e+1)/2-(2e^3+e^2-8e+11)/12=6(e^2-2e+1)-(2e^3+e^2-8e+11)/12=(72e^2-144e+72-2e^3-e^2+8e-11)/12=(-2e^3+71e^2-136e+61)/12。N_3(x)=f[x_0]+f[x_0,x_1](x-x_0)+f[x_0,x_1,x_2](x-x_0)(x-x_1)+f[x_0,x_1,x_2,x_3](x-x_0)(x-x_1)(x-x_2)=1+(e-1)x+[(-e^2+2e-3)/2]x(x-1)+[(-2e^3+71e^2-136e+61)/12]x(x-1)(x-2)=1+(e-1)x+(-e^2+2e-3)/2*(x^2-x)+(-2e^3+71e^2-136e+61)/12*(x^3-3x^2+2x)=(-2e^3+71e^2-136e+61)/12*x^3+[(-e^2+2e-3)/2*(-3)+(e-1)]*x^2+[(-e^2+2e-3)/2*2+(e-1)]*x+1=(-2e^3+71e^2-136e+61)/12*x^3+(-3e^2+6e-9+2e-1)/2*x^2+(-e^2+2e-3+e-1)*x+1=(-2e^3+71e^2-136e+61)/12*x^3+(-3e^2+8e-10)/2*x^2+(-e^2+3e-4)*x+1=(-2e^3+71e^2-136e+61)/12*x^3+(-6e^2+16e-20)/4*x^2+(-e^2+3e-4)*x+1=(-2e^3+71e^2-136e+61)/12*x^3+(-3e^2+8e-10)*x^2+(-e^2+3e-4)*x+1。五、设拟合直线为y=a+bx。根据最小二乘法，系数a,b满足：n*a+(Σx_i)*b=Σy_i(Σx_i^2)*a+(Σx_i^3)*b=Σx_i*y_i代入数据(0,1),(1,3),(2,2),(3,4)：4*a+(0+1+2+3)*b=1+3+2+4=10(0^2+1^2+2^2+3^2)*a+(0*1+1*3+2*2+3*4)*b=0*1+1*3+2*2+3*4=0+3+4+12=19即：4a+6b=1010a+9b=19解此线性方程组：4a+6b=1010a+9b=19乘以2：8a+12b=2010a+9b=19相减：(8a+12b)-(10a+9b)=20-19=>-2a+3b=1=>3b=2a+1=>b=(2a+1)/3代入4a+6b=10：4a+6*[(2a+1)/3]=104a+2(2a+1)=104a+4a+2=108a=8=>a=1b=(2*1+1)/3=3/3=1。拟合直线为y=1+1x，即y=x+1。系数a=1,b=1。六、被积函数f(x)=e^x。步长h=(1-0)/4=0.25。节点x_i=i*h，i=0,1,2,3,4。f(x_i)=e^x_i。f(x_0)=e^0=1f(x_1)=e^0.25≈1.284025f(x_2)=e^0.5≈1.648721f(x_3)=e^0.75≈2.117000f(x_4)=e^1=e≈2.718282梯形公式T(h)=h/2*[f(x_0)+2(f(x_1)+f(x_2)+f(x_3))+f(x_4)]=0.25/2*[1+2(1.284025+1.648721+2.117000)+2.718282]=0.125*[1+2(5.049746)+2.718282]=0.125*[1+10.099492+2.718282]=0.125*(13.817774)≈1.727223。积分真值I=e-1≈2.718282-1=1.718282。梯形公式误差|T(h)-I|≈|1.727223-1.718282|=0.008941。辛普森公式S(h)=h/3*[f(x_0)+4f(x_1)+2f(x_2)+4f(x_3)+f(x_4)]=0.25/3*[1+4(1.284025)+2(1.648721)+4(2.117000)+2.718282]=(1/12)*[1+5.136100+3.297442+8.468000+2.718282]=(1/12)*(20.619824)≈1.718319。辛普森公式误差|S(h)-I|≈|1.718319-1.718282|=0.000037。比较：h=0.25时，梯形公式的误差约为0.008941，辛普森公式的误差约为0.000037。辛普森公式的精度更高。七、初值问题y'=x+y,y(0)=1。步长h=0.1。求y(0.3)。欧拉法公式：y_(k+1)=y_k+h*f(x_k,y_k)。这里f(x,y)=x+y。k=0,x_0=0,y_0=1：y_1=y_0+h*f(x_0,y_0)=1+0.1*(0+1)=1+0.1=1.1。k=1,x_1=0.1,y_1=1.1：y_2=y_1+h*f(x_1,y_1)=1.1+0.1*(0.1+1.1)=1.1+0.1*1.2=1.1+0.12=1.22。k=2,x_2=0.2,y_2=1.22：y_3=y_2+h*f(x_2,y_2)=1.22+0.1*(0.2+1.22)=1.22+0.1*1.42=1.22+0.142=1.362。k=3,x_3=0.3,y_3=1.362：y_4=y_3+h*f(x_3,y_3)=1.362+0.1*(0.3+1.362)=1.362+0.1*1.662=1.362+0.1662=1.5282。所以y(0.3)的近似值为1.5282。八、证明：线性方程组Ax=