C#函数之实现Lagrange插值算法

最新推荐文章于 2024-07-25 23:26:35 发布
最新推荐文章于 2024-07-25 23:26:35 发布
阅读量1.9k 收藏 13
点赞数 3
分类专栏： C#函数 文章标签： 编程语言 线性代数 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/absll/article/details/116828796
版权

博主最近想自己独立地完成一个精密单点定位的应用程序,先实现最简单的功能,那么基础工作就是要对igs最终的精密星历产品进行插值,以满足结算的要求,详细阐述请看前述文章:卫星位置插值方法简介(一)

博主在看过原理后，便使用C#进行了Lagrange函数的编写，话不多说先上代码：

 class interp
    {
        public double lagrange(List<double> T, List<double> f,double t)
        {
            double L = 0;//最后返回的多项式值
            double[] l = new double[T.Count];//定义的每一个多项式系数列表
            for (int i = 0; i < l.Length; i++)
            {
                double l1 = 1;//多项式系数分子
                double l2 = 1;//多项式系数分母
                for (int j = 0; j < T.Count; j++)
                {    
                    if (i != j)
                    {
                        l1 = l1 * (t - T[j]);
                        l2 = l2 * (T[i] - T[j]);

                    }
                    else
                    {
                        continue;
                    
                    }

                }
                l[i] = l1 / l2;
            }
            for (int i = 0; i < l.Length; i++)
            {
                L = L + l[i] * f[i];
            
            }
            return L;
        }
    
    
    }

