2025年大学《数理基础科学》专业题库-偏微分方程的理论与实践探索

2025年大学《数理基础科学》专业题库——偏微分方程的理论与实践探索考试时间：______分钟总分：______分姓名：______一、1.指出下列偏微分方程的阶、线性/非线性、齐次/非齐次，并判断其方程类型（椭圆型、抛物型、双曲型）。(1)u_xx+u_yy-6u_xy+u=sin(x)(2)t(u_t)+x(u_x)=u_xx(3)u_tt-u_xx=0(4)u_xx+u_yy=02.设u(x,y)满足方程u_xx+u_yy=0在区域D内，其中D为不包含原点的上半平面。证明：如果u在D内有界，则u恒为常数。二、3.用分离变量法求解下列定解问题：(1)u_xx-u_yy=0,0<x<π,0<y<h,u(0,y)=0,u(π,y)=0,u(x,0)=f(x),u(x,h)=g(x)。(2)u_tt=u_xx,0<x<1,t>0,u(0,t)=0,u(1,t)=0,u(x,0)=x(1-x),u_t(x,0)=sin(πx)。4.求解下列初值问题（热传导方程）：u_t=u_xx,0<x<L,t>0,u(0,t)=0,u(L,t)=0,u(x,0)=f(x)。三、5.利用傅里叶变换法求解下列定解问题：u_t=k(u_xx),0<x<∞,t>0,u(0,t)=g(t),u(x,0)=f(x)。（提示：考虑使用傅里叶正弦变换或余弦变换，根据边界条件选择）6.求解下列定解问题（波动方程的达朗贝尔公式）：u_tt=a^2u_xx,0<x<L,t>0,u(0,t)=0,u(L,t)=0,u(x,0)=f(x),u_t(x,0)=g(x)。四、7.用特征线法求解下列一阶线性偏微分方程初值问题：yu_x+xu_y=x^2+y^2,x>0,y>0,u(x,1)=x+1。8.设u(x,y)满足拉普拉斯方程u_xx+u_yy=0在单位圆盘D={(x,y)|x^2+y^2<1}内，边界条件为u(x,y)=1/x（圆周上，x≠0）。证明u在D内无解。五、9.一根长为L的均匀细杆，初始温度分布为f(x)，两端分别维持在不同温度T_1和T_2(T_1≠T_2)。若考虑热传导过程，其定解问题为：u_t=ku_xx,0<x<L,t>0,u(0,t)=T_1,u(L,t)=T_2,t>0,u(x,0)=f(x),0<x<L。(1)证明此问题的解可以表示为u(x,t)=T_1+(T_2-T_1)x+∑_{n=1}^∞b_ne^{-k(πn/L)^2t}sin(πnx/L)，其中系数b_n由初始条件f(x)决定。(2)若T_1=T_2=0，即问题变为齐次边界条件的热传导，解的形式如何？10.简述有限差分法求解偏微分方程的基本思想，并以一阶热传导方程u_t=ku_xx的一维初值问题为例，说明如何离散化并建立差分格式。六、11.考虑二维拉普拉斯方程u_xx+u_yy=0在圆盘区域内的定解问题，边界条件为u(r,θ)=f(r,θ)（圆周上）。解释为何在此问题中，使用分离变量法（假设解为R(r)Θ(θ)）会遇到困难，并简述一种可能的处理方法（如极坐标下的分离变量法或直接使用特殊函数解）。12.评述偏微分方程在理解物理现象（如波的传播、热扩散）和解决工程问题（如结构分析、流体力学）中的核心作用，并举例说明如何从实际问题中导出相应的偏微分方程模型。试卷答案一、1.(1)二阶，非线性，非齐次，双曲型。（提示：考察高阶、非线性、非齐次项及混合导数系数）(2)一阶，线性，非齐次，双曲型。（提示：考察一阶方程，线性项和非齐次项）(3)二阶，线性，齐次，双曲型。（提示：标准波动方程）(4)二阶，线性，齐次，椭圆型。（提示：标准拉普拉斯方程）2.证明：令v(x,y)=u(ax+by,cx+dy)，其中a,b,c,d满足a^2+b^2=c^2+d^2≠0。则v_xx=a^2u_xx+2abu_xy+b^2u_yyv_yy=c^2u_xx+2cdu_xy+d^2u_yyv_xy=abu_xx+(a^2+b^2)u_xy+cdu_yy由于u满足u_xx+u_yy=0，代入上式v_xx+v_yy=(a^2+b^2)u_xx+2(ab+cd)u_xy+(c^2+d^2)u_yy=0。又因为a^2+b^2=c^2+d^2，所以v_xx+v_yy=0。即v满足拉普拉斯方程。另一方面，若D'为v变换后的区域，当(x,y)∈D时，(ax+by,cx+dy)∈D'。如果u在D内有界，则v在D'内也有界。由调和函数的性质知，在有界区域上满足拉普拉斯方程且有界的调和函数必为常数。因此，v恒为常数，即u(x,y)=v(1/x,1/y)恒为常数。（提示：利用拉普拉斯方程的调和函数性质，通过坐标变换证明）二、3.(1)解：令u(x,y)=X(x)Y(y)。代入方程得X''Y+XY''=0，即X''/X=-Y''/Y=λ。解特征方程Y''+λY=0，得Y=Asin(√λy)+Bcos(√λy)。解X''-λX=0，得X=Ce^(√λx)+De^(-√λx)。考虑边界条件：u(0,y)=X(0)Y(y)=Ce^(0)+De^(0)Y(y)=(C+D)Y(y)=0。由于Y(y)不恒为0，需C+D=0，即D=-C。所以X(x)=C(e^(√λx)-e^(-√λx))=2Ce^(√λx)sinh(√λx)。u(π,y)=X(π)Y(y)=2Csinh(√λπ)Y(y)=0。由于Y(y)不恒为0，需sinh(√λπ)=0。λ=n^2/(π^2),n=1,2,3,...。此时X_n(x)=2C_nsinh(nx/π)。Y_n(y)=A_nsin(ny/π)+B_ncos(ny/π)。由u(x,0)=f(x)=∑_{n=1}^∞A_nsin(nx/π)，得A_n=(2/π)∫_0^πf(x)sin(nx/π)dx。B_n由u(x,h)=g(x)=∑_{n=1}^∞(B_nsin(nh/π)+A_nsin(nx/π))=∑_{n=1}^∞B_nsin(nh/π)确定，即B_nsin(nh/π)=g(x)-A_nsin(nx/π)=0(积分后)。(需要补充说明B_n的确定方式，但题目未要求)最终解为u(x,y)=∑_{n=1}^∞2C_nsinh(nx/π)(A_nsin(ny/π)+B_nsin(nh/π))。（提示：标准分离变量法，注意边界条件的处理导致本征值和本征函数的确定）(2)解：此为齐次边界条件的波动方程初值问题。令u(x,t)=X(x)T(t)。分离变量得X''/X=T''/T=-ω^2。解得X=Acos(ωx)+Bsin(ωx)，T=Ccos(ωt)+Dsin(ωt)。u(x,0)=X(x)T(0)=X(x)(C+D)=(Acos(ωx)+Bsin(ωx))(C+D)=A(C+D)cos(ωx)+B(C+D)sin(ωx)=x(1-x)。由正交性或直接比较系数，得A(C+D)=0,B(C+D)=0。由于x(1-x)不是正弦或余弦函数，ω不能取简单值。需要将初始位移展开为正弦余弦级数。u(x,0)=x(1-x)=(1/2)-(1/2)x^2。展开为Fourier级数：x(1-x)=∑_{n=1}^∞b_nsin(nπx/π)=∑_{n=1}^∞b_nsin(nx)。b_n=2/π∫_0^πx(1-x)sin(nx)dx=(2/π)[(-1)^nn/(n^2-1)]。所以A(C+D)=0,B(C+D)=2/π[(-1)^nn/(n^2-1)]。为满足非零解，需C+D≠0。取C=0,D=1。u(x,t)=∑_{n=1}^∞B_nsin(nx)sin(ωt)，其中B_n=2/π[(-1)^nn/(n^2-1)],ω=n。最终解为u(x,t)=∑_{n=1}^∞[2/π(-1)^nn/(n^2-1)]sin(nx)sin(nxt)。（提示：齐次边界条件的波动方程，初始位移非正弦余弦形式，需先展开初始位移）4.解：令u(x,t)=X(x)T(t)。分离变量得X''/X=T''/T=-λ。解得X=Acos(√λx)+Bsin(√λx)，T=Ce^(-kλt)。u(0,t)=X(0)T(t)=(Acos(0)+Bsin(0))Ce^(-kλt)=ACe^(-kλt)=0。需A=0。u(L,t)=X(L)T(t)=Bsin(√λL)Ce^(-kλt)=0。由于Ce^(-kλt)≠0，需sin(√λL)=0。√λ=nπ/L,n=1,2,3,...。λ_n=(nπ/L)^2。对应的本征函数为X_n(x)=B_nsin(nπx/L)。将初始条件u(x,0)=f(x)代入u(x,t)=∑_{n=1}^∞B_nsin(nπx/L)e^(-k(nπ/L)^2t)：f(x)=∑_{n=1}^∞B_nsin(nπx/L)。由傅里叶正弦级数展开，系数B_n=(2/L)∫_0^Lf(x)sin(nπx/L)dx。最终解为u(x,t)=∑_{n=1}^∞[2/L∫_0^Lf(x)sin(nπx/L)dx]sin(nπx/L)e^(-k(nπ/L)^2t)。（提示：标准分离变量法，热传导方程，齐次边界条件，初始条件展开为正弦级数）三、5.解：考虑热传导方程的傅里叶余弦变换。设U(ρ,t)=F_c{u(x,t)}(ρ)，ρ=|x|。∂u/∂t=∂u/∂x∂x/∂ρ+∂u/∂x∂x/∂ρ=u_x(1)+u_x(-x/ρ)=u_x-(x/ρ)u_x=(1-x/ρ)u_x(若ρ>0)。但更规范的处理是直接用热传导方程u_t=ku_xx，在0<x<∞上。傅里叶余弦变换定义：F_c{u(x,t)}(ρ)=∫_0^∞u(x,t)cos(ρx)dx。F_c{u_t}=∂/∂t∫_0^∞u(x,t)cos(ρx)dx=∫_0^∞u_t(x,t)cos(ρx)dx。F_c{u_xx}=∫_0^∞u_xx(x,t)cos(ρx)dx=[u_x(x,t)cos(ρx)]_0^∞-∫_0^∞u_x(x,t)(-ρsin(ρx))dx=0-0-∫_0^∞u_x(x,t)ρsin(ρx)dx=-ρ∫_0^∞u_x(x,t)sin(ρx)dx。由分部积分，得∫_0^∞u_x(x,t)sin(ρx)dx=[u(x,t)sin(ρx)]_0^∞-∫_0^∞u(x,t)(ρcos(ρx))dx=0-ρ∫_0^∞u(x,t)cos(ρx)dx=-ρU(ρ,t)。所以F_c{u_xx}=ρ^2U(ρ,t)。代入方程F_c{u_t}=kF_c{u_xx}，得∂U/∂t=kρ^2U。解此ODE得U(ρ,t)=U(ρ,0)e^(kρ^2t)。F_c{u(x,t)}(ρ)=∫_0^∞u(x,t)cos(ρx)dx=U(ρ,0)e^(kρ^2t)。由傅里叶余弦逆变换，得U(ρ,0)=F_c{u(x,0)}(ρ)=F_c{f(x)}(ρ)=∫_0^∞f(x)cos(ρx)dx。所以∫_0^∞u(x,t)cos(ρx)dx=∫_0^∞f(x)cos(ρx)dx*e^(kρ^2t)。令v(x,t)=u(x,t)-(1/π)∫_0^∞(∫_0^∞f(ξ)cos(ρx-ρξ)dξ)e^(kρ^2t)dρ=u(x,t)-(1/π)∫_0^∞f(ξ)(∫_0^∞cos(ρ(x-ξ))e^(kρ^2t)dρ)dξ。令w(ρ,t)=∫_0^∞cos(ρ(x-ξ))e^(kρ^2t)dρ。利用已知积分结果或换元，可得w(ρ,t)=(1/2)e^(kρ^2t)(erf(ρ(x-ξ)/√(2kt))+erf(ρ(x+ξ)/√(2kt)))。因此v(x,t)=u(x,t)-(1/π)∫_0^∞f(ξ)(1/2)e^(kρ^2t)(erf(ρ(x-ξ)/√(2kt))+erf(ρ(x+ξ)/√(2kt)))dρ。最终解为u(x,t)=(1/π)∫_0^∞f(ξ)(erf(ρ(x-ξ)/√(2kt))+erf(ρ(x+ξ)/√(2kt)))dρ+(1/π)∫_0^∞(∫_0^∞f(ξ)cos(ρx-ρξ)dξ)e^(kρ^2t)dρ。（提示：热传导方程的傅里叶余弦变换，需要逆变换技巧，积分可能较复杂）6.解：此为齐次边界条件的波动方程初值问题。令u(x,t)=X(x)T(t)。分离变量得X''/X=T''/T=-ω^2。解得X=Acos(ωx)+Bsin(ωx)，T=Ccos(ωt)+Dsin(ωt)。u(x,0)=X(x)T(0)=X(x)(C+D)=(Acos(ωx)+Bsin(ωx))(C+D)=A(C+D)cos(ωx)+B(C+D)sin(ωx)=f(x)。u_t(x,0)=X(x)T'(0)=X(x)(-ωC+ωD)=ωDX(x)=g(x)。由正交性或直接比较系数，得A(C+D)=∫_0^Lf(x)cos(ωx)dx,B(C+D)=∫_0^Lf(x)sin(ωx)dx。D=g(x)/(ωX(x))。由边界条件u(0,t)=0,u(L,t)=0，得X(0)=0,X(L)=0。即B(C+D)sin(ωL)=0,A(C+D)cos(ωL)=0。由于A,B,C,D不全为0，需sin(ωL)=0,cos(ωL)=0。即ωL=nπ,n=1,2,3,...。ω_n=nπ/L。此时X_n(x)=B_nsin(nπx/L)。对应的D_n=g_n/(ω_nX_n(x))=g_n/(nπ/L)sin(nπx/L)。由正弦级数展开，系数B_n=(2/L)∫_0^Lf(x)sin(nπx/L)dx。最终解为u(x,t)=∑_{n=1}^∞B_nsin(nπx/L)(D_ncos(ω_nt)+E_nsin(ω_nt))，其中E_n=D_n/ω_n。u(x,t)=∑_{n=1}^∞[2/L∫_0^Lf(x)sin(nπx/L)dx]sin(nπx/L)(g_n/(n^2π^2)cos(nπt/L)+g_n/(n^3π^3)sin(nπt/L))。u(x,t)=(1/L)∑_{n=1}^∞(g_n/n^2π^2)sin(nπx/L)(cos(nπt/L)+(n/ω_n)sin(nπt/L))。（提示：齐次边界条件的波动方程，初始条件非齐次，需先展开初始位移和初始速度为正弦级数）四、7.解：此为一阶线性偏微分方程的初值问题。方程可写为(yu_x+xu_y)=x^2+y^2。求特征方程：dx/dy=y/x，即xdx=ydy。积分得x^2/2=y^2/2+C，即x^2-y^2=2C。特征线族为x^2-y^2=C。初值线为x=1,y=1，代入得C=0。所以通过点(1,1)的特征线为x^2-y^2=0，即y=x。沿特征线，设v=u(x,y)，则v满足v_t=(x^2+y^2)/(yu_x+xu_y)沿特征线。在特征线y=x上，u=v，u_x=v_x,u_y=v_y。沿此线，x=y，dx=dy。代入方程得v_x+v_y=x^2+y^2=2x^2。沿y=x，此为一阶常微分方程v'=2x。初始条件为当x=y时，x=1,v=u(1,1)=1+1=2。积分得v=x^2+C。当x=1时，v=2，得C=1。所以v=x^2+1。因为沿特征线y=x，v=u，所以u(x,y)=x^2+1。（提示：标准特征线法，求出特征线族，确定通过初值点的特征线，将问题转化为沿特征线的常微分方程求解）8.证明：令u(x,y)=u(ρ,θ)，其中ρ=√(x^2+y^2),θ=arctan(y/x)。拉普拉斯方程在极坐标下为u_rr+(1/ρ)u_r+(1/(ρ^2))u_θθ=0。边界条件为ρ=1时，u(ρ,θ)=1/x=cos(θ)/1=cos(θ)。若解存在，则u在圆盘D内连续可微。设ρ=1时，u(1,θ)=cos(θ)。由于圆盘不包含原点(0,0)，u在D内无奇点。考虑调和函数的性质：在无源区域（无奇点）内，调和函数的值由其边界值唯一确定（调和函数的保号性或唯一性定理）。因此，如果u在D内存在解，则u在D内必须恒等于边界上的值，即u(ρ,θ)=cos(θ)/ρ。但ρ=1时，u(1,θ)=cos(θ)，而题目给定的边界条件是u(ρ,θ)=1/x=cos(θ)/1=cos(θ)。这与调和函数的唯一性矛盾（除非解处处为常数）。另一种思考方式：如果u在D内存在非零解，则u(ρ,θ)=C/ρ(由调和函数性质)。边界上u(1,θ)=cos(θ)。这意味着C=cos(θ)，C是θ的函数，与C为常数矛盾。因此，u在D内无解。（提示：利用调和函数在无奇点区域内的唯一性定理，或直接计算极坐标下的拉普拉斯方程，结合边界条件分析）五、9.(1)证明：令u(x,t)=X(x)T(t)。分离变量得X''/X=T''/T=-λ。解得X=Acos(√λx)+Bsin(√λx)，T=Ce^(-λkt)。边界条件u(0,t)=0,u(L,t)=0。X(0)T(t)=0,X(L)T(t)=0。需X(0)=0,X(L)=0。X(0)=Acos(0)+Bsin(0)=A=0。X(L)=Bsin(√λL)=0。由于B≠0(否则解恒为0)，需sin(√λL)=0。√λ=nπ/L,n=1,2,3,...。λ_n=(nπ/L)^2。对应的本征函数为X_n(x)=B_nsin(nπx/L)。将初始条件u(x,0)=f(x)代入u(x,t)=∑_{n=1}^∞B_nsin(nπx/L)e^(-λ_nkt)：f(x)=∑_{n=1}^∞B_nsin(nπx/L)。由傅里叶正弦级数展开，系数B_n=(2/L)∫_0^Lf(x)sin(nπx/L)dx。最终解为u(x,t)=∑_{n=1}^∞[2/L∫_0^Lf(x)sin(nπx/L)dx]sin(nπx/L)e^(-k(nπ/L)^2t)。这个解可以分解为两部分：u(x,t)=∑_{n=1}^∞[2/L∫_0^Lf(x)sin(nπx/L)dx]sin(nπx/L)e^(-k(nπ/L)^2t)+T(x)。其中T(x)=T_1+(T_2-T_1)x。T(x)显然满足齐次边界条件u(0,t)=T_1,u(L,t)=T_2，并且是时间的常数函数，满足齐次热传导方程u_t=ku_xx(因为T_t=0,T_xx=0)。所以u(x,t)=T(x)+∑_{n=1}^∞[2/L∫_0^Lf(x)sin(nπx/L)dx]sin(nπx/L)e^(-k(nπ/L)^2t)。（2）若T_1=T_2=0，即问题变为齐次边界条件的热传导，解为u(x,t)=∑_{n=1}^∞[2/L∫_0^Lf(x)sin(nπx/L)dx]sin(nπx/L)e^(-k(nπ/L)^2t)。其中T(x)=0。（提示：标准分离变量法，处理齐次边界条件得到本征值本征函数。对于非齐次边界条件，解可以分解为满足边界条件的齐次解（常数项）和满足热传导方程的齐次边界条件解（傅里叶级数部分）。）10.解：有限差分法是数值求解偏微分方程的一种方法。其基本思想是将连续的偏微分方程和连续的空间、时间变量，用离散的网格点（如矩形网格）和有限差分格式来近似。以一维热传导方程u_t=ku_xx,0<x<L,t>0,u(0,t)=0,u(L,t)=0,u(x,0)=f(x)为例。(1)离散化：将空间[0,L]分成N等份，步长Δx=L/N。令x_i=iΔx,i=0,1,...,N。将时间[0,T]分成M等份，步长Δt=T/M。令t_n=nΔt,n=0,1,...,M。(2)差分格式：在网格点(x_i,t_n)处近似u的值。对u_t使用向前差分：u_t(x_i,t_n)≈(u(x_i,t_{n+1})-u(x_i,t_n))/Δt。对u_xx使用中心差分：u_xx(x_i,t_n)≈(u(x_{i+1},t_n)-2u(x_i,t_n)+u(x_{i-1},t_n))/Δx^2。(3)建立方程：将差分格式代入原方程，得到：(u(x_i,t_{n+1})-u(x_i,t_n))/Δt=k(u(x_{i+1},t_n)-2u(x_i,t_n)+u(x_{i-1},t_n))/Δx^2。(4)迭代求解：整理上式，得到隐式差分格式：u(x_i,t_{n+1})=[1-(2kΔt/Δx^2)]u(x_i,t_n)+(kΔt/Δx^2)[u(x_{i+1},t_n)+u(x_{i-1},t_n)]。这是一个关于u在时间层t_{n+1}的隐式方程。对于每个i=1,...,N-1，得到N-1个方程，包含N-1个未知数u(x_1,t_{n+1}),...,u(x_{N-1},t_{n+1})。需要结合边界条件u(0,t)=0,u(L,t)=0，即u(x_0,t)=0,u(x_N,t)=0。对于第一个时间步t_1(n=0)，需要初始条件u(x_i,0)=f(x_i)=f(x_iΔx)。将这些值作为初始猜测，代入隐式格式，通过迭代方法（如高斯消元法）求解整个时间层t_1的解。对于后续时间步t_2,t_3,...,可以利用已知的解u(x_i,t_{n}),u(x_{i+1},t_{n}),u(x_{i-1},t_{n})，直接求解当前时间步t_{n+1}的值。(5)稳定性与收敛性：有限差分格式的使用需要考虑其稳定性和收敛性。例如，对于热传导方程的显式格式和隐式格式，有严格的稳定性条件（如CFL条件）。（提示：概述有限差分法的基本步骤：离散化、差分格式、建立方程、迭代求解，并以热传导方程一维定解问题为例进行说明，提及稳定性问题）六、11.证明：在直角坐标系下，拉普拉斯方程u_xx+u_yy=0具有分离变量解。令u(x,y)=X(x)Y(y)。代入方程得X''Y+XY''=0，即X''/X=-Y''/Y=λ。解得Y=Asin(√λy)+Bcos(√λy)。在极坐标系(r,θ)中，拉普拉斯方程为u_rr+(1/r)u_r+(1/r^2)u_θθ=0。令u(r,θ)=R(r)Θ(θ)。代入方程得R''Θ+(1/r)R'Θ+(1/r^2)RΘ''=0。分离变量R''/R+(1/r)R'/R=-Θ''/Θ=λ。即R''/R+(1/r)R'/R=λ，Θ''/Θ=-λ。解Θ：若λ≥0，设λ=μ^2。则Θ''-μ^2Θ=0。解为Θ=Ae^(μθ)+Be^(-μθ)。由于解在θ的全平面上（除去原点）需有界，需A=B=0，或考虑μ=inπ(n为整数)。得Θ=Acos(nπθ)+Bsin(nπθ)。解R：R''+(1/r)R'+n^2π^2R=0。这是一个贝塞尔方程（或可转化为贝塞尔方程）。解为R=CJ_n(μr)+DY_n(μr)=CJ_n(nπr)+DY_n(nπr)。其中J_n为第一类贝塞尔函数，Y_n为第二类贝塞尔函数。Y_n在原点发散，通常在物理问题中不考虑。所以解为u(r,θ)=∑_{n=1}^∞[A_nJ_n(nπr)+B_nsin(nπθ)]。初始条件或边界条件通常在极坐标下给出。例如，若边界条件为u(r=1,θ)=f(θ)，则u(r,θ)=∑_{n=1}^∞[C_nJ_n(nπr)+D_nsin(nπθ)]，其中C_n,D_n由边界条件确定。如果要求解在r=1时为u(1,θ)=f(θ)=1/x=cos(θ)/1=cos(θ)。即C_nJ_n(nπ)+D_nsin(nπθ)=cos(θ)。这要求C_n=0(因为J_n(nπ)≠0)，D_nsin(nπθ)=cos(θ)。这意味着D_nsin(nπθ)是cos(θ)的傅里叶级数展开。这通常需要使用欧拉公式或引入复变量方法。例如，令u(r,θ)=R(r)Θ(θ)。若考虑u(r,θ)=R(r)e^(iθ)=R(r)(因为e^(iθ)=cos(θ)+isin(θ)，实部解为R(r)cos(θ)。但题目要求解为u(r,θ)=∑[C_nJ_n(nπr)+D_nsin(nπθ)]，这更像是针对r=1时u(1,θ)=f(θ)=cos(θ)的解法推导，而非直接求解u(r,θ)满足u(r=1,θ)=f(θ)=cos(θ)。更准确的极坐标下拉普拉斯方程的解法通常需要引入特殊函数（如贝塞尔函数）。（提示：分离变量法在极坐标下的处理比直角坐标复杂，需要引入贝塞尔函数等特殊函数。题目要求的解形式u(r,θ)=∑[C_nJ_n(nπr)+D_nsin(nπθ)]，可能是在r=1时u(1,θ)=f(θ)=cos(θ)的解法推导或形式表达。标准的极坐标下拉普拉斯方程的解法通常涉及贝塞尔函数。）12.偏微分方程是描述自然现象和工程系统中的基本数学模型工具，在理解物理现象（如波的传播、热扩散、电磁场分布、流体运动等）和解决工程问题（如结构力学、量子力学、热传导问题、波动问题）中起着核心作用。例如：*波动方程可用于描述弦