能量不等式证明

news/2024/4/19 13:06:57/文章来源:https://blog.csdn.net/coldasice342/article/details/136476029

波动方程初值问题能量不等式的证明

Gronwall 不等式

若非负函数 G ( τ ) G(\tau) G(τ) [ 0 , T ] [0,T] [0,T] 上连续可微, G ( 0 ) = 0 G(0)=0 G(0)=0,且对 τ ∈ [ 0 , T ] \tau\in[0,T] τ[0,T]满足 d G ( τ ) d τ ≤ C G ( τ ) + F ( τ ) \frac{dG(\tau)}{d\tau}\leq CG(\tau)+F(\tau) dτdG(τ)CG(τ)+F(τ) 其中 C C C 为常数且 C > 0 C>0 C>0 F ( τ ) F(\tau) F(τ) [ 0 , T ] [0,T] [0,T] 上不减的非负可积函数,

那么有:

d G ( τ ) d τ ≤ e C τ F ( τ ) \frac{dG(\tau)}{d\tau}\leq e^{C\tau}F(\tau) dτdG(τ)eCτF(τ) G ( τ ) ≤ C − 1 ( e C τ − 1 ) F ( τ ) G(\tau)\leq C^{-1}(e^{C\tau}-1)F(\tau) G(τ)C1(eCτ1)F(τ)

常用的替换技巧

u t u t t = 1 2 ⋅ 2 u t ⋅ u t t = 1 2 ( u t 2 ) t u_{t}u_{tt}=\frac{1}{2}\cdot 2u_{t}\cdot u_{tt}=\frac{1}{2}(u_{t}^2)_t ututt=212ututt=21(ut2)t u u t = 1 2 ⋅ 2 u ⋅ u t = 1 2 ( u 2 ) t uu_t=\frac{1}{2}\cdot 2u\cdot u_t =\frac{1}{2}(u^2)_t uut=212uut=21(u2)t u u x x = ( u u x ) x − u x 2 = u x 2 + u u x x − u x 2 uu_{xx}=(uu_x)_x-u_x^2=u_x^2+uu_{xx}-u_x^2 uuxx=(uux)xux2=ux2+uuxxux2 ∵ ( u t u x ) x = u t x u x + u t u x x \because(u_{t}u_{x})_x=u_{tx}u_{x}+u_{t}u_{xx} (utux)x=utxux+utuxx ∴ u t u x x = ( u t u x ) x − u t x u x \therefore u_{t}u_{xx}=(u_{t}u_{x})_x-u_{tx}u_{x} utuxx=(utux)xutxux

u u u 二阶导函数连续时, u x t = u t x u_{xt}=u_{tx} uxt=utx ∴ u t x u x = 1 2 ( u x 2 ) t = 1 2 ⋅ 2 u x ⋅ u x t = u x u t x \therefore u_{tx}u_{x} = \frac{1}{2}(u_{x}^2)_t=\frac{1}{2}\cdot 2u_{x}\cdot u_{xt}=u_{x}u_{tx} utxux=21(ux2)t=212uxuxt=uxutx ∴ u t u x x = ( u t u x ) x − 1 2 ( u x 2 ) t \therefore u_{t}u_{xx}=(u_{t}u_{x})_x-\frac{1}{2}(u_{x}^2)_t utuxx=(utux)x21(ux2)t

能量不等式证明过程

规定 x x x 增长的方向为正方向,而 G r e e n Green Green 公式曲线积分时,组成梯形区域 K τ K_{\tau} Kτ 边界 ∂ K τ \partial K_{\tau} Kτ 的四条线段按照右手法则,只有 Ω 0 \Omega_0 Ω0 这条线段的方向是与右手法则一致,所以有:

∂ K τ = Ω 0 ∪ ( − Γ τ 2 ) ∪ ( − Ω τ ) ∪ ( − Γ τ 1 ) \partial K_{\tau} =\Omega_0\cup(-\Gamma_{\tau_{2}})\cup(-\Omega_{\tau})\cup(-\Gamma_{\tau_{1}}) Kτ=Ω0(Γτ2)(Ωτ)(Γτ1)
在这里插入图片描述
此外,在 ∂ K τ \partial K_{\tau} Kτ 计算曲线积分时,关于 x , t x,t x,t 都积分,但是这两条边 Ω 0 , Ω τ \Omega_0,\Omega_{\tau} Ω0,Ωτ 只在 x x x 方向有增长,而在 t t t 方向无增长,那么这两条线段在计算曲线积分时关于 d t dt dt 的积分项为 0 0 0

