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

2024-02-03 03:30

本文主要是介绍C#函数之实现Lagrange插值算法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

博主最近想自己独立地完成一个精密单点定位的应用程序,先实现最简单的功能,那么基础工作就是要对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文件,敬请期待哦!
在这里插入图片描述

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

这篇关于C#函数之实现Lagrange插值算法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/672865

相关文章

pandas中位数填充空值的实现示例

《pandas中位数填充空值的实现示例》中位数填充是一种简单而有效的方法,用于填充数据集中缺失的值,本文就来介绍一下pandas中位数填充空值的实现,具有一定的参考价值,感兴趣的可以了解一下... 目录什么是中位数填充?为什么选择中位数填充?示例数据结果分析完整代码总结在数据分析和机器学习过程中,处理缺失数

Golang HashMap实现原理解析

《GolangHashMap实现原理解析》HashMap是一种基于哈希表实现的键值对存储结构,它通过哈希函数将键映射到数组的索引位置,支持高效的插入、查找和删除操作,:本文主要介绍GolangH... 目录HashMap是一种基于哈希表实现的键值对存储结构,它通过哈希函数将键映射到数组的索引位置,支持

Pandas使用AdaBoost进行分类的实现

《Pandas使用AdaBoost进行分类的实现》Pandas和AdaBoost分类算法,可以高效地进行数据预处理和分类任务,本文主要介绍了Pandas使用AdaBoost进行分类的实现,具有一定的参... 目录什么是 AdaBoost?使用 AdaBoost 的步骤安装必要的库步骤一:数据准备步骤二:模型

使用Pandas进行均值填充的实现

《使用Pandas进行均值填充的实现》缺失数据(NaN值)是一个常见的问题,我们可以通过多种方法来处理缺失数据,其中一种常用的方法是均值填充,本文主要介绍了使用Pandas进行均值填充的实现,感兴趣的... 目录什么是均值填充?为什么选择均值填充?均值填充的步骤实际代码示例总结在数据分析和处理过程中,缺失数

Java对象转换的实现方式汇总

《Java对象转换的实现方式汇总》:本文主要介绍Java对象转换的多种实现方式,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录Java对象转换的多种实现方式1. 手动映射(Manual Mapping)2. Builder模式3. 工具类辅助映

Go语言开发实现查询IP信息的MCP服务器

《Go语言开发实现查询IP信息的MCP服务器》随着MCP的快速普及和广泛应用,MCP服务器也层出不穷,本文将详细介绍如何在Go语言中使用go-mcp库来开发一个查询IP信息的MCP... 目录前言mcp-ip-geo 服务器目录结构说明查询 IP 信息功能实现工具实现工具管理查询单个 IP 信息工具的实现服

Python的time模块一些常用功能(各种与时间相关的函数)

《Python的time模块一些常用功能(各种与时间相关的函数)》Python的time模块提供了各种与时间相关的函数,包括获取当前时间、处理时间间隔、执行时间测量等,:本文主要介绍Python的... 目录1. 获取当前时间2. 时间格式化3. 延时执行4. 时间戳运算5. 计算代码执行时间6. 转换为指

SpringBoot基于配置实现短信服务策略的动态切换

《SpringBoot基于配置实现短信服务策略的动态切换》这篇文章主要为大家详细介绍了SpringBoot在接入多个短信服务商(如阿里云、腾讯云、华为云)后,如何根据配置或环境切换使用不同的服务商,需... 目录目标功能示例配置(application.yml)配置类绑定短信发送策略接口示例:阿里云 & 腾

Python正则表达式语法及re模块中的常用函数详解

《Python正则表达式语法及re模块中的常用函数详解》这篇文章主要给大家介绍了关于Python正则表达式语法及re模块中常用函数的相关资料,正则表达式是一种强大的字符串处理工具,可以用于匹配、切分、... 目录概念、作用和步骤语法re模块中的常用函数总结 概念、作用和步骤概念: 本身也是一个字符串,其中

python实现svg图片转换为png和gif

《python实现svg图片转换为png和gif》这篇文章主要为大家详细介绍了python如何实现将svg图片格式转换为png和gif,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录python实现svg图片转换为png和gifpython实现图片格式之间的相互转换延展:基于Py