數(shù)值分析實(shí)驗(yàn)應(yīng)用數(shù)學(xué)研究生

1、數(shù)值分析實(shí)驗(yàn)題 理學(xué)院 應(yīng)用數(shù)學(xué)研究生實(shí)驗(yàn)1 病態(tài)問(wèn)題與數(shù)值穩(wěn)定性實(shí)驗(yàn)1.1 病態(tài)問(wèn)題實(shí)驗(yàn)?zāi)康模核惴ㄓ小皟?yōu)”與“劣”之分，問(wèn)題也有“好”和“壞”之別.對(duì)數(shù)值方法的研究而言，所謂壞問(wèn)題就是問(wèn)題本身對(duì)擾動(dòng)敏感者，反之屬于好問(wèn)題.希望讀者通過(guò)本實(shí)驗(yàn)對(duì)此有一個(gè)初步的體會(huì).數(shù)值分析的大部分研究課題中，如線性代數(shù)方程組，矩陣特征值問(wèn)題，非線性方程及方程組等都存在病態(tài)的問(wèn)題.病態(tài)問(wèn)題要通過(guò)研究和構(gòu)造特殊的算法來(lái)解決，當(dāng)然一般要付出一些代價(jià)（如耗用更多的機(jī)器時(shí)間，占用更多的存儲(chǔ)空間等）問(wèn)題提出：考慮一個(gè)高次的代數(shù)多項(xiàng)式， （1）顯然該多項(xiàng)式的全部根為1,2，20共計(jì)20個(gè)，且每個(gè)根都是單重的（也稱為簡(jiǎn)單的）.

2、現(xiàn)考慮該多項(xiàng)式的一個(gè)擾動(dòng), （2）式中，是一個(gè)非常小的數(shù).這相當(dāng)于是對(duì)（1）中的系數(shù)作一個(gè)小的擾動(dòng).希望比較式（1）和（2）根的差別，從而分析方程（1）的解對(duì)擾動(dòng)的敏感性.實(shí)驗(yàn)內(nèi)容：為了實(shí)現(xiàn)方便，介紹兩個(gè)Matlab函數(shù)："roots"和"poly". u=roots(a)，其中若變量a存儲(chǔ)n+1維的向量，則該函數(shù)的輸出u為一個(gè)n維的向量.設(shè)a的元素依次為，則輸出u的各分量是多項(xiàng)式方程 的全部根；而函數(shù)b=poly(v)的輸出b是一個(gè)維變量，它是以維變量的各分量為根的多項(xiàng)式系數(shù).可見(jiàn)"roots"和"poly"是兩

3、個(gè)互逆的運(yùn)算函數(shù).ve=zeros(1,21); ve(2)=ess; roots(poly(1:20)+ve)上述簡(jiǎn)單的Matlab程序便得到（2）的全部根，程序中的“ess”即是（2）中的.實(shí)驗(yàn)要求：（1）選擇充分小的ess，反復(fù)進(jìn)行上述實(shí)驗(yàn)，記錄結(jié)果的變化并分析它們.如果擾動(dòng)項(xiàng)的系數(shù)很小，我們自然感覺(jué)式（1）和式（2）的解應(yīng)當(dāng)相差很小.計(jì)算中你有什么出乎意料的發(fā)現(xiàn)？表明有些解關(guān)于如此的擾動(dòng)敏感性如何？（2）將方程（2）中的擾動(dòng)項(xiàng)改成或其他形式，實(shí)驗(yàn)中又有怎樣的現(xiàn)象出現(xiàn)？（3）請(qǐng)從理論上分析產(chǎn)生這一問(wèn)題的根源.注意我們可以將方程（2）寫(xiě)成展開(kāi)的形式 . （3）同時(shí)將方程的解看成是系數(shù)的函數(shù)