首先我是先建立了一个插值函数的类，这里面封装之后会写到的所有的插值函数，包括Lagrange插值和Newton插值法，然后具体实现了Lagrange函数
那么写完之后该怎样调用呢？请先看如下代码：

 public int c;
    private void comboBox1_SelectedIndexChanged(object sender, EventArgs e)
    {
            if (comboBox1.Text == "三阶")
            {
                c = 1;
            
            }
            if (comboBox1.Text == "五阶")
            {
                c = 2;

            }
            if (comboBox1.Text == "六阶")
            {
                c = 3;

            }
            if (comboBox1.Text == "八阶")
            {
                c = 4;

            }
            if (comboBox1.Text == "十二阶")
            {
                c = 5;

            }
    }
        private void button1_Click(object sender, EventArgs e)
        {
            
            interp interp = new interp();
            List<double> T = new List<double>();
            List<double> ftx = new List<double>();
            List<double> fty = new List<double>();
            List<double> ftz = new List<double>();
            if (c == 1)
            {
                T.Add(0);
                T.Add(0.25);
                T.Add(0.5);
                T.Add(1);
                ftx.Add(13686.913694);
                ftx.Add(14069.247504);
                ftx.Add(14261.728751);
                ftx.Add(14209.784087);
                fty.Add(-22099.331881);
                fty.Add(-22252.702708);
                fty.Add(-22148.736352);
                fty.Add(-21086.622807);
                ftz.Add(-4728.984963);
                ftz.Add(-1842.522185);
                ftz.Add(1076.236929);
                ftz.Add(6805.670796);
                double satx = interp.lagrange(T, ftx, 0.75);
                double saty = interp.lagrange(T, fty, 0.75);
                double satz = interp.lagrange(T, ftz, 0.75);
                textBox1.Text = satx.ToString();
                textBox2.Text = saty.ToString();
                textBox3.Text = satz.ToString();
                textBox4.Text = (Convert.ToDouble(textBox1.Text) - 14296.367615).ToString();
                textBox5.Text = (Convert.ToDouble(textBox2.Text) + 21764.787764).ToString();
                textBox6.Text = (Convert.ToDouble(textBox3.Text) - 3976.035651).ToString();
               

            }
            if (c == 2)
            {
                T.Add(0);
                T.Add(0.25);
                T.Add(0.5);
                T.Add(1);
                T.Add(1.25);
                T.Add(1.5);
                ftx.Add(13686.913694);
                ftx.Add(14069.247504);
                ftx.Add(14261.728751);
                ftx.Add(14209.784087);
                ftx.Add(14041.598305);
                ftx.Add(13832.692938);
                fty.Add(-22099.331881);
                fty.Add(-22252.702708);
                fty.Add(-22148.736352);
                fty.Add(-21086.622807);
                fty.Add(-20109.159114);
                fty.Add(-18836.848057);
                ftz.Add(-4728.984963);
                ftz.Add(-1842.522185);
                ftz.Add(1076.236929);
                ftz.Add(6805.670796);
                ftz.Add(9514.965215);
                ftz.Add(12055.736360);
                double satx = interp.lagrange(T, ftx, 0.75);
                double saty = interp.lagrange(T, fty, 0.75);
                double satz = interp.lagrange(T, ftz, 0.75);
                textBox1.Text = satx.ToString();
                textBox2.Text = saty.ToString();
                textBox3.Text = satz.ToString();
                textBox4.Text = (Convert.ToDouble(textBox1.Text) - 14296.367615).ToString();
                textBox5.Text = (Convert.ToDouble(textBox2.Text) + 21764.787764).ToString();
                textBox6.Text = (Convert.ToDouble(textBox3.Text) - 3976.035651).ToString();

            }
            if (c == 3)
            {
                T.Add(0);
                T.Add(0.25);
                T.Add(0.5);
                T.Add(1);
                T.Add(1.25);
                T.Add(1.5);
                T.Add(1.75);
                ftx.Add(13686.913694);
                ftx.Add(14069.247504);
                ftx.Add(14261.728751);
                ftx.Add(14209.784087);
                ftx.Add(14041.598305);
                ftx.Add(13832.692938);
                ftx.Add(13623.421198);
                fty.Add(-22099.331881);
                fty.Add(-22252.702708);
                fty.Add(-22148.736352);
                fty.Add(-21086.622807);
                fty.Add(-20109.159114);
                fty.Add(-18836.848057);
                fty.Add(-17283.676462);
                ftz.Add(-4728.984963);
                ftz.Add(-1842.522185);
                ftz.Add(1076.236929);
                ftz.Add(6805.670796);
                ftz.Add(9514.965215);
                ftz.Add(12055.736360);
                ftz.Add(14382.736529);
                double satx = interp.lagrange(T, ftx, 0.75);
                double saty = interp.lagrange(T, fty, 0.75);
                double satz = interp.lagrange(T, ftz, 0.75);
                textBox1.Text = satx.ToString();
                textBox2.Text = saty.ToString();
                textBox3.Text = satz.ToString();
                textBox4.Text = (Convert.ToDouble(textBox1.Text) - 14296.367615).ToString();
                textBox5.Text = (Convert.ToDouble(textBox2.Text) + 21764.787764).ToString();
                textBox6.Text = (Convert.ToDouble(textBox3.Text) - 3976.035651).ToString();

            }
            if (c == 4)
            {
                T.Add(0);
                T.Add(0.25);
                T.Add(0.5);
                T.Add(1);
                T.Add(1.25);
                T.Add(1.5);
                T.Add(1.75);
                T.Add(2);
                T.Add(2.25);
                ftx.Add(13686.913694);
                ftx.Add(14069.247504);
                ftx.Add(14261.728751);
                ftx.Add(14209.784087);
                ftx.Add(14041.598305);
                ftx.Add(13832.692938);
                ftx.Add(13623.421198);
                ftx.Add(13451.837748);
                ftx.Add(13352.029724);
                fty.Add(-22099.331881);
                fty.Add(-22252.702708);
                fty.Add(-22148.736352);
                fty.Add(-21086.622807);
                fty.Add(-20109.159114);
                fty.Add(-18836.848057);
                fty.Add(-17283.676462);
                fty.Add(-15472.783785);
                fty.Add(-13435.708124);
                ftz.Add(-4728.984963);
                ftz.Add(-1842.522185);
                ftz.Add(1076.236929);
                ftz.Add(6805.670796);
                ftz.Add(9514.965215);
                ftz.Add(12055.736360);
                ftz.Add(14382.736529);
                ftz.Add(16454.540357);
                ftz.Add(18234.356133);
                double satx = interp.lagrange(T, ftx, 0.75);
                double saty = interp.lagrange(T, fty, 0.75);
                double satz = interp.lagrange(T, ftz, 0.75);
                textBox1.Text = satx.ToString();
                textBox2.Text = saty.ToString();
                textBox3.Text = satz.ToString();
                textBox4.Text = (Convert.ToDouble(textBox1.Text) - 14296.367615).ToString();
                textBox5.Text = (Convert.ToDouble(textBox2.Text) + 21764.787764).ToString();
                textBox6.Text = (Convert.ToDouble(textBox3.Text) - 3976.035651).ToString();

            }
            if (c == 5)
            {
                T.Add(0);
                T.Add(0.25);
                T.Add(0.5);
                T.Add(1);
                T.Add(1.25);
                T.Add(1.5);
                T.Add(1.75);
                T.Add(2);
                T.Add(2.25);
                T.Add(2.5);
                T.Add(2.75);
                T.Add(3);
                T.Add(3.25);
                ftx.Add(13686.913694);
                ftx.Add(14069.247504);
                ftx.Add(14261.728751);
                ftx.Add(14209.784087);
                ftx.Add(14041.598305);
                ftx.Add(13832.692938);
                ftx.Add(13623.421198);
                ftx.Add(13451.837748);
                ftx.Add(13352.029724);
                ftx.Add(13352.620967);
                ftx.Add(13475.514776);
                ftx.Add(13734.929294);
                ftx.Add(14136.765807);
                fty.Add(-22099.331881);
                fty.Add(-22252.702708);
                fty.Add(-22148.736352);
                fty.Add(-21086.622807);
                fty.Add(-20109.159114);
                fty.Add(-18836.848057);
                fty.Add(-17283.676462);
                fty.Add(-15472.783785);
                fty.Add(-13435.708124);
                fty.Add(-11211.291704);
                fty.Add(-8844.292298);
                fty.Add(-6383.760457);
                fty.Add(-3881.252747);
                ftz.Add(-4728.984963);
                ftz.Add(-1842.522185);
                ftz.Add(1076.236929);
                ftz.Add(6805.670796);
                ftz.Add(9514.965215);
                ftz.Add(12055.736360);
                ftz.Add(14382.736529);
                ftz.Add(16454.540357);
                ftz.Add(18234.356133);
                ftz.Add(19690.739661);
                ftz.Add(20798.192448);
                ftz.Add(21537.629852);
                ftz.Add(21896.709250);
                double satx = interp.lagrange(T, ftx, 0.75);
                double saty = interp.lagrange(T, fty, 0.75);
                double satz = interp.lagrange(T, ftz, 0.75);
                textBox1.Text = satx.ToString();
                textBox2.Text = saty.ToString();
                textBox3.Text = satz.ToString();
                textBox4.Text = (Convert.ToDouble(textBox1.Text) - 14296.367615).ToString();
                textBox5.Text = (Convert.ToDouble(textBox2.Text) + 21764.787764).ToString();
                textBox6.Text = (Convert.ToDouble(textBox3.Text) - 3976.035651).ToString();

            }



        }

