splinetest.m

来自「计算方法中的三弯矩方法」· M 代码 · 共 11 行

x=[5.79997266200454   6.00970053330447   6.56195092941116   7.13117188826633   7.28885670047093   7.89060457126068 8.46123974367823   8.75186063074590]
y=[-5.69051253227686  -5.86969339062953  -6.16597113194670  -6.28264375084247  -6.28759349831078  -6.19354829641297 -5.91565533140666  -5.67311770545967]
f1dao=tan(atan(0.01087)-pi/4);
f8dao=tan(atan(100)-pi/4);
cs = spline(x,[f1dao y f8dao]);
xx = linspace(5.79997266200454,8.75186063074590,101);
plot(x,y,'o',xx,ppval(cs,xx),'-')
title('Matlab自带函数计算（2）中三次样条插值函数')
xlabel ('X轴')
ylabel ('Y轴')
grid