數值分析實驗報告二

1、精選優(yōu)質文檔-傾情為你奉上貴州師范大學數學與計算機科學學院學生實驗報告課程名稱： 數值分析 班級 數學與應用數學2班 實驗日期：2013年9月19日學 號： 姓名 指導教師： 實驗成績： 一、實驗名稱 實驗一:遞推法的穩(wěn)定性,秦九韶算法二、實驗目的及要求 1.熟悉數值穩(wěn)定的概念, 通過上機計算，了解舍入誤差所引起的數值不穩(wěn)定性. 2. 培養(yǎng)Matlab編程與上機調試能力.三、實驗環(huán)境每人一臺計算機，要求安裝Windows XP操作系統，Microsoft office2003、MATLAB6.5(或7.0).四、實驗內容 1.教材例1.13中,取15位數字計算,并分析、比較計算結果. 2.設,

2、用秦九韶算法編程計算在上的值.五、算法描述及實驗步驟針對內容一：（1）先用積分求出（2）根據遞推公式依次求出的值 （n=13,12,11,.）針對內容二：（1） 首先在Matlab里編制函數式M文件，函數名為Qinjiu(a,n,x).function y=Qinjiu(a,n,x)v=a(1);for k=1:n v=v.*x+a(k+1);endy=v;（2）然后調用函數Qinjiu分別計算出x=1，2，3，4上的值.6、 調試過程及實驗結果針對問題一：（1） 在命令窗口運行程序: int(x.14/(x+5),0,1)ans = *log(2)+*log(3)-*log(5)-047/7

3、2072 format long ans = 1.8359e-002所以 = 0.5684 （取15位小數后的結果）（2） 運用for循環(huán)算出n=13,12,.0的值 int n; s= 1.8359e-002 int n;for(n=14;0)s=1./(5*n)-s./5printf(%d,s)n-;（3）計算的結果為： （均保留15位小數）0.7058 0.05000.0567 0.37050.5259 0.11700.7766 0.18750.4958 0.50080.8998 0.88670.0226 0.9396針對問題二：（1）在命令窗口運行程序：a=101:-1:1;（2） 調

4、用函數Qinjiu的計算結果為： Qinjiu(a,100,1)ans = 5151 Qinjiu(a,100,2)ans = 2.5353e+032 Qinjiu(a,100,3)ans =7.7693e+049 Qinjiu(a,100,4)ans = 2.1569e+062所以x=1,2,3,4時的值分別為： 5151 2.5353e+032 7.7693e+049 2.1569e+0627、 總結 在整個實踐過程中學會了用秦九韶算法編程計算多項式的值并能夠采用穩(wěn)定的數值計算方法比較精確的計算定積分。學會自己發(fā)現錯誤解釋錯誤糾正錯誤。同時也培養(yǎng)了自己的編程能力以及上機調試能力。8、 附錄（源程序清單）附錄1：M文件 function y=Qinjiu(a,n,x)v=a(1);for k=1:n v=v.*x+a(k+1);endy=v; int n;附錄2：遞推公式求值 s=