4、，考察方程的某個(gè)解關(guān)于的擾動(dòng)是否敏感，與研究它關(guān)于的導(dǎo)數(shù)的大小有何關(guān)系？為什么？你發(fā)現(xiàn)了什么現(xiàn)象，哪些根關(guān)于的變化更敏感？實(shí)驗(yàn)1.2 誤差傳播與算法穩(wěn)定性實(shí)驗(yàn)?zāi)康模后w會(huì)穩(wěn)定性在選擇算法中的地位，誤差擴(kuò)張的算法是不穩(wěn)定的，是我們所不期望的；誤差衰竭的算法是穩(wěn)定的，是我們努力尋求的，這是貫穿本課程的目標(biāo).問(wèn)題提出：考慮一個(gè)簡(jiǎn)單的由積分定義的序列.顯然利用分部積分易得,. （1）又有.實(shí)驗(yàn)內(nèi)容：由遞推關(guān)系（1），可以得到計(jì)算積分序列的兩種算法 算法1： (2)算法2： (3)實(shí)驗(yàn)要求：（1）分別用算法1、算法2并在計(jì)算中分別采用5位、6位和7位有效數(shù)字，請(qǐng)判斷哪種算法能給出更精確的結(jié)果.（2）兩種算

5、法的優(yōu)劣，與你的第一感覺(jué)是否吻合.請(qǐng)從理論上證明你實(shí)驗(yàn)得出的結(jié)果，解釋實(shí)驗(yàn)的結(jié)果.設(shè)算法1中的計(jì)算誤差為,由遞推計(jì)算的誤差為；算法2中的誤差為，由向前遞推計(jì)算誤差為.如果在上述兩算法中都假定后面的計(jì)算不再引入其他誤差，試給出與的關(guān)系和的關(guān)系.（3）算法1中通常會(huì)很小，當(dāng)n增大時(shí)，的變化趨勢(shì)如何？算法2中通常相對(duì)比較大，當(dāng)n減小時(shí)，誤差又是如何傳播的？也就是說(shuō)比較一下上述兩個(gè)算法，當(dāng)某一步產(chǎn)生誤差后，該誤差對(duì)后面的影響是衰減還是擴(kuò)張的.（4）通過(guò)理論分析與計(jì)算實(shí)驗(yàn)，針對(duì)算法1和算法2的穩(wěn)定性，給出你的結(jié)論.實(shí)驗(yàn)2 非線性方程（組）實(shí)驗(yàn)2.1 迭代法、初始值與收斂性實(shí)驗(yàn)?zāi)康模撼醪秸J(rèn)識(shí)非線性問(wèn)題的迭

6、代法與線性問(wèn)題迭代法的差別，探討迭代法及初始值與迭代收斂性的關(guān)系.問(wèn)題提出：迭代法是求解非線性方程及方程組的基本思想方法，與線性方程的情況一樣，其構(gòu)造方法可以有多種多樣，但關(guān)鍵是怎樣才能使迭代收斂且有較快的收斂速度.實(shí)驗(yàn)內(nèi)容：考慮一個(gè)簡(jiǎn)單的代數(shù)方程可以構(gòu)造多種迭代法，如 (1) (2) (3)在實(shí)軸上取初始值,請(qǐng)分別用迭代（1）（3）做實(shí)驗(yàn)，記錄各算法的迭代過(guò)程.實(shí)驗(yàn)要求：（1）取定某個(gè)初始值，分別計(jì)算（1）（3）迭代結(jié)果，收斂性如何？重復(fù)選取不同的初始值，反復(fù)實(shí)驗(yàn)，請(qǐng)讀者自行設(shè)計(jì)一種比較形象的記錄方式（如利用Matlab的圖形功能），分析三種迭代法的收斂性與初值選取的關(guān)系？（2）對(duì)三個(gè)迭代法

7、中的某個(gè)，取不同的初始值進(jìn)行迭代，結(jié)果如何？試分析迭代法對(duì)不同的初值是否有差異？（3）線性方程組迭代法的收斂性是不依賴初始值選取的，比較線性與非線性問(wèn)題迭代的差異，有何結(jié)論與問(wèn)題？實(shí)驗(yàn)2.2 迭代與混沌實(shí)驗(yàn)?zāi)康模貉芯恳话愕降膹?fù)雜行為，初步看到混沌現(xiàn)象.實(shí)驗(yàn)內(nèi)容：考慮迭代公式.實(shí)驗(yàn)要求：取中不同的值，并取進(jìn)行迭代，畫(huà)出不同情況下的的圖形，并分析取值與圖形的關(guān)系.你將對(duì)于迭代法有更深刻的理解.實(shí)驗(yàn)3 線性方程組實(shí)驗(yàn)3.1主元的選取與算法的穩(wěn)定性問(wèn)題提出：Gauss消去法是我們?cè)诰€性代數(shù)中已經(jīng)熟悉的，但由于計(jì)算機(jī)的數(shù)值運(yùn)算是在一個(gè)有限的浮點(diǎn)數(shù)集合上進(jìn)行的，如何才能確保Gauss消去法作為數(shù)值算

