x = linspace(-1,1,10);
x1 = linspace(-1,1,100);
y =1./(1+25*x.^2);
y1 = interp1(x,y,x1,'cubic');
y2 = interp1(x,y,x1,'spline');
figure();
plot(x,y,'ro');
hold on;
plot(x1,y1,'b',x1,y2,'g');
legend('插值点','多项式','三次样条');
m程序如下:
x=-1:0.0001:1;
y=1./(1+25*x.^2);
x1=-1:0.2:1;
y1=interp1(x,y,x1,'spline');
plot(x1,y1,'o',x,y)
grid on
xlabel('x')
ylabel('y')
y1
结果如图:
m程序如下:
x=-1:0.0001:1;
y=1./(1+25*x.^2);
x1=-1:0.2:1;
y1=interp1(x,y,x1,'spline');
plot(x1,y1,'o',x,y)
grid on
xlabel('x')
ylabel('y')
y1