通過matlab建模,用分段線性和三次樣條兩種插值方法計算x座標軸每改變0.1時得出y座標的變化。
x |
0 |
3 |
5 |
7 |
9 |
11 |
12 |
13 |
14 |
15 |
y |
0 |
1.2 |
1.7 |
2.0 |
2.1 |
2.0 |
1.8 |
1.2 |
1.0 |
1.6 |
通過對上述數據的分析得出步長很小的(x,y)座標
其中主要步驟如下:
(1)分別寫出線性插值與三次樣條插值的代碼
(2)繪畫其圖形
(3)分析並總結
(4)得出結論並檢驗。
運行結果如下:
代碼如下:
附錄
x0=[0 3 5 7 9 11 12 13 14 15];
y0=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6];
x=0:0.1:15;
y1=interp1(x0,y0,x);
y2=interp1(x0,y0,x,'spline');
pp1=csape(x0,y0);
y3=ppval(pp1,x);
pp2=csape(x0,y0,'second');
y4=ppval(pp2,x);
[x',y1',y2',y3',y4']
subplot(1,3,1)
plot(x0,y0,'+',x,y1)
title('Piecewise linear')
subplot(1,3,2)
plot(x0,y0,'+',x,y2)
title('Spline1')
subplot(1,3,3)
plot(x0,y0,'+',x,y3)
title('Spline2')
dx=diff(x);
dy=diff(y3);
dy_dx=dy./dx;
dy_dx0=dy_dx(1)
ytemp=y3(131:151);
ymin=min(ytemp);
index=find(y3==ymin);
xmin=x(index);
[xmin,ymin]