不定期备考小tips[数模][0] #20210529

---

本专栏主要作个人笔记，有相关知识预备的同学也可作复习用。不保证无相应基础的人士能看明白。
万一考试考到了，或者对你的学习有较大帮助，一键三连不过分吧（斜眼笑）

边值问题

S-L微分方程中两种内积的理解

一种内积：$(u,v)=\int uvdx$

S-L微分方程$\mathcal{L}\varphi +\lambda \sigma \varphi =0$，其中$\mathcal{L}\varphi =\frac{d}{dx}(p\frac{d\varphi }{dx})+q\varphi ,p>0,\sigma >0$
此处定义内积$(u,v)={\int }_a^{b}uvdx$，那么直接计算，得到格林公式
$$(u,\mathcal{L}v)-(v,\mathcal{L}u)=p(u{v}^{\prime }-v{u}^{\prime }){|}_a^b$$
对于Regular S-L问题，加入边界条件
$${\beta }_1\varphi (a)+{\beta }_2{\varphi }^{\prime }(a)=0,{\beta }_3\varphi (a)+{\beta }_4{\varphi }^{\prime }(a)=0$$
此时容易验证格林公式右侧变为0。则这时$\mathcal{L}$是自伴算子（可以类比实对称矩阵相应的线性变换）。

另一种内积：$(u,v{)}_1=\int u\sigma vdx$

实际上，我们要求对$\mathcal{L}\varphi +\lambda \sigma \varphi =0$的解进行考察，得到正交等性质。此处${\varphi }_1$与${\varphi }_2$正交的定义是$({\varphi }_1,{\varphi }_2{)}_1:={\int }_a^b{\varphi }_1\sigma {\varphi }_{2}dx=0$，内积定义和上一节所述不同。那么上一节仅仅证明了$(\cdot ,\cdot )$上的自伴性，又怎么能考察新的内积下性质呢？
这个问题可以类比线性代数中相应的问题。设$A$实对称，则$Ax=\lambda x$的各个解（特征向量）正交。此时设一个正定对角矩阵$D$，考察$Ax=\lambda Dx$的解。这等价于考察$S^{T}AS(S^{-1}x)=\lambda S^{T}DS(S^{-1}x)$的解，其中$S$是可逆矩阵（但不一定是正交矩阵）。
根据线性代数容易证明存在$S=D_{1}T$使得$S^{T}DS(=I),S^{T}AS$都是对角矩阵。其中$D_1$是正定对角矩阵，$T$是正交矩阵，$D_{1}D{D}_1=I$。
此时设两个解$x_1,x_2$满足$S^{T}AS(S^{-1}{x}_1)={\lambda }_1(S^{-1}{x}_1),S^{T}AS(S^{-1}{x}_2)={\lambda }_2(S^{-1}{x}_2),{\lambda }_1\ne {\lambda }_2$，则根据$S^{T}AS$对角容易证明$(S^{-1}{x}_1,S^{-1}{x}_2):=x_1^T(S^{-1}{)}^T{S}^{-1}{x}_2=0$（第一种内积）
此时则$(x_1,x_2{)}_1:=x_1^{T}D{x}_2=x_1^T(S^{-1}{)}^T{S}^{T}DS{S}^{-1}{x}_2=x_1^T(S^{-1}{)}^T{T}^T{D}_1^{T}D{D}_{1}T{S}^{-1}{x}_2=x_1^T(S^{-1}{)}^T{S}^{-1}{x}_2=0$（第二种内积）
这可以帮助我们理解为什么$(\cdot ,\cdot )$这个旧内积上的自伴性可以决定$(\cdot ,\cdot {)}_1$这个新内积对应的性质。

F.A.

F.A. (Fredholm Alternative, 弗雷德霍姆二择一)实际上是两次判断，三种情况，千万注意不是简单的”二选一“。此处的整个判断过程可以类比解线性方程组时的过程。

线性方程组

为了类比S-L微分方程，可以考察实数域上$Ax=b$，其中$A$是对角元依次严格递增的对角矩阵或者和这样的对角矩阵正交相似的实对称矩阵。
注：为什么要考察这样的矩阵？因为Regular S-L特征值问题可以证明具有类似的性质。具体证明暂时不用管。
若$A$可逆（等价于$Ax=0$没有非平凡解）则方程有唯一解。
否则

