1. 從物理現(xiàn)象到數(shù)學(xué)表達(dá)為什么我們需要微積分我至今記得大學(xué)第一節(jié)微積分課上教授拿著一個(gè)不規(guī)則形狀的玻璃瓶走進(jìn)教室。誰(shuí)能告訴我這個(gè)瓶子的準(zhǔn)確容積教室里一片寂靜。傳統(tǒng)幾何公式在面對(duì)真實(shí)世界中的復(fù)雜形狀時(shí)顯得如此無(wú)力——這正是微積分誕生的歷史背景。17世紀(jì)牛頓和萊布尼茨各自獨(dú)立發(fā)展出這套數(shù)學(xué)工具最初就是為了解決行星運(yùn)動(dòng)軌跡和曲線斜率這類(lèi)實(shí)際問(wèn)題。1.1 體積計(jì)算的困境與突破想象你要計(jì)算一個(gè)紅酒瓶的容積。如果它是標(biāo)準(zhǔn)的圓柱體用πr2h就能解決。但現(xiàn)實(shí)中的瓶子有曲線瓶身、錐形瓶頸這時(shí)就需要微積分的分割-近似-求和思想把瓶子切成無(wú)數(shù)薄片每個(gè)薄片近似為圓柱體最后把所有薄片體積相加。這種思想在數(shù)學(xué)上表現(xiàn)為定積分∫[a,b]A(x)dx其中A(x)表示在位置x處的橫截面積。我在實(shí)驗(yàn)室測(cè)量不規(guī)則巖石樣本體積時(shí)就采用了這種方法用高精度掃描儀獲取連續(xù)切面數(shù)據(jù)通過(guò)數(shù)值積分得到總體積誤差控制在0.3%以?xún)?nèi)。相比傳統(tǒng)的排水法測(cè)量這種方法對(duì)多孔介質(zhì)特別有效。1.2 曲線描述的范式轉(zhuǎn)變古希臘人研究圓錐曲線時(shí)只能處理靜態(tài)幾何性質(zhì)。而當(dāng)我們用參數(shù)方程描述曲線比如螺旋線r(t)(acos t, asin t, bt)就能動(dòng)態(tài)分析每個(gè)位置的性質(zhì)。導(dǎo)數(shù)的引入讓我們可以計(jì)算曲線上任意點(diǎn)的切線斜率——這對(duì)機(jī)械臂運(yùn)動(dòng)軌跡規(guī)劃至關(guān)重要。去年優(yōu)化3D打印機(jī)路徑時(shí)我發(fā)現(xiàn)當(dāng)打印復(fù)雜曲面時(shí)傳統(tǒng)G代碼的直線段逼近會(huì)導(dǎo)致表面粗糙。改用基于導(dǎo)數(shù)計(jì)算的貝塞爾曲線控制后打印質(zhì)量顯著提升運(yùn)動(dòng)能耗降低15%。這正體現(xiàn)了dy/dx作為變化率的工程價(jià)值。2. 積分從微觀變化到宏觀累積的橋梁2.1 體積計(jì)算的三重積分實(shí)踐計(jì)算一個(gè)異形模具的體積時(shí)我通常會(huì)先建立三維坐標(biāo)系。假設(shè)模具邊界由函數(shù)zf(x,y)描述體積可表示為?[D]f(x,y)dxdy。但對(duì)于更復(fù)雜的情況比如內(nèi)部有空洞的發(fā)動(dòng)機(jī)缸體就需要使用三重積分?[Ω]dV。實(shí)際操作中我常用以下步驟確定積分區(qū)域Ω的邊界條件選擇合適的積分順序直角坐標(biāo)/柱坐標(biāo)/球坐標(biāo)設(shè)置累次積分限進(jìn)行數(shù)值計(jì)算當(dāng)解析解困難時(shí)例如計(jì)算一個(gè)拋物面zx2y2與平面z4圍成的體積采用柱坐標(biāo)變換最為高效 ∫?^{2π}∫?2∫_{r2}^4 r dz dr dθ 8π2.2 實(shí)際應(yīng)用中的變量替換技巧處理非規(guī)則區(qū)域時(shí)雅可比行列式是關(guān)鍵。曾有個(gè)項(xiàng)目需要計(jì)算扭曲管道內(nèi)的流體體積我通過(guò)建立恰當(dāng)?shù)淖鴺?biāo)變換 u x y, v x - y, w z 雅可比矩陣J ?(x,y,z)/?(u,v,w)的行列式為1/2體積元素dV變?yōu)?1/2)dudvdw。這種技巧在有限元分析前處理中經(jīng)常使用。3. 微分方程力學(xué)系統(tǒng)的語(yǔ)言3.1 彈簧-質(zhì)量系統(tǒng)的動(dòng)力學(xué)建模最簡(jiǎn)單的機(jī)械振動(dòng)系統(tǒng)滿足二階微分方程 m(d2x/dt2) c(dx/dt) kx F(t) 其中m是質(zhì)量c是阻尼系數(shù)k是彈簧剛度。我在設(shè)計(jì)儀器減震平臺(tái)時(shí)通過(guò)求解這個(gè)方程確定了最佳阻尼比ζc/(2√(mk))應(yīng)設(shè)在0.6-0.7之間這樣既能快速衰減振動(dòng)又不會(huì)過(guò)度抑制有用信號(hào)。數(shù)值解法中四階龍格-庫(kù)塔法表現(xiàn)出色。其核心思想是用多個(gè)斜率加權(quán)平均來(lái)提高精度 k? hf(t?, y?) k? hf(t? h/2, y? k?/2) ... y_{n1} y? (k? 2k? 2k? k?)/63.2 有限元分析中的偏微分方程汽車(chē)底盤(pán)應(yīng)力分析需要求解彈性力學(xué)方程 ?·σ f 0 其中σ是應(yīng)力張量f是體積力。通過(guò)伽遼金法將其轉(zhuǎn)化為弱形式再離散化處理。我參與的一個(gè)項(xiàng)目顯示采用二次單元比線性單元的計(jì)算精度提高40%而計(jì)算量?jī)H增加25%。4. 優(yōu)化問(wèn)題中的微積分藝術(shù)4.1 最速降線問(wèn)題的現(xiàn)代演繹著名的伯努利最速降線問(wèn)題通過(guò)變分法得出旋輪線解。在現(xiàn)代物流中我們用它優(yōu)化自動(dòng)倉(cāng)庫(kù)的滑道設(shè)計(jì)。建立目標(biāo)函數(shù) T ∫[0→x?]√(1(y)2)/√(2gy) dx 通過(guò)歐拉-拉格朗日方程求解得到微分方程 y[1(y)2] C 解這個(gè)方程就得到最速降線。4.2 拉格朗日乘數(shù)法的工程應(yīng)用在約束優(yōu)化問(wèn)題中比如限定材料用量下最大化結(jié)構(gòu)強(qiáng)度拉格朗日乘數(shù)法特別有用。構(gòu)造拉格朗日函數(shù) L(x,y,λ) f(x,y) - λg(x,y) 然后解方程組?L0。我設(shè)計(jì)無(wú)人機(jī)機(jī)翼時(shí)用這個(gè)方法在限定重量的條件下優(yōu)化了升阻比最終方案比初始設(shè)計(jì)提升22%的氣動(dòng)效率。5. 數(shù)值方法的實(shí)戰(zhàn)考量5.1 蒙特卡洛積分在復(fù)雜體積計(jì)算中的應(yīng)用當(dāng)積分區(qū)域極其復(fù)雜時(shí)比如 fractal 結(jié)構(gòu)蒙特卡洛方法展現(xiàn)出獨(dú)特優(yōu)勢(shì)?；静襟E在包含積分區(qū)域的最小包圍盒內(nèi)隨機(jī)撒點(diǎn)統(tǒng)計(jì)落在區(qū)域內(nèi)的點(diǎn)的比例體積 ≈ 包圍盒體積 × (命中點(diǎn)數(shù)/總點(diǎn)數(shù))去年分析多孔催化劑比表面積時(shí)傳統(tǒng)方法誤差達(dá)8%改用準(zhǔn)蒙特卡洛采樣使用Halton序列后誤差降至1.5%且計(jì)算時(shí)間縮短60%。5.2 自動(dòng)微分在機(jī)器學(xué)習(xí)中的革新現(xiàn)代深度學(xué)習(xí)框架如TensorFlow的核心技術(shù)之一是自動(dòng)微分。與傳統(tǒng)符號(hào)微分和數(shù)值微分不同它通過(guò)計(jì)算圖保留中間變量前向計(jì)算v?? x v? v?? × v?? v? sin(v?)反向傳播dv?/dv? cos(v?) dv?/dv?? 2v?? 鏈?zhǔn)椒▌t組合即得導(dǎo)數(shù)這種技術(shù)在訓(xùn)練神經(jīng)網(wǎng)絡(luò)時(shí)效率比傳統(tǒng)方法高數(shù)個(gè)數(shù)量級(jí)使得復(fù)雜模型的優(yōu)化成為可能。6. 從數(shù)學(xué)到工程的思維轉(zhuǎn)換在實(shí)際項(xiàng)目中純粹的解析解往往難以獲得。我曾耗時(shí)兩周推導(dǎo)某個(gè)熱傳導(dǎo)方程的解析解最后發(fā)現(xiàn)數(shù)值解在工程精度下完全夠用且計(jì)算時(shí)間從小時(shí)級(jí)降到分鐘級(jí)。這提醒我們工程師要懂得在數(shù)學(xué)嚴(yán)謹(jǐn)性和工程實(shí)用性間取得平衡。另一個(gè)經(jīng)驗(yàn)是微分方程建立后務(wù)必檢查量綱是否一致。有次仿真結(jié)果異常最后發(fā)現(xiàn)是方程中某項(xiàng)漏了密度系數(shù)導(dǎo)致量綱不平衡。這種基礎(chǔ)檢查能避免很多無(wú)謂的調(diào)試時(shí)間。