實(shí)驗(yàn)四數(shù)據(jù)插值與擬合

實(shí)驗(yàn)四數(shù)據(jù)插值與擬合第1頁，共49頁，2023年，2月20日，星期六4.1MATLAB中的插值函數(shù)函數(shù)插值來源于函數(shù)的以下問題：只知道函數(shù)在某區(qū)間有定義且已得到區(qū)間內(nèi)一些離散點(diǎn)的值，希望用簡單的表達(dá)式近似給出函數(shù)在此區(qū)間上的整體描述，并能與已知離散點(diǎn)上的值相等。插值法按插值函數(shù)的形式主要分為以下幾種形式：（1）代數(shù)多項(xiàng)式插值；（2）三角多項(xiàng)式插值；（3）有理分式插值。第2頁，共49頁，2023年，2月20日，星期六代數(shù)多項(xiàng)式插值是最常用的插值方式，其內(nèi)容也是最豐富的，它又可分為以下幾種插值方式：（1）非等距節(jié)點(diǎn)插值，包括拉格朗日插值、利用均差的牛頓插值和埃特金插值；（2）非等距節(jié)點(diǎn)插值，包括利用差分的牛頓插值和高斯插值等；（3）在插值中增加了導(dǎo)數(shù)的Hermite（埃爾米特）插值；（4）分段插值，包括分段線性插值、分段Hermite（埃爾米特）插值和樣條函數(shù)插值；（5）反插值。按被插值函數(shù)的變量個(gè)數(shù)還可把插值法分為一元插值和多元插值。第3頁，共49頁，2023年，2月20日，星期六4.1.1一元插值函數(shù)MATLAB中的一元插值函數(shù)為interp1()，它的功能是一維數(shù)據(jù)插值（表格查找）。該命令對(duì)數(shù)據(jù)點(diǎn)之間進(jìn)行計(jì)算內(nèi)插值，它出一元函數(shù)f(x)在中間點(diǎn)的數(shù)值，其中函數(shù)f(x)由所給數(shù)據(jù)決定。一元插值函數(shù)interp1()的幾種調(diào)用格式如表4-1所示。表4-1一維插值插值函數(shù)interp1的語法格式語法形式說明y=interp1(x,Y,xi)由已知點(diǎn)集(x,Y）插值計(jì)算xi上的函數(shù)值y=interp1(x,Y,xi)相當(dāng)于x=1:length(Y)的interp(x,Y,xi)y=interp1(x,Y,xi,method)用指定插值方法計(jì)算插值點(diǎn)xi上的函數(shù)值y=interp1(x,Y,xi,method,’extrap’)對(duì)xi中超出已知點(diǎn)集的插值點(diǎn)用指定插值方法計(jì)算函數(shù)值y=interp1(x,Y,xi,method,’extrap’,extrapval)用指定方法插值xi上的函數(shù)值，超出已知點(diǎn)集處函數(shù)值取extrapvaly=interp1(x,Y,xi,method,’pp’)用指定方法插值，但返回結(jié)果為分段多項(xiàng)式第4頁，共49頁，2023年，2月20日，星期六method方法描述‘nearest’最鄰近插值：插值點(diǎn)處函數(shù)值取與插值點(diǎn)最鄰近的已知點(diǎn)的函數(shù)值‘liner’分段線性插值：插值點(diǎn)處函數(shù)值由連接其最鄰近的兩側(cè)點(diǎn)的線性函數(shù)預(yù)測，MATLAB中interp1的默認(rèn)方法‘spline’樣條插值：默認(rèn)為三次樣條插值?？捎胹pline函數(shù)代替‘pchip’三次Hermite多項(xiàng)式插值?？捎胮chip函數(shù)代替‘cubic’同‘pchip’，三次Hermite多項(xiàng)式插值MATLAB中一維插值有多種算法，由interp1函數(shù)中的method指定。MATLAB中一維插值的各種算法如表4-2所示。表4-2一維插值算法(method)第5頁，共49頁，2023年，2月20日，星期六1.Linear（分段線性插值）它的算法是在每個(gè)小區(qū)間[xi,xi+1]上采用簡單的線性插值。在區(qū)間[xi,xi+1]上的子插值多項(xiàng)式為：由此整個(gè)區(qū)間[xi,xi+1]上的插值函數(shù)為：

其中定義如下：第6頁，共49頁，2023年，2月20日，星期六分段線性插值方法在速度和誤差之間取得了比較好的均衡，其插值函數(shù)具有連續(xù)性，但在已知數(shù)據(jù)點(diǎn)處的斜率一般不會(huì)改變，因此不是光滑的。分段線性插值方法是MATLAB一維插值默認(rèn)的方法。第7頁，共49頁，2023年，2月20日，星期六2.Spline（樣條插值）樣條插值是用分段低次多項(xiàng)式去逼近函數(shù)。樣條函數(shù)可以給出光滑的插值曲線，只要在插值區(qū)間端點(diǎn)提供某些導(dǎo)數(shù)信息，樣條插值可以適應(yīng)不同光滑需求。三次樣條是使用最為廣泛的樣條插值，它在每個(gè)子區(qū)間[xi,xi+1]上都是有二階連續(xù)導(dǎo)數(shù)的三次多項(xiàng)式，即

其中都是三次多項(xiàng)式。第8頁，共49頁，2023年，2月20日，星期六對(duì)于給定的離散的測量數(shù)據(jù)經(jīng)x，y（稱為斷點(diǎn)），要尋找一個(gè)三次多項(xiàng)式y(tǒng)=p(x)，以逼近每對(duì)數(shù)據(jù)(xi,yi)點(diǎn)間曲線。過兩點(diǎn)(xi,yi)和(xi+1,yi+1)只能確定一條直線，而通過一點(diǎn)的三次多項(xiàng)式曲線有無窮多條。為使通過中間斷點(diǎn)的三次多項(xiàng)式曲線具有唯一性，要增加以下的連續(xù)條件和邊界條件（因?yàn)槿味囗?xiàng)式有4個(gè)系數(shù)）：（1）三次多項(xiàng)式在點(diǎn)(xi,yi)處有：；（2）三次多項(xiàng)式在點(diǎn)(xi,yi)處有：；（3）三次多項(xiàng)式在點(diǎn)(xi,yi)處有：；（4）邊界條件：。第9頁，共49頁，2023年，2月20日，星期六表4-2中各種方法中：（1）nearest方法速度最快，占用內(nèi)存最小，但一般來說誤差最大，插值結(jié)果最不光滑；（2）spline三次樣條插值是所有插值方法中運(yùn)行耗時(shí)最長的，其插值函數(shù)以及插值函數(shù)的一階、二階導(dǎo)函數(shù)都連續(xù)，因此是最光滑的插值方法，占用內(nèi)存上比cubic方法小，但當(dāng)已知數(shù)據(jù)點(diǎn)不均勻分布時(shí)可能出現(xiàn)異常結(jié)果。（3）cubic三次多項(xiàng)式插值法中插值函數(shù)及其一階導(dǎo)數(shù)都是連續(xù)的，因此其插值結(jié)果也比較光滑，運(yùn)算速度比spline方法略快，但占用內(nèi)存最多。在實(shí)際的使用中，應(yīng)根據(jù)實(shí)際需求和運(yùn)算條件選擇合適的算法。第10頁，共49頁，2023年，2月20日，星期六例4-1用interp1對(duì)sin函數(shù)進(jìn)行分段線性插值。解：在MATLAB命令窗口中輸入以下命令：