首先说明一下，本次编写只为完成基本功能，因此已知数值等信息是直接传入到列表中。
下面展示一下程序具体界面：
在这里插入图片描述三阶（x,y米级，z厘米级）

在这里插入图片描述五阶（x毫米级,y 厘米级，z亚毫米级）

在这里插入图片描述六阶（x毫米级,y 亚毫米级，z亚毫米级）
在这里插入图片描述八阶（x米级,y 米级，z米级）
在这里插入图片描述
十二阶（x米级,y 米级，z米级）
我们可以看到，插值精度并不会随着插值阶数的不断增加而不断提高，到达一定阶数后，精度反而会降低，产生龙格现象。

以上就是C#实现Lagrange插值算法的全部内容啦！

博主后面会将整个精密星历文件都进行Lagrange插值，并输出为.sp3文件，敬请期待哦！
在这里插入图片描述

感谢关注！感谢支持！一起学习，一起进步！文中若有疏漏之处，还请指正！

十八与她
关注
  • 3
    点赞
  • 踩
  • 13
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 6
    评论
专栏目录
C#:Lagrange插值法​(附完整源码)
希望我的博客，能帮上你解决学习中工作中所遇到的问题
09-15 273
C#:Lagrange插值法​(附完整源码)
用python实现分段函数_python分段函数
weixin_39630466的博客
11-20 5752
广告关闭腾讯云11.11云上盛惠 ，精选热门产品助力上云，云服务器首年88元起，买的越多返的越多，最高返5000元！segmentation fault (core dumped) python脚本: #! usrbinenv python import ctypesimport osimport sys # linux machine codebuf +=xbdxfcxa1x5dx63xd9xe...
等值线C#算法
01-07
C# 实现的等值线 分为两个版本 供大家学习 http://blog.csdn.net/hzgjf/article/details/4230109
C# 实现Excel的线性插值Forecast函数
better me
09-16 2438
最近做了一个项目涉及到插值算法，有源数据和插值后的数据，但不知道使用的具体是什么插值算法，尝试了5，6种不同的算法，但均和目标数据不匹配。后决定试试Excel里面的线性插值，一对比数据，结果发现使用的就是这个算法。 其中做excel实验的方法是模仿的这个网址： http://forums.codeguru.com/showthread.php?461020-Forecast-function
C#常用数学插值法
panjinliang066333的博客
07-31 1799
【代码】C#常用数学插值法。
C#算法完整教程专栏完整目录
希望我的博客，能帮上你解决学习中工作中所遇到的问题
02-15 1242
C#算法完整教程专栏完整目录
C#版本拉格朗日插值算法 — 用于MsChart 绘图控件进行插值算法绘图。
宋同学的Solution
03-20 3694
亲测整体可用，但是感觉结尾部分或有数据剧烈波动的地方，不是很准确。       From:http://www.cnblogs.com/technology/archive/2011/03/12/1982025.html      实际中有时需要对数据进行分析, 最近就遇到了这样一个情况: 有一系列横坐标是时间, 纵坐标是记录值的一些数据, 但横坐标却不是等距记录的. 就是说在第一分钟记录一
C#版本拉格朗日插值算法
weixin_33976072的博客
01-13 945
        实际中有时需要对数据进行分析, 最近就遇到了这样一个情况: 有一系列横坐标是时间, 纵坐标是记录值的一些数据, 但横坐标却不是等距记录的. 就是说在第一分钟记录一次, 第二分钟记录一次, 第四分钟记录一次...不等距. 需求是根据现在的这些数据进行计算, 获得等距时间所对应的记录值, 也就是说第三分钟没有记录, 得通过计算来获得第三分钟的值.         原本想把这些数据显示...