待证明的两个不等式:
∫ Ω τ [ u t 2 ( x , τ ) + a 2 u x 2 ( x , τ ) ] d x ≤ M [ ∫ Ω 0 ( ψ 2 + a 2 φ x 2 ) d x + ∫ K τ f 2 ( x , t ) d x d t ] \int_{\Omega_\tau} \left[u_t^2(x, \tau) + a^2 u_x^2(x, \tau)\right] dx \leq M \left[ \int_{\Omega_0} (\psi^2 + a^2 \varphi_x^2) dx + \int_{K_\tau} f^2(x,t) dxdt \right] Ωτ[ut2(x,τ)+a2ux2(x,τ)]dxM[Ω0(ψ2+a2φx2)dx+Kτf2(x,t)dxdt] ∫ K τ [ u t 2 ( x , t ) + a 2 u x 2 ( x , t ) ] d x d t ≤ M [ ∫ Ω 0 ( ψ 2 + a 2 φ x 2 ) d x + ∫ K τ f 2 ( x , t ) d x d t ] \int_{K_\tau} \left[u_t^2(x,t) + a^2 u_x^2(x,t)\right] dxdt \leq M \left[ \int_{\Omega_0} (\psi^2 + a^2 \varphi_x^2) dx + \int_{K_\tau} f^2(x,t) dxdt \right] Kτ[ut2(x,t)+a2ux2(x,t)]dxdtM[Ω0(ψ2+a2φx2)dx+Kτf2(x,t)dxdt]

  1. 在波动方程 ∂ 2 u ∂ t 2 − a 2 ∂ 2 u ∂ x 2 = f \frac{\partial^2 u}{\partial t^2}-a^2\frac{\partial^2 u}{\partial x^2}=f t22ua2x22u=f 两端同乘以 ∂ u ∂ t \frac{\partial u}{\partial t} tu 并在区域 K τ K_\tau Kτ 上积分并按照前面的常用替换得:

∬ K τ u t u t t − a 2 u t u x x d x d t = ∬ K τ u f d x d t \iint_{K_\tau}u_{t}u_{tt}-a^2u_{t}u_{xx}dxdt=\iint_{K_\tau}ufdxdt Kτututta2utuxxdxdt=Kτufdxdt$

∬ K τ 1 2 ( u t 2 + a 2 u x 2 ) t − a 2 ( u t u x ) x d x d t = ∬ K τ u f d x d t \iint_{K_\tau}\frac{1}{2}(u_{t}^2+a^2u_{x}^2)_t-a^2(u_tu_{x})_{x}dxdt=\iint_{K_\tau}ufdxdt Kτ21(ut2+a2ux2)ta2(utux)xdxdt=Kτufdxdt

  1. 应用下面这个 G r e e n Green Green 公式将上面等式左边替换
    ∬ Ω ∂ P ∂ t + ∂ Q ∂ x d σ = ∮ ∂ Ω − P d x + Q d t \iint_{\Omega}\frac{\partial P}{\partial t}+\frac{\partial Q}{\partial x}d\sigma=\oint_{\partial\Omega}-Pdx+Qdt ΩtP+xQdσ=ΩPdx+Qdt其中 P = 1 2 ( u t 2 + a 2 u x 2 ) P=\frac{1}{2}(u_{t}^2+a^2u_{x}^2) P=21(ut2+a2ux2) Q = − a 2 ( u t u x ) Q=-a^2(u_tu_{x}) Q=a2(utux),得闭合曲线积分

∮ ∂ K τ [ − 1 2 ( u t 2 + a 2 u x 2 ) d x − a 2 ( u t u x ) d t ] \oint_{\partial{K_{\tau}}}[-\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx-a^2(u_tu_{x})dt] Kτ[21(ut2+a2ux2)dxa2(utux)dt]

= − ∮ ∂ K τ a 2 ( u t u x ) d t + 1 2 ( u t 2 + a 2 u x 2 ) d x =-\oint_{\partial{K_{\tau}}}a^2(u_tu_{x})dt+\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx =Kτa2(utux)dt+21(ut2+a2ux2)dx

a 2 ( u t u x ) d t + 1 2 ( u t 2 + a 2 u x 2 ) d x = □ a^2(u_tu_{x})dt+\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx = \square a2(utux)dt+21(ut2+a2ux2)dx=

∵ ∂ K τ = Ω 0 ∪ ( − Γ τ 2 ) ∪ ( − Ω τ ) ∪ ( − Γ τ 1 ) \because \partial K_{\tau} =\Omega_0\cup(-\Gamma_{\tau_{2}})\cup(-\Omega_{\tau})\cup(-\Gamma_{\tau_{1}}) Kτ=Ω0(Γτ2)(Ωτ)(Γτ1)

∴ − ∮ ∂ K τ a 2 ( u t u x ) d t + 1 2 ( u t 2 + a 2 u x 2 ) d x \therefore -\oint_{\partial{K_{\tau}}}a^2(u_tu_{x})dt+\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx Kτa2(utux)dt+21(ut2+a2ux2)dx

= − [ ∫ Ω 0 □ − ∫ Ω τ □ − ∫ Γ τ 1 □ − ∫ Γ τ 2 □ ] =-[\int_{\Omega_0}\square-\int_{\Omega_{\tau}}\square-\int_{\Gamma_{\tau_{1}}}\square-\int_{\Gamma_{\tau_{2}}}\square] =[Ω0ΩτΓτ1Γτ2]

Ω 0 , Ω τ \Omega_0,\Omega_{\tau} Ω0,Ωτ 只在 x x x 方向有增长,而在 t t t 方向无增长,那么这两条线段在计算曲线积分时关于 d t dt dt 的积分项为 0 0 0,可得

− ∮ ∂ K τ a 2 ( u t u x ) d t + 1 2 ( u t 2 + a 2 u x 2 ) d x -\oint_{\partial{K_{\tau}}}a^2(u_tu_{x})dt+\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx Kτa2(utux)dt+21(ut2+a2ux2)dx

