引言

熟悉林丹和李宗偉的都知道，羽毛球是一項激（基）情四射的體育運動。2020 年 7 月 4 日，世界羽壇名將林丹正式宣布退役。

林丹是羽毛球史上第一位集奧運會、世錦賽、世界杯、蘇迪曼杯、湯姆斯杯、亞運會、亞錦賽、全英賽、全運會及多座世界羽聯(lián)超級系列賽冠軍于一身的雙圈全滿貫選手，被譽為中國羽球一哥，是世界羽毛球最知名的明星球員。盡管如此，林丹可能仍然看不懂本文。

要是說“林丹是地球上最會打羽毛球的男人”，應(yīng)該沒有人反對擅長物理的韓同學可能會第一個不服！

圖 3: 擅長裝逼的學霸：打羽毛球基本靠心算

已經(jīng)被保入清華物理系的韓同學表示：羽毛球運動是一個比較簡單的拋體模型，當對手擊球后，通過觀察羽毛球的速度和角度，并用手去感受風速，再考慮空氣動力學的修正，通過快速的心算就可以大概判斷出羽毛球的落點，再選擇合適的動作回擊。盡管我感覺到韓同學可能在裝逼，但機器人確實得靠芯算才能打羽毛球。

圖 4: 確實得靠芯算才能打羽毛球的機器人

實際上，羽毛球的運動并不簡單。由于羽毛球獨特的結(jié)構(gòu)（圓錐形，頭重尾輕），羽毛球在飛行過程中會出現(xiàn)翻轉(zhuǎn)、振蕩和旋轉(zhuǎn)等獨特的行為。例如，運動員在使用球拍拍擊球頭后，羽毛球都要翻轉(zhuǎn)后繼續(xù)向前飛行（圖 5）。本文將建立羽毛球飛行的動力學模型 [1]，研究羽毛球運動的特征，為廣大羽毛球友科學地裝逼提供指南。

圖 5: 羽毛球被球拍擊中后迅速翻轉(zhuǎn)

簡史

在建立模型之前，本文先介紹一下羽毛球的歷史。早在兩千多年前，一種類似羽毛球運動的游戲就在中國、印度等國出現(xiàn)。中國叫打手毽，印度叫浦那，西歐等國則叫做毽子板球。十九世紀七十年代，英國軍人將在印度學到的浦那游戲帶回國，作為茶余飯后的消遣娛樂活動。在 14 世紀末，日本出現(xiàn)了把櫻桃插上美麗的羽毛當球，兩人用木板來回對打的運動。這便是羽毛球運動的原形。

圖 6: 古代中外與毽子和羽毛球相關(guān)畫作 [2-4]

現(xiàn)代羽毛球運動誕生在英國。1873年，在英國格拉斯哥郡的伯明頓鎮(zhèn)有一位叫鮑弗特的公爵（dukes of Beaufort），在他的領(lǐng)地開游園會，有幾個從印度回來的退役軍官就向大家介紹了一種隔網(wǎng)用拍子來回擊打毽球的游戲，人們對此產(chǎn)生了很大的興趣。因這項活動極富趣味性，很快就在上層社會社交場上風行開來。“伯明頓”（Badminton）即成為羽毛球的英文名字。羽毛球于 1992 年巴塞羅那奧運會上被列為正式比賽項目。裝逼指南：要熟知以上羽毛球的起源、歷史和文化，這可以作為裝逼時的談資。比如當見到心儀的女球友拿著羽毛球時，你可以這樣夸贊她：“你拿羽毛球的樣子真美，如果法國的夏爾丹在世，一定會重繪一幅《拿羽毛球的姑娘》（圖 6 右）”。記住，在說“夏爾丹（Chardin）”時一定要用純正的法語口音。

圖 7: 現(xiàn)代羽毛球運動場地、球拍和羽毛球

