matlab探究龙格现象

在数值分析领域中，龙格现象是在一组等间插值点上使用具有高次多项式的多项式插值时出现的区间边缘处的振荡问题。 它是由卡尔·龙格（Runge）在探索使用多项式插值逼近某些函数时的错误行为时发现的。这一发现非常重要，因为它表明使用高次多项式插值并不总能提高准确性。 该现象与傅里叶级数近似中的吉布斯现象相似。
–维基百科

龙格函数使用各种插值方法的结果

龙格函数：
$$f(x)=\frac{1}{1+25x^2}$$
运行结果：
多项式插值现象：

图像上方的数字表示插值数据点的个数。 该图像表明，龙格函数当使用高次的插值多项式时造成的误差会随着n的增大而增大。 当n增大时，部分区间上多项式截断误差偏大的现象更严重，这种现象称为龙格现象。

与之对比的是分段插值，这里提供分段线性插值和三次样条插值：

分段线型插值现象：

三次样条插值现象：

matlab 代码：

多项式插值代码：

len=[3:1:20];% len存储插值点的数目
pic_num = 1
for i=1:length(len)
    x=linspace(-1,1,len(i));
    y = 1.0./(1+25.*x.*x); % 龙格函数
    a=polyfit(x,y,length(x)-1); % 获取多项式插值的系数矩阵
    
    
    xi=[-1:0.001:1];
    yi=polyval(a,xi);% 插值多项式的图像
    cla reset % 重新绘图
    plot(x, y,'*r');% 绘制插值数据点
    hold on
    plot(xi, yi,'r');%绘制插值多项式函数
    hold on
    plot(xi, 1.0./(1+25.*xi.*xi));%绘制龙格函数
    ylim([-0.5, 1.2]);
    title(len(i));%在图像顶部显示插值数据点数目
    
    % 生成GIF动图并保存
    F=getframe(gcf);
    I=frame2im(F);
    [I,map]=rgb2ind(I,256);
    if pic_num == 1
        imwrite(I,map,'test.gif','gif', 'Loopcount',inf,'DelayTime',0.3);
    else
        imwrite(I,map,'test.gif','gif','WriteMode','append','DelayTime',0.3);
    end
    pic_num = pic_num + 1;
end

分段插值代码：

len=[3:1:20];% len存储插值点的数目
pic_num = 1
for i=1:length(len)
    x=linspace(-1,1,len(i));
    y = 1.0./(1+25.*x.*x); % 龙格函数
    
    xi=[-1:0.001:1];
    yi=interp1(x,y,xi, 'linear'); % 第四个参数linear表示线性插值，改为spline则成为三次样条插值
    cla reset% 重新绘图
    plot(x, y,'*r');% 绘制插值数据点
    hold on
    plot(xi, yi,'r');%绘制插值多项式函数
    hold on
    plot(xi, 1.0./(1+25.*xi.*xi));%绘制龙格函数
    ylim([-0.5, 1.2]);
    % xlim([-1.2, 1.2]);
    % saveas(gcf,['./','block_',num2str(i),'.jpg']);
    title(len(i));%在图像顶部显示插值数据点数目
    
    % 生成GIF动图并保存
    F=getframe(gcf);
    I=frame2im(F);
    [I,map]=rgb2ind(I,256);
    if pic_num == 1
        imwrite(I,map,'test.gif','gif', 'Loopcount',inf,'DelayTime',0.3);
    else
        imwrite(I,map,'test.gif','gif','WriteMode','append','DelayTime',0.3);
    end
    pic_num = pic_num + 1; 
end

---

参考资料：

1. Matlab如何制作和保存gif动图 https://blog.csdn.net/lusongno1/article/details/78632457
2. matlab多项式、插值、数据拟合 https://wenku.baidu.com/view/6c7a6619227916888486d7a6.html