>>x=0:2*pi;>>y=sin(x);>>xx=0:0.5:2*pi>>yy=interp1(x,y,xx);>>plot(x,y,'s',xx,yy)注：例4-1中用默認(rèn)的(分段線性插值的linear)對(duì)已知的7個(gè)sin函數(shù)的數(shù)據(jù)點(diǎn)進(jìn)行插值，用plot畫出插值結(jié)果。從圖中可以看出分段線性就是聯(lián)結(jié)兩個(gè)鄰近的已知點(diǎn)的線性函數(shù)插值計(jì)算該區(qū)間內(nèi)插值點(diǎn)上的函數(shù)值。第11頁，共49頁，2023年，2月20日，星期六例4-2用其他一維插值方法對(duì)以下7個(gè)離散數(shù)據(jù)點(diǎn)

(1,3.5)、(2,2.1)、(3,1.3)、(4.0.8)、(5,2.9)、(6,4.2)、(7,5.7)進(jìn)行一維插值方法。解：在MATLAB命令窗口中輸入以下命令：

>>x=[1234567];>>y=[3.52.11.30.82.94.25.7];>>xx=1:0.5:7;>>y1=interp1(x,y,xx,'nearest');>>y2=interp1(x,y,xx,'spline');>>y3=interp1(x,y,xx,'cubic');>>plot(x,y,'o',xx,y1,'-',xx,y2,'-.',xx,y3,':')第12頁，共49頁，2023年，2月20日，星期六第13頁，共49頁，2023年，2月20日，星期六4.2拉格朗日插值法拉格朗日插值法是基于基函數(shù)的插值方法，插值多項(xiàng)式可表示為其中稱為i次基函數(shù)：第14頁，共49頁，2023年，2月20日，星期六在MATLAB中編程實(shí)現(xiàn)拉格朗日插值法函數(shù)為：Language。功能：求已知數(shù)據(jù)點(diǎn)的拉格朗日多項(xiàng)式；調(diào)用格式：f=

