-
2021-04-22 18:11:53
二阶非线性常微分方程的打靶法
二阶非线性常微分方程的打靶法 计算思路 主要分为以下五步: 给定容许误差ε,迭代初始值γ1,对k=1,2,...做: (1)用四阶Runge-Kutta 方法求解初值问,得出u1之后取其在γk的值,从而得到F(γk)=u1(b,γk)-β (2)若|F(γk)|<ε,则u1 即为所求,跳出循环。 否则: (3)用四阶Runge-Kutta 方法求解初值问题 由此得到F'(γk)=v1(b,γk) (4)用牛顿迭代计算γk+1,即 (5)置k+1→k,转(1)直到误差在范围内。 综上,实现了打靶法对非线性方程的拟合。 接下来是每一步的代码: 程序开头各变量的设置 function ys=ndbf(f,g,a,b,alfa,beta,n,eps,s0) %f 为二阶导数,y''=f(x,y,y'),g 为f 对y 求偏导后的 %(a,b)为自变量迭代区间 %alfa,beta 为给定的边值条件 %eps 题目规定的精度 %n 为迭代次数 %选取适当的s0 的初值 循环第一步 x0=[alfa,s0]; %选取合适的迭代初值 y0=RK4(f,a,b,h,x0); %龙格库塔算出u1(γk) constant=y0(n,1)-beta; 这里的y0(n,1)即是u1(γk,b),constant即为F(γk) 循环第二步 检验误差以跳出 if abs(constant)
更多相关内容 -
二阶非线性偏微分方程之间Backlund变换的分类 (2008年)
2021-05-30 09:03:16讨论二阶非线性偏微分方程zxy=F(z)与φxy=G(φ)(其中函数F和G未必相同)之间由可积系统φ=P(φ,z,zx,zy),φy=Q(φ,z,zx,zy)定义的Backlund变换z→φ的分类问题,得到函数F和G的具体表达式,并且给出... -
二阶非线性偏微分方程Dirichlet问题粘性解的存在性和唯一性 (2001年)
2021-05-20 17:30:17运用粘性解理论研究了二阶非线性椭圆偏微分方程Dirichlet问题解的存在性和唯一性。首先建立比较定理,保证了解的唯一性;然后运用Perron方法构造解,保证了解的存在性。通过以上结果,解的存在唯一性和存在性得以解决。 -
完整word版matlab解非线性规划例题_非线性偏微分方程例题
2020-08-14 17:14:16关于非线性规划问题 背景 线性规划问题即目标函数和约束条件都是线性函数的规划问题但在实际工作中还常常会遇到另一类更一般的规划问题即目标函数和约束条件中至少有一个是非线性函数问题即非线性规划问题 求解方法... -
一类二阶非线性偏微分方程Baecklund变换的分类 (2011年)
2021-04-26 12:05:50"讨论形如 ut = F(u,ux,uxx)的非线性偏微分方程由可积系统 vx = P(v,u,ux),vt = Q( v,u,ux)定义的 -
二阶非线性偏微分椭圆方程Neurnann问题的粘性解 (2001年)
2021-05-13 23:13:00本文运用粘性解理论解决了一类非线性椭回方程的Neumann问题。首先建立比较定理,然后用Perron方法构造解,从而得到解的存在唯一性。 -
Matlab求解非线性偏微分方程组 - 数学 - 小木虫 - 学术 科研 互动社区
2021-04-23 23:34:33赡芄酵频加写砦螅绻凑漳愕汲龅墓郊扑悖峒扑愠鲂橹担缓蟛煌5乃姥贰H缃玴u3/pt的最后k*Ac前面的负号改为正号则有解。从你的边界条件来看,似乎这个分量也应该跟其他两个应该是不一样的。...赡芄酵频加写砦螅绻凑漳愕汲龅墓郊扑悖峒扑愠鲂橹担缓蟛煌5乃姥贰H缃玴u3/pt的最后k*Ac前面的负号改为正号则有解。从你的边界条件来看,似乎这个分量也应该跟其他两个应该是不一样的。
function [ output_args ] = PDEs1( input_args )
m=0;
x=linspace(0,0.02,3);
t=linspace(0,1,3);
sol=pdepe(m,@pdes1pde,@pdes1ic,@pdes1bc,x,t);
u1=sol(:,:,1);
u2=sol(:,:,2);
u3=sol(:,:,3);
figure
surf(x,t,u1)
title('u1(x,t)')
xlabel('Distance x')
ylabel('Time t')
figure
surf(x,t,u2)
title('u2(x,t)')
xlabel('Distance x')
ylabel('Time t')
figure
surf(x,t,u3)
title('u3(x,t)')
xlabel('Distance x')
ylabel('Time t')
function[c,f,s]=pdes1pde(x,t,u,DuDx)
c=[1;1;1];
D1=1.25*10^(-9);
D2=1.838*10^(-9);
D3=1.856*10^(-9);
V=0.59;
k=1.03*10^(-4)/60;
Ac=1.9;
Ksp=5.51*10^(-14);
f=[D1/(1+V);D2/(1+V);D3/(1+V)].*DuDx;
s=[-k*Ac/(1+V)*((u(1)*u(2)*u(3)/Ksp)^(1/3)-1);-k*Ac/(1+V)*((u(1)*u(2)*u(3)/Ksp)^(1/3)-1);k*Ac/(1+V)*((u(1)*u(2)*u(3)/Ksp)^(1/3)-1);];
end
function u0 = pdes1ic(x)
u0 = [0;0;0];
end
function [pl,ql,pr,qr] = pdes1bc(xl,ul,xr,ur,t)
pl = [ul(1)-0.0493; ul(2)-0.07142;ul(3)-0.0387];
ql = [0;0;0];
pr = [ur(1)-0.0106;ur(2)-0.03272;ur(3)];
qr = [0;0;0];
end
end,
-
二阶线性偏微分方程的分类和标准式 | 椭圆型、抛物线形、双曲线型 | 偏微分方程(十一)
2020-05-07 17:36:45一般地,n个自变量的二阶线性偏微分方程可表示为 ∑i,j=1naij(x1,⋅⋅⋅,xn)∂2u∂xi∂xj+∑j=1nbj(x1,⋅⋅⋅,xn)∂u∂xj+c(x1,⋅⋅⋅,xn)u=f(x1,⋅⋅⋅,xn)(1) \sum_{i,j=1}^na_{ij}(x_1,···,x_n)\frac{\partial...一般地,n个自变量的二阶线性偏微分方程可表示为
∑ i , j = 1 n a i j ( x 1 , ⋅ ⋅ ⋅ , x n ) ∂ 2 u ∂ x i ∂ x j + ∑ j = 1 n b j ( x 1 , ⋅ ⋅ ⋅ , x n ) ∂ u ∂ x j + c ( x 1 , ⋅ ⋅ ⋅ , x n ) u = f ( x 1 , ⋅ ⋅ ⋅ , x n ) (1) \sum_{i,j=1}^na_{ij}(x_1,···,x_n)\frac{\partial^2 u}{\partial x_i\partial x_j}+\sum_{j=1}^nb_j(x_1,···,x_n)\frac{\partial u}{\partial x_j}+c(x_1,···,x_n)u=f(x_1,···,x_n) \tag{1} i,j=1∑naij(x1,⋅⋅⋅,xn)∂xi∂xj∂2u+j=1∑nbj(x1,⋅⋅⋅,xn)∂xj∂u+c(x1,⋅⋅⋅,xn)u=f(x1,⋅⋅⋅,xn)(1)
当系数 a i j , b j , c a_{ij},b_j,c aij,bj,c均为常数时,称为常系数线性偏微分方程,否则为变系数的,且总可假定 a i j = a j i a_{ij}=a_{ji} aij=aji。方程中不含未知函数的项 f ( x 1 , ⋅ ⋅ ⋅ , x n ) f(x_1,···,x_n) f(x1,⋅⋅⋅,xn)称为非齐次项,当 f ≡ 0 f\equiv 0 f≡0时方程为二阶齐次线性偏微分方程,否则为非齐次的。作为(1)式的特殊情况,一维波动方程能分解为两个一阶线性偏微分方程,从而利用特征线求出其通解。对于一般的二阶线性偏微分方程(1),接下来将通过自变量的变量代换化简方程中的二阶偏导数部分,并用方程在变量代换下的不变性对方程进行分类。这里只讨论两个自变量的情形。
1. 特征方程和特征线
两个自变量的二阶齐次线性偏微分方程为
a 11 ∂ 2 u ∂ x 2 + 2 a 12 ∂ 2 u ∂ x ∂ y + a 22 ∂ 2 u ∂ y 2 + b 1 ∂ u ∂ y + b 2 ∂ u ∂ y + c u = 0 (2) a_{11}\frac{\partial^2u}{\partial x^2}+2a_{12}\frac{\partial^2u}{\partial x\partial y}+a_{22}\frac{\partial^2u}{\partial y^2}+b_1\frac{\partial u}{\partial y}+b_2\frac{\partial u}{\partial y}+cu=0 \tag{2} a11∂x2∂2u+2a12∂x∂y∂2u+a22∂y2∂2u+b1∂y∂u+b2∂y∂u+cu=0(2)
其中, a 11 , a 12 , a 22 a_{11},a_{12},a_{22} a11,a12,a22不同时为0,这里略写了已知函数 a i j , b j ( i , j = 1 , 2 ) a_{ij},b_j(i,j=1,2) aij,bj(i,j=1,2)和c的自变量。设有自变量的变量代换
{ ξ = φ ( x , y ) , η = ψ ( x , y ) , \begin{cases} \xi=\varphi(x,y), \\ \eta=\psi(x,y), \end{cases} {ξ=φ(x,y),η=ψ(x,y),
其中,函数 φ , ψ \varphi,\psi φ,ψ有二阶连续偏导,雅可比(Jacobi)行列式
J = ∂ ( φ , ψ ) ∂ ( x , y ) = ∣ ∂ φ ∂ x ∂ φ ∂ y ∂ ψ ∂ x ∂ ψ ∂ y ∣ ≠ 0 J=\frac{\partial (\varphi,\psi)}{\partial(x,y)}= \begin{vmatrix} \frac{\partial \varphi}{\partial x} & \frac{\partial \varphi}{\partial y} \\ \frac{\partial \psi}{\partial x} & \frac{\partial \psi}{\partial y} \end{vmatrix} \neq 0 J=∂(x,y)∂(φ,ψ)=∣∣∣∣∣∂x∂φ∂x∂ψ∂y∂φ∂y∂ψ∣∣∣∣∣=0
以保证函数 φ , ψ \varphi,\psi φ,ψ相互独立,有反函数存在。未知函数u作为新的自变量 ξ , η \xi,\eta ξ,η的函数仍记为 u ( ξ , η ) u(\xi,\eta) u(ξ,η),利用复合函数求导的链式法则,可将方程(2)变成一个新的二阶线性偏微分方程
A 11 ∂ 2 u ∂ ξ 2 + 2 A 12 ∂ 2 u ∂ ξ ∂ η + A 22 ∂ 2 u ∂ η 2 + B 1 ∂ u ∂ ξ + B 2 ∂ u ∂ η + C u = 0 (a) A_{11}\frac{\partial^2u}{\partial \xi^2}+2A_{12}\frac{\partial^2u}{\partial \xi\partial \eta}+A_{22}\frac{\partial^2u}{\partial \eta^2}+B_1\frac{\partial u}{\partial \xi}+B_2\frac{\partial u}{\partial \eta}+Cu=0 \tag{a} A11∂ξ2∂2u+2A12∂ξ∂η∂2u+A22∂η2∂2u+B1∂ξ∂u+B2∂η∂u+Cu=0(a)
其中,二阶偏导数项的系数
A 11 = a 11 ( ∂ φ ∂ x ) 2 + 2 a 12 ∂ φ ∂ x ∂ φ ∂ y + a 22 ( ∂ φ ∂ y ) 2 , A 12 = a 11 ∂ φ ∂ x ∂ ψ ∂ x + 2 a 12 ( ∂ φ ∂ x ∂ φ ∂ y + ∂ φ ∂ y ∂ φ ∂ x ) + a 22 ∂ φ ∂ y ∂ ψ ∂ y ( 3 ) A 11 = a 11 ( ∂ ψ ∂ x ) 2 + 2 a 12 ∂ ψ ∂ x ∂ ψ ∂ y + a 22 ( ∂ ψ ∂ y ) 2 A_{11}=a_{11}(\frac{\partial \varphi}{\partial x})^2+2a_{12}\frac{\partial \varphi}{\partial x}\frac{\partial \varphi}{\partial y}+a_{22}(\frac{\partial \varphi}{\partial y})^2 ,\\ A_{12}=a_{11}\frac{\partial \varphi}{\partial x}\frac{\partial \psi}{\partial x}+2a_{12}(\frac{\partial \varphi}{\partial x}\frac{\partial \varphi}{\partial y}+\frac{\partial \varphi}{\partial y}\frac{\partial \varphi}{\partial x})+a_{22}\frac{\partial \varphi}{\partial y}\frac{\partial \psi}{\partial y} \quad \quad (3)\\ A_{11}=a_{11}(\frac{\partial \psi}{\partial x})^2+2a_{12}\frac{\partial \psi}{\partial x}\frac{\partial \psi}{\partial y}+a_{22}(\frac{\partial \psi}{\partial y})^2 A11=a11(∂x∂φ)2+2a12∂x∂φ∂y∂φ+a22(∂y∂φ)2,A12=a11∂x∂φ∂x∂ψ+2a12(∂x∂φ∂y∂φ+∂y∂φ∂x∂φ)+a22∂y∂φ∂y∂ψ(3)A11=a11(∂x∂ψ)2+2a12∂x∂ψ∂y∂ψ+a22(∂y∂ψ)2
其他系数可推导而得。为化简方程中的二阶偏导数部分,由(3)式可见,若选取 φ ( x , y ) \varphi(x,y) φ(x,y)或 ψ ( x , y ) \psi(x,y) ψ(x,y)是一阶非线性偏微分方程
a 11 ( ∂ z ∂ x ) 2 + 2 a 12 ∂ z ∂ x ∂ z ∂ y + a 22 ( ∂ z ∂ y ) 2 = 0 (4) a_{11}(\frac{\partial z}{\partial x})^2+2a_{12}\frac{\partial z}{\partial x}\frac{\partial z}{\partial y}+a_{22}(\frac{\partial z}{\partial y})^2=0 \tag{4} a11(∂x∂z)2+2a12∂x∂z∂y∂z+a22(∂y∂z)2=0(4)
的解,则新方程中 A 11 A_{11} A11与 A 22 A_{22} A22至少有一个为0,方程化简。类似于之前对一阶线性偏微分方程的讨论,一阶偏微分方程(4)的求解归结为一阶常微分方程
a 11 ( d y ) 2 − 2 a 12 d x d y + a 22 ( d x ) 2 = 0 (5) a_{11}(dy)^2-2a_{12}dxdy+a_{22}(dx)^2=0 \tag{5} a11(dy)2−2a12dxdy+a22(dx)2=0(5)
的求解。相应地,有定理1:若 φ ( x , y ) = h \varphi(x,y)=h φ(x,y)=h(h为常数)是常微分方程(5)的隐式通解(积分曲线族),则 z = φ ( x , y ) z=\varphi(x,y) z=φ(x,y)是偏微分方程(4)的解。定理的逆命题也成立。
上述定理揭示了常微分方程(5)与二阶线性偏微分方程(2)之间的关系,提供了化简二阶线性偏微分方程的具体方法。称一阶常微分方程(5)为二阶线性偏微分方程(2)的特征方程,称特征方程的积分曲线为方程(2)的特征曲线。
2. 方程的分类、化简和标准形
特征方程(4)的解取决于它的判别式
Δ = a 12 2 − a 11 a 22 \Delta= a_{12}^2-a_{11}a_{22} Δ=a122−a11a22
由(3)式可推出,二阶线性偏微分方程(2)经过变量代换后得到的新方程式的判别式为
A 12 2 − A 11 A 22 = J 2 Δ A_{12}^2-A_{11}A_{22}=J^2\Delta A122−A11A22=J2Δ
当Jacobi行列式 J ≠ 0 J\neq 0 J=0时,在自变量的变量代换下,判别式 Δ \Delta Δ的符号保持不变,反映了方程的本性,据此可将方程分类。在 ( x 0 , y 0 ) (x_0,y_0) (x0,y0)点,若 Δ > 0 \Delta>0 Δ>0,则称二阶线性偏微分方程(2)式在 ( x 0 , y 0 ) (x_0,y_0) (x0,y0)点为双曲型的;若 Δ = 0 \Delta=0 Δ=0则称(2)式在 ( x 0 , y 0 ) (x_0,y_0) (x0,y0)为抛物型的;若 Δ < 0 \Delta<0 Δ<0,则称(2)式在 ( x 0 , y 0 ) (x_0,y_0) (x0,y0)为椭圆型的。若在平面区域D内有 Δ > 0 \Delta>0 Δ>0或 Δ = 0 \Delta=0 Δ=0,或 Δ < 0 \Delta<0 Δ<0,则相应地称方程(2)在区域D内为双曲型、抛物型或椭圆型方程。若方程在区域D的一部分是双曲型的,另一部分是椭圆型的,而在交界线上是抛物型的,则称该方程在D内是混合型方程。
易见,弦振动方程
∂ 2 u ∂ t 2 = a 2 ∂ 2 u ∂ x 2 \frac{\partial^2u}{\partial t^2}=a^2\frac{\partial^2u}{\partial x^2} ∂t2∂2u=a2∂x2∂2u
的 Δ ≡ a 2 > 0 \Delta\equiv a^2>0 Δ≡a2>0一维热传导方程
∂ u ∂ t = a 2 ∂ 2 u ∂ x 2 \frac{\partial u}{\partial t}=a^2\frac{\partial^2u}{\partial x^2} ∂t∂u=a2∂x2∂2u
的 Δ ≡ 0 \Delta \equiv 0 Δ≡0二维拉普拉斯方程
∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2 = 0 \frac{\partial^2 u}{\partial x^2}+\frac{\partial^2u}{\partial y^2}=0 ∂x2∂2u+∂y2∂2u=0
的 Δ = − 1 < 0 \Delta=-1<0 Δ=−1<0,它们在全平面上依次为双曲型方程、抛物型方程和椭圆型方程。对于三种类型的二阶线性偏微分方程,可按以下步骤化简。
当 Δ > 0 \Delta>0 Δ>0,特征方程(5)可分解为两个一阶常微分方程。不妨设 a 11 ≠ 0 a_{11}\neq 0 a11=0,此时,有
d y d x = a 12 + Δ a 11 , d y d x = a 12 − Δ a 11 \frac{dy}{dx}=\frac{a_{12}+\sqrt{\Delta}}{a_{11}},\quad \frac{dy}{dx}=\frac{a_{12}-\sqrt{\Delta}}{a_{11}} dxdy=a11a12+Δ,dxdy=a11a12−Δ
解得两族特征线
φ ( x , y ) = h 1 , ψ ( x , y ) = h 2 \varphi(x,y)=h_1,\quad \psi(x,y)=h_2 φ(x,y)=h1,ψ(x,y)=h2
令 ξ = φ ( x , y ) , η = ψ ( x , y ) \xi=\varphi(x,y),\eta =\psi(x,y) ξ=φ(x,y),η=ψ(x,y),由定理1和公式(3)知 A 11 = A 22 = 0 A_{11}=A_{22}=0 A11=A22=0,新方程(a)有简单形式
∂ 2 u ∂ ξ ∂ η + 1 2 A 12 ( B 1 u ξ + B 2 u η + C u ) = 0 (7) \frac{\partial^2u}{\partial \xi\partial \eta}+\frac{1}{2A_{12}}(B_1u_{\xi}+B_2u_{\eta}+Cu)=0 \tag{7} ∂ξ∂η∂2u+2A121(B1uξ+B2uη+Cu)=0(7)
因为沿两族特征线分别有
∂ φ ∂ x d x + ∂ φ ∂ y d y = 0 , ∂ ψ ∂ x d x + ∂ ψ ∂ y d y = 0 \frac{\partial \varphi}{\partial x}dx+\frac{\partial \varphi}{\partial y}dy=0, \quad \frac{\partial \psi}{\partial x}dx+\frac{\partial \psi}{\partial y}dy=0 ∂x∂φdx+∂y∂φdy=0,∂x∂ψdx+∂y∂ψdy=0
从而
∂ φ ∂ x / ∂ φ ∂ y = − a 12 + Δ a 11 , ∂ φ ∂ x / ∂ ψ ∂ y = − a 12 − Δ a 11 \frac{\partial \varphi}{\partial x}/\frac{\partial \varphi}{\partial y}=-\frac{a_{12}+\sqrt{\Delta}}{a_{11}},\quad \frac{\partial \varphi}{\partial x}/\frac{\partial \psi}{\partial y}=-\frac{a_{12}-\sqrt{\Delta}}{a_{11}} ∂x∂φ/∂y∂φ=−a11a12+Δ,∂x∂φ/∂y∂ψ=−a11a12−Δ
由 a 11 ≠ 0 a_{11}\neq 0 a11=0可知 ∂ φ ∂ y ∂ ψ ∂ y ≠ 0 \frac{\partial \varphi}{\partial y}\frac{\partial \psi}{\partial y}\neq 0 ∂y∂φ∂y∂ψ=0。从雅克比行列式算得
J = − 2 Δ a 11 ∂ φ ∂ y ∂ ψ ∂ y ≠ 0 J=-\frac{2\sqrt{\Delta}}{a_{11}}\frac{\partial \varphi}{\partial y}\frac{\partial \psi}{\partial y}\neq 0 J=−a112Δ∂y∂φ∂y∂ψ=0
故 ξ , η \xi,\eta ξ,η是两个独立的变量,同理由式(3)可算得
A 12 = − 2 Δ a 11 ∂ φ ∂ y ∂ ψ ∂ y ≠ 0 A_{12}=-\frac{2\Delta}{a_{11}}\frac{\partial \varphi}{\partial y}\frac{\partial \psi}{\partial y}\neq 0 A12=−a112Δ∂y∂φ∂y∂ψ=0
而 B 1 , B 2 , C B_1,B_2,C B1,B2,C也可由推导出的公式算出若对式(7)再作变量代换
s = 1 2 ( ξ + η ) , t = 1 2 ( ξ − η ) s=\frac{1}{2}(\xi+\eta),\quad t=\frac{1}{2}(\xi-\eta) s=21(ξ+η),t=21(ξ−η)
可得方程
∂ 2 u ∂ t 2 − ∂ 2 u ∂ s 2 + 1 A ^ 11 ( B ^ 1 u t + B ^ 2 u s + C ^ u ) = 0 (8) \frac{\partial^2u}{\partial t^2}-\frac{\partial^2u}{\partial s^2}+\frac{1}{\hat A_{11}}(\hat B_1u_t+\hat B_2u_s+ \hat Cu)=0 \tag{8} ∂t2∂2u−∂s2∂2u+A^111(B^1ut+B^2us+C^u)=0(8)
方程(7)和(8)都称为双曲型方程的标准型。当 Δ = 0 \Delta=0 Δ=0时,由特征方程(5)只能得到一个一阶线性常微分方程。不妨设 a 11 ≠ 0 a_{11}\neq 0 a11=0,该方程为
d y d x = a 12 a 11 \frac{dy}{dx}=\frac{a_{12}}{a_{11}} dxdy=a11a12
解得一族特征曲线 φ ( x , y ) = h \varphi(x,y)=h φ(x,y)=h。任取二元函数 ψ ( x , y ) \psi(x,y) ψ(x,y),使得 J = ∂ ( φ , ψ ) ∂ ( x , y ) ≠ 0 J=\frac{\partial (\varphi,\psi)}{\partial(x,y)}\neq 0 J=∂(x,y)∂(φ,ψ)=0。作变量代换
ξ = φ ( x , y ) , η = ψ ( x , y ) \xi=\varphi(x,y),\quad \eta=\psi(x,y) ξ=φ(x,y),η=ψ(x,y)
可得 A 11 = 0 , A 22 ≠ 0 A_{11}=0,A_{22}\neq 0 A11=0,A22=0。因 a 12 2 = a 11 a 22 a_{12}^2=a_{11}a_{22} a122=a11a22, ∂ φ ∂ x / ∂ φ ∂ y = − a 22 a 11 \frac{\partial \varphi}{\partial x}/\frac{\partial \varphi}{\partial y}=-\frac{a_{22}}{a_{11}} ∂x∂φ/∂y∂φ=−a11a22,从式(3)可得
A 12 = 1 a 11 ( a 11 ∂ φ ∂ x + a 12 ∂ φ ∂ y ) ( a 11 ∂ ψ ∂ x + a 12 ∂ ψ ∂ y ) = 0 A_{12}=\frac{1}{a_{11}}(a_{11}\frac{\partial \varphi}{\partial x}+a_{12}\frac{\partial \varphi}{\partial y})(a_{11}\frac{\partial \psi}{\partial x}+a_{12}\frac{\partial \psi}{\partial y})=0 A12=a111(a11∂x∂φ+a12∂y∂φ)(a11∂x∂ψ+a12∂y∂ψ)=0
新方程为
∂ 2 u ∂ η 2 + 1 A 22 ( B 1 ∂ u ∂ ξ + B 2 ∂ u ∂ η + C u ) = 0 \frac{\partial^2 u}{\partial \eta^2}+\frac{1}{A_{22}}(B_1\frac{\partial u}{\partial \xi}+B_2\frac{\partial u}{\partial \eta}+Cu)=0 ∂η2∂2u+A221(B1∂ξ∂u+B2∂η∂u+Cu)=0
称为抛物线方程的标准形。当 Δ < 0 \Delta<0 Δ<0时,特征方程(5)只能在复数域内分解成两个一阶方程。不妨设 a 11 ≠ 0 a_{11}\neq 0 a11=0,相应的一阶方程为
d y d x = a 12 + i − Δ a 11 和 d y d x = a 12 − i − Δ a 11 \frac{dy}{dx}=\frac{a_{12}+i\sqrt{-\Delta}}{a_{11}}\quad 和 \quad \frac{dy}{dx}=\frac{a_{12}-i\sqrt{-\Delta}}{a_{11}} dxdy=a11a12+i−Δ和dxdy=a11a12−i−Δ
此时,不存在实特征线,特征方程(5)的隐式通解为
φ ( x , y ) ± i ψ ( x , y ) = h \varphi(x,y)\pm i\psi(x,y)=h φ(x,y)±iψ(x,y)=h
为避免引入复变量,作变换
{ ξ = φ ( x , y ) η = ψ ( x , y ) \begin{cases} \xi=\varphi(x,y) \\ \eta=\psi(x,y) \end{cases} {ξ=φ(x,y)η=ψ(x,y)
由定理1将 z = φ ( x , y ) ± i ψ ( x , y ) z=\varphi(x,y)\pm i\psi(x,y) z=φ(x,y)±iψ(x,y)代入式(4),分别取实、虚部,可得 A 11 = A 22 , A 12 = 0 A_{11}=A_{22}, A_{12}=0 A11=A22,A12=0,则有新方程
∂ 2 u ∂ ξ 2 + ∂ 2 u ∂ η 2 + 1 A 11 ( B 1 ∂ u ∂ ξ + B 2 ∂ u ∂ η + C u ) = 0 \frac{\partial^2 u}{\partial \xi^2}+\frac{\partial^2u}{\partial\eta^2}+\frac{1}{A_{11}}(B_1\frac{\partial u}{\partial \xi}+B_2\frac{\partial u}{\partial \eta}+Cu)=0 ∂ξ2∂2u+∂η2∂2u+A111(B1∂ξ∂u+B2∂η∂u+Cu)=0
称方程为椭圆型方程的标准型。一维波动方程、热传导方程和二维调和方程正式三类方程的最简单的标准性。
例1:求方程
x 2 ∂ 2 u ∂ x 2 − 2 x y ∂ 2 u ∂ x ∂ y + y 2 ∂ 2 u ∂ y 2 + x ∂ u ∂ x + y ∂ u ∂ y = 0 x^2\frac{\partial^2 u}{\partial x^2}-2xy\frac{\partial^2 u}{\partial x\partial y}+y^2\frac{\partial^2u}{\partial y^2}+x\frac{\partial u}{\partial x}+y\frac{\partial u}{\partial y}=0 x2∂x2∂2u−2xy∂x∂y∂2u+y2∂y2∂2u+x∂x∂u+y∂y∂u=0
的通解。解: Δ = ( x y ) 2 − x 2 y 2 = 0 \Delta=(xy)^2-x^2y^2=0 Δ=(xy)2−x2y2=0,方程为抛物型,特征方程
x 2 ( d y ) 2 + 2 x y d x d y + y 2 ( d x ) 2 = ( x d y + y d x ) 2 = 0 x^2(dy)^2+2xydxdy+y^2(dx)^2=(xdy+ydx)^2=0 x2(dy)2+2xydxdy+y2(dx)2=(xdy+ydx)2=0
即
y d x + x d y = 0 ydx+xdy=0 ydx+xdy=0
积分得一族特征线
x y = c xy=c xy=c
作变量代换
ξ = x y , η = y \xi=xy,\quad \eta=y ξ=xy,η=y
可把方程变为
η u η η + u η = 0 \eta u_{\eta\eta}+u\eta=0 ηuηη+uη=0
令 v = u η v=u_\eta v=uη,得
η v η + v = 0 \eta v_\eta+ v=0 ηvη+v=0
解之得
u η = v = φ ( ξ ) η − 1 u_{\eta}=v=\varphi(\xi)\eta^{-1} uη=v=φ(ξ)η−1
所以
u = φ ( ξ ) l n ∣ η ∣ + ψ ( ξ ) = φ ( x y ) l n ∣ y ∣ + ψ ( x , y ) , y ≠ 0 u=\varphi(\xi)ln|\eta|+\psi(\xi)=\varphi(xy)ln|y|+\psi(x,y),\quad y\neq 0 u=φ(ξ)ln∣η∣+ψ(ξ)=φ(xy)ln∣y∣+ψ(x,y),y=0
其中, φ , ψ \varphi,\psi φ,ψ是任意二次可微函数。 -
用python科学计算解偏微分方程
2020-03-17 11:00:27用python科学计算解偏微分方程,特别是用高斯迭代法计算拉普拉斯方程 -
算子法二阶线性非齐次微分方程的通解
2021-11-18 14:25:24对于一般的二阶线性非齐次微分方程 我们先利用格林函数的思想来求解 固定每一个时刻,假定我们能求解 那么可以证明原方程的解就是 我们引入阶跃函数将上面等式上限变为正无穷 我们将上式改写为卷积...对于一般的二阶线性非齐次微分方程
我们先利用格林函数的思想来求解
固定每一个
时刻,假定我们能求解
那么可以证明原方程的解就是
我们引入阶跃函数将上面等式上限变为正无穷
我们将上式改写为卷积形式
到这里我们已经形式上解决了这个一般的问题。我们只需要求解G(t),通过一个积分计算,即解决了这个一般性问题。当然这样计算未必很好算。
求解G(t)最直接的方法是拉普拉斯变换法,这样即可以解出G(t)。这样对所有f(t), 我们都可以通过相同的方法找到通解,这就是信号与系统中核心思想。
直接运用拉普拉斯变换也可以导出相同的结论,很容易的得到
Y(s)=F(s)(1/as^2 + bs + c)
我们将上式反变换写成卷积的形式,并注意到阶跃函数的存在,可以发现积分范围这里是o到t。
实际上,当存在两个实根时,我们是可以快速写出其一般形式解的。
-
如何用matlab解二维的非线性偏微分方程组,其中每个方程是抛物线型的
2021-04-18 13:58:37二是PDE工具箱,可以求解特殊PDE问题,PDEtool有较大的局限性,比如只能求解二阶PDE问题,并且不能解决偏微分方程组,但是它提供了GUI界面,从繁杂的编程中解脱出来了,同时还可以通过File->Save As直接生成M代... -
高阶完全非线性偏微分方程组的机器学习近似算法(笔记)
2020-11-13 14:56:25将神经网络方法应用于偏微分方程,是近年的一个前沿方向。本文是对论文Machine learning approximation algorithms for high-dimensional fully nonlinear partial differential equations and second-order ... -
同伦摄动法在求解非线性偏微分方程中的应用 (2014年)
2021-04-23 13:42:39文章主要研究了同伦摄动法在求解非线性偏微分方程中的应用问题。简要介绍了同伦摄动法,该法的基本思想是通过行波变换并结合同伦摄动理论,把求解某些非线性偏微分方程的问题转化为求解常微分方程的初值问题,最后得出... -
数学物理方程 第二章 线性偏微分方程的通解
2020-05-12 23:37:36n阶线性偏微分方程的求解对于我目前的问题来说没什么用,我只是解二阶线性偏微分方程。 下面我只讲二阶方程如何解 2.1 一维波动方程的解 2.2常系数线性齐次偏微分方程的通解 Lu=0 (当然2.2.1简单的讨论了非齐... -
微分方程求解二(偏微分方程求解)
2020-11-09 20:36:51文章目录(1)偏微分方程的类型(二阶)(2)抛物线型1.显式法2.Crank-Nicholson隐式算法 (3)双曲线型(4)椭圆型 (1)偏微分方程的类型(二阶) a∂2u∂x2+b∂2u∂y∂x+c∂2u∂x2+d∂u∂x+e∂u∂y+fu+g=0a\frac{\partial^2u}{\... -
一阶线性偏微分方程通解法和特征线法(二)| n个自变量情况 | 偏微分方程(八)
2020-05-03 01:26:14n个自变量的一阶线性偏微分方程(n≥2n\geq2n≥2) n个自变量的一阶线性偏微分方程的一般形式为 ∑j=1nbj∂u∂xj+cu=f(7) \sum_{j=1}^nb_j\frac{\partial u}{\partial x_j}+cu=f \tag{7} j=1∑nbj∂xj∂u+cu... -
线性偏微分方程的通解
2022-01-02 22:40:54之前我们介绍的如波动方程、扩散方程等二阶偏微分方程都是线性偏微分方程,也就是说,在方程中只出现对于末知函数的线性运算。引进线性算符 L^\hat{L}L^ 的记号,则线性偏微分方程都可以写成 L^[u]=f \hat{L}[u]=f L... -
应用Matlab求解二阶偏微分方程的基…
2017-03-14 10:50:01>利用Matlab求解二阶偏微分方程的一般有以下步骤 → 题目定义:由方程(3.4.33)和(3.4.35)可以看出,参量 是二阶偏微 分方程的主要参量,只要这几个参量确定,就可以定下偏微分方程的结构。此外要做的事是确定偏... -
解非线性双曲型偏微分方程的预测-校正法 (1994年)
2021-05-27 12:53:50把求解二阶常微分方程的预测-校正法推广到解非线性双曲型偏微分方程,并给出几种预测-校正格式.并用数值例子证明这些格式是有效的. -
二阶椭圆偏微分方程的Pinching一估计 (2007年)
2021-05-15 12:04:27Pinching一估计是研究解的凸性的一种重要方法,主要给出了半线性二阶椭圆偏微分方程的Pinching-估计,并将其推广到一类完全非线性二阶椭圆偏微分方程。 -
非线性微分方程的线性化
2020-04-14 20:38:45非线性微分方程的求解和控制系统性能研究非常复杂,而线 性化后的模型可借助叠加原理的性质,简化系统分析。 因此,研究非线性微分方程的线性化具有较强的工程实用价 值。 什么是非线性数学模型的线性化? 在一定的... -
MATLAB解含参数方程、矩阵方程、二阶微分方程组
2021-04-18 03:18:361、如下公式: 其中x、y、m、n为参数,a、b为未知数,利用MATLAB求解方程; syms x y m n a b;[a,b]=solve('x=m*cos(a)+n*cos(a+b)','y=m*sin(a)+n*sin(a+b)','a','b');会得到四组解;对a(3)、b(3)、a(4)、b(4)进行... -
非线性微分方程有限差分解法
2020-11-23 14:22:55一般的, 对于线性偏微分方程其有限差分解法比较容易编写. 对于非线性微分方程的数值求解而言, 一般的, 如果采用有限差分的方法进行离散, 那么离散后方程中就会出现非线性项. 此时可以将非线性项看作是新的变量, 利用... -
python:利用多种方式解微分方程(以二阶微分系统零状态响应为例)
2020-12-24 14:04:32引入首先用高数知识求解非齐次常系数微分方程image.png再利用信号与系统中冲激响应求解验证image.png利用MATLAB求解验证y=dsolve('D2y+3*Dy+2*y=exp(-t)','y(0)=1','Dy(0)=2','t')得出结果:y =(t - 2 exp(-t) + 3) ... -
一阶线性偏微分方程特征线解法.pdf
2020-12-19 22:11:16• 偏微分方程理论 主要研究具有实际背景的偏微分方程或偏微分方程组。是数学的基础学科之一。第 1页二、《数学物理方程》课程的特点:1、数学理论、解题方法与物理实际有机结合。可以学到:如何根据物理现象建立偏... -
[数理] 非线性偏微分复方程 PDF
2009-05-24 14:39:031 带弱条件的非线性椭圆型方程与方程组 2 高阶椭圆型方程组 3 可测系数的二阶非线性非散度型抛物型方程 4 可测系数的二阶非线性非散度型抛物型方程组 5 一阶与二阶双曲型复方程 6 一阶与二阶混合型复方程 -
科学网—求解偏微分方程开源有限元软件deal.II学习--Step 3 - 亓欣波的博文
2021-05-23 04:38:47这样显著地跟程序说明需要计算什么,就能加速运算,因为有的软件所有东西不管有用没用都一块计算,比如二阶导数、法向量等。注意最后用了按位或这一运算符,这在c语言中很常见,这里的意思就是我要计算谁“和”谁。1... -
偏微分方程的物理含义
2022-03-22 10:47:59偏微分方程的物理含义1.实际问题采用偏微分方程的原因2.偏微分方程的物理含义(1) 双曲型方程。其代表是 波动方程(2)抛物型方程。其代表是 热传导方程(3) 椭圆型方程。其代表是 拉普拉斯方程(4)薛定锷方程成为量子... -
常微分方程的解析解(方法归纳)以及基于Python的二阶微分方程边值问题的数值算例实现
2021-10-05 13:17:18常微分方程的解析解(方法归纳)...一阶齐次(非齐次)线性微分方程(一阶) 二阶常系数微分方程(二阶) 高阶常系数微分方程(nnn阶) 可分离变量的微分方程(一阶) 这类微分方程可以变形成如下形式: f(x)dx=g(y)dy -
Python小白的数学建模课-11.偏微分方程数值解法
2021-08-17 14:10:06偏微分方程可以描述各种自然和工程现象, 是构建科学、工程学和其他领域的数学模型主要手段。 本文采用有限差分法求解偏微分方程,通过案例讲解一维平流方程、一维热传导方程、二维双曲方程、二维抛物方程和二维椭圆...