= ∫ Γ τ 1 ∪ Γ τ 2 a 2 ( u t u x ) d t + 1 2 ( u t 2 + a 2 u x 2 ) d x − ∫ Ω 0 1 2 ( u t 2 + a 2 u x 2 ) d x + ∫ Ω τ 1 2 ( u t 2 + a 2 u x 2 ) d x =\int_{\Gamma_{\tau_{1}}\cup\Gamma_{\tau_{2}}}a^2(u_tu_{x})dt+\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx-\int_{\Omega_0}\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx+\int_{\Omega_\tau}\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx =Γτ1Γτ2a2(utux)dt+21(ut2+a2ux2)dxΩ021(ut2+a2ux2)dx+Ωτ21(ut2+a2ux2)dx

记上面等式的三项分别为 J 1 , J 2 , J 3 J_1,J_2,J_3 J1,J2,J3,即:

J 1 = ∫ Γ τ 1 ∪ Γ τ 2 a 2 ( u t u x ) d t + 1 2 ( u t 2 + a 2 u x 2 ) d x J_1=\int_{\Gamma_{\tau_{1}}\cup\Gamma_{\tau_{2}}}a^2(u_tu_{x})dt+\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx J1=Γτ1Γτ2a2(utux)dt+21(ut2+a2ux2)dx

J 2 = − ∫ Ω 0 1 2 ( u t 2 + a 2 u x 2 ) d x J_2=-\int_{\Omega_0}\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx J2=Ω021(ut2+a2ux2)dx

J 3 = ∫ Ω τ 1 2 ( u t 2 + a 2 u x 2 ) d x J_3=\int_{\Omega_\tau}\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx J3=Ωτ21(ut2+a2ux2)dx

  1. 应用下面 d x dx dx d t dt dt 的关系,将 J 1 J_1 J1 统一为只关于 d t dt dt 的积分
    Γ τ 1 : d x = a d t \Gamma_{\tau_{1}}:dx=adt Γτ1:dx=adt Γ τ 2 : d x = − a d t \Gamma_{\tau_{2}}:dx=-adt Γτ2:dx=adt

J 1 = ∫ Γ τ 1 a 2 ( u t u x ) d t + 1 2 a ( u t 2 + a 2 u x 2 ) d t + ∫ Γ τ 2 a 2 ( u t u x ) d t − 1 2 a ( u t 2 + a 2 u x 2 ) d t J_1=\int_{\Gamma_{\tau_{1}}}a^2(u_tu_{x})dt+\frac{1}{2}a(u_{t}^2+a^2u_{x}^2)dt+\int_{\Gamma_{\tau_{2}}}a^2(u_tu_{x})dt-\frac{1}{2}a(u_{t}^2+a^2u_{x}^2)dt J1=Γτ1a2(utux)dt+21a(ut2+a2ux2)dt+Γτ2a2(utux)dt21a(ut2+a2ux2)dt
= ∫ Γ τ 1 a 2 ( 2 a u t u x + u t 2 + a 2 u x 2 ) d t + ∫ Γ τ 1 a 2 ( 2 a u t u x − u t 2 − a 2 u x 2 ) d t =\int_{\Gamma_{\tau_{1}}}\frac{a}{2}(2au_tu_{x}+u_{t}^2+a^2u_{x}^2)dt+\int_{\Gamma_{\tau_{1}}}\frac{a}{2}(2au_tu_{x}-u_{t}^2-a^2u_{x}^2)dt =Γτ12a(2autux+ut2+a2ux2)dt+Γτ12a(2autuxut2a2ux2)dt

利用完全平方公式得

J 1 = ∫ Γ τ 1 a 2 ( u t + a u x ) 2 d t − ∫ Γ τ 2 a 2 ( u t − a u x ) 2 d t J_1=\int_{\Gamma_{\tau_{1}}}\frac{a}{2}(u_t+au_x)^2dt-\int_{\Gamma_{\tau_{2}}}\frac{a}{2}(u_t-au_x)^2dt J1=Γτ12a(ut+aux)2dtΓτ22a(utaux)2dt

因为在 Γ τ 1 \Gamma_{\tau_{1}} Γτ1 线段上规定的正方向是 t t t 增长的方向,而 Γ τ 2 \Gamma_{\tau_{2}} Γτ2 上相反,所以 t t t 的积分限分别是 0 0 0 τ \tau τ τ \tau τ 0 0 0 ,也即

J 1 = ∫ 0 τ a 2 ( u t + a u x ) 2 d t − ∫ τ 0 a 2 ( u t − a u x ) 2 d t ≥ 0 J_1=\int_{0}^{\tau}\frac{a}{2}(u_t+au_x)^2dt-\int_{\tau}^0\frac{a}{2}(u_t-au_x)^2dt\geq0 J1=0τ2a(ut+aux)2dtτ02a(utaux)2dt0