8、法的穩(wěn)定性呢？Gauss消去法從理論上算法到數(shù)值算法，其關(guān)鍵是主元的選擇.主元的選擇從數(shù)學(xué)理論上看起來(lái)平凡，它卻是數(shù)值分析中十分典型的問(wèn)題.實(shí)驗(yàn)內(nèi)容：考慮線性方程組.編制一個(gè)能自動(dòng)選取主元，又能手動(dòng)選取主元的求解線性代數(shù)方程組的Gauss消去法過(guò)程.實(shí)驗(yàn)要求： （1）取矩陣則方程有解.取計(jì)算矩陣的條件數(shù).讓程序自動(dòng)選取主元，結(jié)果如何？（2）現(xiàn)選擇程序中手動(dòng)選取主元的功能.每步消去過(guò)程總選取按模最小或按模盡可能小的元素作為主元，觀察并記錄計(jì)算結(jié)果，若每步消去過(guò)程總選取按模最大的元素作為主元，結(jié)果又如何？分析實(shí)驗(yàn)的結(jié)果.（3）取矩陣階數(shù) 或者更大，重復(fù)上述實(shí)驗(yàn)過(guò)程，觀察記錄并分析不同的問(wèn)題及消去過(guò)

9、程中選擇不同的主元時(shí)計(jì)算結(jié)果的差異，說(shuō)明主元素的選取在消去過(guò)程中的作用.（4）選取其他你感興趣的問(wèn)題或者隨機(jī)生成矩陣，計(jì)算其條件數(shù).重復(fù)上述實(shí)驗(yàn)，觀察記錄并分析實(shí)驗(yàn)的結(jié)果.實(shí)驗(yàn)3.2 條件數(shù)與線性方程組的性態(tài)實(shí)驗(yàn)?zāi)康模豪斫鈼l件數(shù)的意義和方程組的性態(tài)對(duì)解向量的影響.實(shí)驗(yàn)內(nèi)容：設(shè)和，其中，.由相應(yīng)矩陣的元素計(jì)算，其計(jì)算公式為.實(shí)驗(yàn)要求： 對(duì)取，下面均用Matlab函數(shù)求解方程組.（1）取=4,6,8，分別計(jì)算和，判斷和是否為病態(tài)矩陣？隨的增大，矩陣性態(tài)的變化如何？（2）取=6，分別求出兩個(gè)方程組的解向量；（3）取=6，不變，對(duì)的元素和分別加一個(gè)攝動(dòng)10-6，分別求出的解向量；（4）取=6，不變，對(duì)

10、的元素和分別加一個(gè)攝動(dòng)10-7，分別求出的解向量；不變，對(duì)的元素加一個(gè)攝動(dòng)10-4，求出的解向量；（5）觀察和分析和的微小擾動(dòng)對(duì)解向量的影響，得出你的結(jié)論；（6）求，和的計(jì)算結(jié)果和理論估計(jì).實(shí)驗(yàn)3.3 迭代法初值、收斂性與收斂速度實(shí)驗(yàn)?zāi)康模豪斫獾ㄊ諗康暮x以及迭代初值和方程組系數(shù)矩陣性質(zhì)對(duì)收斂速度的影響.實(shí)驗(yàn)內(nèi)容：用迭代法解方程組，其中實(shí)驗(yàn)要求：（1）選取不同的初始向量和不同的方程組右端向量，給定迭代誤差要求，用雅可比迭代法與高斯塞德?tīng)柕ㄓ?jì)算，觀測(cè)得到的迭代向量序列是否收斂？若收斂，記錄迭代次數(shù)，分析計(jì)算結(jié)果并得出你的結(jié)論；（2）取定初始向量和右端向量，將的主對(duì)角線元素成倍增長(zhǎng)若干次，

