Skip to content

连续方程的存在性终点任务 ​

返回核心路线

本组任务从两种近似对象出发:Euler 折线与有限维凸包中的近似不动点。目标是证明它们留下的极限真正满足方程,并准确说明结论停在哪里。第一至三题完成 ODE 核心任务,第四至五题完成非压缩积分方程,第六至七题检查换条件后的迁移与边界。

需要补课时,先读一致收敛、连续模和Ascoli 定理。比较方法时,可按条件选读压缩映射、ODE 唯一性、有限维不动点与变分直接法。

任务 ​

  1. 对 y′=t+|y|、y(0)=0,在 [0,1/4]×[−1,1] 上选共同高度界,并证明任意 Euler 分割都不逃出矩形。使用步长 1/10,把最后不足整步的一段也写出来,求全部网点值
  2. 推导最大步长为 h 时的一致积分残差界。证明 h≤1/1000 足以使证书小于 1/100,并说明证书是否同时给出了到指定解的距离
  3. 从任意 hn→0 的分割序列取得经典解,写明通过积分的极限估计。把方程改为二维系统 x1′=t+|x2|、x2′=|x1|、x(0)=0,在欧氏单位球中重新计算共同界。若时间向左,怎样改写
  4. 在实直线的紧集 K=[0,2] 上取网点 0,1,2,半径 ε=3/2,对 y=3/4 算出 Schauder 证明中的连续权重与加权点。解释为何这三个网点不必线性无关,并补齐从近似不动点到真不动点的子列论证
  5. 在 ‖u‖∞≤2 的连续函数球上认证
(Fu)(x)=x+14∫01(1+x)extsin⁡(8u(t))dt.

分别检查自映射、连续性和紧像。先求零函数沿常函数方向在 x=0 的导数,再用 u≡1/64 与 v≡0 给出一个有限的非压缩见证。你究竟得到了哪几项存在、唯一、迭代收敛或停止证书 6. 将上一题换成

(Tu)(x)=x2+14∫01(1+|x−t|)|u(t)|dt,‖u‖∞≤2.

核不可处处求导、非线性也不 Lipschitz。重新给出三项适用条件,判断是否仍有解 7. 为什么不能把“像相对紧”换成“定义域闭且有界”?验证 ℓ2 单位球上 Sx=(1−‖x‖22,x1,x2,…) 的失败机制,并说明上一题为什么不直接属于变分直接法的能量输入

Euler 折线的答案 ​

一:总时间控制位移,末段使用实际长度 ​

右端绝对值至多为 M=5/4。每段斜率都至多为 M,因此把之前所有段及当前段的位移相加,得到 |yh(t)|≤Mt≤5/16<1。这个估计可以逐段归纳使用,保证每次求右端都仍在给定矩形内,避免先假定全部折线已存在。

步长 1/10 给网格

0,110,15,14,

末段长 1/20。节点值依次为

y0=0,y1=0,y2=110⋅110=1100,y3=1100+120(15+1100)=140.

在节点之间按当前左端斜率作直线。不能为了凑整步把终点推进到 3/10,那已经换了任务区间。

二:小的是方程残差 ​

一段内 |s−τj|≤h,|yh(s)−yh(τj)|≤5h/4。平方根不等式给

|f(s,yh(s))−f(τj,yh(τj))|≤h+5h/4.

逐段积分,统一得到

‖Rh‖∞≤h+5h/44.

h=1/1000 时,1/800<39/1000,故证书小于 (1+39)/4000=1/100。更小的 h 也满足。这里没有选出特定解,也没有证明残差到解误差的稳定性,因此不能把这个 1/100 改成解误差。

三:极限论证与向量迁移 ​

共同高度和共同 Lipschitz 界让 Ascoli 给出 yhnk→y 一致。紧矩形使右端一致连续。若 ek=‖yhnk−y‖∞,则在这个例子中

sups|f(s,yhnk(s))−f(s,y(s))|≤ek.

因此积分方程的极限误差至多为

ek+14(hnk+5hnk/4)+14ek,

它趋于零。极限满足积分方程,连续被积函数再由微积分基本定理给经典解。全列收敛、唯一与速率均没有由这一步得到。

二维单位球内,两个右端分量至多分别为 5/4 与 1,故欧氏界可取 M=41/4;M/4=41/16<1,不逃出。先后对两个坐标抽取子列可保留两坐标的一致收敛。重新组合分量误差,得到

‖Rh‖∞≤14(h+Mh)2+Mh.

