資源描述：
《matlab上機(jī)實(shí)驗(yàn)報(bào)告》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫。
1、上絳電力嗲紇數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)報(bào)告院系專業(yè)年級(jí)學(xué)生姓名學(xué)號(hào)指導(dǎo)教師能源與機(jī)械丁.程學(xué)院動(dòng)力機(jī)械及丁?程2012級(jí)張亞杰VS1210124014黃建雄2012年12月26日數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)題(1)由遞推公式八=-5/^+1，從/0的幾個(gè)近似值出發(fā)，計(jì)算/2();nr1解：1{)=<=0.1823Jo5+x計(jì)算120編輯matlab命令如下：1=0.1823forn=l:1:20,I=-5*I+l/n；fprintf('%.Id%.4f',n,I);end10.088520.057530.045840.020850.09586-0.312571.70548-8.401
2、8942.120010-210.5002111052.591912-5262.87621326314.457814-131572.217515657861.154216-3289305.70831716446528.600218-82232642.945519411163214.780320-2055816073.8513(2)粗糙估計(jì)/，。，用—丄八+丄，計(jì)算/。；55/7解：l20=—5+x使川復(fù)合中點(diǎn)公式進(jìn)行積分，相應(yīng)的matlab程序如下:1=0;forh=0:0.001:1,m=h+0.0005;1=1+0.001*m^20/(5+m)；fprintf(1%.Id%
3、.4f1,m,I)；enddisp(I)；fork=l:20,n=21-k；1=0.2*(1/n-I)；fprintf('%.Id%.4f'zn,I)；enddisp(I)結(jié)果：ae-uiU.UU&l8e-0l0.00528e-0l0.00548e-0l0.0055200.00848e-0l0.0056190.00898e-0l0.0057180,00938e-0l0.0058170.00999e-0l0.0060160.01059e-0l0.0061150.01129e-0l0.0062140.01209e-0l0.0064130?01309e-0l0.006512
4、0.01419e-0l0.0066110.01549e-0l0.0068100.01699e-0l0.006990.01889e-0l0.007180.02129e-0l0.0072—0.02430e+000.007460.02850e+000.007550.03430e+000.007740.04310e+000.007830.05800e+000.008020.08840e+000.008210.18230.00820.1823程序結(jié)朿時(shí)輸出兩個(gè)I值，第一個(gè)表示120,笫二個(gè)表示分別為12()=0.00821()=0.1823(1)分析結(jié)果的可靠性及產(chǎn)生此現(xiàn)象的原因(重點(diǎn)
5、分析原因)。從上述計(jì)算屮分析得到如果先得到10，再?gòu)腎Q由遞推公式得到I2o,I2e結(jié)果跟精確值相比誤差很人；如果先估算I2C,在從ho冇遞推公式彳?到Ie，Ie的結(jié)果跟精確值相比近似相等。原因分析：如果從10推I2Q的近似值，需要用到遞推公式In=-5Ind+l/n,10本身結(jié)果是有誤差的：經(jīng)過遞推公式計(jì)算20次，就等于誤差被認(rèn)為的放大5的20次方倍，所以得到的I2Q與其精確值相差其遠(yuǎn)。如果從I2o推Io的近似值,需要用到Ind=0.2(l/n-In),盡管I2Q本身有誤差，但是經(jīng)過20次運(yùn)算，其誤差縮小到原來的0.2的20次方估，所以得到的Io與其精確值比較相近。2.求方
6、程eA+10x—2=0的近似根，要求
7、人+1-a
8、<5x10'并比較計(jì)算量。(1)在［0,1］上用二分法；Matlab程序如下:a=0;b=l；c=b-a；n=0whilec>0.0005,x=(a+b)/2；f=exp(x)+10*x-2；iff>0,b=x；c=b-a;elseiff<0,a=x；c=b-a；elsex=x；c=0；endn=n+l；fprintf(1%.Id%.4fend%.4f1,n,x,c);結(jié)果如下：10.50000.500020.25000.250030.12500.125040.06250.062550.09380.031360.07810
9、.0156i0.08590.007880.08980.003990.09180.0020100.09080.0010110.09030.0005解得到;x=0.0903(2)取初值x0=0,2-ex10采川matlab進(jìn)行迭代的程序如卜‘：x=0;c=l;n=0；whilec>0.0005,m=x；m=(2-exp(m))/10；c=abs(m-x);x=m；n=n+l；fprintf(1%.Id%.4f%.4f',n,x,c)；end(3)10.10000.100020.08950.010530.