Language(x,y)或f=

Language(x,y,x0)。其中，x為已知數(shù)據(jù)點(diǎn)的x坐標(biāo)向量；

y為已知數(shù)據(jù)點(diǎn)的y坐標(biāo)向量；

x0為插值點(diǎn)的x坐標(biāo)；

f為求得的拉格朗日多項(xiàng)式或x0處的插值。第15頁，共49頁，2023年，2月20日，星期六

functionf=Language(x,y,x0)%求已知數(shù)據(jù)點(diǎn)的拉格朗日多項(xiàng)式%已知數(shù)據(jù)點(diǎn)的x坐標(biāo)向量:x%已知數(shù)據(jù)點(diǎn)的y坐標(biāo)向量:y%為插值點(diǎn)的x坐標(biāo):x0%求得的拉格朗日多項(xiàng)式或x0處的插值:fsymst;if(length(x)==length(y))n=length(x);elsedisp('x和y的維數(shù)不相等！');return;end%檢錯(cuò)f=0.0;for(i=1:n)l=y(i);for(j=1:i-1)l=l*(t-x(j))/(x(i)-x(j));end;for(j=i+1:n)l=l*(t-x(j))/(x(i)-x(j));%計(jì)算拉格朗日基函數(shù)

end;f=f+l;%計(jì)算拉格朗日插值函數(shù)

simplify(f);%化簡

if(i==n)if(nargin==3)f=subs(f,'t',x0);%計(jì)算插值點(diǎn)的函數(shù)值

elsef=collect(f);%將插值多項(xiàng)式展開

f=vpa(f,6);%將插值多項(xiàng)式的系數(shù)化成6位精度的小數(shù)

endendend第16頁，共49頁，2023年，2月20日，星期六例4-3根據(jù)下表的數(shù)據(jù)點(diǎn)求出其拉格朗日插值多項(xiàng)式，并計(jì)算當(dāng)x=1.6時(shí)y的值。解：>>

x=[11.21.82.54];>>y=[0.84150.93200.97380.5985-0.7568];>>f=language(x,y)f=1.05427*t-.145485e-1*t^2-.204917*t^3+.328112e-1*t^4-.261189e-1>>f=language(x,y,1.6)f=0.9992x11.21.82.54y0.84150.93200.97380.5985-0.7568第17頁，共49頁，2023年，2月20日，星期六4.3利用均差的牛頓插值法函數(shù)f的零階均差定義為，一階定義均差為一般地，函數(shù)f的k階均差定義為：利用均差的牛頓插值法多項(xiàng)式為

