8.3 脈衝函數與 Heaviside#
本章導言#
真實世界的輸入不溫柔:開關合上(階躍)、錘子敲擊(脈衝)、閃電擊中(尖峰)、雷射脈衝。處理它們需要兩件數學武器:Heaviside 階躍函數 \(u(t)\) 與Dirac 脈衝 \(\delta(t)\) 。後者在 1930 年代被物理學家 Paul Dirac(狄拉克,1902–1984) 隨手使用(「在 \(t = 0\) 處無窮高、積分為 1、其他地方為零」的怪獸),直到 1945 年 Laurent Schwartz(施瓦茨,1915–2002) 的分布理論才賦予嚴格身份——Schwartz 因此獲 1950 年菲爾茲獎。本節講述這對工具的變換性質與應用。
Heaviside 階躍:開關的數學#
時移律(8.1 節性質 5):
「延遲 = 乘 \(e^{-as}\) 」——一切時間延遲系統(傳輸線、加工流水線、神經傳導延遲、控制系統的死區時間)的統一記號。方波: \(u(t) - u(t - T)\) ;脈衝串: \(\sum u(t - nT)\) ——複雜輸入由階躍積木搭成。
Dirac 脈衝:極限的幽靈#
直觀構造:寬 \(\varepsilon\) 、高 \(\frac{1}{\varepsilon}\) 、面積 1 的矩形 \(d_\varepsilon(t)\) , \(\varepsilon \to 0\) :
操作性質(定義的現代形式):
** \(\delta\) 是「取值機」不是「函數」:它作用於函數、吐出一個數——Schwartz 的定義:分布 = 檢驗函數空間上的連續線性泛函**。
變換與微分的禮物:
更驚人的: \(u' = \delta\) (階躍的導數是脈衝——「從 0 跳到 1 的速度是無窮尖峰」)。在分布意義下微分是無障礙的:處處不可微的階躍、甚至更野的對象都能求導——分布理論把微分學的疆域推到極限之外。Schwartz 的口號:「微分不再需要連續。」
捲積身份: \(\delta\) 是乘法單位元#
** \(\delta * f = f\) ——脈衝是捲積代數的單位元(如同 1 是乘法的單位元)。8.2 節場景二現在獲得完備表述:系統對 \(\delta\) 的響應 \(h\) = 傳遞函數的逆變換;任意輸入的響應 = \(h * u\) 。「敲一下,聽全聲」(模態測試)、「打一針,看血藥濃度」(藥動學 bolus 注射)——以單位輸入探測系統**是實驗科學的通用策略, \(\delta\) 是它的數學簽名。
一個具體例題#
例題 8.3.1:系統 \(\ddot{y} + 3\dot{y} + 2y = \delta(t)\) 、 \(y(0^-) = \dot{y}(0^-) = 0\) (錘擊靜止質量—彈簧系統)。求 \(y(t)\) ,並解釋 \(t = 0\) 處發生了什麼。
解:變換(初值為零):
檢查 \(t = 0\) 的邊界行為: \(y(0^+) = 0 = y(0^-)\) ——位置連續; \(y'(t) = -e^{-t} + 2e^{-2t}\) , \(y'(0^+) = 1 \neq 0 = y'(0^-)\) ——速度跳變 1!
物理解讀:錘擊(衝量 \(J = \int F\, dt = 1\) )瞬間改變動量 \(m\Delta v = J\) ,位置來不及動。 \(\delta\) 輸入的效果: \(y\) 連續、 \(y'\) 跳躍(跳幅 = 衝量)。一般地,方程右端含 \(\delta\) 時,最高階導數以下的項保持連續、最高階導數跳變——這是求解脈衝問題的「跳變條件」,等價於變換法自動完成的事。✓(驗證: \(y'(0^+) - y'(0^-) = 1 =\) 衝量。)
後續: \(t > 0\) 後輸入歸零,系統以齊次模態衰減( \(e^{-t}, e^{-2t}\) ——特徵根 \(-1, -2\) 的自由振鈴)。脈衝響應 = 「被打一拳之後的自由振盪」——這就是為什麼敲擊測試能讀出系統的全部模態。
歷史啟示#
- 物理先行的數學對象: \(\delta\) 在 Heaviside(1890s 運算微積分中 \(\frac{d}{dt}u\) )與 Dirac(1930 量子力學《原理》)手中是「合法的非法居民」;Schwartz 1945 年發給它戶口。數學的戶籍制度常比數學實踐晚五十年(呼應 1.1 節無窮小的平反史)。
- 單位元的力量: \(\delta\) 作為捲積單位元,使「系統 = 對單位輸入的響應」成為精確定理——抽象代數的「單位元」概念在最意外的角落獲得物理身體。
- 微分疆域的擴張:階躍可微(= \(\delta\) )、脈衝可再微( \(= \delta'\) )、甚至白噪聲(廣義隨機過程)都進入了微積分——「可微」從函數的屬性變成對象的屬性,這是二十世紀分析學最大的地盤擴張之一。
與其他章節及領域的關聯#
- 8.1、8.2 節:時移律與脈衝響應(本節是它們的主角化)。
- 7.1 節:敲擊測試與模態分析; \(Q\) 從脈衝響應的振鈴讀出。
- 9.2 節:d'Alembert 解與波傳播——脈衝在弦上分為兩個行波;Huygens 原理(3D 波的清晰尾跡)。
- 1.1 節:Berkeley 攻擊的「幽靈」與 Robinson 非標準分析、Schwartz 分布——無窮小/無窮尖峰的三百年平反史在此合流。
- 量子力學:位置本徵態 \(|x\rangle\) 是 \(\delta(x)\) 型分布;動量表象是傅立葉變換——量子力學的整個語法建立在分布上。
- 訊號處理:採樣定理(Shannon 1949)把離散樣本寫成 \(\delta\) 串乘連續信號;數位世界的每一次採樣都是一次 \(\delta\) 應用。
- 偏微分方程:基本解(fundamental solution)——熱核、波傳播的 Green 函數都以 \(\delta\) 為源;9 章的 \(G\) 函數是 \(\delta\) 響應的空間版。
練習#
- 用時移律表示方波 \(f(t) = 1\) ( \(0 < t < 2\) )、 \(0\) (其他)的 Laplace 變換。
- 求 \(y' + y = \delta(t - 2)\) 、 \(y(0) = 0\) 的解,並驗證 \(y\) 在 \(t = 2\) 處的跳變行為(本例是一階: \(y\) 本身跳變!)。
- 求 \(\ddot{y} + \omega_0^2 y = \delta(t)\) (無阻尼)的脈衝響應,並討論其永續振鈴(無衰減 = \(Q \to \infty\) 的極端,7.1 節)。
- 證明 \(u' = \delta\) :對任意檢驗函數 \(\phi\) ,計算 \(\int u \cdot (-\phi')\, dt\) (分布導數的定義 \(\langle u', \phi \rangle = -\langle u, \phi' \rangle\) )並等於 \(\phi(0)\) 。
- 思考題:跳變條件說「 \(\delta\) 使最高階導數跳、其餘連續」。請論證:對 \(n\) 階系統 \(\delta\) 輸入導致 \(y^{(n-1)}\) 跳變 \(1\) (歸一化後),而更低階導數連續——並從牛頓力學(衝量—動量)解釋為什麼這在物理上必然。
進階研讀:分布世界的五個深化場景#
場景一:Schwartz 空間與 Fourier 的完備化。分布理論使每一個局部可積函數都有 Fourier 變換(廣義意義):常數 1 的 Fourier 變換是 \(2\pi\delta(\omega)\) ——直流是零頻率的脈衝; \(\cos\omega_0 t\) 的譜是 \(\pm\omega_0\) 的兩根 \(\delta\) ——純音是頻譜上的兩根針。頻譜分析儀螢幕上的尖峰,數學上是分布。沒有分布理論,沒有現代訊號處理。
場景二:Green 函數與基本解。 \(L G = \delta\) 的解 \(G\) 是算子 \(L\) 的「單位脈衝響應」:熱方程的基本解 \(K(x, t) = \frac{1}{\sqrt{4\pi t}} e^{-x^2/4t}\) (高斯核,9.1 節);波方程的基本解(3D 是球面殼——Huygens 原理,9.2 節)。「點源 → 場」是物理建模的原子操作, \(\delta\) 是點源的數學名字。
場景三:採樣與重建。Shannon 採樣定理:帶限信號 \(f\) ( \(|\hat f|\) 支撐於 \([-B, B]\) )由樣本 \(f(n/2B)\) 完全決定,重建 \(f(t) = \sum f(n/2B) \operatorname{sinc}(2Bt - n)\) 。證明的核心:樣本 = \(f \cdot \delta\) 串,其譜 = 原譜的週期複製。CD 音質 44.1 kHz 的由來(人耳 20 kHz 帶限 × 採樣定理 × 工程餘量)—— \(\delta\) 串支撐了整個數位音訊工業。
場景四:隨機脈衝與白噪聲。白噪聲 = \(\delta\) 相關的廣義隨機過程( \(\mathbb{E}[\dot W(t)\dot W(s)] = \delta(t - s)\) ——維納過程的導數,Itô 微積分 8.1 節場景五的對象)。濾波器輸出「著色」噪聲的計算全靠分布語言。通訊系統的誤碼率分析、金融的波動率建模,都在 \(\delta\) 相關函數上工作。
場景五:微局部分析——現代的深入。奇異性的傳播(波前集,Hörmander 1970s)追蹤「分布在哪裡不光滑、沿哪個方向不光滑」——微分算子的代數(象徵演算)與波的幾何(幾何光學)在微局部層面統一。從 Schwartz 1945 到 Hörmander 1970s,分布理論長成了 PDE 理論的骨架(Atiyah–Singer 指標定理 1963 也以其為語言)。
補充:Schwartz 與其時代#
Laurent Schwartz 是托洛茨基主義者、反殖民主義者(公開反對阿爾及利亞戰爭)、布爾巴基學派成員。1944–45 年在 Grenoble 孤立環境中建立分布理論,1950 年菲爾茲獎(首屆)。他晚年回憶:分布的動機部分來自解決 Heaviside/Dirac 的「非法」操作——「數學家的工作之一,是把物理學家的成功合法化。」1950 年代的教科書把 \(\delta\) 稱為「奇異函數」並小心規避;今天的工程教材隨手使用——數學合法化的傳播需要一代人。
補充:為什麼 \(\delta\) 不是「無窮高的函數」#
三個層次的失敗證明它是泛函而非函數:(1)若 \(\delta(0) = \infty\) ,則不是實值函數;(2)若取「 \(a \cdot \delta_{a>0}\) 峰 + 處處零」的初等函數 \(d\) ,其 Lebesgue 積分 = 0(單點測度零)≠ 1——在正統積分下「面積 1」與「其他處為零」不可兼得;(3)唯一出路:改變「函數」的定義——讓對象以對其他函數的作用被認定(泛函),而非以逐點值被認定。身份的轉換(從「值」到「作用」)是二十世紀數學的普遍運動(亦見:測度、層、概形——Grothendieck 的「以映射定義對象」)。
名詞速查表(本節出現的核心術語)#
- Heaviside 函數: \(u(t)\) ;開關與延遲的積木。
- 時移律:延遲 \(a\) ⇔ 乘 \(e^{-as}\) 。
- Dirac 脈衝: \(\delta\) ;捲積單位元、取值機。
- 分布(廣義函數):檢驗函數空間上的連續線性泛函(Schwartz 1945)。
- 跳變條件: \(\delta\) 輸入下最高階導數跳變其餘連續。
- 脈衝響應: \(h = \mathcal{L}^{-1}\{G\}\) ;系統的完整指紋。
- 基本解: \(LG = \delta\) 的解;Green 函數的核心。
- 採樣定理:帶限信號的樣本完備性(Shannon 1949)。
- 白噪聲: \(\delta\) 相關的廣義過程。
- 波前集:奇異性位置與方向的精細記錄(Hörmander)。
年表速查(本節涉及的關鍵年份)#
- 1880s:Heaviside 運算中 \(\frac{d}{dt}u\) 型對象。
- 1926:Dirac 量子力學中系統使用 \(\delta\) 。
- 1932:Kirchhoff/波傳播基本解的經典理論整合。
- 1945:Schwartz 分布理論。
- 1949:Shannon 採樣定理。
- 1950:Schwartz 獲首屆菲爾茲獎。
- 1960s:Hörmander 微局部分析的先驅工作。
- 1970s:波前集理論完成。
深化問答三則#
問: \(\delta\) 輸入下「 \(y\) 連續、 \(y'\) 跳」,這是巧合嗎? 答:是力學的必然。二階方程 \(\ddot y + \cdots = \delta\) 中, \(\ddot y\) 必須含 \(\delta\) 才能平衡右端; \(y'\) 是 \(\ddot y\) 的積分,積分 \(\delta\) 得階躍(跳); \(y\) 是 \(y'\) 的積分,積分階躍得連續折線。每一次積分都把奇異性「洗淡一級」—— \(\delta\) → 跳 → 連續。這個「奇異性階梯」是分布理論的正則性階梯(Sobolev 空間)的初等形式。
問:脈衝響應與頻率響應(7.1 節)是什麼關係? 答:一對 Fourier/Laplace 變換: \(h(t) \leftrightarrow G(i\omega)\) 。時域敲擊(脈衝)與頻域掃頻(正弦穩態)是同一系統資訊的兩種讀取——測量學中分別叫時域法(modal impact test)與頻域法(sweep test),結果經 FFT 互換。理論上等價的方法,實務上因噪聲與儀器限制而各有擅場——數學等價性與工程最優性是兩個問題。
問:分布理論的「以作用定義對象」思想還出現在哪裡? 答:無處不在。(1)測度論:測度由其對可積函數的積分定義;(2)代數幾何:Grothendieck 概形由其點上的「函數環」定義;(3)量子力學:量子態由其對觀測量的期望值定義;(4)機器學習:核方法以核函數(相似度作用)定義隱空間。「對象 = 它如何作用於其他對象」——這是二十世紀數學最深的範式遷移之一, \(\delta\) 函數是其最直觀的教學案例。
(本節完)