[工作總結(jié)]學(xué)生實(shí)習(xí)實(shí)訓(xùn)總結(jié)報(bào)告2011版

1、.重慶科技學(xué)院學(xué)生實(shí)習(xí)（實(shí)訓(xùn)）總結(jié)報(bào)告 學(xué) 院:_ 數(shù)理學(xué)院_ 專業(yè)班級(jí):_ 力學(xué)2009-01_學(xué)生姓名:_ 學(xué) 號(hào):_ _實(shí)習(xí)(實(shí)訓(xùn))地點(diǎn):_ 校內(nèi)數(shù)理學(xué)院I217機(jī)房_ _報(bào)告題目:_ 工程數(shù)值方法上機(jī)實(shí)習(xí)報(bào)告 _ 報(bào)告日期： 年 月 日 指導(dǎo)教師評(píng)語(yǔ): _ _ _ _ 成績(jī)（五級(jí)記分制）:_ _ 指導(dǎo)教師（簽字）:_正文第一周 熟悉Fortran/C/C+語(yǔ)言和Visual Fortan/C+ 編輯編譯環(huán)境第二周實(shí)驗(yàn)名稱：迭代法求解線性方程組實(shí)驗(yàn)?zāi)康模赫莆涨蠼饩€性方程組的常見(jiàn)直接解法，培養(yǎng)編程與上機(jī)調(diào)試能力?；疽螅簯?yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程

2、，并上機(jī)調(diào)試通過(guò)；3學(xué)時(shí)。 實(shí)驗(yàn)結(jié)果：（1）全選主元高斯消去法： x1=1.172493x2=0.209292x3=-0.045339x4=0.798758（2）LU分解法： x1=1.172493x2=0.209292x3=-0.045339x4=0.798758第三周實(shí)驗(yàn)名稱：迭代法求解線性方程組實(shí)驗(yàn)?zāi)康模赫莆涨蠼饩€性方程組的常見(jiàn)迭代方法，培養(yǎng)編程與上機(jī)調(diào)試能力。基本要求：應(yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)；3學(xué)時(shí)。 實(shí)驗(yàn)結(jié)果（給出迭代的初值，控制的精度；迭代步數(shù)和迭代結(jié)果）：（1）簡(jiǎn)單迭代法（Jacobi迭代法）迭代的初值x（0）=0x（

3、1）=0x（2）=0x（3）=0精度0。000001迭代步數(shù)42和迭代結(jié)果x（0）=1.1724927x（1）=0.2092920x（2）=-0.0453392x（3）=0.7987577（2）Gauss-Seidel迭代法 迭代精度=0.000001第三周實(shí)驗(yàn)名稱：矩陣特征值的數(shù)值計(jì)算實(shí)驗(yàn)?zāi)康模赫莆涨缶仃嚨奶卣髦岛椭魈卣飨蛄康膬绶?，培養(yǎng)編程與上機(jī)調(diào)試能力?；疽螅簯?yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)；3學(xué)時(shí)。已知矩陣 ，（1）應(yīng)用冪法和反冪法計(jì)算該矩陣的絕對(duì)值最大和最小特征值. （2）選擇不同的初值，觀察所需的迭代次數(shù)和迭代結(jié)果. （3）給出

4、該矩陣的所有特征值和相應(yīng)的特征向量.實(shí)驗(yàn)結(jié)果： 矩陣求逆輸出結(jié)果：相應(yīng)的特征向量=（0,0,1）.所有特征值=(1,2,3,4,5,6).乘冪法計(jì)算矩陣絕對(duì)值最大特征值和最小特征值：實(shí)驗(yàn)名稱：非線性方程求根實(shí)驗(yàn)?zāi)康模赫莆辗蔷€性方程求根的常見(jiàn)迭代方法，培養(yǎng)編程與上機(jī)調(diào)試能力?；疽螅簯?yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)；3學(xué)時(shí)。已知方程在x1.5附近有一個(gè)實(shí)根。（1）分別對(duì)下面的三種迭代格式應(yīng)用Aitken迭代法求出方程在x1.5附近的實(shí)根. （2）應(yīng)用Newton迭代法求出方程在x1.5附近的實(shí)根. （3）應(yīng)用插值法（割線法）求出方程在x1.5

5、附近的實(shí)根.實(shí)驗(yàn)結(jié)果（給出迭代控制的精度；迭代步數(shù)和迭代結(jié)果）：（1）梯度法 迭代精度=0.000001（2）擬牛頓法 精度：0.000001；迭代步數(shù)：11;迭代結(jié)果：x=2.000000；y=-3.000000；z=1.000000第四周實(shí)驗(yàn)名稱：非線性方程組求根實(shí)驗(yàn)?zāi)康模赫莆辗蔷€性方程組求根的梯度法、擬牛頓法，培養(yǎng)編程與上機(jī)調(diào)試能力。基本要求：應(yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)；3學(xué)時(shí)。 （初值取x=1, y= -1, z=0）實(shí)驗(yàn)結(jié)果（給出迭代控制的精度；迭代步數(shù)和迭代結(jié)果）：（1）梯度法 迭代精度=0.000001（2）擬牛頓法 精度

