讀博后降維打擊數(shù)學建模！

文章來自知乎提問：如何評價 2021 年美國大學生數(shù)學建模競賽？今年賽情如何？

來湊個熱鬧降維打擊美賽了，寫一寫前兩題的思路。2018年讀大二的時候拿了一次美賽H，現(xiàn)在已經(jīng)讀PhD了，再看當時做的跟屎一樣。看了幾個回答發(fā)現(xiàn)很多人甚至對題目的理解都是錯的，私以為比賽最重要的是解決問題，而不是堆砌一堆花里胡哨的很Fancy的方法。

A題

第一問：描述 breakdown of ground litter and woody fibers 在 ?fungal activity in multiple species of fungi 下的情況

要解決這個問題首先要搞明白用什么來量化分解情況？有什么影響分解的情況？fungal activity真菌活動是什么？多種真菌到底是哪幾種？

介于我們不是養(yǎng)蘑菇專家，我們能獲取的信息基本都來源于A trait-based understanding of wood decomposition by fungi 這篇文章，下載下來看看有哪些信息。雖然我們不能抄襲論文里面的方法，但可以幫助我們了解一些基本概念。

在Introduction里面我們發(fā)現(xiàn)一句話：

Models traditionally used microbial biomass as a proxy of decomposer activity (3, 6), treating the microbial community as a single homogeneous group or a small number of functionally distinct pools (7).

傳統(tǒng)方法一般使用 microbial biomass 生物量來代表分解者活動，這個問題里面顯然fungi就是decomposer， 所以fungi的量應該就代表這個問題里面的真菌活動狀況。閱讀文章后，我們發(fā)現(xiàn)實驗中的hyphal density(菌絲密度) 應該就是對應這個生物量。

To characterize the growth of each isolate, we quantified its hyphal extension rate and hyphal density (presented in ref. 34)

那么用什么來量化分解情況呢？論文里面是這樣說的：

Here, we explore which fungal characteristics predict wood decomposition rate across a range of common wood decomposer fungi.

就是wood decomposition rate。具體的，這個根據(jù)木塊被真菌吃掉的數(shù)量來計算的。

we measured the mass loss of wood blocks when colonized by each fungus to estimate a standardized wood decomposition rate

下一個問題是多種真菌到底有哪幾種? 從第二題要求往模型里面加入Fungi不同種類之間的相互作用，我們可以大體推測第一問不需要考慮多種fungi的內(nèi)部作用。

OK 第一問我們的思路已經(jīng)清楚了，題目想讓我們建立一個數(shù)學模型（回歸模型）來刻畫自變量（真菌的量fungi biomass）和因變量(分解率wood decomposition rate)的關系，這就是傳統(tǒng)方法的做法。我覺得在下面兩個文章有可能找到關于兩個變量之間關系的一些基本假設。

K. L. McGuire, K. K. Treseder, Microbial communities and their relevance for ecosystem models: Decomposition as a case study. Soil Biol. Biochem. 42, 529–535 (2010)

S. Manzoni, A. Porporato, Soil carbon and nitrogen mineralization: Theory and models across scales. Soil Biol. Biochem. 41, 1355–1379 (2009).

第二問：考慮不同真菌間的相互作用，在模型里面加入Fig1和Fig2的信息

首先搞明白上一問中的遺留問題，真菌到底有哪些類型? 在文章中我們發(fā)現(xiàn)：

First, we use a database of 22 fungal traits previously measured in each of 34 wood rot fungi collected from a wide geographic range across North America (25, 34–36) to identify potential drivers of wood decomposition.

這里作者用了22種真菌特性進行研究，那真菌特性是什么東西？為什么里面沒有出現(xiàn)真菌的類型？我們繼續(xù)從論文中探索，發(fā)現(xiàn)這么一句話：

However, in recent years, the development of traitbased approaches has begun to transform our understanding of broad-scale functional patterns. By linking traits to ecosystem functioning, these approaches have been used in plant and animal ecology to infer the functioning of novel communities without prior knowledge of the taxa that are present.

大概就是現(xiàn)在一種新的 traitbased 研究方法，在僅知道一群真菌的特性的情況下（不需要考慮他們的種類），就可以代替知道真菌的種類進行研究。所以我們這一題的核心或許是不考慮具體的species，而是考慮一堆具有相同traits的不同類真菌（例如長得快的菌，抗凍的菌……等等），用他們的common trait作為種類的依據(jù)。

那么這些traits又是啥呢？從論文里我們發(fā)現(xiàn)這22種(2+11+9)traits可以分為三大類:

Hyphal extension rate and hyphal density reflect hyphal morphology and growth strategy,? ecological performance traits (11 in total) relate to combative ability.tolerance of a range of temperature and moisture conditions (25),? and finally, the production of oxidative and hydrolytic enzymes (9 traits) promotes nutrient acquisition from organic resources (26).

但要注意問題里面說了，我們只關注growth rate和moisture tolerance就好：

For this MCM Problem you should focus on just two traits of a fungus: the growth rate of the fungus and the fungus’ tolerance to moisture.

OK 那么問題明確了，這里我們需要加入growth rate和moisture tolerance到我們在第一問中建立的一元回歸模型中，把它變成一個多元回歸模型，相當于是在原來的傳統(tǒng)方法上，引入這篇論文里的方法，因此需要使用這篇文章的數(shù)據(jù)：

