一、目的
1、掌握拉格朗日插值基础原理
2、掌握拉格朗日插值组的步骤
3、能用程序语言对拉格朗日插值进行编程实现
二、代码
1、Form1.cs
using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;
namespace Lagrange
{
public partial class Form1 : Form
{
List<double> X = new List<double>();
List<double> Y = new List<double>();
int n;
public double Lix(double x, int i)
{
double s = 1;
for (int j = 0; j <= n; j++)
{
if (j != i)
{
s = s * (x - X[j]) / (X[i] - X[j]);
}
}
return s;
}
public double Lnx(double x)
{
double s = 0;
for (int i = 0; i <= n; i++)
{
s = s + Lix(x, i) * Y[i];
}
return s;
}
public Form1()
{
InitializeComponent();
}
private void Form1_Load(object sender, EventArgs e)
{
}
private void Form1_MouseDown(object sender, MouseEventArgs e)
{
if (e.Button == MouseButtons.Left)
{
X.Add(e.X);
Y.Add(e.Y);
n = X.Count - 1;
Graphics g = this.CreateGraphics();
g.DrawEllipse(new Pen(Color.Red), e.X - 3, e.Y - 3, 6, 6);
}
else if (e.Button == MouseButtons.Right)
{
Graphics g = this.CreateGraphics();
for (double xx = X[0]; xx <= X[n]; xx++)
{
double y = Lnx(xx);
g.DrawEllipse(new Pen(Color.Blue), (int)(xx - 1), (int)(y - 1), 2, 2);
}
}
}
}
}
2、Program.cs
using System;
using System.Collections.Generic;
using System.Linq;
using System.Threading.Tasks;
using System.Windows.Forms;
namespace Lagrange
{
static class Program
{
/// <summary>
/// 应用程序的主入口点。
/// </summary>
[STAThread]
static void Main()
{
Application.EnableVisualStyles();
Application.SetCompatibleTextRenderingDefault(false);
Application.Run(new Form1());
}
}
}
三、程序运行截图
先点击鼠标左键生成若干点然后点击鼠标右键绘制图像。