第一项平方项且积分下限小于上限故其积分非负,而第二项是平方项在积分下限大于上限加上前面的负号,故也是非负,所以 J 1 ≥ 0 J_1\geq0 J10

  1. J 1 , J 2 , J 3 J_1,J_2,J_3 J1,J2,J3 写回来,根据 J 1 ≥ 0 J_1\geq0 J10 构造不等关系
    J 1 + J 2 + J 3 = ∬ K τ u t f d x d t J_1+J_2+J_3=\iint_{K_{\tau}}u_{t}fdxdt J1+J2+J3=Kτutfdxdt
    ∵ J 1 ≥ 0 \because J_1\geq0 J10
    ∴ J 3 ≤ J 1 + J 3 = J 1 + J 2 + J 3 − J 2 = ∬ K τ u t f d x d t − J 2 \therefore J_3\leq J_1+J_3=J_1+J_2+J_3-J_2=\iint_{K_{\tau}}u_{t}fdxdt-J_2 J3J1+J3=J1+J2+J3J2=KτutfdxdtJ2
    也即 ∫ Ω τ 1 2 ( u t 2 + a 2 u x 2 ) d x ≤ ∬ K τ u t f d x d t + ∫ Ω 0 1 2 ( u t 2 + a 2 u x 2 ) d x \int_{\Omega_\tau}\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx\leq \iint_{K_{\tau}}u_{t}fdxdt+\int_{\Omega_0}\frac{1}{2}(u_{t}^2+a^2u_{x}^2)dx Ωτ21(ut2+a2ux2)dxKτutfdxdt+Ω021(ut2+a2ux2)dx
    两边同乘 2 2 2 Ω 0 \Omega_0 Ω0 t = 0 t=0 t=0 的线段,而 t = 0 t=0 t=0 时, u t ( x , 0 ) = ψ , u x ( x , 0 ) = φ x u_t(x,0)=\psi,u_x(x,0)=\varphi_{x} ut(x,0)=ψ,ux(x,0)=φx
    ∫ Ω τ ( u t 2 + a 2 u x 2 ) d x ≤ ∬ K τ 2 u t f d x d t + ∫ Ω 0 ( ψ 2 + a 2 φ x 2 ) d x \int_{\Omega_\tau}(u_{t}^2+a^2u_{x}^2)dx\leq \iint_{K_{\tau}}2u_{t}fdxdt+\int_{\Omega_0}(\psi^2+a^2\varphi_{x}^2)dx Ωτ(ut2+a2ux2)dxKτ2utfdxdt+Ω0(ψ2+a2φx2)dx

对上面不等式右边 ∬ K τ 2 u t f d x d t \iint_{K_{\tau}}2u_{t}fdxdt Kτ2utfdxdt 这一项应用 C a u c h y Cauchy Cauchy 不等式 2 a b ≤ a 2 + b 2 2ab\leq a^2+b^2 2aba2+b2
∫ Ω τ ( u t 2 + a 2 u x 2 ) d x ≤ ∬ K τ u t 2 + f 2 d x d t + ∫ Ω 0 ( ψ 2 + a 2 φ x 2 ) d x \int_{\Omega_\tau}(u_{t}^2+a^2u_{x}^2)dx\leq \iint_{K_{\tau}}u_{t}^2+f^2dxdt+\int_{\Omega_0}(\psi^2+a^2\varphi_{x}^2)dx Ωτ(ut2+a2ux2)dxKτut2+f2dxdt+Ω0(ψ2+a2φx2)dx

此时,对比要证明的不等式
∫ Ω τ ( u t 2 + a 2 u x 2 ) d x ≤ ∬ K τ f 2 d x d t + ∫ Ω 0 ( ψ 2 + a 2 φ x 2 ) d x \int_{\Omega_\tau}(u_{t}^2+a^2u_{x}^2)dx\leq \iint_{K_{\tau}}f^2dxdt+\int_{\Omega_0}(\psi^2+a^2\varphi_{x}^2)dx Ωτ(ut2+a2ux2)dxKτf2dxdt+Ω0(ψ2+a2φx2)dx
发现只是右端多了一项 ∬ K τ u t 2 d x d t \iint_{K_{\tau}}u_{t}^2dxdt Kτut2dxdt

  1. 利用 G r o n w a l l Gronwall Gronwall 不等式将该项消去

G ( τ ) = ∬ K τ ( u t 2 + a 2 u x 2 ) d x d t = ∫ 0 τ ∫ x 0 − a ( t 0 − t ) x 0 + a ( t 0 − t ) ( u t 2 + a 2 u x 2 ) d x d t G(\tau)=\iint_{K_{\tau}}(u_{t}^2+a^2u_{x}^2)dxdt=\int_{0}^{\tau}\int_{x_0-a(t_0-t)}^{x_0+a(t_0-t)}(u_{t}^2+a^2u_{x}^2)dxdt G(τ)=Kτ(ut2+a2ux2)dxdt=0τx0a(t0t)x0+a(t0t)(ut2+a2ux2)dxdt

d G ( τ ) d τ = ∫ Ω τ ( u t 2 + a 2 u x 2 ) d x \frac{dG(\tau)}{d\tau}=\int_{\Omega_\tau}(u_{t}^2+a^2u_{x}^2)dx dτdG(τ)=Ωτ(ut2+a2ux2)dx 是待证明不等式的左端

