能量不等式证明

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

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.mzph.cn/news/732282.shtml

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

相关文章

Flutter图片内存占用过大问题

图片(Image)加载原理: Image : 显示图⽚的Widget,通过ImageState管理ImageProvider的⽣命周期。 ImageProvider: 图⽚的抽象概念。 根据Image创建实例时调用的工厂方法的不同(Image.network或者Image.assetImage&#…

【FPGA/IC】CRC电路的Verilog实现

前言 在通信过程中由于存在各种各样的干扰因素,可能会导致发送的信息与接收的信息不一致,比如发送数据为 1010_1010,传输过程中由于某些干扰,导致接收方接收的数据却成了0110_1010。为了保证数据传输的正确性,工程师们…

禁止使用搜索引擎,你了解吗?

员工A:“我今天想搜索的时候,用不了浏览器了,你能用么?” 员工B:“不知道啊我试一下啊” “也不行” 员工C:“为什么啊?” 针对上述对话,我们不禁思考: 公司为什么禁…

python基础9_序列类型

回顾: 什么是变量?,有什么用? 可以变化的量, 就是个容器,多次变化,方便后续使用, 前面介绍了哪些数据类型? bool, str, int, float 用什么函数查看数据的类型? a "hello" print(type(a)) 到了这一步,,我们认识了哪些数据类型呢? int 整型(整数), float…

office下常见问题总结——(持续更新学习记录中......)

目录 Wordword2019中, 当给选定的汉字设置格式后,其他相同汉字也会自动应用相同的格式?在Word中,当输入数字后加上句点(.)时会自动被识别为标题,如何关闭功能?如何让当前的word中的样式 ,匹配全局模版中的样式?在word中,为什么…

一、NLP中的文本分类

目录 1.0 文本分类的应用场景 1.1 文本分类流程 ​编辑 1.2 判别式模型 1.3 生成式模型 1.4 评估 1.5 参考文献 1.0 文本分类的应用场景 (1)情感分析:中性,正向评论,负向评论,黄色言论,暴…

Java基础 - 8 - 算法、正则表达式、异常

一. 算法 什么是算法? 解决某个实际问题的过程和方法 学习算法的技巧? 先搞清楚算法的流程,再直接去推敲如何写算法 1.1 排序算法 1.1.1 冒泡排序 每次从数组中找出最大值放在数组的后面去 public class demo {public static void main(S…

Matlab偏微分方程拟合 | 完整源码 | 视频教程

专栏导读 作者简介:工学博士,高级工程师,专注于工业软件算法研究本文已收录于专栏:《复杂函数拟合案例分享》本专栏旨在提供 1.以案例的形式讲解各类复杂函数拟合的程序实现方法,并提供所有案例完整源码;2.…

力扣hot100题解(python版63-68题)

63、搜索插入位置 给定一个排序数组和一个目标值,在数组中找到目标值,并返回其索引。如果目标值不存在于数组中,返回它将会被按顺序插入的位置。 请必须使用时间复杂度为 O(log n) 的算法。 示例 1: 输入: nums [1,3,5,6], target 5 输…

【npm】前端工程项目配置文件package.json详解

简言 详细介绍了package.json中每个字段的作用。 package.json 本文档将为您介绍 package.json 文件的所有要求。它必须是实际的 JSON,而不仅仅是 JavaScript 对象文字。 如果你要发布你的项目,这是一个特别重要的文件,其中name和version是…

王升:Audio电感对车载功放EMC的影响 | 演讲嘉宾公布

一、智能车载音频 II 专题论坛 智能车载音频 II 专题论坛将于3月28日同期举办! 我们正站在一个前所未有的科技革新的交汇点上,重塑我们出行体验的变革正在悄然发生。当人工智能的磅礴力量与车载音频相交融,智慧、便捷与未来的探索之旅正式扬帆…

若依/RuoYi-Vue使用docker-compose部署

系统需求 JDK > 1.8 MySQL > 5.7 Maven > 3.0 Node > 12 Redis > 3 思路 前端服务器 nginx 后端服务器代码打包 java、maven、node 数据库/缓存 mysql、redis 开始 创建目录ruoyi并进入 克隆若依代码 git clone RuoYi-Vue: 🎉 基于Spring…

试用期自我总结报告10篇

试用期自我总结报告(篇1) 一转眼试用期的时间飞快就过去了,在这段时间里我学习到了很多,也把自己在过去学习的东西得已融会贯通。能够来到幼儿园里成为一名老师是我一直以来的目标,而我也终于完成了自己的目标&#x…

Winform窗体随着屏幕的DPI缩放,会引起窗体变形及字体变形,superTabControl标签字体大小不匹配

一、前言 superTabControl做的浏览器标签(cefsharp)在缩放比例(125%,150%时字体不协调) 物联网浏览器,定制浏览器,多媒体浏览器(支持H264)参考栏目文章即可 二、配置参数 app.manifest参数 dpiAware =true <application xmlns="urn:schemas-microsoft-c…

用 Axios 提升前端异步请求的效率

&#x1f90d; 前端开发工程师、技术日更博主、已过CET6 &#x1f368; 阿珊和她的猫_CSDN博客专家、23年度博客之星前端领域TOP1 &#x1f560; 牛客高级专题作者、打造专栏《前端面试必备》 、《2024面试高频手撕题》 &#x1f35a; 蓝桥云课签约作者、上架课程《Vue.js 和 E…

ubuntu22.01安装及配置

前言 本次安装基于VMware Pro 16进行安装。 ubuntu版本&#xff1a;ubuntu-22.04.3-live-server-amd64.iso 1、下载 1.1官网下载 https://ubuntu.com/download 1.2、清华大学镜像网站下载 https://mirrors.tuna.tsinghua.edu.cn/ 进入网站后搜索ubuntu&#xff0c;选择ubu…

Claude3超大杯发布,将取代ChatGPT4?

演示站点&#xff1a; https://ai.uaai.cn 创作模块 &#xff08;Claude3已接入&#xff0c;欢迎体验&#xff09;官方论坛&#xff1a; www.jingyuai.com 前言 今天【超越GPT-4 接近人类水平&#xff01;Claude 3系列大模型惊艳问世】登上了全网热搜&#xff0c;那么【超…

微信小程序实现上拉加载更多

一、前情提要 微信小程序中实现上拉加载更多&#xff0c;其实就是pc端项目的分页。使用的是scroll-view&#xff0c;scroll-view详情在微信开发文档/开发/组件/视图容器中。每次上拉&#xff0c;就是在原有数据基础上&#xff0c;拼接/合并上本次上拉请求得到的数据。这里采用…

电动机工作原理图

电机&#xff08;俗称“马达”&#xff09;是指依据电磁感应定律实现电能转换或传递的一种电磁装置。它的主要作用是产生驱动转矩&#xff0c;作为用电器或各种机械的动力源。 电动机&#xff0c;转子置于旋转磁场中&#xff0c;在旋转磁场的作用下&#xff0c;获得一个转动力…

SSRF漏洞基础原理(浅层面解释 + 靶场演示)

一、SSRF漏洞的基本概念&#xff1a; SSRF--全名&#xff1a;Server-Side Request Forgery&#xff0c;汉译&#xff1a;服务端请求伪造&#xff0c;漏洞别名“借刀杀人”。 想象以下&#xff0c;现存在一个 Web应用&#xff0c;这个Web应用可以帮助我们能爬取互联网上的其他…