三次樣條插值函數(shù)逼近車門曲線

1、精選優(yōu)質(zhì)文檔-傾情為你奉上精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)三次樣條插值函數(shù)逼近車門曲線。插值節(jié)點數(shù)據(jù)：x012345678910y2.513.304.044.705.225.545.785.405.575.705.80y0.80.2要求1）寫出程序清單，生成M文件2）畫出車門曲線圖，給出插值結(jié)果，生成DOC文件。1、function yi=cubic_spline2(x,y,ydot,xi)% 三次樣條插值公式（三彎矩方程，第一類邊界條件）% x為向量，全部的插值節(jié)點；% y為向量，插值節(jié)點處的函數(shù)值；% ydot為向量，端點

2、處的導(dǎo)數(shù)值；若缺省，則用均差代替導(dǎo)數(shù)；% xi為標(biāo)量，自變量x；% yi為xi處的函數(shù)估計值。n=length(x);ny=length(y);% 輸入的插值點與它的函數(shù)值應(yīng)有相同的個數(shù)if n=ny error(The lengths of X and Y must be equal); return;end% 如果沒有給出端點處的導(dǎo)數(shù)值，則用均差代替導(dǎo)數(shù)if isempty(ydot)=1 ydot=(y(2)-y(1)/(x(2)-x(1) (y(n)-y(n-1)/(x(n)-x(n-1);endh=zeros(1,n);lambda=ones(1,n);mu=ones(1,n);m=

3、zeros(n,1);d=zeros(n,1);for k=2:nh(k)=x(k)-x(k-1); % 插入節(jié)點必須互異 if abs(h(k)eps error(the DATA is error!); return; endendfor k=2:n-1 lambda(k)=h(k+1)/(h(k)+h(k+1); mu(k)=1-lambda(k);d(k)=3*(mu(k)*(y(k+1)-y(k)/h(k+1)+lambda(k)*(y(k)-y(k-1)/h(k);endd(2)=d(2)-lambda(2)*ydot(1);d(n-1)=d(n-1)-mu(n-1)*ydot(2

4、);d(n)=;d(1)=;lambda(n)=;lambda(1)=;mu(n)=;mu(1)=;A=diag(2*ones(1,n-2);for i=1:n-3 A(i,i+1)=mu(i); A(i+1,i)=lambda(i+1);endm=Ad;m=ydot(1);m;ydot(2);for k=2:n if x(k-1)=xi & xi x=0 1 2 3 4 5 6 7 8 9 10;y=2.51 3.30 4.04 4.70 5.22 5.54 5.78 5.40 5.57 5.70 5.80;ydot=0.8 0.2;xx=0.0:0.01:10.0;yy=cubic_spline2(x,y,ydot,xx);plot(x,y,o,xx,yy,r)% 求插值函數(shù)系數(shù)：for k=2:n a(k-1)=y(k-1)/h(k)3; b(k-1)=y(k)/h(k)3; c(k-1)=2*(y(k-1)/h(k)3)+m(k-