第18頁，共49頁，2023年，2月20日，星期六系數(shù)的計(jì)算過程如表4-3所示表4-3均差計(jì)算表格一階均差二階均差三階均差……n階均差……………第19頁，共49頁，2023年，2月20日，星期六在MATLAB中編程實(shí)現(xiàn)利用均差牛頓插值法函數(shù)為：Newton。功能：求已知數(shù)據(jù)點(diǎn)的均差形式的牛頓插值多項(xiàng)式；調(diào)用格式：f=

Newton(x,y)或f=

Newton(x,y,x0)。其中，x為已知數(shù)據(jù)點(diǎn)的x坐標(biāo)向量；

y為已知數(shù)據(jù)點(diǎn)的y坐標(biāo)向量；

x0為插值點(diǎn)的x坐標(biāo)；

f為求得的牛頓插值法多項(xiàng)式或x0處的插值。第20頁，共49頁，2023年，2月20日，星期六functionf=Newton(x,y,x0)%求已知數(shù)據(jù)點(diǎn)的均差形式牛頓插值多項(xiàng)式%已知數(shù)據(jù)點(diǎn)的x坐標(biāo)向量:x%已知數(shù)據(jù)點(diǎn)的y坐標(biāo)向量:y%為插值點(diǎn)的x坐標(biāo):x0%求得的均差形式牛頓插值多項(xiàng)式或x0處的插值:fsymst;if(length(x)==length(y))n=length(x);c(1:n)=0.0;elsedisp('x和y的維數(shù)不相等！');return;end第21頁，共49頁，2023年，2月20日，星期六f=y(1);y1=0;l=1;

for(i=1:n-1)for(j=i+1:n)y1(j)=(y(j)-y(i))/(x(j)-x(i));endc(i)=y1(i+1);l=l*(t-x(i));f=f+c(i)*l;simplify(f);y=y1;

if(i==n-1)if(nargin==3)f=subs(f,'t',x0);elsef=collect(f);%將插值多項(xiàng)式展開

f=vpa(f,6);endendend第22頁，共49頁，2023年，2月20日，星期六例4-4根據(jù)下表的數(shù)據(jù)點(diǎn)求出其均差形式牛頓插值多項(xiàng)式，并計(jì)算當(dāng)x=2.0時(shí)y的值。解：>>x=[11.21.82.54];>>y=[11.443.246.2516];>>f=Newton(x,y)f=.182711e-14-.482154e-14*t+1.00000*t^2-.169177e-14*t^3+.211471e-15*t^4>>f=Newton(x,y,2.0)f=4x11.21.82.54y11.443.246.2516第23頁，共49頁，2023年，2月20日，星期六4.4利用差分的牛頓插值法差分分為向前差分、后向差分和中心差分三種，它們的記法及定義如下為：

其中：代表向前差分；代表向后差分；代表向后差分。第24頁，共49頁，2023年，2月20日，星期六

假設(shè)。為了方便，可構(gòu)造如表4-4所示的差分表（）。表4-4差分計(jì)算表格………………第25頁，共49頁，2023年，2月20日，星期六4.4.1向前牛頓插值向前牛頓插值多項(xiàng)式可表示如下：其中叫做步長，，且的取值范圍為。第26頁，共49頁，2023年，2月20日，星期六在MATLAB中編程實(shí)現(xiàn)向前牛頓插值法函數(shù)為：Newtonforward。功能：求已知數(shù)據(jù)點(diǎn)的向前牛頓插值法多項(xiàng)式；調(diào)用格式：f=

Newtonforward(x,y)或

f=

Newtonforward

(x,y,x0)。其中，x為已知數(shù)據(jù)點(diǎn)的x坐標(biāo)向量；

y為已知數(shù)據(jù)點(diǎn)的y坐標(biāo)向量；

x0為插值點(diǎn)的x坐標(biāo)；

