




下載本文檔
版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
1、實用標準數學建模作業姓名:學院:計算機科學與技術班級:學號1.在區域xe-2,2,y -2,3內繪制函數z=expA(-x 2-y 2)曲面圖及等值線圖。解:曲面圖如下:回 S3Figure- 1> > x=-2:0.5:2;> > y=-2:0.5:3;> > X,Y=meshgrid(x,y);>> Z=exp(-X.A2-''Y.A2);>> mesh(X,Y,Z)>>等值線圖如下:文案大全>> x=-2:0.5:2;>> y=-2:0.5:3;>> X,Y=mes
2、hgrid(x,y);>> Z=exp(-X.A2-Y.A2);>> mesh(X,Y,Z)>> surf(X,Y,Z)>> surf(X,Y,Z)>> contour(X,Y,Z)>>2.已知一組觀測數據,如表 1所示.1: NNM-1.7*1.4-1.1*08*0.5-0 20.10 4擊0.102B90.117410.I3I580.144830.156560.166220.173320.177500.17853工10 71L1.61,9T22.52.83 1Vi0.176350 1710P0.163020.1S255
3、0.14020.126550 112190.097680.08353鼻3,43-744.34.64.9*0 070150 056760.046870 037290.029140.0223(1) 試用差值方法繪制出 xW -2,4.9區間內的光滑曲線,并比較各種差值算法的優 劣.(2) 試用最小二乘多項式擬合的方法擬合表中的數據,選擇一個能較好擬合數據點的多項式的階次,給出相應多項式的系數和偏差平方和(3) 若表中數據滿足正態分布函數y(x)-e,xW/2Q2 .試用最小二乘非線性擬2 合的方法求出分布參數,二值,并利用鎖求參數值繪制擬合曲線,觀察擬合效果解:(1)分別用最領近插值,分段線性插值
4、(缺省值),分段三次樣條插值,保形分段三 次插值方法繪制在 xW -2,4.9的光滑曲線,圖形如下:回就樣條插值效果最好,其次線性插值,最近點插值效果最差,在這里效果好像不太明顯。最近點插值優點就是速度快,線性插值速度稍微慢一點, 但效果好不少。所以線性插值是個不錯的折中方法。樣條插值,它的目的是試圖讓插值的曲線顯得更平滑,為了這個目的,它們不得不利用到周圍若干范圍內的點,不過計算顯然要比前兩種大許多。MATLA改件如下:>> x0=-2:0.3:4.9;>> y0=0.10289 0.11741 0.13158 0.14483 0.15656 0.16622 0.17
5、332 0.17750 0.17853 .0.17635 0.17109 0.16302 0.15255 0.1402 0.12655 0.11219 0.09768 0.08353 .0.07015 0.05876 0.04687 0.03729 0.02914 0.02236;>> cx=-2:0.3:4.9;>> y1=interp1(cx,y0,cx,'nearest');>> y2=interp1(cx,y0,cx,'linear');>> y3=interp1(cx,y0,cx,'spline&
6、#39;);>> y4=interp1(cx,y0,cx,'cubic');> > subplot(2,2,1),plot(cx,y0,'o',cx,y1,'-r'),title('Nearest Interpolant');> > subplot(2,2,2),plot(cx,y0,'o',cx,y1,'-k'),title('Linear Interpolant');> > subplot(2,2,3),plot(cx,y0,
7、9;o',cx,y1,'-b'),title('Spline Interpolant');> > subplot(2,2,4),plot(cx,y0,'o',cx,y1,'-k'),title('Cubic Interpolant');> > subplot(2,2,1),plot(cx,y0,'o',cx,y1,'-r'),title('Nearest Interpolant');(2),從圖形可以看出曲線函數遵從募函數的形式,設備函
8、數形式為:y =axP可化為ln y =ln a + B ln x.即把非線性函數轉化為線性函數,原線性函數形式為 p(x) = a1x + a0由此我們可以得出p(x)等價于lny;x等價于lnx; P = a1 , lna =a0我們可以先求出a1,a0。求一個線性多項式 p(x) = ax + a0使之在最小二乘準則下擬合這些觀測值,問題即化為匹 =0,畦 =0 ,用MATLAB工具我們可以求得最后的結-a 朋求a°a1使E(a0 al尸min £ y - (ax +a。)利用多元函數極值原理可知,若目標函數E( a0 a1)的極小值存在,一定有果。>>
9、log(x0);> > 10g(y0);> > x0=log(x0);> > y0=log(y0);> > n=length(x0);> > a=sum(x0);> > b=sum(y0);> > c=sum(x0.*y0);> > d=sum(x0.A2);> > a0=(d*b-c*a)*(n*d-aA2);> > a1=(n*c-a*b)/(n*d-aA2);> > a0,a1a0 =-2.5891e+05 - 1.7515e+06i al =0.1045
10、- 0.3558i即系數 a0 為-2.5891e+05 - 1.7515e+06i,al 為 0.1045 - 0.3558i其相應多項式的系數和偏差平方和.我們可以求出 E= -7.2019e+13 + 2.1767e+13i其MATLABC件如下:> > Y=a1*x0+a0;> > e=Y-y0;> > E=sum(e.A2)E =-7.2019e+13 + 2.1767e+13i即其相應多項式的系數和偏差平方和為-7.2019e+13 + 2.1767e+13i? 3.將某物體放置在空氣中, 在t=0時刻測得其溫度 U0=150度,10min后測得
11、溫度 3=87度, 假設空氣白溫度為24度。試建立數學模型給出物體的溫度 u與時間t的關系, 并計算20min后物體的溫度。解:為了解決上述問題,我們首先需要了解有關熱力學的一些基本規律:比如:熱量總是從溫度高的物體向溫度低的物體傳導的;在一定的溫度范圍(其中包括了上述問題的溫度在內),一個物體的溫度與這物體的溫度和其所在介質 的溫度的差值成正比例。這是已為實驗證明了的牛頓冷卻定律。設空氣的溫度為ua,物體在時刻t的溫度為u=u(t),則溫度的變化速度為業。注意熱量總是從溫度高的物體向溫度低的物體傳導的,因而初始溫 dt度大于空氣溫度,即(u°>“),所以溫差u7a恒正;又因為
12、物體的溫度將隨時間而逐漸冷卻,故溫度變化速度 曲包負。因此,由牛頓冷卻定律得到dtdu二一K(u -ua) 出這里的K>0是比例常數。此(1)方程就是冷卻過程的數學模型。為了確定溫度u與時間t的關系,我們需要從上面(1)的方程中解出u。又因為ua是常數,并且u-ua>0,所以我們可以將上述式子改寫成d(u-ua) = -Kdt將此式積分可得到如下式子u -ualn( u -ua) = -Kt c1u -ua = eA (一Kt c1) = ceA( 一Kt)即 u=U3+ceA(-Kt)根據初始條件:t=0時,U=Uo代入上式得C=UO-U a于是 u=uo+(uo-u a)eA(
13、-Kt)又根據條件,當t=10時,u=ui代入上式得ui = ua + (u 0-u a)eA(-10K) 一 1 .K = In (u o-ua)/(u 1-u a) 10根據題意我們可知u°=150, u=87,ua=24,代入得到 “ 1150 -241 .K=一 In= ln 2 =0.0691087 -2410從而 u=24+126eA(-0.069t)這就是物體冷卻時溫度u隨著時間t的變化規律。用t=20代入得u=55.7度4.假設在某商場中,某種商品在t時刻的價格為P(t),若假定其變化率與商 品的需求量D和供給量S之差成正比(比例系數為k),若D = a -bP,S
14、= -c dP其中a,b,c,d均為正常數,若已知初始價格為Po,求任意時刻t時該商品的價格解:一般情況下,某種商品的價格主要服從市場供求關系,由題意我們可知 商品需求量D是價格P的單調遞減函數,商品供給量S是價格P的單調遞增函數, 即D =a -bP,S = -c dP (1)其中a,b,c,d均為常數,且b>0,d>0.當需求量與供給量相等時,由(1)可得供求平衡時的價格 Pe=a上,并稱R b d為均衡價格。由題意得:dp =kD(p)-S(p)dt其中比例系數k>0,用來反應價格的調整進度。將(1)式代入方程可得其中常數=k(b+d),>0,所以此方程的通解為P
15、(t)=P e+CeA(- - t)由于初始價格P(0) = Po代入上式,得C = RPe于是我們可以求出任意時刻價 格P與時刻t之間的函數為:P(t) = P e+(Po-Pe)A(- Zt),并且我們可以得出,因為九>0知,tT +8時P(t) T Pe,說明隨著時間的不 斷推延,實際價格P(t)將逐漸趨近均衡價格Pe。5.農場種植計劃問題 三等,相應的耕地面積 分別為100、300和200km,計劃種植水稻、大豆和玉米.要求三種作物的最低收 獲量分別為190、130和350噸.I 、II 、III等耕地種植三種作物的單產如 表所小.去2:不同等級耕施八同作物單產(t/km?) 類
16、別 | IIIIlTK 米 1412 10若三種作物的售價分別為水稻1.2元/kg,大豆1.50元/kg ,玉米0.80元/kg. 那么(1)如何制訂種植計劃,才能使總產量最大?(2)如何制訂種植計劃,才能使總產值最大? 解:(1):問題分析:確定種植最佳土地分配,即每種等級耕地分別種植水稻、大豆、玉米的 面積模型建立:1, ,決策變量:令x1,x2,x3分別為I II III 三等耕地上種植的水稻面積,令x4,x5,x6分別為I II III三等耕地上種植的大豆面積,令x7,x8,x9 三等耕地上種植的玉米面積。且令為xi (1<=i<=9)面積的耕 地上的產量為Ci.2, 目標
17、函數:總產量最大,即 max- 9T cxi3, 約束條件:最低產量限制:最低水稻產量190噸,最低大豆產量130噸,最低 玉米產量350噸11x1+9.5x 2+9x3 1908x4+6.8X5+6X6 呈 13014x7+12x8+10X9 呈 350耕地面積恒定:Xi +X4+X7=100X2+X5+X8=300X3+X6+X9=200上負條件: Xi, X2, X3, X4, X5, X6, X7, X8, X9= 0數學模型:maX=1lX+9.5X2+9X3+8X4+6.8X 5+6X6+14X7+12X8+10X9'-11x1-9.5x2-9x3190-8x4-6.8x5
18、-6x6-130-14x7-12x8-10x9 <-350x1 + x4 + x7 = 100Jx2 + x5+x8 = 300x3 + x6 + x9= 200x1, x2, x3, x4, x5, x6, x7, x8, x9 之 0用MATLABt解,用命令格式III ,文件如下: >>c=11 9.5 9 8 6.8 6 14 12 10;>> A=-11 -9.5 -9 0 0 0 0 0 00 0 0 -8 -6.8 -6 0 0 0 0 0 0 0 0 0 -14 -12 -10;>> b=-190;-130;-350;>>
19、Aeq=1 0 0 1 0 0 1 0 00 1 0 0 1 0 0 1 00 0 1 0 0 1 0 0 1;>> beq=100;300;200;>> vlb=0;0;0;0;0;0;0;0;0;>> vub=;>> x,fval=linprog(c,A,b,Aeq,beq,vlb,vub) Optimization terminated.x =17.27270.00000.000082.7273300.0000165.00000.00000.000035.0000fval =4.2318e+03即,模型的最優解為(17.2727 0.0 0.0 82.7273 300.0 165.0 0.00.0 35.0) T,目標函數最優值為4.231 1
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
- 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- JJF(紡織)106-2023織物通用磨損性能測試儀校準規范
- JJF(紡織)022-2023纖維切斷器校準規范
- 房地產代理合同匯編(17篇)
- 公司廣告贊助合同書(7篇)
- 2025年初升高暑期數學講義專題05 分式方程與無理方程的重難點突破(含答案)
- 考研復習-風景園林基礎考研試題帶答案詳解(研優卷)
- 2025-2026年高校教師資格證之《高等教育法規》通關題庫附參考答案詳解(完整版)
- 2024年山東華興機械集團有限責任公司人員招聘筆試備考題庫附答案詳解(滿分必刷)
- 2025年黑龍江省五常市輔警招聘考試試題題庫及一套參考答案詳解
- 2025年河北省定州市輔警招聘考試試題題庫附答案詳解(黃金題型)
- 全面指南:2024年醫學整形美容醫院員工手冊
- 2025年中國經濟信息社福建分公司招聘筆試參考題庫含答案解析
- 2025年度食用菌產業園區公共設施運營管理合同3篇
- 《費孝通-鄉土中國》差序格局
- 2023-2024學年天津市和平區八年級(下)期末數學試卷(含答案)
- 2021去遠方上海研學旅行方案申請及綜合反思表
- 2025屆北京市人大學附屬中學高三下學期一模考試英語試題含解析
- 二級等保標準
- 藥棒穴位按摩技術
- 行政案例分析-終結性考核-國開(SC)-參考資料
- 操作系統-001-國開機考復習資料
評論
0/150
提交評論