數(shù)值分析第二次程序題-資料下載頁
2025-03-25 02:50本頁面
【正文】 y(m(i)+1)y(m(i)))/(x(m(i)+1)x(m(i)))*(z(i)x(m(i)))。endl(2001)=y(201)。f(2001)=y(201)。v=linspace(1,0,100)。u=sin(pi*v)。plot(v,u,39。k39。)hold onv=linspace(0,50)。u=cos(pi*v)。plot(v,u,39。k39。)hold onv=linspace(,1,50)。u=0。plot(v,u,39。k39。)hold onplot(z,l,39。r39。)三次樣條插值x=1::1。n=length(x)。for i=1:21 if x(i)0 y(i)=sin(pi*x(i))。 elseif x(i) y(i)=0。 else y(i)=cos(pi*x(i))。 endend for i=1:n1h(i)=x(i+1)x(i)。endfor i=1:n2u(i)=h(i)/(h(i+1)+h(i))。r(i)=1u(i)。endG=zeros(n2,n2)。for i=1:n2 G(i,i)=2。endfor i=2:n2 G(i,i1)=u(i1)。 G(i,i+1)=r(i1)。endd=zeros(1,n2)。for i=1:n2 d(i)=6*((y(i+2)y(i+1))/h(i+1)(y(i+1)y(i))/h(i))/(h(i+1)+h(i))。endd=d39。M=G\d。M=[0。M。0]。for i=1:n1 z=[x(i)::x(i+1)]。 m=length(z)。 for j=1:m s(j)=M(i)*(x(i+1)z(j))^3/+M(i+1)*(z(j)x(i))^3/+(y(i)M(i)*)*(x(i+1)z(j))/+(y(i+1)M(i+1)*)*(z(j)x(i))/。 end plot(z,s,39。* r39。,39。MarkerSize39。,3) hold onendv=linspace(1,0,100)。u=sin(pi*v)。plot(v,u,39。k39。)hold onv=linspace(0,50)。u=cos(pi*v)。plot(v,u,39。k39。)hold onv=linspace(,1,50)。u=0。plot(v,u,39。k39。)hold on
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號-1