f為求得的向前牛頓插值法多項(xiàng)式或x0處的插值。第27頁，共49頁，2023年，2月20日，星期六functionf=Newtonforward(x,y,x0)%求已知數(shù)據(jù)點(diǎn)的向前差分牛頓插值多項(xiàng)式%已知數(shù)據(jù)點(diǎn)的x坐標(biāo)向量:x%已知數(shù)據(jù)點(diǎn)的y坐標(biāo)向量:y%為插值點(diǎn)的x坐標(biāo):x0%求得的向前差分牛頓插值多項(xiàng)式或x0處的插值:fsymst;if(length(x)==length(y))n=length(x);c(1:n)=0.0;elsedisp('x和y的維數(shù)不相等！');return;end第28頁，共49頁，2023年，2月20日，星期六f=y(1);y1=0;xx=linspace(x(1),x(n),(x(2)-x(1)));if(xx~=x)disp('節(jié)點(diǎn)之間不是等距的！');return;endfor(i=1:n-1)for(j=1:n-i)y1(j)=y(j+1)-y(j);endc(i)=y1(1);l=t;for(k=1:i-1)l=l*(t-k);end;f=f+c(i)*l/factorial(i);simplify(f);y=y1;if(i==n-1)if(nargin==3)f=subs(f,'t',(x0-x(1))/(x(2)-x(1)));elsef=collect(f);f=vpa(f,6);endendend第29頁，共49頁，2023年，2月20日，星期六4.4.2向前牛頓插值向后牛頓插值多項(xiàng)式可表示如下：其中叫做步長，，且的取值范圍為。第30頁，共49頁，2023年，2月20日，星期六在MATLAB中編程實(shí)現(xiàn)向后牛頓插值法函數(shù)為：Newtonback。功能：求已知數(shù)據(jù)點(diǎn)的向前牛頓插值法多項(xiàng)式；調(diào)用格式：f=

Newtonback

(x,y)或

f=

Newtonback(x,y,x0)。其中，x為已知數(shù)據(jù)點(diǎn)的x坐標(biāo)向量；

y為已知數(shù)據(jù)點(diǎn)的y坐標(biāo)向量；

x0為插值點(diǎn)的x坐標(biāo)；

f為求得的向前牛頓插值法多項(xiàng)式或x0處的插值。第31頁，共49頁，2023年，2月20日，星期六functionf=Newtonback(x,y,x0)%求已知數(shù)據(jù)點(diǎn)的向后差分牛頓插值多項(xiàng)式%已知數(shù)據(jù)點(diǎn)的x坐標(biāo)向量:x%已知數(shù)據(jù)點(diǎn)的y坐標(biāo)向量:y%為插值點(diǎn)的x坐標(biāo):x0%求得的向前差分牛頓插值多項(xiàng)式或x0處的插值:fsymst;if(length(x)==length(y))n=length(x);c(1:n)=0.0;elsedisp('x和y的維數(shù)不相等！');return;end第32頁，共49頁，2023年，2月20日，星期六f=y(n);y1=0;xx=linspace(x(1),x(n),(x(2)-x(1)));if(xx~=x)disp('節(jié)點(diǎn)之間不是等距的！');return;endfor(i=1:n-1)for(j=i+1:n)y1(j)=y(j)-y(j-1);endc(i)=y1(n);l=t;for(k=1:i-1)l=l*(t+k);end;

