數(shù)值分析實驗報告二

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

2、用秦九韶算法編程計算在上的值.五、算法描述及實驗步驟針對內(nèi)容一：（1）先用積分求出（2）根據(jù)遞推公式依次求出的值 （n=13,12,11,.）針對內(nèi)容二：（1） 首先在Matlab里編制函數(shù)式M文件，函數(shù)名為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）然后調(diào)用函數(shù)Qinjiu分別計算出x=1，2，3，4上的值.6、 調(diào)試過程及實驗結(jié)果針對問題一：（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位小數(shù)后的結(jié)果）（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）計算的結(jié)果為： （均保留15位小數(shù)）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） 調(diào)

4、用函數(shù)Qinjiu的計算結(jié)果為： 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、 總結(jié) 在整個實踐過程中學(xué)會了用秦九韶算法編程計算多項式的值并能夠采用穩(wěn)定的數(shù)值計算方法比較精確的計算定積分。學(xué)會自己發(fā)現(xiàn)錯誤解釋錯誤糾正錯誤。同時也培養(yǎng)了自己的編程能力以及上機(jī)調(diào)試能力。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=