MATLAB三次樣條插值之三轉(zhuǎn)角法

1、精選優(yōu)質(zhì)文檔-傾情為你奉上非常類似前面的三彎矩法，這里的sanzhj函數(shù)和intersanzhj作用相當(dāng)于前面的sanwanj和intersanwj，追趕法程序通用，代碼如下。%function newu,w,newv,d=sanzhj(x,y,x0,y0,y1a,y1b)% 三轉(zhuǎn)角樣條插值% 將插值點(diǎn)分兩次輸入，x0 y0 單獨(dú)輸入% 邊值條件a的一階導(dǎo)數(shù) y1a 和b的一階導(dǎo)數(shù) y1bn=length(x);m=length(y);if m=nerror('x or y 輸入有誤，再來(lái)');endv=ones(n-1,1);u=ones(n-1,1);d=zeros(n-1

2、,1);w=2*ones(n-1,1);h0=x(1)-x0;h=zeros(n-1,1);for k=1:n-1h(k)=x(k+1)-x(k);endv(1)=h0/(h0+h(1);u(1)=1-v(1);d(1)=3*(v(1)*(y(2)-y(1)/h(1)+u(1)*(y(1)-y0)/h0);%for k=2:n-1v(k)=h(k-1)/(h(k-1)+h(k);u(k)=1-v(k);d(k)=3*(v(k)*(y(k+1)-y(k)/h(k)+u(k)*(y(k)-y(k-1)/h(k-1);endd(1)=d(1)-u(1)*y1a;d(n-1)=d(n-1)-v(n-1

3、)*y1b;newv=v(1:n-2,:);newu=u(2:n-1,:);%function intersanzhj(x,y,x0,y0,y1a,y1b)% 三轉(zhuǎn)角樣條插值%第一部分n=length(x);m=length(y);if m=nerror('x or y 輸入有誤，再來(lái)');end%重新定義hh=zeros(n,1);h(1)=x(1)-x0;for k=2:nh(k)=x(k)-x(k-1);end% 調(diào)用三轉(zhuǎn)角函數(shù)a,b,c,d=sanzhj(x,y,x0,y0,y1a,y1b);% 三對(duì)角方程m=chase(a,b,c,d);% 求MM=1;m;0;% 求

4、插值函數(shù)fprintf('三次樣條（三轉(zhuǎn)角）插值的函數(shù)表達(dá)式n');syms X ;fprintf('S0-1:n');S(1)=collect(y0/h(1).3)*(X-x(1).2*(h(1)+2*(X-x0)+(y(1)/h(1).3)*(X-x0).2*(h(1)+2*(x(1)-X)+(M(1)/h(1).2)*(X-x0)*(X-x(1).2+(M(2)/h(1).2)*(X-x(1)*(X-x0).2);for k=2:nfprintf('S%d-%d:n',k-1,k);S(k)=collect(y(k-1)/h(k).3)*(

5、X-x(k).2*(h(k)+2*(X-x(k-1)+(y(k)/h(k).3)*(X-x(k-1).2*(h(k)+2*(x(k)-X)+(M(k)/h(k).2)*(X-x(k-1)*(X-x(k).2+(M(k+1)/h(k).2)*(X-x(k)*(X-x(k-1).2);endS=S.'disp(S);fprintf('以上為樣條函數(shù)（三轉(zhuǎn)角）解析式，顯示為手寫(xiě)習(xí)慣如下：n');pretty(S);%第二部分%是否繼續(xù)運(yùn)行程序myloop=input('繼續(xù)運(yùn)行程序輸入“1”，否則輸入“0”n');if myloopwhile myloopxi

6、=input('輸入需要計(jì)算的點(diǎn)的值，并按回車(chē)鍵n');if xi>x0|xi<x(n)fprintf('現(xiàn)在開(kāi)始計(jì)算輸入點(diǎn)的插值函數(shù)值n');elsefprintf('輸入數(shù)值不在插值范圍內(nèi)，請(qǐng)重新輸入n');xi=input('輸入需要計(jì)算的點(diǎn)的值，并按回車(chē)鍵n');end% 確定輸入的數(shù)值應(yīng)該使用哪個(gè)解析式newx=x0;x;r,suoy=min(abs(newx-xi);fprintf('輸入點(diǎn)的插值函數(shù)值為:nn');fprintf('t');if xi<=newx(suoy)f=subs(S(suoy-1),X,xi);elsef=subs