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

2024-02-03 03:30
文章标签 c# 算法 实现 函数 .net lagrange 插值 netcore

本文主要是介绍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。

相关文章

Nginx 配置跨域的实现及常见问题解决

Nginx 配置跨域的实现及常见问题解决

《Nginx配置跨域的实现及常见问题解决》本文主要介绍了Nginx配置跨域的实现及常见问题解决,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来... 目录1. 跨域1.1 同源策略1.2 跨域资源共享（CORS）2. Nginx 配置跨域的场景2.1
阅读更多...
Python中提取文件名扩展名的多种方法实现

Python中提取文件名扩展名的多种方法实现

《Python中提取文件名扩展名的多种方法实现》在Python编程中,经常会遇到需要从文件名中提取扩展名的场景,Python提供了多种方法来实现这一功能,不同方法适用于不同的场景和需求,包括os.pa... 目录技术背景实现步骤方法一：使用os.path.splitext方法二：使用pathlib模块方法三
阅读更多...
CSS实现元素撑满剩余空间的五种方法

CSS实现元素撑满剩余空间的五种方法

《CSS实现元素撑满剩余空间的五种方法》在日常开发中，我们经常需要让某个元素占据容器的剩余空间，本文将介绍5种不同的方法来实现这个需求，并分析各种方法的优缺点，感兴趣的朋友一起看看吧... css实现元素撑满剩余空间的5种方法 在日常开发中，我们经常需要让某个元素占据容器的剩余空间。这是一个常见的布局需求
阅读更多...
HTML5 getUserMedia API网页录音实现指南示例小结

HTML5 getUserMedia API网页录音实现指南示例小结

《HTML5getUserMediaAPI网页录音实现指南示例小结》本教程将指导你如何利用这一API，结合WebAudioAPI，实现网页录音功能,从获取音频流到处理和保存录音，整个过程将逐步... 目录1. html5 getUserMedia API简介1.1 API概念与历史1.2 功能与优势1.3
阅读更多...
Java实现删除文件中的指定内容

Java实现删除文件中的指定内容

《Java实现删除文件中的指定内容》在日常开发中,经常需要对文本文件进行批量处理,其中,删除文件中指定内容是最常见的需求之一,下面我们就来看看如何使用java实现删除文件中的指定内容吧... 目录1. 项目背景详细介绍2. 项目需求详细介绍2.1 功能需求2.2 非功能需求3. 相关技术详细介绍3.1 Ja
阅读更多...
使用Python和OpenCV库实现实时颜色识别系统

使用Python和OpenCV库实现实时颜色识别系统

《使用Python和OpenCV库实现实时颜色识别系统》：本文主要介绍使用Python和OpenCV库实现的实时颜色识别系统,这个系统能够通过摄像头捕捉视频流,并在视频中指定区域内识别主要颜色（红... 目录一、引言二、系统概述三、代码解析1. 导入库2. 颜色识别函数3. 主程序循环四、HSV色彩空间详解
阅读更多...
PostgreSQL中MVCC 机制的实现

PostgreSQL中MVCC 机制的实现

《PostgreSQL中MVCC机制的实现》本文主要介绍了PostgreSQL中MVCC机制的实现,通过多版本数据存储、快照隔离和事务ID管理实现高并发读写,具有一定的参考价值,感兴趣的可以了解一下... 目录一 MVCC 基本原理python1.1 MVCC 核心概念1.2 与传统锁机制对比二 Postg
阅读更多...
SpringBoot整合Flowable实现工作流的详细流程

SpringBoot整合Flowable实现工作流的详细流程

《SpringBoot整合Flowable实现工作流的详细流程》Flowable是一个使用Java编写的轻量级业务流程引擎,Flowable流程引擎可用于部署BPMN2.0流程定义,创建这些流程定义的... 目录1、流程引擎介绍2、创建项目3、画流程图4、开发接口4.1 Java 类梳理4.2 查看流程图4
阅读更多...
C++中零拷贝的多种实现方式

C++中零拷贝的多种实现方式

《C++中零拷贝的多种实现方式》本文主要介绍了C++中零拷贝的实现示例,旨在在减少数据在内存中的不必要复制,从而提高程序性能、降低内存使用并减少CPU消耗,零拷贝技术通过多种方式实现,下面就来了解一下... 目录一、C++中零拷贝技术的核心概念二、std::string_view 简介三、std::stri
阅读更多...
C++高效内存池实现减少动态分配开销的解决方案

C++高效内存池实现减少动态分配开销的解决方案

《C++高效内存池实现减少动态分配开销的解决方案》C++动态内存分配存在系统调用开销、碎片化和锁竞争等性能问题,内存池通过预分配、分块管理和缓存复用解决这些问题,下面就来了解一下... 目录一、C++内存分配的性能挑战二、内存池技术的核心原理三、主流内存池实现：TCMalloc与Jemalloc1. TCM
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2