現(xiàn)代羽毛球運動通常是兩人（單打）或四人（雙打）的球場運動。標準球場長 13.4 m，寬 5.2 m，被高 1.55 m 的網(wǎng)分隔成兩部分（圖 7 左）。比賽時，運動員站在各自的球場上，通過球拍擊打羽毛球，從而使羽毛球越過網(wǎng)并落入對手的半場。羽毛球過網(wǎng)之前，每側(cè)只能擊打一次。一旦羽毛球落地或運動員犯規(guī)，則該球結(jié)束并計分。羽毛球比賽通常采用“21 分制，3局2勝”，一場完整的羽毛球賽持續(xù) 1 小時左右，每局持續(xù) 20 分鐘左右，每球持續(xù) 10 秒左右。如圖 7 右所示，羽毛球由 16 根羽毛構(gòu)成的球裙和軟木球頭組成，其質(zhì)量 m = 5.3 g 左右，長度 L = 10 cm 左右，直徑 D = 68 mm 左右，羽毛夾角 $Λ=45 度左右 [5]。裝逼指南：要記住羽毛球場和羽毛球的這些參數(shù)，在適當?shù)臅r候準確的報出這些數(shù)字，能瞬間提高逼格，引起異性球友的關(guān)注。

實驗

大量文獻從實驗、理論和數(shù)值角度對羽毛球的運動軌跡進行了廣泛的研究 [6]。圖 8 是一次實驗中拍攝到的羽毛球運動軌跡，快照間隔為 50 ms，羽毛的初始速度為 58 m/s，出球角度（初速度與水平方向的夾角）為 52 。

圖 8: 羽毛球運動軌跡 [6]

從圖中可以看出，開始羽毛球幾乎以直線飛行，并快速減速，在軌跡頂部速度最小，然后在重力作用重新加速，最終以恒定速度 6.7 m/s 垂直下落。顯然這并不是一條伽利略拋物線：羽毛球的軌跡并不具有對稱性，其射程（圖中顯示約為 9 m）遠小于預期（以相同速度和出球角度，按照理想拋物線計算約為 240 m）。這種不對稱性和射程的減小是由空氣阻力引起的。

也有少量文獻研究了羽毛球在飛行過程中的翻轉(zhuǎn) [1]。圖 9 顯示的是一次實驗中羽毛球與球拍碰撞后，球軸與速度方向的夾角 φ 隨時間的演變。快照間隔為 5 ms，初始速度 v0 = 18.6 m/s，初始翻轉(zhuǎn)角速度為 w = 206 rad/s。

圖 9: 羽毛球離拍后的角度與速度夾角 [1]

從圖中可以看出，與球拍接觸后（約 1 ms 的時間）, 羽毛球經(jīng)過 tf = 15 ms 的時間， 球軸方向發(fā)生 180 的翻轉(zhuǎn)。然后，球軸以速度方向為平衡位置進行阻尼振蕩，直到球軸與速度方向一致為止。振蕩時間 to = 85 ms。 經(jīng)過約 ts = 130 ms 后, 球軸與速度方向一致。

模型

為了分析羽毛球的運動特性，需要對作用在羽毛球上的力加以分析。羽毛球被球拍擊中后，主要受到重力 mg 和空氣阻力的作用。在引入空氣阻力之前，我們首先介紹一下流體力學最重要的概念之一，雷諾數(shù)：

其中 ρ表示流體（空氣）密度，μ表示流體動力粘度， v是物體（羽毛球）相對流體的速度，D 是物體的特征長度。雷諾數(shù)實際上是流體的慣性力與粘性力比值的量度，雷諾數(shù)較小時，粘性力對流場的影響大于慣性力，流場中流速的擾動會因粘性力而衰減，流體流動穩(wěn)定，為層流；反之，若雷諾數(shù)較大時，慣性力對流場的影響大于粘性力，流體流動較不穩(wěn)定，流速的微小變化容易發(fā)展、增強，形成紊亂、不規(guī)則的紊流流場。在羽毛球運動中，雷諾數(shù)大約在 10000 到 200000 之間，屬于典型的湍流區(qū)。若假設(shè)羽毛球在運動的過程中，迎風截面積不變，則空氣的阻力可以表示為