Lagrange插值算法
10-14
Lagrange插值算法属于缺失值补全算法，觉得很有意思，分享给大家，如果觉得我写的不太好可以看这位大佬写的，非常喜欢https://www.jianshu.com/p/2238ef6f3a59
拉格朗日Lagrange插值法的matlab实现
05-02
拉格朗日Lagrange插值法的matlab实现。给出一批离样点，做出一条通过这些点的光滑曲线，构造一个简单函数来近似。 本篇为Lagrange插值法，构造插值多项式。
MFC示例：Lagrange插值算法的实现
08-20
本文将详细介绍如何使用C++编程语言，结合Microsoft Foundation Classes (MFC) 框架，在Visual Studio 2008环境下实现Lagrange插值算法。 首先，理解Lagrange插值的基本概念至关重要。Lagrange插值通过构造一组...
数值计算使用Lagrange插值的示例
05-13
拉格朗日插值是一种在数值分析中广泛使用的数学方法，用于通过一组已知的数据点构建一个多项式函数，使得这个函数在每个已知数据点上都与实际的y值相匹配。这种方法对于预测和估算未知数据点的值非常有用，尤其是在...
Lagrange_Newton_Aitken插值算法.rar_Aitken_aitken插值_lagrange_newton_插
09-14
在DOS界面中实现这些插值算法，意味着用户可能需要输入数据点，然后程序会输出对应的插值结果。这通常涉及到编程，如C语言或者Fortran，它们在科学计算中被广泛使用，因为它们能够高效地处理数学运算。 在实际应用...
【数模笔记03】插值算法
2301_79843311的博客
07-17 143
也就是说，其实这个多项式是人为规定的，且通过矩阵可以证明在n阶时候的多项式唯一。
微积分-微分应用2（平均值定理）
weixin_45911156的博客
07-20 1318
平均值定理
C语言100道基础拔高题（1）
2303_78660417的博客
07-23 3325
最后按最小值，中间值，最大值的顺序打印即可。解题思路：首先输出由这几个数字所组成的所有三位数，接着再设置条件，使其输出的三位数不重复，下面我们来看下源代码。值得注意的是：所以题目的代码都是作者自行编写，如有更好的思路或者代码的优化，还请指出更正。解题思路：首先找到能整除这个正整数的最小正数，接着重置循环初值，让它继续查找最小正数，最后将得到的所有正数放到一个数组里，最后设置好格式之后循环打印出来。解题思路：首先确定好区间，接着判断是否是素数，素数就是一个大于1的数，除了1和它本身之外，不能被其他数整除。
算法 —— 暴力枚举
最新发布
fen_0108的博客
07-25 3911
第二看纵向宽发生改变，长重置为1，长为1，宽为2的时候，第一行个数为7个，总共有5 x 7 个，长为2，宽为2的时候 第一行个数为6个，共有5 x 6个……首先是横着的长方形，宽始终为1，长不断发生改变，可以看出长为2的时候，第一行个数为6个，总共有6 x 6个，长为3的时候，总共有6 x 5个……这是一道简单的模拟题，枚举出所有可能情况，不会超过规定时间的，以下附上k
2024.7.24 作业
qq_64434282的博客
07-24 2687
【代码】2024.7.24 作业。
评论 6
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

十八与她

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值