F ( τ ) = ∬ K τ f 2 d x d t + ∫ Ω 0 ( ψ 2 + a 2 φ x 2 ) d x F(\tau)=\iint_{K_{\tau}}f^2dxdt+\int_{\Omega_0}(\psi^2+a^2\varphi_{x}^2)dx F(τ)=Kτf2dxdt+Ω0(ψ2+a2φx2)dx

∫ Ω τ ( u t 2 + a 2 u x 2 ) d x ≤ ∬ K τ u t 2 + f 2 d x d t + ∫ Ω 0 ( ψ 2 + a 2 φ x 2 ) d x \int_{\Omega_\tau}(u_{t}^2+a^2u_{x}^2)dx\leq \iint_{K_{\tau}}u_{t}^2+f^2dxdt+\int_{\Omega_0}(\psi^2+a^2\varphi_{x}^2)dx Ωτ(ut2+a2ux2)dxKτut2+f2dxdt+Ω0(ψ2+a2φx2)dx 右边添加 ∬ K τ a 2 u x 2 d x d t \iint_{K_{\tau}}a^2u_{x}^2dxdt Kτa2ux2dxdt,则
∫ Ω τ ( u t 2 + a 2 u x 2 ) d x ≤ ∬ K τ u t 2 + a 2 u x 2 + f 2 d x d t + ∫ Ω 0 ( ψ 2 + a 2 φ x 2 ) d x \int_{\Omega_\tau}(u_{t}^2+a^2u_{x}^2)dx\leq \iint_{K_{\tau}}u_{t}^2+a^2u_{x}^2+f^2dxdt+\int_{\Omega_0}(\psi^2+a^2\varphi_{x}^2)dx Ωτ(ut2+a2ux2)dxKτut2+a2ux2+f2dxdt+Ω0(ψ2+a2φx2)dx

那么此时该不等式满足 G r o n w a l l Gronwall Gronwall 不等式的前提条件
d G ( τ ) d τ ≤ C G ( τ ) + F ( τ ) \frac{dG(\tau)}{d\tau}\leq CG(\tau)+F(\tau) dτdG(τ)CG(τ)+F(τ)
此时 C = 1 C=1 C=1,那么由 G r o n w a l l Gronwall Gronwall 不等式的两个结论,可分别得到如下两个能量不等式,其中 M = e τ M=e^\tau M=eτ
∫ Ω τ [ u t 2 ( x , τ ) + a 2 u x 2 ( x , τ ) ] d x ≤ M [ ∫ Ω 0 ( ψ 2 + a 2 φ x 2 ) d x + ∫ K τ f 2 ( x , t ) d x d t ] \int_{\Omega_\tau} \left[u_t^2(x, \tau) + a^2 u_x^2(x, \tau)\right] dx \leq M \left[ \int_{\Omega_0} (\psi^2 + a^2 \varphi_x^2) dx + \int_{K_\tau} f^2(x,t) dxdt \right] Ωτ[ut2(x,τ)+a2ux2(x,τ)]dxM[Ω0(ψ2+a2φx2)dx+Kτf2(x,t)dxdt]
∫ K τ [ u t 2 ( x , t ) + a 2 u x 2 ( x , t ) ] d x d t ≤ M [ ∫ Ω 0 ( ψ 2 + a 2 φ x 2 ) d x + ∫ K τ f 2 ( x , t ) d x d t ] \int_{K_\tau} \left[u_t^2(x,t) + a^2 u_x^2(x,t)\right] dxdt \leq M \left[ \int_{\Omega_0} (\psi^2 + a^2 \varphi_x^2) dx + \int_{K_\tau} f^2(x,t) dxdt \right] Kτ[ut2(x,t)+a2ux2(x,t)]dxdtM[Ω0(ψ2+a2φx2)dx+Kτf2(x,t)dxdt]

波动方程混合问题能量不等式的证明

待证明的两个不等式:
∫ 0 l [ u t 2 ( x , τ ) + a 2 u x 2 ( x , τ ) ] d x ≤ M [ ∫ 0 l ( ψ 2 + a 2 φ x 2 ) d x + ∫ Q τ f 2 d x d t ] \int_{0}^{l} \left[u_t^2(x, \tau) + a^2 u_x^2(x, \tau)\right] dx \leq M \left[ \int_{0}^{l} (\psi^2 + a^2 \varphi_x^2) dx + \int_{Q_\tau} f^2 dxdt \right] 0l[ut2(x,τ)+a2ux2(x,τ)]dxM[0l(ψ2+a2φx2)dx+Qτf2dxdt] ∬ Q τ ( u t 2 + a 2 u x 2 ) d x d t ≤ M [ ∫ 0 l ( ψ 2 + a 2 φ x 2 ) d x + ∫ Q τ f 2 d x d t ] \iint_{Q_\tau} \left(u_t^2 + a^2 u_x^2\right) dxdt \leq M \left[ \int_{0}^{l} (\psi^2 + a^2 \varphi_x^2) dx + \int_{Q_\tau} f^2 dxdt \right] Qτ(ut2+a2ux2)dxdtM[0l(ψ2+a2φx2)dx+Qτf2dxdt]

