數(shù)值分析上機實習報告(編輯修改稿)
2025-02-14 21:49 本頁面
【文章內容簡介】 Q=zeros(2,1)。x=[1。2。3。4。5。6。7。8。9。10]。y=[0。。]。h=zeros(10,1)。d=zeros(10,1)。u=zeros(10,1)。v=zeros(10,1)。r=zeros(10,1)。l=zeros(10,1)。z=zeros(10,1)。m=zeros(10,1)。for t=1:1:9。h(t)=x(t+1)x(t)。endd(1)=6/h(1)*((y(2)y(1))/h(1)1)。d(10)=6/h(9)*((y(10)y(9))/h(9))。for t=1:1:8u(t+1)=h(t)/(h(t)+h(t+1))。v(t+1)=1u(t+1)。d(t+1)=6/(h(t)+h(t+1))*((y(t+2)y(t+1))/(x(t+2)x(t+1))(y(t+1)y(t))/(x(t+1)x(t)))。endu(10)=1。v(1)=1。r(1)=d(1)。for t=2:1:10l(t)=u(t)/r(t1)。r(t)=d(t)l(t)*v(t1)。endz(1)=d(1)。for t=2:1:10z(t)=d(t)l(t)*z(t1)。endm(10)=z(10)/r(10)。for t=9:1:1m(t)=(z(t)v(t)*m(t+1))/r(t)。endfor t=1:1:10if p=tamp。amp。p(t+1)Q(:,1)=feval(ssss,p,t,x,m,h,y)。breakendendfunction Q=ssss(p,t,x,m,h,y)Q=zeros(2,1)。Q(1,1)=((power((x(t+1)p),3)*m(t)+power((px(t)),3)*m(t+1))/6+(y(t)m(t)*h(t)*h(t)/6)*(x(t+1)p)+(y(t+1)m(t+1)*h(t)*h(t)/6)*(px(t)))/h(t)。Q(2,1)=((power((x(t+1)p),2)*m(t)+power((px(t)),2)*m(t+1))/2+(y(t)m(t)*h(t)*h(t)/6)+(y(t+1)m(t+1)*h(t)*h(t)/6))/h(t)。end，三對角陣需要滿足：(i=1,2,…,n)均非奇異，保證有唯一的Doolittle分解?！?。, 三次樣條插值的解存在且唯一。 第三題（允許誤差ε=）。 數(shù)值積分的Romberg算法計算步驟如下: 數(shù)值分析上機實習報告 2012年12月 當時，就停機。includeincludedefine N 9float f(float x) /*定義函數(shù)f(x)*/{ float y。 y=pow(3,x)*pow(x,)*(5*x+7)*sin(x*x)。 return(y)。} main(){ float T[N+1][N+1],h[N+1],a=1,b=3,m[N+1]。 int i,l。 T[1][0]=(ba)*(f(a)+f(b))/2。 l=1。 while(l=N) { m[l]=0。 for(i=1。i=(pow(2,l1))。i++) m[l]+=f(a+(2*i1)*(ba)/pow(2,l))。 T[1][l]=(T[1][l1]+(ba)*m[l]/pow(2,l1))/2。 l++。 } i=1。 while(i=N) { for(l=1。l=Ni+1。l++) T[i+1][l1]=(pow(4,i)*T[i][l]T[i][l1])/(pow(4,i)1)。 h[i]=T[i][0]T[i+1][0]。 if(fabs(h[i])=1e5) break。 i++。 } printf(The answer is:%f,T[i+1][0])。} MATLAB上機程序function [T,n]=mromb(f,a,b,eps)if nargin4,eps=1e6。endh=ba。R(1,1)=(h/2)*(feval(f,a)+feval(f,b))。n=1。J=0。err=1。while (erreps) J=J+1。h=h/2。S=0。 for i=1:n x=a+h*(2*i1)。 S=S+feval(f,x)
教學課件相關推薦
鄂ICP備17016276號-1