數(shù)學建模-物理模型

-.z.2015年學院第十屆數(shù)學建模競賽承諾書我們仔細閱讀了中國大學生數(shù)學建模競賽的競賽規(guī)則.我們完全明白，在競賽開場后參賽隊員不能以任何方式〔包括、電子、網(wǎng)上咨詢等〕與隊外的任何人研究、討論與賽題有關的問題。我們知道，抄襲別人的成果是違反競賽規(guī)則的,如果引用別人的成果或其他公開的資料〔包括網(wǎng)上查到的資料〕，必須按照規(guī)定的參考文獻的表述方式在正文引用處和參考文獻中明確列出。我們重承諾，嚴格遵守競賽規(guī)則，以保證競賽的公正、公平性。如有違反競賽規(guī)則的行為，我們將受到嚴肅處理。我們參賽的題目是：B題物理實驗模型參賽的組〔甲組或乙組〕：甲組所屬系：土木工程參賽隊員(打印并簽名)：序號班級簽名聯(lián)系1王紅超13土本一班213土本一班3日期：年月日評閱編號〔由組委會評閱前進展編號〕：物理實驗模型摘要首先對數(shù)據(jù)進展綜合處理，給出函數(shù)的定義、性質(zhì)建立模型并用數(shù)學軟件求解物理模型。在MATLAB中運用插值法與最小二乘法進展模型擬合，利用函數(shù)轉(zhuǎn)化為一般的非線性規(guī)劃問題從而有效地進展模型建立，模型求解，模型改良。模型改良中根據(jù)擬合曲線的擬合程度進展確立模型，求解模型之后比擬前后模型的置信區(qū)間、誤差項平方和、AdjustedR-ssquare、R-square之間的大小已確立模型改良的是否成功。關鍵字：插值法最小二乘法擬合曲線擬合程度置信區(qū)間誤差項平方和AdjustedR-squareR-square問題重述：學院物理實驗室一組實驗數(shù)據(jù)〔不考慮物理量的含義〕如表建立*與y之間的數(shù)學模型。模型假設：在e*cel中輸入數(shù)據(jù)，利用插入圖表作出散點圖〔如下列圖〕根據(jù)散點圖走向大致為指數(shù)型函數(shù)，建立模型為其中a>0模型建立：〔1〕在MATLAB中運用插值法輸入以下數(shù)據(jù)：*data=[-1.78

-1.1

-0.86

-0.78

-0.7

-0.67

-0.66

-0.65

-0.64

-0.59

-0.56

-0.49

-0.4

-0.26

-0.12

-0.08

0

0.07

0.25

0.4

0.65];ydata=[-17.5

-15

-10

-7

-3

-0.1

0

0.5

1

4

6

10

20

30

40

45

54

60

80

90

150];

plot(*data,ydata,'bo');title('original

data

points'

);*=linspace(min(*data),ma*(*data),100);

y=interp1(*data,ydata,*,'spline');

hold

on

plot(*,y,'r-');

legend('數(shù)據(jù)點','插值曲線')>>

cftool選擇模型得出擬合曲線〔2〕在MATLAB中運用多項式法輸出數(shù)據(jù)：*data=[-1.78

-1.1

-0.86

-0.78

-0.7

-0.67

-0.66

-0.65

-0.64

-0.59

-0.56

-0.49

-0.4

-0.26

-0.12

-0.08

0

0.07

0.25

0.4

0.65];

ydata=[-17.5

-15

-10

-7

-3

-0.1

0

0.5

1

4

6

10

20

30

40

45

54

60

80

90

150];

plot(*data,ydata,'bo');

title('original

data

points'

);p=polyfit(*data,ydata,2);

*=linspace(min(*data),ma*(*data),100);

y=polyval(p,*);

hold

on

plot(*,y,'r-');

legend('數(shù)據(jù)點','插值曲線')>>

cftool選擇模型得出擬合曲線可見兩種數(shù)據(jù)處理的方法求出最終的擬合曲線一樣，建立模型為。模型求解：對上述建立的模型進展軟件求解得出結(jié)果為：所以得到模型為置信區(qū)間為95%誤差項平方和2676均方根誤差為11.87R-square為0.9246AdjustedR-square為0.9207模型改良：確定另一模型得出擬合曲線如下：根據(jù)圖形的擬合程度可以判定后者模型更加準確。對以上模型求解得：其中a=85.38(72.08,98.69)b=-1.189(-1.731,-1.007)c=-36.7(-47.76,-25.63