矩陣與數(shù)值分析上機實習(xí)設(shè),.,,,:=%N=input(anumber(N>=2)')if<number')S=j=S=S+'S:')aN=aN=aN=:=%eps('single')N=input(anumber(N>=2)')if<number')S=j=S=S+'S:')aN==aN==aN=n.():A=input('請輸入系數(shù),'n=input('請輸入計算變量的值'====i==n+結(jié)果:請輸入系數(shù)，由高次冪開始[73-511]請輸入計算變量的值23ans=85169，.:0000000000000000000000000000000000000000000000000000BB=n==x=solution')i=j=i+n=k=*+;i=in=disp('solution')結(jié)果:-:%0000000000000000000000000000000000000000000000000000BB=n==x=solution')i=%%==j=if>==>ij====j=i+n=k=*+;i=in=disp('solution'):-3和AA::0000000000000000000000000000000000000000000000000000BB=n==x=L=solution')i=j=i+n=k=*+;=-*i=in=disp('C')disp('L'):L0000000000000000000000000000000000000000000000000-0-00000000U000-000000-00000000-000000000-00-0000000000000000000000L=000000000000000000000000000000000000000000000000-0000000000U=000-000000-00000000-000000000-00000000000000000000000000P=100000010000001000000100000010000001000000000000000000000000000000100100000000001求A=YY和Xi==i====A的,1110702116B==L=x=i=i+=-=-=i=i1=/=i=1i==-%disp('x')結(jié)果:x，:A=2003100253006B=2B==c=a=b=f=y=%LU==i======-%LUf=f=f=i=-disp('y')=f=i==-disp('f')=f=y:f1111.A的XY);結(jié)果:C=1.0000-1.0000-2.0000-2.00001.00003.00002.00002.00000-0.50001.500000.50000.50002.5000-0.5000>>[Q,R]=qrhouseholder(C)Q=-0.31620.31620.63250.6325-0.7071-0.7071-0.0000-0.0000-0.25820.2582-0.77460.51640.5774-0.5774-0.00000.5774R=-3.16230.00000.0000-0.00003.1623-2.82840.00000.47430.35362.0555-0.35361.03281.1547-1.5492-0.0000-0.00008.分別應(yīng)用Jacobi迭代法和Gauss-Seidel迭代法求解如下方程組8>?x+x=7<123++x=?>123?+x+=:123Jacobi代碼A=1148121B=B==x==disp('x')n=input('請輸入迭代次數(shù)'n1i=j===x='')disp('x')請輸入迭代次數(shù):10x0.66231.99492.3284Gauss-seidel代碼A=1148121B=B==x=disp('x')n=input(iterations:'i=j===disp('x')結(jié)果x0.66672.00002.33339.分別應(yīng)用Newton12x3?5x?1=0在[1,2]2）xsx=在[-4,-3]上的一個根。Newton迭代法（1）=-6-==1=if=')break;ifdisp('x1')結(jié)果x11.6730（2）=-6-=+==1=%==-ifdisp('x1')結(jié)果：x1-3.1416Newton截斷法（1）代碼=+=-6-===1=-=-break;==disp('x2'結(jié)果x21.6730（2）=+=1=-break;==disp('x2'結(jié)果x2-3.141610.采用二分法計算非線性方程xcosx?2=，查找區(qū)間為[-4,-2]代碼:-2=%a=b===abnumber')===b*>0)=0m=<==<0b=a=disp('root')結(jié)果:root-2.4988e+00011.已知函數(shù)，在[-5,5]上分別取2，插值節(jié)點，用Lagrange方法插值，并把原函數(shù)圖與插值函數(shù)圖比較，觀察插值效果代碼:%%+=input(step:'a=b=n==j==a+*='x')i=v=j=j%%j==*=+disp('result')結(jié)果x135x0123451919x1-0y11--::212.用三次樣條插值上題中的插值節(jié)點，并畫圖比較插值效果。（提示：原函數(shù)在兩個端點-5,5的導(dǎo)數(shù)值可作為邊界條件）代碼:21]=input(step:'a=b=n=disp('n')====i==a+=+='complete'j==='x''y'on;on;disp('pp.coefs')結(jié)果n50nn13.已知，分別用復(fù)合梯形法和復(fù)合Simpson公式計算積分，區(qū)間分為20，40，80，200個小區(qū)間，并計算其精確值，比較計算精度情況。精確結(jié)果為：0復(fù)合梯形法%t13--?%204080200個小=%=T=a=b=i=h=j=-=a+h=+=++'len')'T')len20T7.9936e-016len40T7.9936e-016len80T9.5479e-016len100T5.3291e-017復(fù)合simpon代碼=T=a=b=i=h==j==a+h*=a+h==+=+'xk')'xk1')'xk12')=++a++a+=++'len')'T')結(jié)果len20T6.8094e-016len40T8.2897e-016len80T7.1794e-016len100T5.3291e-017>>14.用2點，3點和5點Gauss積分法分別計算定積分，并與真值作比較。xfun=fcnchk(fun,'vectorize'disp('ql')結(jié)果::=2:>>fun=@(x)(x.^2.*cos(x));>>clear>>fun=@(x)(x.^2.*cos(x));>>a=0;>>b=pi/2;>>n=1;>>tol=0.00005;>>[ql,Ak,xk]=guasslegendre(fun,a,b,n,tol)ql4.7464e-001ql=4.7464e-0013點:>>fun=@(x)(x.^2.*cos(x));clearfun=@(x)(x.^2.*cos(x));a=0;b=pi/2;n=2;tol=0.00005;>>[ql,Ak,xk]=guasslegendre(fun,a,b,n,tol)ql4.6724e-001ql=4.6724e-0015點:>>fun=@(x)(x.^2.*cos(x));clearfun=@(x)(x.^2.*cos(x));a=0;b=pi/2;n=4;tol=0.00005;>>[ql,Ak,xk]=guasslegendre(fun,a,b,n,tol)ql4.6740e-001ql=4.6740e-001，畫法=+=0%a=b===step')n=anumber')N=n<0n>Nnumber')j==a+*=++=len(i)')u')步長：0.1lengthofstep0.1len=1.0000e-001pleaseanumber10n=10len(i)1.0000e-001u1.2260e+004>>clearlengthofstep0.05len=5.0000e-002pleaseanumber20n=20len(i)5.0000e-002u5.8700e+007lengthofstep0.01len=1.0000e-002pleaseanumber90n=90len(i)1.0000e-002u4.3690e+034>>法=+=0%a=b===step')n=a