11、非主對(duì)角線元素不變，每次用雅可比迭代法計(jì)算.要求迭代誤差滿足，比較收斂速度，分析現(xiàn)象并得出你的結(jié)論.實(shí)驗(yàn)3.4 病態(tài)線性方程組求解問(wèn)題提出：理論的分析表明，求病態(tài)的線性方程組是困難的.實(shí)際情況是否如此，會(huì)出現(xiàn)怎樣的現(xiàn)象呢？實(shí)驗(yàn)內(nèi)容：考慮方程組=的求解，其中系數(shù)矩陣為Hilbert矩陣，.這是個(gè)著名的病態(tài)問(wèn)題.通過(guò)首先給定解（例如取為各個(gè)分量均為1）再計(jì)算出右端的辦法給出確定的問(wèn)題.實(shí)驗(yàn)要求：（1）選擇問(wèn)題的維數(shù)為6，分別用Gauss消去法（即LU分解）、J迭代方法，GS迭代和SOR迭代求解方程組，其各自的結(jié)果如何？將計(jì)算結(jié)果與問(wèn)題的解比較，結(jié)論如何.（2）逐步增大問(wèn)題的維數(shù)，仍然用上述的方法來(lái)

12、解它們，計(jì)算的結(jié)果如何？計(jì)算的結(jié)果說(shuō)明了什么？（3）討論病態(tài)問(wèn)題求解的算法.實(shí)驗(yàn)4 插值法實(shí)驗(yàn)4.1 多項(xiàng)式插值的振蕩現(xiàn)象問(wèn)題提出: 考慮在一個(gè)固定的區(qū)間上用插值逼近一個(gè)函數(shù).顯然拉格朗日插值中使用的節(jié)點(diǎn)越多，插值多項(xiàng)式的次數(shù)就越高.我們自然關(guān)心插值多項(xiàng)式的次數(shù)增加時(shí)，是否也更加靠近被逼近的函數(shù).龍格給出的一個(gè)例子是極著名并富有啟發(fā)性的.設(shè)區(qū)間上函數(shù).實(shí)驗(yàn)內(nèi)容：考慮區(qū)間的一個(gè)等距劃分，分點(diǎn)為,則拉格朗日插值多項(xiàng)式為，其中是次拉格朗日插值基函數(shù).實(shí)驗(yàn)要求：（1）選擇不斷增大的分點(diǎn)數(shù)目畫(huà)出原函數(shù)及插值多項(xiàng)式函數(shù)在上的圖像，比較并分析實(shí)驗(yàn)結(jié)果.（2）選擇其他的函數(shù)，例如定義在區(qū)間上的函數(shù) ， 重復(fù)上

13、述的實(shí)驗(yàn)看其結(jié)果如何.（3）區(qū)間上切比雪夫點(diǎn)的定義為，以為插值節(jié)點(diǎn)構(gòu)造上述各函數(shù)的拉格朗日插值多項(xiàng)式，比較其結(jié)果.實(shí)驗(yàn)4.2 樣條插值的收斂性問(wèn)題提出：多項(xiàng)式插值是不收斂的，即插值的節(jié)點(diǎn)多，效果不一定就好，對(duì)樣條函數(shù)插值又如何呢？理論上證明樣條插值的收斂性是比較困難的，也超出了本課程的內(nèi)容.通過(guò)本實(shí)驗(yàn)可以驗(yàn)證這一理論結(jié)果.實(shí)驗(yàn)內(nèi)容：請(qǐng)按一定的規(guī)劃分別選擇等距離或者非等距離的插值節(jié)點(diǎn)，并不斷增加插值節(jié)點(diǎn)的個(gè)數(shù).考慮實(shí)驗(yàn)4.1中的函數(shù)或者選擇其它你有興趣的函數(shù)，可以用matlab的函數(shù)“spline”作此函數(shù)的三次樣條插值. matlab還提供有spline工具箱（toolbox）.你可以找到極豐

