數(shù)值分析實(shí)驗(yàn)報(bào)告二

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

2、用秦九韶算法編程計(jì)算在上的值.五、算法描述及實(shí)驗(yàn)步驟針對(duì)內(nèi)容一：（1）先用積分求出（2）根據(jù)遞推公式依次求出的值 （n=13,12,11,.）針對(duì)內(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分別計(jì)算出x=1，2，3，4上的值.6、 調(diào)試過(guò)程及實(shí)驗(yàn)結(jié)果針對(duì)問(wèn)題一：（1） 在命令窗口運(yùn)行程序: 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） 運(yùn)用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）計(jì)算的結(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針對(duì)問(wèn)題二：（1）在命令窗口運(yùn)行程序：a=101:-1:1;（2） 調(diào)

4、用函數(shù)Qinjiu的計(jì)算結(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時(shí)的值分別為： 5151 2.5353e+032 7.7693e+049 2.1569e+0627、 總結(jié) 在整個(gè)實(shí)踐過(guò)程中學(xué)會(huì)了用秦九韶算法編程計(jì)算多項(xiàng)式的值并能夠采用穩(wěn)定的數(shù)值計(jì)算方法比較精確的計(jì)算定積分。學(xué)會(huì)自己發(fā)現(xiàn)錯(cuò)誤解釋錯(cuò)誤糾正錯(cuò)誤。同時(shí)也培養(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=