資源描述：
《數(shù)值分析實(shí)驗(yàn)new》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、《數(shù)值分析》計算實(shí)習(xí)報告冊專業(yè)信息與計算科學(xué)學(xué)號20112984姓名張妃___2012~2013年第一學(xué)期實(shí)驗(yàn)一數(shù)值計算的工具M(jìn)atlab1.解釋下MATLAB程序的輸出結(jié)果程序：t=0.1n=1:10e=n/10-n*t答案：t=0.1000小數(shù)點(diǎn)后輸出四位n=12345678910n取1到10e=1.0e-015*00-0.055500-0.1110-0.1110000e等于-0.0555×2.下面MATLAB程序的的功能是什么？程序：x=1;while1+x>1,x=x/2,pause(0.02),endx=1;whilex+x>x,x=2*x,pause(
2、0.02),endx=1;whilex+x>x,x=x/2,pause(0.02),end答案：當(dāng)x=1時，對x無限二分，到1+x>1結(jié)束，且循環(huán)之間間隔0.02秒當(dāng)x=1時，對無限循環(huán)，到x+x>x結(jié)束，且循環(huán)之間間隔0.02秒當(dāng)x=1時，對x無限二分，到x+x>x結(jié)束，且循環(huán)之間間隔0.02秒3.考慮下面二次代數(shù)方程的求解問題公式是熟知的，與之等價地有，對于，應(yīng)當(dāng)如何選擇算法。答案：clc,cleara=1;b=100000000;c=1;x1=(-b+sqrt(b*b-4*a*c))/2*ax2=(-b+sqrt(b*b-4*a*c))/2*ax3=2*c/
3、(-b+sqrt(b*b-4*a*c))x4=2*c/(-b+sqrt(b*b-4*a*c))d=x1*x1+100000000*x1+1e=x2*x2+100000000*x2+1f=x3*x3+100000000*x3+1g=x4*x4+100000000*x4+1x1=-7.4506e-009=-7.4506×x2=-7.4506e-009x3=-134217728x4=-134217728d=0.25494194030762e=0.25494194030762f=4.592625709481985e+015g=4.592625709481985e+015因
4、為d，e比f，g更接近0所以公式更好。4.函數(shù)有冪級數(shù)展開利用冪級數(shù)計算的MATLAB程序?yàn)閒unctions=powersin(x)s=0;t=x;n=1;whiles+t~=s;s=s+t；t=-x^2/（(n+1)*(n+2)）*t；n=n+2；end修改后程序functions=powersin(x)s=0;x=pi/2;t=x;k=0;n=1;whiles+t~=s;s=s+t;t=-x^2/((n+1)*(n+2))*t;n=n+2;k=k+1;endK結(jié)果X=21*pi/2k=60ans=0.9999X=11*pi/2k=37ans=-1.0000X
5、=pi/2k=11ans=1.0000(a)解釋上述程序的終止準(zhǔn)則。t=0時，終止(b)對于計算的進(jìn)度是多少？分別計算多少項(xiàng)？ans=1.0000ans=-1.0000ans=0.999911,37,605.考慮調(diào)和級數(shù)，它是微積分中的發(fā)散級數(shù)，在計算機(jī)上計算該級數(shù)的部分和，會得到怎么樣的結(jié)果，為什么？clc;clears=0;n=20;fori=n-1:n+10000000000t=1/i;s=s+t;endS結(jié)果s=18.5697fori=n-1:n+100000s=8.5952fori=n-1:n+10000000000000s=18.5697S會隨著i的項(xiàng)
6、數(shù)的增大而增大，到一定程度后不變。因?yàn)樵谟嬎銠C(jī)中如果一個數(shù)很小，計算機(jī)會默認(rèn)為0.6.指數(shù)函數(shù)的級數(shù)展開，如果對于，用上述的級數(shù)近似計算指數(shù)函數(shù)的值，這樣的算法結(jié)果是否會好，為什么？1.clc,clears=1;t=1;x=-1;fori=1:100000t=t*i;tt=x^i/t;s=s+tt;endsX=1;s=2.7183x=-1;s=0.36792.functions=powersin(x)s=0;m=1;t=1;n=1;whiles+t~=s;s=s+t;m=m*n;t=-x^(n)/m;n=n+2;endx=0,f=1;x=1,y=-0.4107;x
7、=-1,f=2,41077.考慮數(shù)列，它的統(tǒng)計平均定義為，標(biāo)準(zhǔn)差數(shù)學(xué)上等價于作為標(biāo)準(zhǔn)差的兩種算法，你將如何評價他們的得與失。1.clc,cleara=input('請輸入數(shù)據(jù)')n=length(a);sm=sum(a);ave=sm/n;b1=0;fori=1:nb1=(a(i)-ave)^2+b1;endb1b2=((1/(n-1))*b1)^(1/2)結(jié)果b1=133.9368b2=11.57312.clc,cleara=input('請輸入數(shù)據(jù)')n=length(a);sm=sum(a);ave=sm/n;b3=0;fori=1:nb3=(a(i))^2
8、+b3;e