6、：0.000001；迭代步數(shù)：11;迭代結(jié)果：x=2.000000；y=-3.000000；z=1.000000第五周實(shí)驗(yàn)名稱：多項(xiàng)式插值實(shí)驗(yàn)一實(shí)驗(yàn)?zāi)康模赫莆斩囗?xiàng)式插值的拉格朗日方法和逐步埃特肯方法，培養(yǎng)編程與上機(jī)調(diào)試能力?；疽螅簯?yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)；3學(xué)時(shí)。 已知函數(shù)f(x)的值表如下表所示，試分別用拉格朗日插值方法和逐步埃特肯方法線性插值方法計(jì)算 f(1.813)的近似值。 x1.6151.6341.7021.8281.921f(x)2.414502.464592.652713.030353.34066實(shí)驗(yàn)結(jié)果（給出插值多

7、項(xiàng)式的次數(shù)，控制的精度和插值結(jié)果）：（1）拉格朗日插值方法 插值多項(xiàng)式的次數(shù)：4，控制的精度：0.00001，插值結(jié)果：2.98281。（2）逐步埃特肯方法線性插值方法 插值多項(xiàng)式的次數(shù)：4，控制的精度：0.000001，插值結(jié)果：2.983322。第六周實(shí)驗(yàn)名稱：多項(xiàng)式插值實(shí)驗(yàn)二實(shí)驗(yàn)?zāi)康模毫私獠煌逯捣椒ǖ恼`差，掌握三次樣條插值方法，培養(yǎng)編程與上機(jī)調(diào)試能力?；疽螅簯?yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)；3學(xué)時(shí)。 已知函數(shù)f(x)的值表如下表所示： x00.7851.572.3563.1423.9274.7125.4986.283f(x)00.

8、70710.7070-0.707-1-0.7070實(shí)驗(yàn)結(jié)果（給出插值多項(xiàng)式的次數(shù)，控制的精度和插值結(jié)果）：（1）三點(diǎn)拋物線插值方法 x= 0.524, f(x)=0.517875（2）九點(diǎn)多項(xiàng)式插值方法 t= 0.524, z=0.500420（3）采用自然邊界條件的三次樣條函數(shù)插值方法 t(i)= 0.524000 ，z (i)=0.500230（4）采用周期邊界條件的三次樣條函數(shù)插值方法 （5）四種插值方法的絕對(duì)誤差分析 第一種絕對(duì)誤差6.47%第一種絕對(duì)誤差2.66%第一種絕對(duì)誤差3.15%第七周實(shí)驗(yàn)名稱：離散數(shù)據(jù)的最小二乘擬合實(shí)驗(yàn)?zāi)康模赫莆兆钚《嗽?，?huì)求離散數(shù)據(jù)的最小二乘擬合曲線并

9、畫(huà)圖。基本要求：應(yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)；應(yīng)用Excel或者Origin軟件繪制離散數(shù)據(jù)和擬合曲線的對(duì)比圖；3學(xué)時(shí)。1數(shù)據(jù)擬合Malthus人口指數(shù)增長(zhǎng)模型中參數(shù)從17901980年間美國(guó)每隔10年的人口記錄如下表：年 份1790180018101820183018401850人口(×106)3.95.37.29.612.917.123.2年 份1860187018801890190019101920人口(×106)31.438.650.262.976.092.0106.5年 份193019401950196019

10、701980  人口(×106)123.2131.7150.7179.3204.0226.5  用以上數(shù)據(jù)檢驗(yàn)馬爾薩斯(Malthus)人口指數(shù)增長(zhǎng)模型，根據(jù)檢驗(yàn)結(jié)果進(jìn)一步討論馬爾薩斯人口模型的改進(jìn)。提示：Malthus 模型的基本假設(shè)是：人口的增長(zhǎng)率為常數(shù)，記為 r。記時(shí)刻t的人口為x(t），（即x(t）為模型的狀態(tài)變量）且初始時(shí)刻的人口為x0，于是得到如下微分方程：微分方程的解為，再用數(shù)據(jù)擬合模型中的參數(shù)。2舊車價(jià)格預(yù)測(cè)  某年美國(guó)舊車價(jià)格的調(diào)查資料如下表，其中xi表示轎車的使用年數(shù)，yi表示相應(yīng)的平均價(jià)格。試分析用什么形式的曲線來(lái)擬合上述的數(shù)據(jù)（用

11、線性和二次多項(xiàng)式），并預(yù)測(cè)使用4.5年后轎車的平均價(jià)格大致為多少？xi12345678910yi2615194314941087765538484290226204實(shí)驗(yàn)結(jié)果（要求繪出離散數(shù)據(jù)和擬合曲線的對(duì)比圖）：1數(shù)據(jù)擬合Malthus人口指數(shù)增長(zhǎng)模型中參數(shù) 2舊車價(jià)格預(yù)測(cè) 第八、九、十周實(shí)驗(yàn)名稱：數(shù)值積分實(shí)驗(yàn)一實(shí)驗(yàn)?zāi)康模赫莆仗菪喂健impson算法、復(fù)化梯形公式、復(fù)化Simpson算法?；疽螅簯?yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)； 3學(xué)時(shí)。已知飛機(jī)在高度H的上升速度如下表所示：先用二次拋物線擬合函數(shù)關(guān)系，然后分別用變步長(zhǎng)梯形求積法和變步

12、長(zhǎng)Simpson求積法計(jì)算飛機(jī)從地面（）上升到千米高空所需要的時(shí)間實(shí)驗(yàn)結(jié)果：1變步長(zhǎng)梯形求積法 t=0.3479612步長(zhǎng)Simpson求積法 t=0.347961 第十一、二周實(shí)驗(yàn)名稱：數(shù)值積分實(shí)驗(yàn)二實(shí)驗(yàn)?zāi)康模赫莆正堌惛穹e分法和高斯求積法?；疽螅簯?yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)； 3學(xué)時(shí)。1. 分別用變步長(zhǎng)梯形求積法和龍貝格求積法計(jì)算積分（要求計(jì)算精度為） 2. 下圖所示的均質(zhì)橢圓鋼板密度為7.8×103kg/m3，長(zhǎng)軸為0.5m，短軸為0.3m，厚度為0.4m。 分別用高斯求積法和龍貝格求積法計(jì)算均質(zhì)橢圓鋼板的轉(zhuǎn)動(dòng)慣量Jx，

13、要求計(jì)算精度為 。實(shí)驗(yàn)結(jié)果：1變步長(zhǎng)梯形求積法和龍貝格求積法計(jì)算積分 變步長(zhǎng)梯形：t=0.842709，龍貝格：t=0.8427092均質(zhì)橢圓鋼板的轉(zhuǎn)動(dòng)慣量高斯法：g=33.082106，龍貝格：t=33.076978第十三、四周實(shí)驗(yàn)名稱：常微分方程實(shí)驗(yàn)一實(shí)驗(yàn)?zāi)康模赫莆障蚯皻W拉公式、向后歐拉公式、改進(jìn)的歐拉公式、龍格-庫(kù)塔法?；疽螅簯?yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)； 3學(xué)時(shí)。至少采用兩種不同的方法、不同的步長(zhǎng)，求解初值問(wèn)題 ，比較不同方法、不同步長(zhǎng)的計(jì)算y(0.5), 并比較誤差。實(shí)驗(yàn)結(jié)果： 前歐拉公式h=0.1Y1=1Y2=1.1Y3=

14、1.191818Y4=1.277438Y5=1.358213Y6=1.435133Y7=1.508966Y8=1.580338Y9=1.649783Y10=1.717779龍格-庫(kù)塔法y0.000000=1.001003y0.150000=1.132129y0.300000=1.246727y0.450000=1.348094y0.510000=1.385342y0.600000=1.437864y0.750000=1.516604y0.900000=1.584032y0.960000=1.607615第十五、六周實(shí)驗(yàn)名稱：常微分方程實(shí)驗(yàn)二（導(dǎo)彈跟蹤問(wèn)題仿真實(shí)驗(yàn)）實(shí)驗(yàn)?zāi)康模赫莆粘Ｎ⒎址匠探M的數(shù)

15、值解法?；疽螅簯?yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)； 3學(xué)時(shí)。某軍的一導(dǎo)彈基地發(fā)現(xiàn)正北方向120 km處海面上有敵艇一艘以90 km/h的速度向正東方向行駛. 該基地立即發(fā)射導(dǎo)彈跟蹤追擊敵艇, 導(dǎo)彈速度為450 km/h，自動(dòng)導(dǎo)航系統(tǒng)使導(dǎo)彈在任一時(shí)刻都能對(duì)準(zhǔn)敵艇.試問(wèn)導(dǎo)彈在何時(shí)何處擊中敵艇？ （導(dǎo)彈和敵艇距離為1米時(shí)即可認(rèn)為擊中）當(dāng)t =0時(shí),導(dǎo)彈位于原點(diǎn)O,敵艇位于(0,H)點(diǎn), H= 120 (km). P(x, y)當(dāng)時(shí)刻t ,導(dǎo)彈位于P(x(t),y(t)敵艇位于(90t,H)點(diǎn)O120120x東y北實(shí)驗(yàn)結(jié)果： t = 0.2778 x = 25.0058 y = 120.0031第十七周綜合訓(xùn)練實(shí)驗(yàn)?zāi)康模簯?yīng)用數(shù)值方法求解力學(xué)問(wèn)題?；疽螅簯?yīng)用C語(yǔ)言(或者C+、Fortran、Matlab、Maple)編程，并上機(jī)調(diào)試通過(guò)； 3學(xué)時(shí)。1. 運(yùn)用節(jié)點(diǎn)法，自主設(shè)計(jì)編程數(shù)值計(jì)算下圖所示桁架結(jié)構(gòu)中各桿的內(nèi)力（P