上式中右側(cè)負號表示阻力方向與相對速度方向相反。 ρ為空氣密度， S 為羽毛球的截面積，v為羽毛球的速度矢量， |v| 為速度大小， Cd 是阻力系數(shù)。實驗表明，在上述雷諾數(shù)范圍內(nèi)， Cd 幾乎是常數(shù)，大概在 0.6 到 0.7 之間 [7,8]。上式表明，物體在空氣（流體）中受到的阻力大小正比于迎風橫截面積和相對速度的平方。這很容易理解，生活中有很多例子，例如騎自行車速度越大感受到的風阻也越大，專業(yè)自行車選手會通過壓低身體來減小迎風截面積。

圖 10: 騎自行車受到空氣的阻力與車速，迎風橫截面積相關(guān)

軌跡

為了像學霸一樣能夠判斷出羽毛球的落點，我們需要討論羽毛球的宏觀軌跡。從圖 9 的實驗中可知，羽毛球在被擊出之后，其朝向會迅速掉頭，經(jīng)過很短的震蕩過程后朝向會與速度方向保持一致。因此，其迎風的截面積在絕大部分時間內(nèi)是不變的。通過上文空氣阻力的分析，可知羽毛球的質(zhì)心運動方程為

由 dv/dt =0 可以確定出羽毛球的最終下落速度

如果將羽毛球的運動限定在 x-z 平面中，則羽毛球的質(zhì)心運動方程可改寫為

其中點點x 、點 x 分別表示水平方向坐標 x 對時間的一階導數(shù)（速度）和二階導數(shù)（加速度）。應(yīng)用 MATLAB 中的 ode45 函數(shù)對以上微分方程組進行數(shù)值積分可求得羽毛球的軌跡。在本文的計算中，我們?nèi)∮鹈虻馁|(zhì)量 m = 5.3 g，截面積 S = 36 cm^2，阻力系數(shù) Cd = 0.63，空氣密度ρ=1.2 kg/m^3。圖 11 為實驗（點）和計算（線）得到的羽毛球運動軌跡。

圖 11: 實驗和計算得到的羽毛球運動軌跡

結(jié)果表明，計算給出的不同初速度和出球角度的軌跡曲線與實驗拍攝到的軌跡點非常吻合。這也證實了在研究羽毛球的軌跡時，可以忽略羽毛球在運動的過程中迎風截面積的變化。計算所使用的程序見附錄。

此外，我們通過軌跡的計算，還可以獲得不同速度大小和出球角度的羽毛球射程。計算結(jié)果如圖 12 所示，所使用的程序見附錄。

圖 12: 羽毛球射程與初速度和角度的關(guān)系

圖中每條曲線上的最高點表示該初速度下最大射程對應(yīng)的出球角度（最大射程角）?？梢宰C明，在忽略空氣阻力的情況下，出球角度為 45 度時射程最大。而有空氣阻力時，初速度越小，羽毛球受空氣阻力影響越小，最大射程角越接近于 45 度。初速度越大，最大射程角越小。羽毛球的最高飛行速度可達 170 m/s [9]，容易算出該球速下的最大射程約為 15.5 m。因此，羽毛球不太容易飛出長度為 13.4 m 的羽毛球場。這也是羽毛球平均每回合擊球數(shù)高達 13.5 次的原因。而網(wǎng)球平均每回合擊球數(shù)只有 3.5 次 [10]。網(wǎng)球的最高射程高達 67 m，這遠大于只有 24 m 長的網(wǎng)球場。裝逼指南：這提醒我們裝逼要注意科學性，千萬不要說出類似“我一拍能打出 16 米”這樣有失逼準的話。羽毛球的質(zhì)心運動方程具有解析解 [11]，該解析解可以給出近似的射程

其中L= 2m/(ρS Cd) 為空氣動力學特征長度。上式表明，羽毛球的射程取決于初速度 v0 和角度θ 的大小。我猜測學霸應(yīng)該就是通過對類似公式進行心算來估計羽毛球落點的。從質(zhì)心運動方程可以看出，羽毛球的軌跡和射程還與空氣的密度相關(guān)。而空氣的密度又取決于溫度 T 和相對濕度 δ [12]：

