數學建模：最小二乘擬合實驗

1、 數學建模期末實驗作業(yè) 院 系：數學學院 專 業(yè)：信息與計算科學 年 級：2014級 試題編號：37胡克定律的綜合評價分析背景摘要：利用一個打蛋器和一個物理學公式，毀掉一面六英寸厚的承重墻，這么天方夜譚的事你能相信嗎？但它卻真的發(fā)生了！越獄這一電視劇相信很多人都耳熟，即使沒看過里面的內容，但應該都曾經聽過它的大名。在越獄第一季第六集中，Michael要通過地下管道爬到醫(yī)務室的下面，但是一條重要通道是被封死的，因此必須要把這個封死的墻破壞掉，由于是混凝土結構，因此破壞起來很難，Michael從紋身上拓下魔鬼的畫像，投影在掩住管道入口的墻上，用“胡克定律”計算出最佳位置，再用小巧的打蛋器在承重墻上

2、鉆出了幾個小洞，最后借助這幾個小洞毀掉了這堵承重墻。相信大多數人都覺的很夢幻很不科學，但事實就是這樣的令人驚訝。搜狐娛樂曾經報道過，有越獄粉絲不相信這一情節(jié)，在現實生活中進行實驗，結果真的重現了“胡克定律”鑿墻這一情節(jié)。胡克定律的表達式為F=k·x或F=k·x，其中k是常數，是物體的勁度（倔強）系數。在國際單位制中，F的單位是牛，x的單位是米，它是形變量（彈性形變），k的單位是牛/米。倔強系數在數值上等于彈簧伸長（或縮短）單位長度時的彈力。彈性定律是胡克最重要的發(fā)現之一，也是力學最重要基本定律之一。在現代，仍然是物理學的重要基本理論。胡克的彈性定律指出：彈簧在發(fā)生彈性形變時

3、，彈簧的彈力Ff和彈簧的伸長量（或壓縮量）x成正比，即F= -k·x。k是物質的彈性系數，它由材料的性質所決定，負號表示彈簧所產生的彈力與其伸長（或壓縮）的方向相反。但當我們進行多次實驗，便會發(fā)現隨著F的逐步增大，便不再服從胡克定律。為此我們應當運用插值與擬合的內容，探索更加準確的公式。一、建模問題1.問題提出1.1 問題背景彈簧在壓力F 的作用下伸長x，一定范圍內服從胡克定理：F與x成正比，即F=kx?，F在得到下面一組F，x數據，并在（x，F）坐標下作圖，可以看到當F大到一定數據值后，就不服從這個定律了。表1-1X1247912131517F1.53.96.611.715.618.

4、819.620.621.11.2 問題提出試根據上述所給出的數據及已知的胡克公式，解決一下問題：（1）試由數據確定k（2）給出不服從胡克定理時的近似公式1.3 問題分析這是一道關于彈簧勁度系數的問題，對于此類建模有實際的價值，而且也可以讓我們拓寬物理學習的視野，很有價值。二、模型假設通過閱讀題目與查閱資料，我們可以發(fā)現，F的值是隨著X的改變而改變的，當X小于某一值時，F遵循胡克定律，而當X大于某一值時，F便不再遵循胡克定律，故我們可以提出以下假設。假設1：當X<9時，F遵循胡克定律。假設2：當X>9時，F不遵循胡克定律。三、模型建立 已知胡克定律為：F=KX，但通過簡單的計算題目中

5、所給的數據，便會發(fā)現K的值并非固定值，我們可假設F=KX中還有第三個未知量S。故建立模型公式：F=KX+S運用數學建模與數學實驗（第四版）7.4.1線性最小二乘擬合內容，在matlab程序上可進行求解。四、符號說明表4-1 符號名稱符號含義f彈簧受到的壓力數值k彈簧的彈勁系數x彈簧增加的距離S胡克定律中的未知數值五、模型求解解：（1）試由數據確定k輸入以下程序：x=0 1 2 4 7 9 12 13 15 17;f=0 1.5 3.9 6.6 11.7 15.6 18.8 19.6 20.6 21.1;k=f/x可得結果：k =1.4377（2）給出不服從胡克定理時的近似公式：輸入以下程序：x

