Lagrnge插值程序1.docVIP

在Matlab中，可以编写如下程序来利用Lagrange插值公式进行计算： function f=Lagrange(x,fx,inx) n=length(x);m=length(inx); for i=1:m; z=inx(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x(j))/(x(k)-x(j)); end end s=p*fx(k)+s; end f(i)=s; end plot(x,fx,O,inx,f) 例：已知数据如下，求1—12之间以0.2为间隔的所有值，并画出图形： 1 2 3 4 5 6 7 8 9 10 11 12 12 234 34 -1 34 2 5 23 34 9 45 23 解：执行命令 x=[1:12] fx=[12 234 34 -1 34 2 5 23 34 9 45 23] xi=[1:0.2:12] Lagrange(x,fx,xi) 得出结果： 12.0000 -60.5937 18.2765 124.9778 202.5952 234.0000 223.3757 184.1249 131.4738 78.4253 34.0000 2.9467 -13.6885 -17.5810 -12.0379 -1.0000 11.7556 23.1624 31.1611 34.7730 34.0000 29.6054 22.8332 15.1153 7.8099 2.0000 -1.6307 -2.8397 -1.7907 1.0404 5.0000 9.4024 13.6643 17.4033 20.4834 23.0000 25.2037 27.3769 29.6858 32.0400 34.0000 34.7742 33.3426 28.7320 20.4439 9.0000 -3.4848 -12.8605 -12.8873 4.0592 45.0000 112.3788 197.1817 267.9699 254.3439 23.0000 拉格朗日插值法理论介绍 在 HYPERLINK /wiki/%E6%95%B0%E5%80%BC%E5%88%86%E6%9E%90 \o 数值分析 数值分析中，拉格朗日插值法是以 HYPERLINK /wiki/%E6%B3%95%E5%9B%BD \o 法国 法国十八世纪 HYPERLINK /wiki/%E6%95%B0%E5%AD%A6%E5%AE%B6 \o 数学家 数学家 HYPERLINK /wiki/%E7%BA%A6%E7%91%9F%E5%A4%AB%C2%B7%E8%B7%AF%E6%98%93%E6%96%AF%C2%B7%E6%8B%89%E6%A0%BC%E6%9C%97%E6%97%A5 \o 约瑟夫·路易斯·拉格朗日 约瑟夫·路易斯·拉格朗日命名的一种 HYPERLINK /wiki/%E5%A4%9A%E9%A1%B9%E5%BC%8F%E6%8F%92%E5%80%BC \o 多项式插值 多项式插值方法。许多实际问题中都用 HYPERLINK /wiki/%E5%87%BD%E6%95%B0 \o 函数 函数来表示某种内在联系或规律，而不少函数都只能通过实验和观测来了解。如对实践中的某个 HYPERLINK /wiki/%E7%89%A9%E7%90%86 \o 物理 物理量进行观测，在若干个不同的地方得到相应的观测值，拉格朗日插值法可以找到一个 HYPERLINK /wiki/%E5%A4%9A%E9%A1%B9%E5%BC%8F \o 多项式 多项式，其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日（插值）多项式。 HYPERLINK /wiki/%E6%95%B0%E5%AD%A6 \o 数学 数学上来说，拉格朗日插值法可以给出一个恰好穿过二维 HYPERLINK /wiki/%E5%B9%B3%E9%9D%A2 \o 平面 平面上若干个已知点的多项式函数。拉格朗日插值法最早被 HYPERLINK /wiki/%E8%8B%B1%E5%9B%BD \o 英国 英国数学家 HYPERLINK /wiki/%E7%88%B1%E5%BE%B7%E5%8D%8E%C2%B7%E5%8D%8E%E