C#,数值计算——插值和外推,多项式插值与外推插值(Poly_interp)的计算方法与源程序

1 文本格式

using System;

namespace Legalsoft.Truffer
{
    /// <summary>
    /// 多项式插值与外推插值
    /// Polynomial Interpolation and 
    /// Extrapolation interpolation routines for one dimension
    /// </summary>
    public class Poly_interp : Base_interp
    {
        public double dy { get; set; }

        public Poly_interp(double[] xv, double[] yv, int m) : base(xv, yv[0], m)
        {
            this.dy = 0.0;
        }

        /// <summary>
        /// Given a value x, and using pointers to data xx and yy, this routine returns
        /// an interpolated value y, and stores an error estimate dy. The returned
        /// value is obtained by mm-point polynomial interpolation on the subrange
        /// xx[jl..jl + mm - 1].
        /// </summary>
        /// <param name="jl"></param>
        /// <param name="x"></param>
        /// <returns></returns>
        /// <exception cref="Exception"></exception>
        public override double rawinterp(int jl, double x)
        {
            int ns = 0;
            //double[] xa = xx.range(jl);
            //double[] ya = yy.range(jl);
            double[] c = new double[mm];
            double[] d = new double[mm];
            //double dif = Math.Abs(x - xa[0]);
            double dif = Math.Abs(x - xx[jl + 0]);
            for (int i = 0; i < mm; i++)
            {
                double dift = Math.Abs(x - xx[jl + i]);
                if ((dift) < dif)
                {
                    ns = i;
                    dif = dift;
                }
                c[i] = yy[jl + i];
                d[i] = yy[jl + i];
            }
            double y = yy[jl + ns--];
            for (int m = 1; m < mm; m++)
            {
                for (int i = 0; i < mm - m; i++)
                {
                    double ho = xx[jl + i] - x;
                    double hp = xx[jl + i + m] - x;
                    double w = c[i + 1] - d[i];
                    double den = ho - hp;
                    //if ((den) == 0.0)
                    if (Math.Abs(den) <= float.Epsilon)
                    {
                        throw new Exception("Poly_interp error");
                    }
                    den = w / den;
                    d[i] = hp * den;
                    c[i] = ho * den;
                }
                y += (dy = (2 * (ns + 1) < (mm - m) ? c[ns + 1] : d[ns--]));
            }
            return y;
        }
    }
}
 

2 代码格式