在这里插入图片描述
二维区域的情况下,这里可以不用 G r e e n Green Green 公式将面积分转为线积分,可以直接按照二重积分展开。

  1. 在波动方程 u t t − a 2 u x x = f u_{tt}-a^2u_{xx}=f utta2uxx=f 两端同乘 u t u_t ut 并在区域 Q τ = ( 0 , l ) × ( 0 , τ ) Q_\tau=(0,l)\times(0,\tau) Qτ=(0,l)×(0,τ) 上积分

∬ Q τ u t u t t − a 2 u t u x x d x d t = ∬ Q τ u t f d x d t \iint_{Q_\tau}u_{t}u_{tt}-a^2u_{t}u_{xx}dxdt=\iint_{Q_\tau}u_tf dxdt Qτututta2utuxxdxdt=Qτutfdxdt

利用常用的替换技巧
u t u t t = 1 2 ⋅ 2 u t ⋅ u t t = 1 2 ( u t 2 ) t u_{t}u_{tt}=\frac{1}{2}\cdot 2u_{t}\cdot u_{tt}=\frac{1}{2}(u_{t}^2)_t ututt=212ututt=21(ut2)t u t u x x = ( u t u x ) x − 1 2 ( u x 2 ) t u_{t}u_{xx}=(u_{t}u_{x})_x-\frac{1}{2}(u_{x}^2)_t utuxx=(utux)x21(ux2)t ∬ Q τ 1 2 ( u t 2 ) t + 1 2 a 2 ( u x 2 ) t − ( a 2 u t u x ) x d x d t \iint_{Q_\tau}\frac{1}{2}(u_{t}^2)_t+\frac{1}{2}a^2(u_{x}^2)_t-(a^2u_tu_x)_xdxdt Qτ21(ut2)t+21a2(ux2)t(a2utux)xdxdt,这里用二重积分展开

∬ Q τ 1 2 ( u t 2 ) t + 1 2 a 2 ( u x 2 ) t − ( a 2 u t u x ) x d x d t \iint_{Q_\tau}\frac{1}{2}(u_{t}^2)_t+\frac{1}{2}a^2(u_{x}^2)_t-(a^2u_tu_x)_xdxdt Qτ21(ut2)t+21a2(ux2)t(a2utux)xdxdt

= ∫ 0 τ ∫ 0 l [ ( 1 2 u t 2 + 1 2 a 2 u x 2 ) t − ( a 2 u t u x ) x ] d x d t =\int_{0}^{\tau}\int_{0}^{l}[(\frac{1}{2}u_{t}^2+\frac{1}{2}a^2u_{x}^2)_t-(a^2u_tu_x)_x]dxdt =0τ0l[(21ut2+21a2ux2)t(a2utux)x]dxdt

= ∫ 0 l 1 2 u t 2 ( x , τ ) + 1 2 a 2 u x 2 ( x , τ ) d x − ∫ 0 l 1 2 u t 2 ( x , 0 ) + 1 2 a 2 u x 2 ( x , 0 ) d x − ∫ 0 τ a 2 u t ( l , t ) u x ( l , t ) d t + ∫ 0 τ a 2 u t ( 0 , t ) u x ( 0 , t ) d t =\int_{0}^{l}\frac{1}{2}u_{t}^{2}(x,\tau)+\frac{1}{2}a^2u_{x}^{2}(x,\tau)dx-\int_{0}^{l}\frac{1}{2}u_{t}^{2}(x,0)+\frac{1}{2}a^2u_{x}^{2}(x,0)dx-\int_{0}^{\tau}a^2u_{t}(l,t)u_{x}(l,t)dt+\int_{0}^{\tau}a^2u_{t}(0,t)u_{x}(0,t)dt =0l21ut2(x,τ)+21a2ux2(x,τ)dx0l21ut2(x,0)+21a2ux2(x,0)dx0τa2ut(l,t)ux(l,t)dt+0τa2ut(0,t)ux(0,t)dt

  1. 根据初始条件和齐次边界条件替换上面表达式

初始条件:
u t ( x , 0 ) = ψ ( x ) u_t(x,0)=\psi(x) ut(x,0)=ψ(x) u x ( x , 0 ) = φ x ( x ) u_x(x,0)=\varphi_x(x) ux(x,0)=φx(x)
边界条件:
u ( 0 , t ) = 0 u(0,t)=0 u(0,t)=0 u ( l , t ) = 0 u(l,t)=0 u(l,t)=0即杆的两端是固定的,也即速度为 0 0 0
u t ( l , t ) = u t ( 0 , t ) = 0 u_t(l,t)=u_t(0,t)=0 ut(l,t)=ut(0,t)=0 带入 1. 1. 1. 得到的等式并两边同乘 2 2 2

= ∫ 0 l u t 2 ( x , τ ) + a 2 u x 2 ( x , τ ) d x = ∫ 0 l ψ 2 + a 2 φ x 2 d x + 2 ∫ 0 τ ∫ 0 l u t f d x d t =\int_{0}^{l}u_{t}^{2}(x,\tau)+a^2u_{x}^{2}(x,\tau)dx=\int_{0}^{l}\psi^2+a^2\varphi_{x}^{2}dx+2\int_{0}^{\tau}\int_{0}^{l}u_tfdxdt =0lut2(x,τ)+a2ux2(x,τ)dx=0lψ2+a2φx2dx+20τ0lutfdxdt
对等式右边应用 C a u c h y Cauchy Cauchy 不等式以及添加非负项 ∫ 0 τ ∫ 0 l a 2 u x 2 d x d t \int_{0}^{\tau}\int_{0}^{l}a^2u_{x}^2dxdt 0τ0la2ux2dxdt

