實(shí)驗(yàn)次樣條插值實(shí)驗(yàn)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上數(shù)值分析實(shí)驗(yàn)報(bào)告姓 名忘川學(xué) 號系 別數(shù)學(xué)系班級12級 主講教師指導(dǎo)教師實(shí)驗(yàn)日期2014/6/25專業(yè)信息與計(jì)算科學(xué)專業(yè) 課程名稱數(shù)值分析同組實(shí)驗(yàn)者無一、實(shí)驗(yàn)名稱： 實(shí)驗(yàn)四、三次樣條插值實(shí)驗(yàn)二、 實(shí)驗(yàn)?zāi)康模?掌握三次樣條插值的運(yùn)用；2了解拉格朗日插值在高次上的誤差。三、實(shí)驗(yàn)內(nèi)容及要求：給定函數(shù)節(jié)點(diǎn)(1)編寫三次樣條程序求，取 (2)作原函數(shù)、Langrage插值函數(shù)和三次樣條插值函數(shù)的圖像，并比較它們的區(qū)別。附：算法描述：三次樣條插值（Cubic Spline Interpolation）PURPOSE: To find a piecewise cubic splin

2、e function . (, ).where , and ,. INPUT: interpolated points; clamped boundary condition.OUTPUT: Step 1 For Step 2 By clamped boundary condition, Set ,Step 3 Solve tri-diagonal systemStep 4 Output . Stop四、實(shí)驗(yàn)步驟（或記錄）(1)編寫三次樣條程序求，取 程序如下：function s,m=selfspline(x0,y0,df,x,conds)n=length(x0);h=diff(x0);b=

3、ones(1,n)*2;for j=2:n-1 a(j)=h(j-1)/(h(j-1)+h(j); c(j)=h(j)/(h(j-1)+h(j); d(j)=6*(h(j-1)*(y0(j+1)-y0(j)-h(j)*(y0(j)-y0(j-1)/(h(j)*h(j-1)*(h(j)+h(j-1);enda(1:n-2)=a(2:n-1);switch condscase 1 a(n-1)=1; c(1)=1; d(1)=6*(y0(2)-y0(1)/(x0(2)-x0(1)-df(1)/(x0(2)-x0(1); d(n)=6*(df(2)-(y0(n)-y0(n-1)/(x0(n)-x0(

4、n-1)/(x0(n)-x0(n-1); case 2 a(n-1)=0; c(1)=0; d(1)=2*df(1); d(n)=2*df(2);otherwise error('conds值錯(cuò)誤， conds只能為1或2')endfor i=2:n r=a(i-1)/b(i-1); b(i)=b(i)-r*c(i-1); d(i)=d(i)-r*d(i-1);endm(n)=d(n)/b(n);for i=n-1:-1:1 m(i)=(d(i)-c(i)*m(i+1)/b(i);endfor j=1:length(x) for i=1:n-1 if x(j)>=x0(i

5、)&x(j)<=x0(i+1) s(j)=m(i)*(x0(i+1)-x(j)3/(6*h(i)+m(i+1)*(x(j)-x0(i)3/(6*h(i)+(y0(i)-m(i)*h(i)2/6)*(x0(i+1)-x(j)/h(i)+(y0(i+1)-m(i+1)*h(i)2/6)*(x(j)-x0(i)/h(i); end endend在matlab的命令窗口輸入：>> k=0:10;>> x0=-5+k;>> y0=1./(x.2+1);>> x=-5：0.5:5;>> y=1./(x.2+1);>> d

6、f=diff(y)得到如下數(shù)據(jù)：df = 0.0086 0.0118 0.0166 0.0245 0.0379 0.0621 0.1077 0.1923 0.3000 0.2000 -0.2000 -0.3000 -0.1923 -0.1077 -0.0621 -0.0379 -0.0245 -0.0166 -0.0118 -0.0086(2)作原函數(shù)、Langrage插值函數(shù)和三次樣條插值函數(shù)的圖像，并比較它們的區(qū)別。先寫拉格朗日插值的程序，如下：function y=lagr1(x0,y0,x)n=length(x0); m=length(x);for i=1:m z=x(i); s=0.

7、0; for k=1:n p=1.0; for j=1:n if j=k p=p*(z-x0(j)/(x0(k)-x0(j); end end s=p*y0(k)+s; end y(i)=s;endend然后直接在matlab的命令窗口輸入：x0=-5 -4 -3 -2 -1 0 1 2 3 4 5 ;y0=1./(x0.2+1);x=-5:0.01:5;df=diff(y0);conds=1;y=selfspline(x0,y0,df,x,conds);y1=1./(x.2+1);y2=lagr1(x0,y0,x);plot(x,y1,'b',x,y,'y-',x,y2,'r:')legend('原函數(shù)','三次樣條插值','Lagrange插值')即可得到以下圖形：(3) 分析結(jié)果：根據(jù)圖像可知，用拉格朗日插值函數(shù)進(jìn)