• $A$不可逆，$Ax=0$有非平凡解。
• 根据$A$正交标准形中对角元至多一个是0，得到$dimKer(A)=1$。
• 由实对称矩阵的性质，$Ker(A)$与$Im(A)$正交。
  于是容易看出，$b\in Im(A)$等价于存在一个非零的$\alpha \in KerA$与$b$正交。并且注意$b\in Im(A)$即$\exists \beta ,A\beta =b$，此时$\beta +k\alpha$对于一切$k$都是$Ax=b$的解。
• $b\notin Im(A)$时，$Ax=b$无解。

总之，当$A$不可逆，如果存在非零$\alpha \in KerA$与$b$正交，则$Ax=b$有无穷多解，否则$Ax=b$无解。（特别地，若$b=0$则$Ax=0$在$A$不可逆时总是有无穷多解）

S-L微分方程

若$\mathcal{L}\varphi =0$没有非平凡解，则$\mathcal{L}\varphi =f$有唯一解。
否则，设$\mathcal{L}\varphi =0$一个非平凡解${\varphi }_0$。此时：

1. 当$({\varphi }_0,f)=\int {\varphi }_{0}fdx=0$时，方程有无穷多解。
2. 当$({\varphi }_0,f)=\int {\varphi }_{0}fdx\ne 0$时，方程无解。

无量纲化

换元的一般表示

比如有一些有量纲物理量$M,L,T$，则一般的做法是设$M=M_{c}m,L=L_{c}l,T=T_{c}t$，其中$M_c$是”特征“物理量，$m=\frac{M}{M_c}$无量纲，是实际物理量和特征物理量的比。
而涉及导数时，由$x=\frac{X}{X_c},t=\frac{T}{T_c}$容易知道$\frac{dx}{dt}=\frac{T_c}{X_c}\frac{dX}{dT}$，显然看到$\frac{dx}{dt}$无量纲。
涉及高阶导是一个坑。$\frac{d^{2}x}{d{t}^2}=\frac{d}{dt}\frac{dx}{dt}=\frac{T_c^2}{X_c}\frac{d^{2}X}{d{T}^2}$，分母$X_c$不是二次！

如何选择特征物理量

在有些物理量趋于0或趋于无穷时，如何适当选择物理量。

• 无量纲化后的方程中的参数不能发散。
• 不能出现”不定式“。（详细意义之后阐述）

例如$\frac{d^{2}Y}{d{T}^2}=-\frac{4/3\cdot \pi R^3\rho G}{(R+Y{)}^2}$，$Y^{\prime }(0)=-V_0,Y(0)=2m(2米)$，终态$Y_{final}=1m$，求$R$在极限情况（且$R$趋于极限过程中$\rho ,G,V_0$都是恒量）时，初态到终态运动时间。
首先，先设待定的$L_c,T_c$为特征长度和特征时间，则
$$\ddot{y}=-\frac{4/3\cdot \pi \rho G{T}_c^2\cdot (R/L_c{)}^3}{y+(R/L_c{)}^2},\dot{y}(0)=-V_0{T}_c/L_c,y(0)=2m/L_c$$

• 错误示例1：在$R\to 0$时，如果令$L_c=R$，那么$y(0)$发散，这不可取。
• 错误示例2：在$R\to \infty$时，如果令$L_c=R$，对于非零有限数$V_0$，如果取特征时间$T_c$使得$V_0{T}_c/L_c=1$那么$T_c\to \infty$，无量纲化后初速$\dot{y}(0)=-1$且$y(0)=y_{final}=0$，则从初态到终态需要0个特征时间。
  0个特征时间，每个特征时间是无穷大，这是$0\cdot \infty$不定式，不可取。
• 参考示例1：在$R\to \infty$时，如果令$L_c=R$，对于非零有限数$V_0$，如果取特征时间$T_c$是非零有限数（不趋于无穷），那么$\dot{y}(0)=0$但是$\ddot{y}$始终非零有限，这没有出现不定式，则可以确定地下结论：物体几乎瞬间从初态运动至终态。这符合实际物理情况。
• 注1：虽然可以对物体运动所需时间下结论，但这时考察物体在任意时刻运动了多少距离就很困难，这是因为物体在所考察$Y(0)=2m,Y_{final}=1m$初终态范围内，始终运动了0个特征长度，每个特征长度又是无穷大，那么对于长度出现了$0\cdot \infty$不定式，考察长度相关的东西自然产生困难。
• 正确示例1：$L_c=1m,\frac{4}{3}\pi \rho G{T}_c^{2}R/L_c=1$（可以看出$T_c\to 0$），此时$\ddot{y}=-1,\dot{y}=0,y(0)=2$，那么物体任意时刻运动了多少距离可以看出，并且初态到终态运动所需时间为$\sqrt{2}$个特征时间，也就是0也符合实际情况！
• 注2：这些结论的前提都是$\rho$保持不变。如果质量$\frac{4}{3}\pi \rho R^3$保持不变，情况又有所不同