∫ 0 l u t 2 ( x , τ ) + a 2 u x 2 ( x , τ ) d x ≤ ∫ 0 l ψ 2 + a 2 φ x 2 d x + ∫ 0 τ ∫ 0 l f 2 d x d t + ∫ 0 τ ∫ 0 l u t 2 + a 2 u x 2 d x d t \int_{0}^{l}u_{t}^{2}(x,\tau)+a^2u_{x}^{2}(x,\tau)dx\leq \int_{0}^{l}\psi^2+a^2\varphi_{x}^{2}dx+\int_{0}^{\tau}\int_{0}^{l}f^2dxdt+\int_{0}^{\tau}\int_{0}^{l}u_{t}^2+a^2u_{x}^2dxdt 0lut2(x,τ)+a2ux2(x,τ)dx0lψ2+a2φx2dx+0τ0lf2dxdt+0τ0lut2+a2ux2dxdt

  1. 应用 G r o n w a l l Gronwall Gronwall 不等式

G ( τ ) = ∫ 0 τ ∫ 0 l u t 2 + a 2 u x 2 d x d t G(\tau) = \int_{0}^{\tau}\int_{0}^{l}u_{t}^2+a^2u_{x}^2dxdt G(τ)=0τ0lut2+a2ux2dxdt

F ( τ ) = ∫ 0 l ψ 2 + a 2 φ x 2 d x + ∫ 0 τ ∫ 0 l f 2 d x d t F(\tau)=\int_{0}^{l}\psi^2+a^2\varphi_{x}^{2}dx+\int_{0}^{\tau}\int_{0}^{l}f^2dxdt F(τ)=0lψ2+a2φx2dx+0τ0lf2dxdt

那么 d G ( τ ) d τ = ∫ 0 l u t 2 ( x , τ ) + a 2 u x 2 ( x , τ ) d x \frac{dG(\tau)}{d\tau}=\int_{0}^{l}u_{t}^{2}(x,\tau)+a^2u_{x}^{2}(x,\tau)dx dτdG(τ)=0lut2(x,τ)+a2ux2(x,τ)dx

G ( τ ) , F ( τ ) , d G ( τ ) d τ G(\tau),F(\tau),\frac{dG(\tau)}{d\tau} G(τ),F(τ),dτdG(τ) 带入 G r o n w a l l Gronwall Gronwall 不等式的两个结论即得到待证不等式,其中 C = 1 C=1 C=1

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.luyixian.cn/news_show_997912.aspx

如若内容造成侵权/违法违规/事实不符,请联系dt猫网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【PyQt】16-剪切板的使用

文章目录 前言一、代码疑惑快捷键 二、现象2.1 复制粘贴文本复制粘贴 2.2 复制粘贴图片复制粘贴 2.3 复制粘贴网页 总结 前言 1、剪切板的使用 2、pycharm的编译快捷键 3、类的属性和普通变量的关系 4、pyqt应该养成的编程习惯-体现在代码里了,自己看看。 一、代码…

Bert Encoder和Transformer Encoder有什么不同

前言:本篇文章主要从代码实现角度研究 Bert Encoder和Transformer Encoder 有什么不同?应该可以帮助你: 深入了解Bert Encoder 的结构实现深入了解Transformer Encoder的结构实现 本篇文章不涉及对注意力机制实现的代码研究。 注:…

CentOS7 利用remi yum源安装php8.1

目录 前言remi yum源remi yum源 支持的操作系统remi yum源 支持的php版本 安装epel源安装remi源安装 php8.1查看php版本查看php-fpm服务启动php-fpm服务查看php-fpm服务运行状态查看php-fpm服务占用的端口查看 php8.1 相关的应用 前言 CentOS Linux release 7.9.2009 (Core) …

maven项目引入私有jar,并打包到java.jar中

私有jar存放位置 maven依赖 <dependency><groupId>com.hikvision.ga</groupId><artifactId>artemis-http-client</artifactId><version>1.1.10</version><scope>system</scope><systemPath>${project.basedir}/s…

【ETCD】简介安装常用操作---图文并茂详细讲解

目录 一 简介 1.1 etcd是什么 1.2. 特点 1.3. 使用场景 1.4 关键字 1.5 工作原理 二 安装 2.1 etcd安装前介绍 2.2 安装 2.3 启动 2.4 创建一个etcd服务 三 常用操作 一 简介 1.1 etcd是什么 etcd是CoreOS团队于2013年6月发起的开源项目&#xff0c;它的目标是构建…

面向对象中类与对象

思考系统1000个对象逻辑结构 理解系统1000个对象物理结构 对象this 引用 类的静态变量和静态函数 静态变量和静态函数属于类本身&#xff0c;而不是类的实例。它们可以在不创建类的实例的情况下直接通过类名访问。静态变量在内存中只有一份拷贝&#xff0c;被所有实例共享&…

从一个问题开始聊聊clickhouse的物化视图

