2025年大学《数学与应用数学》专业题库- 计算数学在地球物理学中的应用

2025年大学《数学与应用数学》专业题库——计算数学在地球物理学中的应用考试时间：______分钟总分：______分姓名：______一、选择题（请将正确选项的字母填在题后的括号内。每小题3分，共15分。）1.在求解线性方程组Ax=b时，若系数矩阵A的行范数、列范数和条件数均大于1，则迭代法（）收敛。A.雅可比B.高斯-赛德尔C.SORD.以上都可能2.已知函数f(x)在[a,b]上连续，下列哪种方法得到的数值积分结果一定精确等于∫[a,b]f(x)dx？（不考虑舍入误差）A.梯形法则B.辛普森法则C.牛顿-柯特斯公式（n=2）D.高斯求积公式3.对于初值问题y'=f(t,y),y(t0)=y0，欧拉方法的局部截断误差阶为（）。A.O(1)B.O(h)C.O(h^2)D.O(h^3)4.在地球物理数据处理中，对地震记录进行滤波以去除噪声，这通常涉及到（）的应用。A.多项式插值B.数值积分C.数值微分D.信号卷积5.若地球物理场满足二维拉普拉斯方程∇²u=0，在矩形区域上采用有限差分法离散时，点(u(i+1,j),u(i-1,j))处的差分格式为（）。A.u(i+1,j)+u(i-1,j)-4u(i,j)=0B.u(i+1,j)+u(i-1,j)-2u(i,j)=0C.u(i+1,j)-u(i-1,j)-2u(i,j)=0D.u(i+1,j)-u(i-1,j)=2u(i,j)二、填空题（请将答案填在题后的横线上。每小题4分，共20分。）6.用迭代法求解Ax=b时，若矩阵A的严格对角占优，则该迭代法（）收敛。7.给定数据点(x0,y0),(x1,y1)，拉格朗日插值多项式L1(x)满足L1(x0)=()，L1(x1)=()。8.数值求积公式∫[a,b]f(x)dx≈Σ[i=0ton]A[i]f(x[i])的代数精度是指该公式能精确积分多项式degreek的函数，当k=()时，该公式具有最高的代数精度。9.用四阶龙格-库塔法求解常微分方程初值问题时，步长h减小，局部截断误差将（）。10.在地球物理正演中，通过数值方法模拟地震波在介质中的传播，目的是计算（）。三、计算题（请写出详细的计算过程。每小题10分，共30分。）11.给定线性方程组：4x1+x2=5x1+3x2=7试用雅可比迭代法求其近似解，迭代两次（初始值取零向量）。12.已知函数f(x)在[0,π/2]上通过点(0,1),(π/4,√2/2),(π/2,1)。试构造二次拉格朗日插值多项式P2(x)，并计算P2(π/6)的值。13.用改进的欧拉法（即梯形法则）求解初值问题y'=x+y,y(0)=1，在区间[0,0.5]上，取步长h=0.25。四、综合应用题（请结合所学知识进行分析和解答。每小题12分，共24分。）14.地震波在均匀介质中传播时，一维波动方程为∂²u/∂t²=c²∂²u/∂x²，其中c为波速。试用显式有限差分法推导时间步长为τ，空间步长为h的差分格式，并分析其稳定性条件（稳定性数）。15.在地球电阻率测量的数值模拟中，常需要求解拉普拉斯方程∇²φ=0以确定电位分布。若在一个二维区域内，边界条件给定，内部包含一个理想导体（电位为常数）。试简述如何利用有限差分法或有限元法的基本思想对该问题进行离散和求解，并说明在处理理想导体边界时可能采用的方法。试卷答案一、选择题1.B2.C3.B4.D5.A二、填空题6.收敛7.y0,y18.n+19.减小10.地震波到达地面的时间与震源位置三、计算题11.解：雅可比迭代格式：x^(k+1)=D⁻¹(b-(L+U)x^(k))其中A=[41;13],b=[5;7]D=[40;03],L=[0-1;-10],U=[01;00]D⁻¹=[1/40;01/3],b=[5;7]x^(k+1)=[1/4*5;1/3*7]+[1/4*1*x2^(k);1/3*(-1)*x1^(k)]x^(k+1)=[5/4+1/4*x2^(k);7/3-1/3*x1^(k)]令x^(k+1)=x^(k)=x，得迭代矩阵：x=[5/4+1/4*x;7/3-1/3*x]迭代方程为：x1^(k+1)=5/4+1/4*x2^(k)x2^(k+1)=7/3-1/3*x1^(k)初始值x^(0)=[0;0]第一次迭代：x1^(1)=5/4+1/4*0=5/4=1.25x2^(1)=7/3-1/3*0=7/3≈2.3333第二次迭代：x1^(2)=5/4+1/4*2.3333=1.25+0.5833=1.8333x2^(2)=7/3-1/3*1.25=2.3333-0.4167≈1.9167近似解为x≈[1.8333;1.9167]12.解：拉格朗日插值多项式P2(x)=L0(x)f(x0)+L1(x)f(x1)+L2(x)f(x2)L0(x)=[(x-x1)(x-x2)]/[(x0-x1)(x0-x2)]=[(x-√2/2)(x-1)]/[(0-√2/2)(0-1)]=x(x-√2/2)/(√2/2)L1(x)=[(x-x0)(x-x2)]/[(x1-x0)(x1-x2)]=[(-√2/2)(x-1)]/[(-√2/2)(-√2/2)]=x-1L2(x)=[(x-x0)(x-x1)]/[(x2-x0)(x2-x1)]=[(x-0)(x-√2/2)]/[(π/2-0)(π/2-√2/2)]=x(x-√2/2)/[π/2*(π/2-√2/2)]P2(x)=[x(x-√2/2)/(√2/2)]*1+(x-1)*(√2/2/2)+[x(x-√2/2)/(π/2*(π/2-√2/2))]*1P2(x)=[2x(x-√2/2)]/√2+[√2/2*(x-1)]/2+[2x(x-√2/2)]/[π(π-√2)]P2(π/6)=[2*(π/6)*(π/6-√2/2)]/√2+[√2/2*(π/6-1)]/2+[2*(π/6)*(π/6-√2/2)]/[π(π-√2)]P2(π/6)=[π(π-3√2)]/(3√2)+[√2(π-6)]/[4π]+[π(π-3√2)]/[3π(π-√2)]P2(π/6)=[π(π-3√2)]/(3√2)+[√2(π-6)]/[4π]+[π-3√2]/[3(π-√2)]P2(π/6)=[π(π-3√2)*(π-√2)+√2(π-6)*3(π-√2)+3√2*4π*(π-3√2)]/[12π√2(π-√2)]P2(π/6)=[π(π^2-4π√2+3√2^2)+3√2(π^2-8π+6√2)+12π√2(π-3√2)]/[12π√2(π-√2)]P2(π/6)=[π^3-4π^2√2+9π+3√2π^2-24π√2+36+12π^2√2-36√2^2]/[12π√2(π-√2)]P2(π/6)=[π^3+(-4+3√2)π^2+(-24+9+12√2)π+(36-36*2)]/[12π√2(π-√2)]P2(π/6)=[π^3+(-4+3√2)π^2+(-15+12√2)π]/[12π√2(π-√2)]13.解：改进的欧拉法（梯形法则）：k₁=f(tⁿ,yⁿ)=tⁿ+yⁿk₂=f(tⁿ+h,yⁿ+hk₁)=(tⁿ+h)+(yⁿ+h(tⁿ+yⁿ))=tⁿ+yⁿ+h(2tⁿ+yⁿ)y^(n+1)=yⁿ+h/2*(k₁+k₂)=yⁿ+h/2*[(tⁿ+yⁿ)+(tⁿ+yⁿ+h(2tⁿ+yⁿ))]=yⁿ+h/2*[2tⁿ+2yⁿ+h(2tⁿ+yⁿ)]=yⁿ+h*(tⁿ+yⁿ)+h²/2*(2tⁿ+yⁿ)在[0,0.5]区间，tⁿ=n*0.25,n=0,1,2;h=0.25;yⁿ=y(0)=1n=0(tⁿ=0,yⁿ=1):k₁=0+1=1k₂=0+1+0.25(2*0+1)=1+0.25*1=1.25y^(1)=1+0.25/2*(1+1.25)=1+0.125*2.25=1+0.28125=1.28125n=1(tⁿ=0.25,yⁿ≈1.28125):k₁=0.25+1.28125=1.53125k₂=0.25+1.28125+0.25(2*0.25+1.28125)=1.53125+0.25(0.5+1.28125)=1.53125+0.25*1.78125=1.53125+0.4453125=1.9765625y^(2)=1.28125+0.25/2*(1.53125+1.9765625)=1.28125+0.125*3.5078125=1.28125+0.4384375=1.7196875四、综合应用题14.解：一维波动方程∂²u/∂t²=c²∂²u/∂x²采用中心差分格式离散时间导数和空间导数：∂u/∂t≈(u(x,t+τ)-u(x,t-τ))/(2τ)∂²u/∂t²≈(u(x,t+τ)-2u(x,t)+u(x,t-τ))/τ²∂u/∂x≈(u(x+h,t)-u(x-h,t))/(2h)∂²u/∂x²≈(u(x+h,t)-2u(x,t)+u(x-h,t))/h²将差分格式代入原方程：(u(x,t+τ)-2u(x,t)+u(x,t-τ))/τ²=c²*(u(x+h,t)-2u(x,t)+u(x-h,t))/h²整理得显式格式：u(x,t+τ)=2u(x,t)-u(x,t-τ)+(cτ/h)²*(u(x+h,t)-2u(x,t)+u(x-h,t))记C=cτ/h，则：u(x,t+τ)=2u(x,t)-u(x,t-τ)+C²*(u(x+h,t)-2u(x,t)+u(x-h,t))稳定性分析：考虑方程u(x,t+τ)-u(x,t)=[2-2cos(ωh)+C²sin²(ωh)]u(x,t)(ω为角频率)对于稳定性，需有|2-2cos(ωh)+C²sin²(ωh)|≤1因为sin²(ωh)=1-cos²(ωh)，代入得：|2-2cos(ωh)+C²(1-cos²(ωh))|=|2-2cos(ωh)+C²-C²cos²(ωh)|≤1令z=cos(ωh)，则|2-2z+C²-C²z²|≤1展开并考虑z²=(cos(ωh))²，得|(2+C²)-2z-C²z²|≤1显式格式的稳定性条件为C≤1，即cτ/h≤1，或h≥cτ。15.解：拉普拉斯方程∇²φ=0在二维区域Ω内，边界条件为Γ。若Ω内包含理想导体，设其区域为Ω_c，电位为常数V_c。原问题可分解为：1.在区域Ω-Ω_c内，求解∇²φ=0，边界条件为Γ上的给定电位值。2.在区域Ω_c内，电位恒定为V_c。数值求解方法：*有限差分法：