向左时令 z(s)=y(t0−s),右端变成 −f(t0−s,z),高度界不变。一般无限维状态不能只靠逐点有界重复这个抽取步骤;Peano 页说明了有限维条件的具体用途。

紧像不动点的答案 ​

四:网点可以相关,权重必须连续 ​

三个未归一化权重分别为

32−34=34,32−14=54,32−54=14.

分母为 9/4,所以归一化权重为 1/3,5/9,1/9,加权点为

Pε(3/4)=0⋅13+1⋅59+2⋅19=79.

误差为 |7/9−3/4|=1/36<3/2。这里 Pε(y)≠y,可见它不必是固定原集每一点的投影。三个网点生成的凸包只是直线区间,实际仿射维数为一;Brouwer 在这个一维空间中应用,不必把网点变成三条独立方向。

一般证明从 Gε=PεF 的不动点得到 ‖xn−Fxn‖<1/n。紧像给 Fxnk→y,同一残差使 xnk→y。闭性使极限留在 C,连续性给 Fxnk→Fy,所以 Fy=y。抽取的是像中的序列,不能凭 C 闭有界就抽取原序列。

五:非压缩积分方程的完整结论 ​

利用 (1+x)ext≤2et 与正弦的绝对值至多一,有

‖Fu‖∞≤1+e−12<2.

正弦的 Lipschitz 常数给

‖Fu−Fv‖∞≤4(e−1)‖u−v‖∞,

所以映射连续。关于 x 的导数满足

|(Fu)′(x)|≤1+14∫01ext[1+(1+x)t]dt≤1+3e4.

共同高度和共同增量界让 Ascoli 给出像相对紧。Schauder 遂给出一条满足原方程的连续函数;它还因属于像而具有连续导数。零函数不是解,因为 F0(x)=x。

常函数 ua≡a 在 x=0 给 (Fua)(0)=sin⁡(8a)/4,相对 a 在零处的导数为 2。要用一个有限输入直接作见证,取 a=1/64,则

‖Fua−F0‖∞‖ua−0‖∞≥16sin⁡(1/8)≥16(18−(1/8)36)=383192>1.

所用 sin⁡z≥z−z3/6 在 0≤z≤1 成立:1−cos⁡s≤s2/2,积分即得 z−sin⁡z≤z3/6。这是真正的下界,足以否定任何 q<1 的压缩常数;仅仅算出一个大于一的上界做不到。

最终结论是存在一条非零解及上述高度、导数约束。唯一性、直接迭代收敛和到指定解的数值停止证书均未得到。具体证明见Schauder 页。

六:不求核的导数也可以取得紧像 ​

||u||≤2,而

supx∫01(1+|x−t|)dt=32,

所以 ‖Tu‖∞≤1/2+32/8<2。平方根差的估计给

‖Tu−Tv‖∞≤38‖u−v‖∞,

仍足以证明连续。由 ||x−t|−|y−t||≤|x−y|,

|(Tu)(x)−(Tu)(y)|≤(12+24)|x−y|.

因此像共同有界且等度连续,三项条件都成立,仍有解。迁移时保留的是统一连续控制,不是上一题的求导公式。

七:准确定位失去的条件 ​

S 在 ℓ2 单位球上连续且 ‖Sx‖2=1。固定点方程迫使所有坐标相等,平方可和性迫使它们全零,又与首坐标方程矛盾。其像包含两两距离 2 的 en+1,所以失去的是相对紧像。

变分直接法需要能量、强制性和弱下半连续性。这里给出的是自映射与它的像控制;这些材料本身并没有产生满足要求的能量。若已有能量模型,可转向原来的弱紧性与极小点路线,其既有一维、三维任务仍各自成立。

可运行的有限检查 ​

下载 Python 标准库检查器。直接运行会用有理数区间传播 Euler 节点,检查末段、自映射常数、有限网权重和有限非压缩见证。默认参数为最大步长 1/1000、残差预算 1/100;可用参数 --step 1/10 查看第三个非整步区间,用 --residual-tol 调整预算,用 --nodes 输出全部节点的有理上下界。

平方根通过整数平方比较得到包围区间;e 通过有限 Taylor 和与正项尾界包围。因此被标记为通过的是这些有限算术检查。若所选步长的上界超过预算,程序报告该预算未获认证;它不声称实际残差必然超过预算。超出预设步数预算则报告受阻,不悄悄缩短时间区间。

程序明确将无限序列的子列抽取和固定点存在证明标为未运行。数值轨道、有限网和算术证书帮助复核正文的具体输入,不会用有限次计算替代两个存在定理的证明。