【问题】 今天有A问我一个问题&#xff0c;我明明创建了一个物化视图&#xff0c;源表是有数据的&#xff0c;为什么查询物化视图就没有数据&#xff1f; 创建物化视图的SQL示意如下&#xff1a; CREATE MATERIALIZED VIEW schema1.test_mvon cluster clusterNameTO schema1…

docker三剑客compose+machine+swarm小结

背景 在容器领域&#xff0c;不少公司会使用docker三剑客composemachineswarm进行容器编排和部署&#xff0c;本文就简单记录下这几个工具的用法 三剑客composemachineswarm compose compose主要是用于容器编排&#xff0c;我们部署容器时&#xff0c;容器之间会有依赖&…

单链表的实现(数据结构)

本篇博客主要是单链表&#xff08;无头单项不循环&#xff09;的实现的代码分享 说明&#xff1a;因为此单链表无头&#xff08;哨兵位&#xff09;&#xff0c;可以说成没有初始化也可以说初始化时没有一个有效地址作为单链表的起始地址 例如下面代码中的plist NULL。 所以在…

神经网络(neural network)

在这一章中我们将进入深度学习算法&#xff0c;学习一些神经网络相关的知识&#xff0c;这些是有更加强大的作用&#xff0c;更加广泛的用途。 神经元和大脑(neurons and the brain): 我们对于我们的编程的进步主要来自我们对于大脑的研究&#xff0c;根据我们对于大脑的研究…

9*9乘法表

单层循环 99乘法表 动手写写画画整个程序的输出结果

Python 读取写入excel文件

使用Python读取和写入excel的xlsx、xls文件 目录 读取xlsx文件 安装三方库 引入三方库 读取数据 打开文件 表名 最大行数 最大列数 读取一张表 读取整个文件 返回xls整体内容 安装三方包 读取内容 写入xls文件 引入三方库 创建文件并写入数据 报错及解决 报错…

Vue.js+SpringBoot开发无代码动态表单系统

目录 一、摘要1.1 项目介绍1.2 项目录屏 二、功能模块2.1 动态类型模块2.2 动态文件模块2.3 动态字段模块2.4 动态值模块 三、系统设计3.1 用例设计3.2 数据库设计3.2.1 动态类型表3.2.2 动态文件表3.2.3 动态字段表3.2.4 动态值表 四、系统展示五、核心代码5.1 查询档案类型5.…

电脑不小心格式化了,怎么恢复?

在这个数字化时代&#xff0c;电脑已经成为我们日常生活和工作中不可或缺的工具。然而&#xff0c;有时我们可能会不小心格式化电脑硬盘&#xff0c;导致重要数据的丢失。那么&#xff0c;电脑不小心格式化了&#xff0c;怎么恢复&#xff1f; 别着急&#xff0c;在本篇攻略中&…

unicloud 云数据库概念及创建一个云数据库表并添加记录(数据)

云数据库概念 uniCloud提供了一个 JSON 格式的文档型数据库。顾名思义&#xff0c;数据库中的每条记录都是一个 JSON 格式的文档。 它是 nosql 非关系型数据库&#xff0c;如果您之前熟悉 sql 关系型数据库&#xff0c;那么两者概念对应关系如下表&#xff1a; 关系型JSON 文…

基于React的低代码开发:探索应用构建的新模式

&#x1f308;个人主页: Aileen_0v0 &#x1f525;热门专栏: 华为鸿蒙系统学习|计算机网络|数据结构与算法|MySQL| ​&#x1f4ab;个人格言:“没有罗马,那就自己创造罗马~” #mermaid-svg-OywB1Epu30PrvOJQ {font-family:"trebuchet ms",verdana,arial,sans-serif;f…

智慧农业5G融合方案(2)

应用场景 农业航空 农业农村部资料显示,植保无人机具有机动灵活、喷施效率高、施药效果好等特点,能够克服复杂地形条件下地面喷雾机具进地难的问题。目前我国有400多家企业从事植保无人机研发、生产、销售等全产业链业务。主要机型以电动多旋翼为主。2018年作业面积约2.67亿…

ElasticSearch之文档的存储

写在前面 本文看下文档的存储相关内容。 1&#xff1a;如何确定文档存储在哪个分片&#xff1f; 我们需要确保文档均匀分布在所有的分片中&#xff0c;避免某些部分机器空闲&#xff0c;部分机器繁忙的情况出现&#xff0c;想要实现均匀分布我们可以考虑如下的几种分片路由算…

Python数据分析库之pandera使用详解

概要 在数据科学和数据分析中,数据的质量至关重要。不良的数据质量可能导致不准确的分析和决策。为了确保数据的质量,Python Pandera 库应运而生。本文将深入介绍 Python Pandera,这是一个用于数据验证和清洗的库,并提供丰富的示例代码,帮助大家充分利用它来提高数据质量…

FPGA AXI4总线操作教程

AXI&#xff08;Advanced Extensible Interface&#xff09;总线是一种高性能、低延迟的片上系统&#xff08;SoC&#xff09;接口标准&#xff0c;广泛应用于现代数字系统设计中。它允许不同的硬件组件以高效、可靠的方式进行数据传输和控制。本教程将介绍AXI总线的基本操作和…