f=f+c(i)*l/factorial(i);simplify(f);y=y1;if(i==n-1)if(nargin==3)f=subs(f,'t',(x0-x(n))/(x(2)-x(1)));elsef=collect(f);f=vpa(f,6);endendend第33頁，共49頁，2023年，2月20日，星期六例4-5根據(jù)下表的數(shù)據(jù)點(diǎn)求出其差分形式的牛頓插值多項(xiàng)式，并計(jì)算當(dāng)x=1.55時(shí)y的值。解：>>x=[11.21.41.61.8];>>y=[0.84150.9320 0.98540.99960.9738];>>f=Newtonforward(x,y)f=.841500+.108025*t-.169042e-1*t^2-.675000e-3*t^3+.541667e-4*t^4>>f=Newtonforward(x,y,1.55)f=0.9998f=Newtonback(x,y)f=.973800-.457417e-1*t-.198042e-1*t^2+.191667e-3*t^3+.541667e-4*t^4>>f=Newtonback(x,y,1.55)f=0.9998x11.21.41.61.8y0.84150.93200.98540.99960.9738第34頁，共49頁，2023年，2月20日，星期六4.5Hermite插值Hermite插值滿足在節(jié)點(diǎn)上等于給定函數(shù)值，而且在節(jié)點(diǎn)上的導(dǎo)數(shù)值也等于給定的導(dǎo)數(shù)值。對(duì)于高階導(dǎo)數(shù)的情況，Hermite插值多項(xiàng)式比較復(fù)雜，在實(shí)際中，常常遇到的是函數(shù)值與一階導(dǎo)數(shù)給定的情況。在此情況下，n個(gè)節(jié)點(diǎn)x1，x2，…，xn的Hermite插值多項(xiàng)式的表達(dá)式如下：其中，，，第35頁，共49頁，2023年，2月20日，星期六在MATLAB中編程實(shí)現(xiàn)Hermite插值法函數(shù)為：Hermite。功能：求已知數(shù)據(jù)點(diǎn)的Hermite插值法多項(xiàng)式；調(diào)用格式：f=

Hermite

(x,y,y_1)或

f=

Hermite

(x,y,y_1,x0)。其中，x為已知數(shù)據(jù)點(diǎn)的x坐標(biāo)向量；

y為已知數(shù)據(jù)點(diǎn)的y坐標(biāo)向量；

y_1為已知數(shù)據(jù)點(diǎn)導(dǎo)數(shù)向量；

x0為插值點(diǎn)的x坐標(biāo)；

f為求得的Hermite插值法多項(xiàng)式或x0處的插值。第36頁，共49頁，2023年，2月20日，星期六functionf=Hermite(x,y,y_1,x0)%求已知數(shù)據(jù)點(diǎn)的向后差分牛頓插值多項(xiàng)式%已知數(shù)據(jù)點(diǎn)的x坐標(biāo)向量:x%已知數(shù)據(jù)點(diǎn)的y坐標(biāo)向量:y%已知數(shù)據(jù)點(diǎn)的導(dǎo)數(shù)向量:y_1%求得的Hermite插值多項(xiàng)式或x0處的插值:fsymst;f=0.0;if(length(x)==length(y))if(length(y)==length(y_1))n=length(x);elsedisp('y和y的導(dǎo)數(shù)的維數(shù)不相等！');return;endelsedisp('x和y的維數(shù)不相等！');return;end第37頁，共49頁，2023年，2月20日，星期六fori=1:nh=1.0;a=0.0;forj=1:nif(j~=i)h=h*(t-x(j))^2/((x(i)-x(j))^2);a=a+1/(x(i)-x(j));endend

f=f+h*((x(i)-t)*(2*a*y(i)-y_1(i))+y(i));

if(i==n)if(nargin==4)f=subs(f,'t',x0);elsef=vpa(f,6);endendend第38頁，共49頁，2023年，2月20日，星期六例4-6根據(jù)下表的數(shù)據(jù)點(diǎn)求出其Hermite插值多項(xiàng)式，并計(jì)算當(dāng)x=1.144時(shí)y的值。解：