6、=0 1 2 4 7 9 12 13 15 17;f=0 1.5 3.9 6.6 11.7 15.6 18.8 19.6 20.6 21.1;a=polyfit(x,f,1)z=polyval(a,x);plot(x,f,'k+',x,z,'r')運行結果：a =1.3340 1.2678可得圖5-2-1：圖5-2-1通過圖5-2-1，可以看到當彈簧伸長10個單位長度后，擬合的情況并不好，偏差較多，且用計算結果得出的公式F=1.3340X+1.2678與胡克定律也相差甚遠，故可以根據圖5-1，將十個數據分為兩組來進行驗算。先取前六個數據的值進行線性擬合：輸入程序

7、：x=0 1 2 4 7 9;f=0 1.5 3.9 6.6 11.7 15.6;a=polyfit(x,f,1)z=polyval(a,x);plot(x,f,'k+',x,z,'r')可得結果：a =1.7085 0.0008可得圖5-2-2：圖5-2-2通過圖5-2-2，可以看到擬合良好，且0.008可以忽略不計，故可以用F=1.7085X來表示力和彈簧伸長的關系，該公式較符合胡克定律。接下來對后面的四個數據進行二次擬合來觀察效果：輸入程序:x=12 13 15 17;f=18.8 19.6 20.6 21.1;a=polyfit(x,f,2)z=poly

8、val(a,x);plot(x,f,'k+',x,z,'r')可得結果：a = -0.0732 2.5790 -1.5834可得圖5-2-3：圖5-2-3通過圖5-2-3可以看出，當才用后四個數據來進行擬合時，擬合情況較為準確，可以接受，近似公式為F=-0.0732+2.5790X-1.5834六、結果分析通過運用matlab分別進行三次擬合，可以發(fā)現三次擬合的結果大不相同。第一次將所有數據進行擬合，擬合的情況并不好，偏差較多，得出的公式與胡克定律也相差甚遠；第二次按照模型假設，只采用0 1 2 4 7 9這六個數據進行擬合，擬合情況較為良好，所得到的公式也極為

9、接近胡克定律；第三次同樣按照模型假設，采用12 13 15 17這四個數據進行擬合，擬合情況同樣良好，且所得公式也符合我們的模型假設。綜合三次擬合，現在可以解答第二步所建立的模型假設。假設1：當X<9時，F遵循胡克定律。假設1結果：當X<9時，F遵循胡克定律，其公式為F=1.7085X假設2：當X>9時，F不遵循胡克定律。假設2結果：當X>9時，F不遵循胡克定律，其近似公式為近似公式為F=-0.0732+2.5790X-1.5834七、實驗心得在進行建模和仿真分析時，人們經常面臨用已知系統實測數據應用數學模型描述對應系統，即對數據進行擬合。擬合的目的是尋找給定的曲線（直

10、線），它在某種準則下最佳的擬合數據。常用的擬合方法之一是多項式的最小二乘擬合，其準則是最小誤差平方和準則，所用的擬合曲線為多項式。在本次建模實驗中，我們所用到的方法就是是線性最小二乘擬合，通過這次實驗，使得我掌握了用線性最小二乘擬合建立回歸數學模型（包括參數估計和模型建立），并通過幾個數據擬合的回歸分析來判斷曲線（直線）擬合的精度，判斷模型建立是否正確。八、參考文獻類別著錄項目專著數學建模與數學實驗/趙靜，但琦主編.4版北京：高等教育出版社，2014.8九、附錄附錄A：其實除我們所熟知的F=KX這一簡單的胡克定律，還有一個廣義的胡克定律，其內容為：在材料的線彈性范圍內（見下圖的材料應力應變曲線的比例極限范圍內），固體的單向拉伸變形與所受的外力成正比；也可表述為：在應力低于比例極限的情況下，固體中的應力與應變成正比，即=，式中E為常數，稱為彈性模量或楊氏模量。把胡克定律推廣應用于三向應力和應變狀態(tài)，則可得到廣義胡克定律。圖