實驗次樣條插值實驗

1、精選優(yōu)質文檔-傾情為你奉上數值分析實驗報告姓 名忘川學 號系 別數學系班級12級 主講教師指導教師實驗日期2014/6/25專業(yè)信息與計算科學專業(yè) 課程名稱數值分析同組實驗者無一、實驗名稱： 實驗四、三次樣條插值實驗二、 實驗目的：1掌握三次樣條插值的運用；2了解拉格朗日插值在高次上的誤差。三、實驗內容及要求：給定函數節(jié)點(1)編寫三次樣條程序求，取 (2)作原函數、Langrage插值函數和三次樣條插值函數的圖像，并比較它們的區(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四、實驗步驟（或記錄）(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值錯誤， 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)得到如下數據：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)作原函數、Langrage插值函數和三次樣條插值函數的圖像，并比較它們的區(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('原函數','三次樣條插值','Lagrange插值')即可得到以下圖形：(3) 分析結果：根據圖像可知，用拉格朗日插值函數進