>>x=1:0.2:1.8;>>y=[11.09541.18321.26491.3416];>>y_1=[0.50.45640.42260.39530.3727];>>f=Hermite(x,y,y_1)f=678.168*(t-1.20000)^2*(t-1.40000)^2*(t-1.60000)^2*(t-1.80000)^2*(-20.3333+21.3333*t)+10850.7*(t-1.)^2*(t-1.40000)^2*(t-1.60000)^2*(t-1.80000)^2*(-10.4063+9.58473*t)+24414.1*(t-1.)^2*(t-1.20000)^2*(t-1.60000)^2*(t-1.80000)^2*(.591560+.422600*t)+10850.7*(t-1.)^2*(t-1.20000)^2*(t-1.40000)^2*(t-1.80000)^2*(17.4978-10.1455*t)+678.168*(t-1.)^2*(t-1.20000)^2*(t-1.40000)^2*(t-1.60000)^2*(50.9807-27.5773*t)>>f=Hermite(x,y,y_1,1.44)f=1.2000x11.21.41.61.8y11.09541.18321.26491.3416y’0.50000.45640.42260.39530.3727第39頁，共49頁，2023年，2月20日，星期六4.6spline三次樣條插值Spline插值為分段三次插值，即：在每一個(gè)小區(qū)間上是不超過三次的多項(xiàng)式且具有二階連續(xù)導(dǎo)數(shù)，利用具有一階導(dǎo)數(shù)邊界條件的源程序?yàn)椋簄aspline.mfunctionm=naspline(x,y,dy0,dyn,xx)%用途：三階樣條插值（一階導(dǎo)數(shù)邊界條件）%格式：x為節(jié)點(diǎn)向量；y為數(shù)據(jù)；dyo，dyn為左右兩端點(diǎn)的一階導(dǎo)數(shù)%如果xx缺省，則輸出各節(jié)點(diǎn)的一階導(dǎo)數(shù)值，否則m為xx的三階樣條插值n=length(x)-1;%計(jì)算小區(qū)間的個(gè)數(shù)h=diff(x);lemda=h(2:n)./(h(1:n-1)+h(2:n));mu=1-lemda;第40頁，共49頁，2023年，2月20日，星期六g=3*(lemda.*diff(y(1:n))./h(1:n-1)+mu.*diff(y(2:n+1))./h(2:n));g(1)=g(1)-lemda(1)*dy0;g(n-1)=g(n-1)-mu(n-1)*dyn;%求解三對(duì)角方程dy=nachase(lemda,2*ones(1:n-1),mu,g);%若給插值節(jié)點(diǎn)，計(jì)算插值m=[dy0;dy;dyn];ifnargin>=5s=zeros(size(xx));fori=1:nifi==1kk=find(xx<=x(2));elseifi==nkk=find(xx>=x(n));elsekk=find(xx>x(i)&xx<=x(i+1));endxbar=(xx(kk)-x(i))/h(i);s(kk)=alpha0(xbar)*y(i)+alpha1(xbar)*y(i+1)+...h(i)*beta0(xbar)*m(i)+h(i)*beta1(xbar)*m(i+1);endm=s;end第41頁，共49頁，2023年，2月20日，星期六%追趕法functionx=nachase(a,b,c,d)n=length(a);fork=2:nb(k)=b(k)-a(k)/b(k-1)*c(k-1);d(k)=d(k)-a(k)/b(k-1)*d(k-1);endx(n)=d(n)/b(n);fork=n-1:-1:1x(k)=(d(k)-c(k)*x(k+1)/b(k));endx=x(:);%基函數(shù)functiony=alpha0(x)y=2*x.^3-3*x.^2+1;functiony=alpha1(x)y=-2*x.^3+3*x.^2;functiony=beta0(x)y=x.^3-2*x.^2+x;functiony=beta1(x)y=x.^3-x.^2;第42頁，共49頁，2023年，2月20日，星期六4.7多項(xiàng)式曲線擬合對(duì)給定的試驗(yàn)數(shù)據(jù)點(diǎn)(xi,yi)(i=1,2,…,N)，可以構(gòu)造m次多項(xiàng)式：由曲線擬合的定義，應(yīng)該使得下式取極小值通過簡單運(yùn)算可以得出系數(shù)是下面線性方程組的解第43頁，共49頁，2023年，2月20日，星期六

其中第44頁，共49頁，2023年，2月20日，星期六在MATLAB中編程實(shí)現(xiàn)