matlab差值后求斜率,用三次样条插值求离散点斜率matlab程序.pdfVIP

matlab差值后求斜率,⽤三次样条插值求离散点斜率matlab程

序

三次样条插值的matlabt程序

:_0{i)^*K,b-i:tfunctionx=followup(a,b,c,d)n=length(d);*i0s$Gt6T$?8a7m

a(1)=0;7Q8c5R.`/N#ba(1)=0;7Q8c5R.`/N#b

L(1)=b(1);7B%~(G/X+G:U3G5z4}

y(1)=d(1)/L(1);

$|#`/f+{)I3eu(1)=c(1)/L(1);8?8I!W8V2Sh%f|.n7I)j6c

fori=2

n-1)

0H;K]6O,y$]$MA$W6rL(i)=b(i)-a(i)*u(i-1);Z#Z8`B.L;B

y(i)=(d(i)-y(i-1)*a(i))/L(i);5dY)f.E/F+c`5H;_;l

u(i)=c(i)/L(i);u(i)=c(i)/L(i);

L(n)=b(n)-a(n)*u(n-1);4G%A+^-O$n:z9O

y(n)=(d(n)-y(n-1)*a(n))/L(n);/B5l*p:p$Y%X

%“赶”的过程0}$C5u-_(b0],t%“赶”的过程0}$C5u-_(b0],t

fori=(n-1):-1:1

5O1K/}9l0h6@(?/AP$mx(i)=y(i)-u(i)*x(i+1);;T3h#w.m8w#L,?5_!o+H!B

endendendend

symst

*}/@8C3H1X*}/@8C3H1X*}/@8C3H1Xg6j;az9En=length(x);g6j;az9En=length(x);g6j;az9En=length(x);

g2K#G;B},t(?h(i)=x(i+1)-x(i);

4i(p;P5_0D+bendo(n0oj3EcX

fori=2:n-1;

8JM)J8[*_+O!Wlamda(i)=h(i)/(h(i-1)+h(i));

0~o9[,o1m1Dmiu(i)=1-lamda(i);

3@-k5_3M!e8[k%~0F~*rg(i)=3*(lamda(i)*((y(i)-y(i-1))/h(i-1))+miu(i)*((y(i+1)-y(i))/h(i)));

.?$b!r0|3t(rend.?$b!r0|3t(rend

g(n)=3*((y(n)-y(n-1))/h(n-1));

4G!W4C2P/z%前边求出lamdamiu和g从⽽可以确定系数矩阵5H6T,D3N+_3R(|(a({4G!W4C2P/z%前边求出lamdamiu和g从⽽可以确定系数矩阵5H6T,D3N+_3R(|(a({4G!W4C2P/z%前边求出lamdamiu和g从⽽可以确定系数矩阵5H6T,D3N+_3R(|(a({

lamda(n)=1;.O/M,V+R8E)Xk

lamda(1)=0;

!?/L0U4P%Z+q8oG$K,B/r%根据第⼆边界条件补充两个lamda和miu的值

#i!p8V!EB/Dfori=1:n#i!p8V!EB/Dfori=1:n

end

)@(]/A%?5X7Sm=followup(lamda,beta,miu,g)r!A8c)^(N6?

%解出m的值从⽽可确定stst为各段的插值多项式%解出m的值从⽽可确定stst