代码怎么写，已知x = [0 2 4 5 8 12 12.8 17.2 19.9 20]，对y=exp(x)sin(x)函数曲线上的数据点进行四种不同插值，并绘制原数据散点图和插值后的曲线图。
% 给定数据点
x = [0 2 4 5 8 12 12.8 17.2 19.9 20];
y = exp(x).*sin(x);
% 用四种不同的插值方法拟合数据
xi = linspace(0,20,200); % 定义插值的新点
yi1 = interp1(x,y,xi,'linear'); % 线性插值
yi2 = interp1(x,y,xi,'nearest'); % 最邻近插值
yi3 = interp1(x,y,xi,'spline'); % 样条插值
yi4 = interp1(x,y,xi,'pchip'); % PCHIP插值
% 绘图
figure;
plot(x,y,'o',xi,yi1,'-',xi,yi2,'--',xi,yi3,':',xi,yi4,'-.');
legend('原始数据','线性插值','最邻近插值','样条插值','PCHIP插值');
xlabel('x');
ylabel('y');
title('四种不同插值方法的比较');