using System;namespace Legalsoft.Truffer
{/// <summary>/// 多项式插值与外推插值/// Polynomial Interpolation and /// Extrapolation interpolation routines for one dimension/// </summary>public class Poly_interp : Base_interp{public double dy { get; set; }public Poly_interp(double[] xv, double[] yv, int m) : base(xv, yv[0], m){this.dy = 0.0;}/// <summary>/// Given a value x, and using pointers to data xx and yy, this routine returns/// an interpolated value y, and stores an error estimate dy. The returned/// value is obtained by mm-point polynomial interpolation on the subrange/// xx[jl..jl + mm - 1]./// </summary>/// <param name="jl"></param>/// <param name="x"></param>/// <returns></returns>/// <exception cref="Exception"></exception>public override double rawinterp(int jl, double x){int ns = 0;//double[] xa = xx.range(jl);//double[] ya = yy.range(jl);double[] c = new double[mm];double[] d = new double[mm];//double dif = Math.Abs(x - xa[0]);double dif = Math.Abs(x - xx[jl + 0]);for (int i = 0; i < mm; i++){double dift = Math.Abs(x - xx[jl + i]);if ((dift) < dif){ns = i;dif = dift;}c[i] = yy[jl + i];d[i] = yy[jl + i];}double y = yy[jl + ns--];for (int m = 1; m < mm; m++){for (int i = 0; i < mm - m; i++){double ho = xx[jl + i] - x;double hp = xx[jl + i + m] - x;double w = c[i + 1] - d[i];double den = ho - hp;//if ((den) == 0.0)if (Math.Abs(den) <= float.Epsilon){throw new Exception("Poly_interp error");}den = w / den;d[i] = hp * den;c[i] = ho * den;}y += (dy = (2 * (ns + 1) < (mm - m) ? c[ns + 1] : d[ns--]));}return y;}}
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/162017.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【ES6.0】- Promise对象

【ES6.0】- Promise对象 文章目录 【ES6.0】- Promise对象一、概述二、Promise状态三、Promise方法3.1 Promise.prototype.then方法&#xff1a;链式操作3.2 Promise.prototype.catch方法&#xff1a;捕捉错误3.3 Promise.race方法&#xff1a;捕捉错误3.4 Promise.any()3.5 Pr…

第三节-Android10.0 Binder通信原理(三)-ServiceManager篇

1、概述 在Android中&#xff0c;系统提供的服务被包装成一个个系统级service&#xff0c;这些service往往会在设备启动之时添加进Android系统&#xff0c;当某个应用想要调用系统某个服务的功能时&#xff0c;往往是向系统发出请求&#xff0c;调用该服务的外部接口。在上一节…

广告机/商业显示屏_基于MT878安卓主板方案

安卓主板在广告机领域扮演着重要的角色。无论是在商场、车站、酒店、电梯、机场还是高铁站&#xff0c;LED广告机广泛应用&#xff0c;并通过不同方式进行播放和管理。 广告机/商业显示屏_基于MT878安卓主板方案 基于MT8788安卓主板方案的广告机采用了联发科MT8788八核芯片方案…

对比两个数组中对应位置的两个元素将每次对比的最大值用于构成新的数组np.maximum()

【小白从小学Python、C、Java】 【计算机等考500强证书考研】 【Python-数据分析】 对比两个数组中对应位置的两个元素 将每次对比的最大值用于构成新的数组 np.maximum() 选择题 以下代码的输出结果为&#xff1f; import numpy as np a1 [1,2,33] a2 [11,2,3] print("…

Axios 默认配置 简化URL 简化代码 多台服务器接口配置

main.js配置 import Axios from axios Axios.defaults.method GET//设置默认的请求类型 Axios.defaults.baseURL https://apis.jxcxin.cn/api//设置接口地址 Axios.defaults.params { token: abc } //每次请求都带上这个参数 Axios.defaults.timeout 5000 //请求的超时时间…

MATLAB - text的两种使用方法

text小技巧 1. 常规使用&#xff08;Method 1&#xff09;2. 在显示画面的相对位置&#xff08;Method 2&#xff09;3. 举个例子 1. 常规使用&#xff08;Method 1&#xff09; text(x,y,txt)2. 在显示画面的相对位置&#xff08;Method 2&#xff09; text(string,‘ABC’,…

使用websocket获取thingsboard设备的实时数据

背景 有一个读者前来咨询,如何实时获取设备的遥测数据。 其实tb是有提供websocket接口来获取设备数据的。而且还支持js跨域调用。下面给大家演示一下。 websocket地址 完整代码 <!DOCTYPE HTML> <html><h

HTTP协议和WebSocket协议之间的区别

HTTP协议和WebSocket协议之间的主要区别在于它们的设计目的和通信方式。 HTTP协议是一种无状态的协议&#xff0c;它的主要设计目的是用于从Web服务器传输超文本到本地浏览器的传输协议。HTTP协议使用请求和响应模型&#xff0c;客户端向服务器发送请求&#xff0c;服务器返回…

【Java并发编程十二】线程池

线程池 用来统一地管理线程&#xff0c;避免线程的重复创建与销毁。使用线程池可以让执行完的线程回到线程池&#xff0c;等待下一次调用。 import jdk.jshell.EvalException; import java.util.concurrent.ExecutorService; import java.util.concurrent.Executors; import j…

Matplotlib颜色条的配置_Python数据分析与可视化

Matplotlib颜色条配置 基本颜色颜色条选择配色方案颜色条刻度的限制与扩展功能的设置离散型颜色条 基本颜色 Matplotlib提供了8种指定颜色的方法&#xff1a; 在[0&#xff0c;1]中的浮点值的RGB或RGBA元组&#xff08;例如 (0.1, 0.2, 0.5) 或&#xff08;0.1&#xff0c; 0.…

C语言中文网 - Shell脚本 - 9

第1章 Shell基础(开胃菜) 9. Shell修改命令提示符 Shell 通过PS1和PS2这两个环境变量来控制提示符的格式,修改PS1和PS2的值就能修改命令提示符的格式。 PS1 控制最外层的命令提示符格式。 PS2 控制第二层的命令提示符格式。 在修改 PS1 和 PS2 之前,我们先用 echo 命令输出…

contos7中mongodb数据库无法备份与还原,数据库工具的安装

由于之前数据库没有卸载干净&#xff0c;导致直接用sudo yum install -y mongodb-org-tools命令无法直接安装&#xff0c;只能选择手动安装了。 先去官网找到mongo-tool工具 MongoDB Database Tools Downloads | MongoDB 然后复制要下载的版本的地址。 然后直接用wget来下载 …

【每日OJ —— 622. 设计循环队列】

每日OJ —— 622. 设计循环队列 1.题目&#xff1a;622. 设计循环队列2.解法2.1.解法讲解2.1.1.算法讲解2.1.2.代码实现2.1.3.提交通过展示 1.题目&#xff1a;622. 设计循环队列 2.解法 1.本题有很多解法&#xff1a;可以使用数组&#xff0c;单链表&#xff0c;双链表&#x…

2023亚太杯数学建模赛题人工精准翻译

大家好&#xff0c;亚太杯今天早上6点已经开赛啦&#xff0c;然后我在这里给大家带来赛题的精准人工翻译&#xff0c;防止大家直接用软件翻译导致某些地方乱码或者翻译不精准&#xff0c;这会导致后续做题过程出现很大偏差。 注意&#xff0c;以下翻译均免费发放word形式的哈&…

【精选】CSS入门必看知识点大合集

CSS简介 CSS概念 CSS&#xff08;Cascading Style Sheets&#xff09;层叠样式表&#xff0c;又叫级联样式表&#xff0c;简称样式表 CSS文件后缀名为.css CSS用于HTML文档中元素样式的定义 为什么需要CSS 使用css的唯一目的就是让网页具有美观一致的页面 语法 CSS 规则…

DB2—03(DB2中常见基础操作)

DB2—03&#xff08;DB2中常见基础操作&#xff09; 1. 前言1.1 oracle和mysql相关 2. db2中的"dual"2.1 SYSIBM.SYSDUMMY12.2 使用VALUES2.3 SYSIBM.SYSDUMMY1 "变" dual 3. db2中常用函数3.1 nvl()、value()、COALESCE()3.2 NULLIF() 函数3.3 LISTAGG() …

论文《Unsupervised Dialog Structure Learning》笔记:详解DD-VRNN

D-VRNN模型和DD-VRNN模型 总体架构 离散-可变循环变分自编码器&#xff08;D-VRNN&#xff09;和直接-离散-可变循环变分自编码器&#xff08;DD-VRNN&#xff09;概述。D-VRNN和DD-VRNN使用不同的先验分布来建模 z t z_t zt​之间的转换&#xff0c;如红色实线所示。 x t x_t…

Rust错误处理:Result

文章目录 简介错误匹配 Rust基础教程&#xff1a; 初步⚙ 所有权⚙ 结构体和枚举类⚙ 函数进阶⚙ 泛型和特征⚙ 并发和线程通信⚙ cargo包管理⚙ 可空类型Option Rust进阶教程&#xff1a; 用宏实现参数可变的函数⚙ 类函数宏 简介 Rust中没有提供类似try…catch之类…

qt双击treeview节点之后,完成编辑,获取完成编辑得信号

Qt 中&#xff0c;要获取 QTreeView 节点完成编辑的信号&#xff0c;您可以使用 QTreeView 的 doubleClicked 信号以及重写相应的编辑触发函数。通常情况下&#xff0c;完成编辑后会触发 closeEditor 信号。 以下是一个示例代码&#xff0c;演示了如何实现这个功能&#xff1a…