數(shù)值分析實(shí)驗(yàn)報告程序?qū)嶒?yàn)

線性方程求while(R>5e-6);if

x=c【NewtonNewton 輸入函f=input('請輸入需要求解函數(shù)%%%求解f(x)的導(dǎo)%%%改進(jìn)常數(shù)或重根%%%初始值x0=input('inputinitialvaluex0>>');max=100;%最大迭代次R=eval(subs(f,'x0','x'));%求解f(x0)，以確定初值x0while(abs(R)>1e-if(eval(subs(f,'x0','x'))<1e-10);

ifk>max;%如果迭代次數(shù)大于給定值，認(rèn)為迭代不收斂，重ss=input('mayberesultiserror,chooseanewx0,y/n?>>','s');ifx0=input('inputinitialvaluex0>>');

【結(jié)果與】1題：實(shí)驗(yàn)二Guass列主元消去A=[afork=1:n-1;%%%是否進(jìn)行主元選ifabs(A(k,k))<yibuseilou;%事先給定的認(rèn)為有必要選主元的小數(shù)

if%%%%選主forr=k+1:n;elsep=k;%%%交換元ifforq=k:n+1;%%%判斷系數(shù)矩陣是否奇異或非常嚴(yán) abs(A(k,k))<1e-forr=k+1:n;forq=k:n+1;

fork=n-1:-1:1;for【結(jié)果與】1題：Guassx--Guassx--x采用Gauss消去法時，如果在消元時對角線上的元素始終較大（假如大于8Runge現(xiàn)象的產(chǎn)生和克本實(shí)驗(yàn)采用的M文件編寫。其中待插值的方程寫成function的方式，functiony=f(x);Lagrange插值源程序n=input('將區(qū)間分為的等份數(shù)輸入s=[-1+2/n*[0:n]];%%%給定的定點(diǎn)，Rf為給定的函x=-forfork=1:n+1;iff=f+Rf(s(q))*l;%求插值函gridonhold分段線性插值源n=input('將區(qū)間分為的等份數(shù)輸入s=[-1+2/n*[0:n]];%%%給定的定點(diǎn)，Rf為給定的函forx=-1:hh:1;fork=1:n+1;%%%switchcaseifl=(x-s(2))./(s(1)-casen+1ififx>=s(k-elseifx>=s(k)&x<=s(k+1);

ff=ff+Rf(s(k))*l;%%求插值函數(shù)gridonhold三次樣條插值源程序（采用第一邊界條件n=input('將區(qū)間分為的等份數(shù)輸入%%%插值區(qū)hh=0.001;%s=[a+(b-a)/n*[0:n]];%%%給定的定點(diǎn)，Rf為給定的函%%%%第一邊界條件Rf"(-1),Rf"(1)fork=1:n;%取出節(jié)點(diǎn)間距h(k)=s(k+1)-fork=1:n-1;%求出系數(shù)向量lamuda,miu%%%%賦值系數(shù)矩陣fork=1:n-forp=1:n-switchpcasekcasek-casek+1%%%%求出dfork=1:n-switchkcase1d(k)=6*f2c([s(k)s(k+1)s(k+2)])-miu(k)*v;casen-1d(k)=6*f2c([s(k)s(k+1)s(k+2)])-la(k)*v;d(k)=6*f2c([s(k)s(k+1)%%%%求解M陣forif

p=ceil((x-s(1))/((b-

ff1=1/h(p)*(s(p+1)-ff2=1/h(p)*(x- holdongrid【結(jié)果與】f(x)

125x

xj=-1+jh,h=0.1，j=0,…，n。下面分別給出Lagrange插值，三次樣條插值，其它曲線沒有標(biāo)出，從數(shù)據(jù)表中可以看出具體的誤差。表中，L10(x)為LagrangeS10(x)S40(x)n=1040X10(x)，X40(x)分別代表n=10，40的線性分段插值函數(shù)。x-1000000-0010000-0010000-0000000-0000000-00-00000-00-00000-00-00000-0000000-0000000-0000000-0000000-0000000-0000000-0000000-0000000-0000000-0000000-0000000-0000000011111100000000000000000000000000000000000000000000000000000000000000000000000000000000000000-0000000-0000000-0