其中空氣相對濕度 0

圖 13: 羽毛的微觀網(wǎng)狀結(jié)構(gòu)

實驗表明，羽毛球質(zhì)量（單位：g）隨空氣相對濕度線性增加：

考慮北方的冬天（假設(shè)溫度為 -10，相對濕度 30%）和南方的夏天（假設(shè)溫度為 30，相對濕度 90 %）羽毛球的飛行軌跡，計算結(jié)果如圖 14，所使用的程序見附錄。

圖 14: 氣溫和濕度對羽毛球軌跡的影響

北方冬天和南方夏天空氣密度可相差 15%，并最終導致羽毛的軌跡和射程有明顯差異。因此，運動員在不同季節(jié)和不同地區(qū)打球，可不單是身體適應(yīng)了溫度和濕度的變化就足夠了的。這個結(jié)果同樣也解釋了為何部分選手在打羽毛球之前需要蒸球。裝逼指南：提前查看天氣預報并記住溫度和濕度，根據(jù)公式計算出空氣密度和羽毛球質(zhì)量。在打球時，就可以“指導”迷妹們?nèi)绾斡煤粑惺軡穸?、皮膚感受溫度，來修正羽毛球的軌跡和落點。

翻轉(zhuǎn)

羽毛球在飛行的過程中，還會出現(xiàn)翻轉(zhuǎn)和振蕩微觀行為。從圖 5 和 9 的實驗中可以看出，迎面飛向運動員的羽毛球都是球頭朝前。當運動員使用球拍擊中羽毛球球頭后，球軸朝向會發(fā)生 180 度的翻轉(zhuǎn)。

圖 15: 不同角度的羽毛球受力示意圖

羽毛球的這種翻轉(zhuǎn)的行為其實是羽毛球的獨特結(jié)構(gòu)導致的。如圖 15 所示，羽毛球的質(zhì)量主要集中于球頭，而球裙的羽毛卻貢獻了大部分截面積。因此，空氣阻力的作用中心和羽毛球質(zhì)心顯然不會重合，這會對羽毛球形成一個凈力矩，使羽毛球頭部朝向其運動方向轉(zhuǎn)動。盡管圖 15d 所示的情況也屬于平衡態(tài)，但這種平衡是不穩(wěn)定平衡，極小的氣流擾動就能夠迅速的破壞這一平衡，使之朝著圖 15b 或 c 的情況轉(zhuǎn)變，并最終變?yōu)榉€(wěn)定平衡狀態(tài)（圖 15a）。

圖 16: 羽毛球受力分析圖

圖 16 左是對羽毛球的受力分析，其飛行速度為v，重心為 G， 空氣阻力的作用中心為 P，球軸與速度夾角為 φ。為了方便分析，我們將羽毛球簡化為圖 16 右所示的兩個圓球。 B 球表示球裙，其截面積為 SB，質(zhì)量為 mB。 C 球表示球頭，其截面積為 SC，質(zhì)量為 mC。由角動量定理可以推導（過程見附錄）出羽毛球以速度方向為平衡位置的振動方程：

考慮到羽毛球重而小的球頭和大而輕的球裙，即 SB*mC >> SC*mB。以及由質(zhì)心的定義可知 mB L_=m_C L_。并假設(shè) vB 和 vC 約等于 v，則上式可化簡為

其中 w0^2 = ρSB Cd v^2/ 2m L_，tau_s = 2m_B(1+m_B/m_c)/\rho S_B C_d v。上式是一個阻尼振蕩方程，描述了羽毛球以運動方向為平衡位置的振動。角頻率的平方項 w0^2 對應(yīng)于空氣阻力產(chǎn)生的穩(wěn)定力矩，而阻尼項 1/\tau_s 則來源于角度變化產(chǎn)生的垂直羽毛球運動方向的速度對應(yīng)的阻力。應(yīng)用 MATLAB 中的 ode45 函數(shù)可對以上振動方程和質(zhì)心運動方程進行耦合求解，圖 17 是通過計算得到的羽毛球自由下落過程，模擬所使用的程序見附錄。