奇异极限

随手举例：$\dot{x}=y,\alpha \dot{y}=y-x$
如果$\alpha \to 0$，方程变为$\dot{x}=x=y$
初始时应该满足$x=y$，否则矛盾。（奇异极限会要求初值条件）
实际上，初始时如果$x>y$，那么$\dot{y}\to \infty$。

• 如果$\alpha <0$，那么$\dot{y}\to +\infty$，也就是$y$快速增大直至$x=y$附近，$y-x$和$\alpha$同一量级。这是合理的奇异极限。
• 如果$\alpha >0$，那么$\dot{y}\to -\infty$，也就是$x-y$会越来越大，最终无法使得系统收敛到奇异极限！

渐近展开

主项平衡原理中“同阶”的理解

$x^2-x+\frac{1}{4}ϵ=0,ϵ\to 0$
现在要把$x$表示成$ϵ$的函数。且我们想考察这个函数$x(ϵ)$的“前若干项近似”。
设$x=x_0{\delta }_0(ϵ)+x_1{\delta }_1(ϵ)+\cdots$，且在$ϵ\to 0$时${\delta }_0\gg {\delta }_1\gg \cdots$，对于代数方程，我们可以认为$x_i$都是常数，和$ϵ$无关。不妨只考察首项$x_0\ne 0,{\delta }_0(ϵ)\ne 0$的情况（非平凡解）。
考察$x(ϵ)$展开的首项$x_0{\delta }_0$，有$x_0^2{\delta }_0^2-x_0{\delta }_0+\frac{1}{4}ϵ=0$.现在对前式的三项分类讨论。

• 第一项是主项或第二项是主项，则${\delta }_0=0$，舍去
• 第三项是主项，无解，舍去
• 三项都是主项，则现在整个方程是一个关于$x_0$的一元二次方程，系数和$ϵ$相关。为了这三项确实都有必要存在，他们应该是同阶的。不妨把同阶条件写成（之后解释）${\delta }_0^2={\delta }_0=ϵ$，而这是不可能的。也就是这三项不可能在实际中同阶。
• 第一项和第二项是主项，则整个方程$x_0^2{\delta }_0^2-x_0{\delta }_0=0$，系数${\delta }_0,{\delta }_0^2$都和$ϵ$相关。为了这两项都有必要存在，它们在$ϵ\to 0$时同阶。故不妨设${\delta }_0=1$。这时，被略去的第三项$\frac{1}{4}ϵ$确实很小，所以这是合理的，则$x_0^2-x_0=0,x_0\ne 0,x_0=1$（注意舍去了$x_0=0$的情况！）
• 对同阶条件的说明：由于$x_0$是非零有限数，所以各项系数同阶也就是各项同阶。至于${\delta }_0$和${\delta }_0^2$同阶如何表示，你可以最简单地取${\delta }_0={\delta }_0^2$，也可以设$ln({\delta }_0+1)=114514{\delta }_0^2$，这都不会影响最后的结果。不信你看：$x_0^2{\delta }_0^2-x_0{\delta }_0=0$，解得$x_0=1/{\delta }_0$.回忆刚刚$x(ϵ)$的第一项$x_0{\delta }_0=1\cdot 1=1$，现在仍然有第一项$x_0{\delta }_0=1$（注意$x_0,{\delta }_0$都是常数），确实对求$x(ϵ)$的渐近展开式没有影响！
• 第一项和第三项是主项，则${\delta }_0=\sqrt{ϵ}$，第二项就太大了。
• 第二项和第三项是主项，则${\delta }_0=ϵ$，第一项${ϵ}^2{x}_0^2$确实很小。所以这也是合理的，且$-ϵx_0+\frac{1}{4}ϵ=0,x_0=\frac{1}{4},x(ϵ)\approx \frac{1}{4}ϵ$