https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.1909166117/-/DCSupplemental

要注意你最后模型的結果，要對比一下原始論文里的Fig.3，看跟這個圖是不是一致。

第三問，分析模型，并刻畫不同種類真菌之間長期與短期的關系

第二問建完模型之后，分析下幾個自變量之間的關系（相關性和獨立性），看一看大概有什么關系，根據(jù)這些關系做一些統(tǒng)計分析。例如是否發(fā)現(xiàn)自變量之間存在此消彼長的情況，可以用不同菌之間的競爭來解釋等等……

重要的是注意這幾點：

• 長期與短期的動態(tài)關系，要考慮他們的變化趨勢（可以參考時間序列中提取趨勢的方法）。考慮他們之間的競爭關系，作者進行了pair wise的競爭對比在這篇文了的Combative ability數(shù)據(jù)里面。
• 檢驗 rapid fluctuations 對各個變量影響的sensitivity （引入噪聲進行假設檢驗）
• 分析氣候變化帶來的影響。此時再加入自變量溫度、濕度等指標。要用到這篇論文里的Temperature and moisture niche數(shù)據(jù)。

……

后面兩問也需要用到前三問建立的回歸模型做一些預測。把基本概念搞清楚以后，思路就有了。

B題

B題看了半天，直接看不懂，先說說我對于題目的理解。

題目中出現(xiàn)了幾種對象：

• EOC：理解為救援指揮中心，應該是一些固定位置的站點，但是Google沒搜到位置，也不知道是一個還是好多個。
• SSA drones：帶攝像頭和傳感器的無人機，用來監(jiān)視救火隊的，理解為在救火隊附近使用的無人機（理解為0距離，攝像頭一般看不了很遠）
• Radio Repeater drones：帶中繼器的無人機，用來幫助EOC和救火隊交流的無人機，可以離救火隊遠一點（20KM）
• Firefighters ：救火隊

我搜了搜SSA是啥，下面有個解釋。情景感知是happening around you, at or near your physical location，應該是觀測很近的那種。我覺得就像是美國特戰(zhàn)隊那種，自己拿著的小型無人機。

Situational Awareness is a concept closely involved with physical security information management (PSIM, see the white papers on this subject). It is usually defined as being aware of what is happening around you, at or near your physical location, or at some other location where you are supposed to be or where there are assets or people that you must protect.?

ref: https://www.titan-vision.com/pages/hot-topics/total-situational-awareness.html

幾個很疑惑的問題：

• 兩種無人機之間有沒有交流？在問題里面沒看到。
• 無人機從哪里開始飛? 是救火隊拿著到時候飛，還是從EOC飛過去？我覺得應該是救火隊拿著飛，而且從下面的飛機懸停也可以大概推測。

大概推測一下情景是這樣的，K個帶中繼器的無人機懸停在距離EOC固定距離(20km)的圓內(nèi)，保證加上救火隊的交流半徑（5km平地，2km市區(qū)）能夠跟EOC交流。因為每個救火隊拿著一個SSA無人機，所以直接把救火隊小組抽象為SSA，SSA是需要跟著救火隊移動的。

你要買很多SSA無人機，以及少量的帶中繼器的無人機。問題說帶中繼器和遙感的無人機一萬美金很貴，應該是那種大型的用來運輸東西的，SSA的小型無人機用來監(jiān)視的帶攝像頭那種，類似大疆的，Amazon搜搜應該很便宜？

Akme Corporation's prototype WileE–15.2X hybrid drone is projected to cost approximately $10,000 (AUD) when equipped with either a radio repeater or video & telemetry capability.

損失函數(shù)可以定義一個 【失聯(lián)的SSA（救火隊）數(shù)量】【沒有及時被救援的火場面積】 與 【買兩種無人機花的錢】的加權，正好體現(xiàn)題目中要求的：

Your model should 【balance capability and safety】 with 【economics】,

你可以改變這三個權重去做trade-off，看看情況是什么樣。

然后我想到的約束有下面幾個：

• 救火隊（SSA）的距離（5km平地，2km市區(qū)）+ ?中繼器的距離(20km) 小于與EOC的距離
• 每個火點的多少距離內(nèi)需要至少1個救火隊
• (SSA)救火隊（SSA）趕到火點的時間需要小于K分鐘

這些約束不需要完全滿足，用拉格朗日松弛一下去求一個下界。

關于是不是SSA和中繼器都用一種無人機，我傾向于不是使用一種無人機，否則第一問balance經(jīng)濟和安全就沒法做了，只需要關注安全就好了。B題的思路是基于我理解正確的情況下，如果我理解的不對，那么思路就是不對的，大家參考下吧。

第二問我的理解是，去隨機生成一些極端的大火情況，然后看第一問種哪些約束被打破了。根據(jù)那些打破的約束，看應當如何進行調(diào)整。

……

我知道美賽對于沒有經(jīng)過系統(tǒng)科研訓練的本科生很難很頭痛，但這也是這類比賽的意義所在，能夠在三天的時間內(nèi)專心獨自解決一個實際問題，是很痛苦但結束后會覺得很享受的事情。貼一句我很喜歡的爆裂鼓手的一句話："If you want the fucking core part, earn it! "