圖 17: 羽毛球自由下落的計算機模擬

根據(jù)振動方程，我們可以估算出羽毛球運動的一些特征時間。實驗表明（圖 9），羽毛球翻轉(zhuǎn)的時間小于穩(wěn)定時間。因此在考慮翻轉(zhuǎn)時間時，可以忽略上式中的阻尼項?？紤]到初始條件：翻轉(zhuǎn)角度和角速度分別為 π 和 點φ0，對上式進行積分有

則羽毛球翻轉(zhuǎn)時間可表示為

羽毛球振蕩時間可表示為

羽毛球的穩(wěn)定時間為

我們?nèi)?shù) mC=3.0 g、mB=2.3 g、SB = 36 cm^2、Cd = 0.63，L_=2 cm 和 ρ=1.2 kg/m^3，從而可以計算出 w0=66.6。代入上式可得羽毛球翻轉(zhuǎn)時間 tf=14 ms，振蕩時間 to=94 ms 和穩(wěn)定時間 ts=127 ms。這和圖 9 實驗給出的翻轉(zhuǎn)時間 15 ms 左右、振蕩時間 85 ms 左右和穩(wěn)定時間 130 mm 左右非常吻合。

此外，初始條件中的初始速度 v0 和初始翻轉(zhuǎn)角速度兩者并不獨立。圖 18 表明羽毛球初始翻轉(zhuǎn)角速度和球裙長度之積 與羽毛球的初速度 之間存在非常強的線性關(guān)系。換句話說，羽毛球的初速度越大，其初始翻轉(zhuǎn)角速度也越大。

圖 18: 初始角速度與初始速度的關(guān)系

前面討論了羽毛球在受到球拍擊打后的翻轉(zhuǎn)。在大多情況下，這種翻轉(zhuǎn)在對方球員接球之前就已經(jīng)完成。只有當羽毛球的穩(wěn)定時間與總飛行時間相當時，羽毛球的翻轉(zhuǎn)才會對對方球員造成影響。假定羽毛球在 1 m 的豎直高度上以 15 度角擊出，圖 19 給出了該種情況下羽毛球的穩(wěn)定時間與總飛行時間之比ts/t0 隨羽毛球射程 x0/Lfield的變化（計算程序見附錄）。

圖 19: 羽毛球的穩(wěn)定時間與射程關(guān)系

其中射程被除了球場長度 Lfield = 13.4 m 來進行歸一化。從圖 19 可以看出，只要水平身程大于 1.7 m（相應(yīng)的初速度大于 3 m/s），羽毛球穩(wěn)定時間與總飛行時間之比就會迅速下降。即對方球員在接球之前，羽毛球的朝向就已經(jīng)基本與飛行方向一致。這同樣也意味著，對于近網(wǎng)球，你不能簡單的假定自己在回擊時，羽毛球的球頭是正對著自己飛過來的。

圖 20: 近網(wǎng)球情況下,羽毛球的翻轉(zhuǎn)

另外，穩(wěn)定時間的計算是建立在羽毛球被正常擊打的前提下的，那么搓球無疑會進一步增大翻轉(zhuǎn)對于羽毛球運動的影響。裝逼指南：可以通過控制球速（通常出球速度小于 3 m/s）來延長羽毛球的翻轉(zhuǎn)、振蕩和穩(wěn)定時間，讓對手無法擊中羽毛球球頭。

旋轉(zhuǎn)

羽毛球除了會翻轉(zhuǎn)之外，也會繞著中心軸旋轉(zhuǎn)。如圖 21 所示，各片羽毛并不是平攤成一個正多邊形，而是以一個小角度 β相互疊在一起形成一個螺旋形。

圖 21: 穿過羽毛球的流線圖及旋轉(zhuǎn)示意圖