14、富的樣條工具，包括B-樣條.實(shí)驗(yàn)要求：（1）隨著節(jié)點(diǎn)個(gè)數(shù)的增加，比較被逼近函數(shù)和樣條插值函數(shù)誤差的變化情況.分析所得結(jié)果并與拉格朗日多項(xiàng)式插值比較.（2）樣條插值的思想是最早產(chǎn)生 于工業(yè)部門(mén).作為工業(yè)應(yīng)用的例子考慮如下問(wèn)題：某汽車制造商用三次樣條插值設(shè)計(jì)車門(mén)的曲線，其中一段的數(shù)據(jù)如下：xk0 1 2 3 4 5 6 7 8 9 10ykyk0.0 0.79 1.53 2.19 2.71 3.03 3.27 2.89 3.06 3.19 3.290.8 0.2（3）對(duì)實(shí)驗(yàn)4.1中的問(wèn)題，計(jì)算他們的B-樣條插值.實(shí)驗(yàn)5 函數(shù)逼近實(shí)驗(yàn)5.1 最小二乘擬合實(shí)驗(yàn)?zāi)康模鹤钚《藬M合經(jīng)驗(yàn)公式.實(shí)驗(yàn)內(nèi)容：某類

15、疾病發(fā)病率為和年齡段(每五年為一段，例如05歲為第一段，610歲為第二段,)之間有形如的經(jīng)驗(yàn)關(guān)系，觀測(cè)得到的數(shù)據(jù)表如下1234567890.8982.383.071.842.021.942.222.774.02101112131415161718194.765.466.5310.916.522.535.750.661.681.8實(shí)驗(yàn)要求：（1）用最小二乘法確定模型中的參數(shù)和.（2）利用MATLAB畫(huà)出離散數(shù)據(jù)及擬合函數(shù)圖形.（3）利用MATLAB畫(huà)出離散點(diǎn)處的誤差圖，并計(jì)算相應(yīng)的均方誤差.（4）談一談你對(duì)最小二乘擬合的理解，并舉出一個(gè)應(yīng)用此方法的例子.實(shí)驗(yàn)5.2 逼近方法比較實(shí)驗(yàn)?zāi)康模豪斫獠煌?/p>

16、的逼近方法.實(shí)驗(yàn)內(nèi)容：對(duì)被逼近函數(shù)，在區(qū)間0,1上分別利用插值方法、最佳一致逼近以及最佳平方逼近三種方法求形如的逼近函數(shù)，并進(jìn)行比較.再畫(huà)出每種方法的逼近函數(shù)曲線以及誤差圖.實(shí)驗(yàn)5.3 最佳平方逼近多項(xiàng)式的收斂性實(shí)驗(yàn)?zāi)康模貉芯孔罴哑椒奖平囗?xiàng)式的收斂性質(zhì).實(shí)驗(yàn)內(nèi)容：取函數(shù)，在-1，1上以勒讓德多項(xiàng)式為基函數(shù)，對(duì)于構(gòu)造最佳平方逼近多項(xiàng)式，令，將的曲線畫(huà)在一個(gè)圖上.令，畫(huà)出的曲線.做出之間的最小二乘曲線，能否提出關(guān)于收斂性的猜測(cè).實(shí)驗(yàn)6 微積分實(shí)驗(yàn)6.1 高斯數(shù)值積分方法用于積分方程求解問(wèn)題提出：線性的積分方程的數(shù)值求解，可以被轉(zhuǎn)化為線性代數(shù)方程組的求解問(wèn)題.而線性代數(shù)方程組所含未知數(shù)的個(gè)數(shù)，與用

17、來(lái)離散積分的數(shù)值方法的節(jié)點(diǎn)個(gè)數(shù)相同.在節(jié)點(diǎn)數(shù)相同的前提下，高斯數(shù)值積分方法有較高的代數(shù)精度，用它通常會(huì)得到較好的結(jié)果.實(shí)驗(yàn)內(nèi)容：求解第二類Fredholm積分方程.首先將積分區(qū)間a,b等分成n份，在每個(gè)子區(qū)間上離散方程中的積分就得到線性代數(shù)方程組.實(shí)驗(yàn)要求：分別使用如下方法，離散積分方程中的積分. 復(fù)化梯形方法； 復(fù)化辛普森方法； 復(fù)化高斯方法.求解如下的積分方程.（1），方程的準(zhǔn)確解為；（2），方程的準(zhǔn)確解為；比較各方法的計(jì)算量和誤差以分析它們的優(yōu)劣.實(shí)驗(yàn)6.2 高斯積分?jǐn)?shù)值計(jì)算的蒙特卡羅方法問(wèn)題提出：高維空間中的積分，如果維數(shù)不很高且積分區(qū)域是規(guī)則的或者能等價(jià)地寫(xiě)成多重積分的形式，可以用一

