版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
1、數(shù)值分析實驗報告0誤差分析誤差分析實驗1.1(問題)實驗?zāi)康模核惴ㄓ小皟?yōu)”與“劣”之分,問題也有“好”與“壞”之別。對數(shù)值方法的研究而言,所謂壞問題就是問題本身對擾動敏感者,反之屬于好問題。通過本實驗可獲得一個初步體會。數(shù)值分析的大部分研究課題中,如線性代數(shù)方程組、矩陣特征值問題、非線性方程及方程組等都存在病態(tài)的問題。病態(tài)問題要通過研究和構(gòu)造特殊的算法來解決,當然一般要付出一些代價(如耗用更多的機器時間、占用更多的存儲空間等)。問題提出
2、:考慮一個高次的代數(shù)多項式)1.1()()20()2)(1()(201????????kkxxxxxp?顯然該多項式的全部根為12…20共計20個,且每個根都是單重的?,F(xiàn)考慮該多項式的一個擾動)2.1(0)(19??xxp?其中是一個非常小的數(shù)。這相當于是對(1.1)中的系數(shù)作一個小的擾動。?19x我們希望比較(1.1)和(1.2)根的差別,從而分析方程(1.1)的解對擾動的敏感性。實驗內(nèi)容:為了實現(xiàn)方便,我們先介紹兩個Matlab函數(shù)
3、:“roots”和“poly”。roots(a)u?其中若變量a存儲n1維的向量,則該函數(shù)的輸出u為一個n維的向量。設(shè)a的元素依次為,則輸出u的各分量是多項式方程121?naaa?01121???????nnnnaxaxaxa?的全部根;而函數(shù)poly(v)b?的輸出b是一個n1維變量,它是以n維變量v的各分量為根的多項式的系數(shù)??梢姟皉oots”和“poly”是兩個互逆的運算函數(shù)。000000001.0?ess)211(zerosve
4、?)2(essve?數(shù)值分析實驗報告2y=poly2sym(a)rr=solve(y)fn=2:21nfm=1:8ess=10^(6m)ve=zeros(121)ve(n)=essa=poly(1:20)vey=poly2sym(a)r=solve(y)6ms=max(abs(rrr))endend數(shù)值實驗結(jié)果及分析:數(shù)值實驗結(jié)果及分析:fmatlong6mn7891022.797226874783311.867536320091581
5、.060527623807480.2527314421904731.693766997674240.923106667069640.084716145697410.4080402640941140.854013934155360.199410220200610.03972935295834050.110311005388710.0429653236284400600007000080000900001000001100001200001
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 眾賞文庫僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 《數(shù)值分析》上機實驗報告
- 數(shù)值分析實驗報告.doc
- 實驗報告范文5篇
- 矩陣與數(shù)值分析實驗報告
- 《數(shù)值分析》實驗報告書
- 清華大學數(shù)值分析實驗報告
- 清華大學數(shù)值分析實驗報告
- 最優(yōu)化數(shù)值實驗報告三
- 數(shù)值計算方法實驗報告
- 數(shù)值計算方法實驗報告
- html5實驗報告
- 數(shù)值分析實驗報告 jacobi迭代和seidel迭代分析資料
- 數(shù)值分析課程設(shè)計實驗報告基本格式.doc
- 需求分析實驗報告
- 實驗報告—差熱分析
- 操作系統(tǒng)實驗報告 (5)
- 整理全站儀閉合導線測量實驗報告(5篇)
- 資本結(jié)構(gòu)分析實驗報告
- 判別分析實驗報告
- 隨機信號分析實驗報告
評論
0/150
提交評論