同風車、電風扇葉片等類似的結(jié)構(gòu)一樣，當空氣相對于羽毛球向后流動時，空氣會在垂直羽片方向施加作用力（圖 21）帶動羽毛球旋轉(zhuǎn)。平衡時，羽毛球以一個穩(wěn)定的角速度旋轉(zhuǎn)。該角速度下空氣施加在羽毛上的推進力矩和摩擦力矩相等，即

其中 Sp 為羽片面積。實驗表明 RΩ

這表明羽毛球的旋轉(zhuǎn)角速度 Ω 與飛行速度 v 存在線性關(guān)系。如圖 22 所示，實驗數(shù)據(jù)的擬合驗證了這種線性關(guān)系的存在。

圖 22:羽毛球旋轉(zhuǎn)角速度與飛行速度的關(guān)系

這種旋轉(zhuǎn)會對羽毛球的運動產(chǎn)生什么樣的影響？是否會因為旋轉(zhuǎn)，離心力將羽毛向外甩開，造成羽毛球有效截面積的增大？實際上，由于風阻的作用，羽毛還會受到一個向內(nèi)的壓力。風洞實驗表明 [1]，無論羽毛球是否旋轉(zhuǎn)，其阻力系數(shù)基本保持不變，即旋轉(zhuǎn)對阻力系數(shù)沒有顯?影響。此外，羽毛球的旋轉(zhuǎn)是否會產(chǎn)生明顯的進動，并且有助于其穩(wěn)定朝向呢？進一步研究證明，羽毛球繞對稱軸的轉(zhuǎn)動慣量太低，以至于轉(zhuǎn)動與平動速度比 ΩR/v 過小而根本起不到穩(wěn)定朝向的作用。這種效應(yīng)只在羽毛球的翻轉(zhuǎn)運動中才會有一定的貢獻。裝逼指南：羽毛球的旋轉(zhuǎn)對羽毛球的運動沒什么明顯的影響，但由于旋轉(zhuǎn)與球速存在正比關(guān)系，對于迎面飛來不太方便觀測速度的球，可以通過球的旋轉(zhuǎn)反推出球速，從而選擇合適的反擊策略。

結(jié)果

前文中，我們建立起了羽毛球的質(zhì)心運動方程、球軸的振動方程以及旋轉(zhuǎn)與飛行速度的關(guān)系。接下來，我們將應(yīng)用 MATLAB 中的 ode45 函數(shù)對質(zhì)心運動方程和球軸振動方程進行耦合求解，由于羽毛球的旋轉(zhuǎn)對其飛行沒有顯?影響，這里不再考慮。圖 23 和 24 是初速度為 50 m/s、初速度角為 50 度的羽毛球運動的模擬結(jié)果。

圖 23: 羽毛球運動的模擬

圖 24: 羽毛球翻轉(zhuǎn)角度隨時間的變化

模擬所使用的參數(shù)為：球頭和球裙質(zhì)量分別為 mC=3.3 g 和 mB=2 g，直徑分別為 dC=3.2 cm 和 dB=6.8 cm，羽毛球質(zhì)心與球頭質(zhì)心的距離 L_=2 cm。模擬所使用的程序見附錄。模擬結(jié)果顯示，羽毛球在飛行的初始階段翻轉(zhuǎn)和振蕩得比較厲害，而后逐漸趨于穩(wěn)定，后期羽毛球朝向基本與飛行速度一致。這和我們之前的分析和實驗觀測一致。裝逼指南：打羽毛球時帶上筆計本電腦，每打完一局休息階段，可用本文提供的模型和程序?qū)Ρ荣愔械拿恳磺蜻M行復盤。必要時，可向心儀的異性球友演示模擬結(jié)果，并專注地看著模擬動畫說出類似這樣的話：果然和我心算的結(jié)果一樣，剛才那球要是出球角再增加 5 度，出球速度再下調(diào) 2.3 m/s 就更完美了。

結(jié)論

