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

2、處的導數值；若缺省，則用均差代替導數；% xi為標量，自變量x；% yi為xi處的函數估計值。n=length(x);ny=length(y);% 輸入的插值點與它的函數值應有相同的個數if n=ny error(The lengths of X and Y must be equal); return;end% 如果沒有給出端點處的導數值，則用均差代替導數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); % 插入節點必須互異 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)% 求插值函數系數：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-