文档介绍：实验次样条插值实验 样条插值实验报告
下面是我整理的试验次样条插值试验 样条插值试验报告，欢迎参阅。
数值分析试验报告 姓 名 忘川 学 号 系 别 数学系 班级 12 级 主讲老师 指导老师 试验日期condition0 0( ) , ( )i n nS x y S x y ? ? ? ? ? ? . OUTPUT: ? ? n i M i , , 1 , 0 ? ? Step 1 For 1 , , 2 , 1 ? ? n i ? 1111 11 1( 1)6( )iii iii i ii ii i i iii i i ihh hhh hy y y yh h h h?? ? ?????? ?? ???? ? ??? ?? ??
Step 2 By clamped boundary condition, Set 01n? ? ? ? ,1 00 01 16( )y yyh h??? ? ?,16( )n nn nn ny yyh h???? ? ? Step 3 Solve tri-diagonal system 0 0 01 1 1 12 2 2 21 1 1 n n n nn n nMMMMM? ?? ? ?? ? ?? ? ?? ?? ? ? ?? ?? ? ? ?? ?? ? ? ?? ?? ? ? ?? ?? ? ? ??? ?? ? ? ?? ?? ? ? ?? ?? ? ? ?? ?? ? ? ?? ?? ? ? ?? ?? ? ? ? Step 4 Output ? ? n i M i , , 1 , 0 ? ? . Stop 四、试验步骤（或记录） (1)编写三次样条程序求10 ( )S x ，取10 10 10( ) ( )( 0,1, ,10), ( 5) ( 5), (5) (5).k kS x f x k S f S fⅱⅱ = = - = - = L 程序如下：
function [s,m]=selfspline(x0,y0,df,x,conds) n=length(x0); h=diff(x0); b=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))); end a(1:n-2)=a(2:n-1); switch conds case 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(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") end for