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

相关文章

不懂推荐算法也能设计推荐系统

本文以商业化应用推荐为例,告诉我们不懂推荐算法的产品,也能从产品侧出发, 设计出一款不错的推荐系统。 相信很多新手产品,看到算法二字,多是懵圈的。 什么排序算法、最短路径等都是相对传统的算法(注:传统是指科班出身的产品都会接触过)。但对于推荐算法,多数产品对着网上搜到的资源,都会无从下手。特别当某些推荐算法 和 “AI”扯上关系后,更是加大了理解的难度。 但,不了解推荐算法,就无法做推荐系

hdu1043(八数码问题,广搜 + hash(实现状态压缩) )

利用康拓展开将一个排列映射成一个自然数,然后就变成了普通的广搜题。 #include<iostream>#include<algorithm>#include<string>#include<stack>#include<queue>#include<map>#include<stdio.h>#include<stdlib.h>#include<ctype.h>#inclu

康拓展开(hash算法中会用到)

康拓展开是一个全排列到一个自然数的双射(也就是某个全排列与某个自然数一一对应) 公式: X=a[n]*(n-1)!+a[n-1]*(n-2)!+...+a[i]*(i-1)!+...+a[1]*0! 其中,a[i]为整数,并且0<=a[i]<i,1<=i<=n。(a[i]在不同应用中的含义不同); 典型应用: 计算当前排列在所有由小到大全排列中的顺序,也就是说求当前排列是第

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

2. c#从不同cs的文件调用函数

1.文件目录如下: 2. Program.cs文件的主函数如下 using System;using System.Collections.Generic;using System.Linq;using System.Threading.Tasks;using System.Windows.Forms;namespace datasAnalysis{internal static

综合安防管理平台LntonAIServer视频监控汇聚抖动检测算法优势

LntonAIServer视频质量诊断功能中的抖动检测是一个专门针对视频稳定性进行分析的功能。抖动通常是指视频帧之间的不必要运动,这种运动可能是由于摄像机的移动、传输中的错误或编解码问题导致的。抖动检测对于确保视频内容的平滑性和观看体验至关重要。 优势 1. 提高图像质量 - 清晰度提升:减少抖动,提高图像的清晰度和细节表现力,使得监控画面更加真实可信。 - 细节增强:在低光条件下,抖

hdu1171(母函数或多重背包)

题意:把物品分成两份,使得价值最接近 可以用背包,或者是母函数来解,母函数(1 + x^v+x^2v+.....+x^num*v)(1 + x^v+x^2v+.....+x^num*v)(1 + x^v+x^2v+.....+x^num*v) 其中指数为价值,每一项的数目为(该物品数+1)个 代码如下: #include<iostream>#include<algorithm>

【数据结构】——原来排序算法搞懂这些就行,轻松拿捏

前言:快速排序的实现最重要的是找基准值,下面让我们来了解如何实现找基准值 基准值的注释:在快排的过程中,每一次我们要取一个元素作为枢纽值,以这个数字来将序列划分为两部分。 在此我们采用三数取中法,也就是取左端、中间、右端三个数,然后进行排序,将中间数作为枢纽值。 快速排序实现主框架: //快速排序 void QuickSort(int* arr, int left, int rig

【C++】_list常用方法解析及模拟实现

相信自己的力量,只要对自己始终保持信心,尽自己最大努力去完成任何事,就算事情最终结果是失败了,努力了也不留遗憾。💓💓💓 目录   ✨说在前面 🍋知识点一:什么是list? •🌰1.list的定义 •🌰2.list的基本特性 •🌰3.常用接口介绍 🍋知识点二:list常用接口 •🌰1.默认成员函数 🔥构造函数(⭐) 🔥析构函数 •🌰2.list对象

【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯,先赞后看哦~🎈🎈 🏆 作者简介:景天科技苑 🏆《头衔》:大厂架构师,华为云开发者社区专家博主,阿里云开发者社区专家博主,CSDN全栈领域优质创作者,掘金优秀博主,51CTO博客专家等。 🏆《博客》:Python全栈,前后端开发,小程序开发,人工智能,js逆向,App逆向,网络系统安全,数据分析,Django,fastapi