18、元函數(shù)積分的數(shù)值方法來(lái)計(jì)算高維空間的積分，蒙特卡羅方法對(duì)計(jì)算復(fù)雜區(qū)域甚至不連通的區(qū)域上的積分并沒(méi)有特殊的困難.實(shí)驗(yàn)內(nèi)容：對(duì)于一般的區(qū)域，計(jì)算其測(cè)度（只要理解為平面上的面積或空間中的體積）的一般方法是：先找一個(gè)規(guī)則的區(qū)域A包含，且A的測(cè)度是已知的.生成區(qū)域A中m個(gè)均勻分布的隨機(jī)點(diǎn)，如果其中有n個(gè)落在區(qū)域中，則區(qū)域的測(cè)度為n/m.函數(shù)在區(qū)域上的積分可以近似為：區(qū)域的測(cè)度與函數(shù)在中n個(gè)隨機(jī)點(diǎn)平均值的乘積，即.實(shí)驗(yàn)要求：假設(shè)冰淇淋的下部為一錐體而上面為一半球，考慮冰淇淋錐體體積問(wèn)題：計(jì)算錐面上方和球面內(nèi)部區(qū)域的體積.如果使用球面坐標(biāo)，該區(qū)域可以表示為如下的積分：.用蒙特卡羅方法可以計(jì)算該積分. 另一方

19、面，顯然這樣的冰淇淋可以裝在如下立方體的盒子里 而該立方體的體積為8.只要生產(chǎn)這個(gè)盒子里均勻分布的隨機(jī)點(diǎn)落入冰淇淋錐點(diǎn)的個(gè)數(shù)與總點(diǎn)數(shù)之比再乘以8就是冰淇淋錐的體積.比較兩種方法所得到的結(jié)果.實(shí)驗(yàn)7 常微分方程實(shí)驗(yàn)7.1 Lorenz問(wèn)題與混沌問(wèn)題提出：考慮著名的Lorenz方程 （1）其中s,r,b為變化區(qū)域有一定限制的實(shí)參數(shù).該方程形式簡(jiǎn)單，表面上看并無(wú)驚人之處，但由該方程揭示出的許多現(xiàn)象，促使“混沌”成為數(shù)學(xué)研究的嶄新領(lǐng)域，在實(shí)際應(yīng)用中也產(chǎn)生了巨大的影響.實(shí)驗(yàn)內(nèi)容：計(jì)算可以直接使用下面的Matlab程序，其中方程的形式與（1）略有差異.% BLZ Plot the orbit around

20、 the Lorenz chaotic attractor.clf, clc% Solve the ordinary differential equation describing the Lorenz chaotic attractor.The equations % are defined in an M-file,blzeq.m. The values of the global parameters are global SIGMA RHO % BETASIGMA=10.; RHO=28.; BETA=8./3.;% The graphics axis limits are set

21、to values known to contain the solution.axis(10 40 -20 20 -20 20), view(3), hold on, title(Lorenz Attractor)y0=0 0 eps; tfinal=100; y=ode23(blzep,0,tfinal,y0)其中的函數(shù)blzeq.m定義為% BLZEQ Equation of the Lorenz chaotic attractor. ydot=lorenzeq(t,y).% The differential equation is written in almost linear fo

22、rm.Global SIGMA RHO BETAA=-BETA 0 y(2)0 -SIGMA SIGMA -y(2) RHO -1 ;Ydot=A*y;實(shí)驗(yàn)要求：（1）請(qǐng)找出（1）與上述程序中使用的方程間的關(guān)系.（2）對(duì)目前取定的參數(shù)值SIGMA、RHO和BETA，選取不同的初始值y0（當(dāng)前的程序中的y0是坐標(biāo)原點(diǎn)），運(yùn)行上述的程序，觀察計(jì)算的結(jié)果有什么特點(diǎn)？解的曲線是否有界？解的曲線是不是周期的或趨于某個(gè)固定的點(diǎn)？（3）在問(wèn)題允許的范圍內(nèi)適當(dāng)改變其中的參數(shù)值SIGMA、RHO和BETA，再選取不同的初始值y0，運(yùn)行上述的程序，觀察并記錄計(jì)算的結(jié)果有什么特點(diǎn)？是否發(fā)現(xiàn)什么不同的現(xiàn)象？實(shí)驗(yàn)7.2 剛性問(wèn)