羽毛球的動力學及其對比賽的影響一直倍受關(guān)注。由于羽毛球獨特的結(jié)構(gòu)（圓錐形，頭重尾輕），使得空氣阻力的作用中心和羽毛球質(zhì)心并不重合。這導致了羽毛球在飛行過程中會出現(xiàn)翻轉(zhuǎn)、振蕩和旋轉(zhuǎn)等獨特的行為。本文通過理論和模擬對羽毛球獨特的行為進行了研究，并為科學的裝逼提供了寶貴的意見：

- 要熟知羽毛球的起源、歷史和文化，這可以作為裝逼時的談資。

- 要記住羽毛球場和羽毛球的參數(shù)，在適當?shù)臅r候準確的報出這些數(shù)字。

- 裝逼要注意科學性，千萬不要說出違反物理規(guī)律、有失逼準的話。

- 提前查好參數(shù)，做好計算。打球時，就可以演示如何用呼吸感受濕度、皮膚感受溫度，來修正羽毛球的軌跡和落點。

- 可以通過控制球速來延長羽毛球的翻轉(zhuǎn)、振蕩和穩(wěn)定時間，來讓對手無法擊中羽毛球球頭。

- 可通過羽毛球的旋轉(zhuǎn)來反推出球速，從而選擇合適的反擊策略。

- 可用本文提供的模型和程序?qū)Ρ荣愔械拿恳磺蜻M行復盤，并向球友演示復盤后優(yōu)化的結(jié)果。

提示：本文的裝逼指南僅適用于羽毛球打的好、數(shù)學物理又不差的人。如果羽毛球打得很爛，又強行裝逼，往往只會顯得畫蛇添足。羽毛球打得好的人，不僅體力好，腦子也快。這樣的人一般很討人異性喜歡。我一陳姓師弟，就因為羽毛球打得好，現(xiàn)已迎娶白富美球友，走向人生巔峰。要是再讓他學會此裝逼指南，必定是如虎添翼。不說了，我要喊他來學裝逼指南了。

參考資料

[1] Caroline Cohen, Baptiste Darbois Texier, David Quere, and Christophe Clanet. The physics of badminton. New Journal of Physics, 17(6):063001, 2015

[2] Wikipedia contributors. Jianzi — Wikipedia, the free encyclopedia, 2020: https://en.wikipedia.org/wiki/Jianzi

[3] Wisual Art Encyclopedia. Three beauties playing battledore and shuttlecock, 2011: https://www.wikiart.org/en/utagawa-toyokuni/three-beauties-playing-battledore-and-shuttlecock

[4] Wikipedia contributors. Girl with a racquet — Wikipedia, the free encyclopedia, 2018: https://en.wikipedia.org/wiki/Girl_with_a_Racquet

[5] Wikipedia contributors. Badminton — Wikipedia, the free encyclopedia, 2020: https://en.wikipedia.org/wiki/Badminton

[6] Baptiste Darbois Texier, Caroline Cohen, Guillaume Dupeux, David Quere, and Christophe Clanet. On the size of sports fields. New Journal of Physics, 16(3):033039, 2014

[7] Alison J Cooke. Shuttlecock aerodynamics. Sports Engineering, 2(2):85–96, 1999

[8] Firoz Alam, Harun Chowdhury, Chavaporn Theppadungporn, and Aleksandar Subic. Measurements of aerodynamic properties of badminton shuttlecocks. Procedia Engineering, 2(2):2487–2492, 2010

[9] Guillaume Laffaye. Comprendre et progresser au badminton: les axes de progression, la technique, les situations. Chiron, 2011

[10] R. Stubbs. The Sports Book: The sports, the rules, the tactics, the techniques. Dorling Kindersley, 2013

[11] Caroline Cohen, Baptiste Darbois-Texier, Guillaume Dupeux, Eric Brunel, David Quéré, and Christophe Clanet. The aerodynamic wall. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 470(2161):20130497, 2014

[12] Wikipedia contributors. Density of air — Wikipedia, the free encyclopedia, 2020: https://en.wikipedia.org/wiki/Density_of_air