C#,数值计算,求解微分方程的预测校正法(修正欧拉法)算法与源代码

Leonhard Euler

1 微分方程

微分方程,是指含有未知函数及其导数的关系式。解微分方程就是找出未知函数。
微分方程是伴随着微积分学一起发展起来的。微积分学的奠基人Newton和Leibniz的著作中都处理过与微分方程有关的问题。微分方程的应用十分广泛,可以解决许多与导数有关的问题。物理中许多涉及变力的运动学、动力学问题,如空气的阻力为速度函数的落体运动等问题,很多可以用微分方程求解。此外,微分方程在化学、工程学、经济学和人口统计等领域都有应用。
数学领域对微分方程的研究着重在几个不同的面向,但大多数都是关心微分方程的解。只有少数简单的微分方程可以求得解析解。不过即使没有找到其解析解,仍然可以确认其解的部分性质。在无法求得解析解时,可以利用数值分析的方式,利用电脑来找到其数值解。 动力系统理论强调对于微分方程系统的量化分析,而许多数值方法可以计算微分方程的数值解,且有一定的准确度。
 

预估-校正法

对于给定的带有初始条件y(x0)=y0的微分方程dy/dx=f(x, y),用预测-校正法求近似解。
预估-校正法: 预估-校正法又称修正-欧拉法。
在欧拉方法中,在一个点绘制切线,并计算给定步长的斜率。
因此,这种方法对于线性函数效果最好,但是对于其他情况,仍然存在截断误差。
为了解决这个问题,引入了修正欧拉法。

3 源程序

using System;
using System.Collections;
using System.Collections.Generic;

namespace Legalsoft.Truffer.Algorithm
{
    /// <summary>
    /// 给定微分方程的一阶偏导方程
    /// </summary>
    /// <param name="x"></param>
    /// <param name="y"></param>
    /// <returns></returns>
    public delegate double SDE_Equation(double x, double y);

    /// <summary>
    /// 求解微分方程的预测校正法或修正欧拉法
    /// </summary>
    public static partial class Algorithm_Gallery
    {
        public static SDE_Equation dydx = null;

        private static double Predict(double x, double y, double h)
        {
            double y1p = y + h * dydx(x, y);
            return y1p;
        }

        private static double Correct(double x, double y, double x1, double y1, double h)
        {
            double e = 0.00001;
            double y1c = y1;

            do
            {
                y1 = y1c;
                y1c = y + 0.5 * h * (dydx(x, y) + dydx(x1, y1));
            }
            while (Math.Abs(y1c - y1) > e);

            return y1c;
        }

        public static double DE_Modified_Euler_Solve(double x, double xn, double y, double h)
        {
            while (x < xn)
            {
                double x1 = x + h;
                double y1p = Predict(x, y, h);
                double y1c = Correct(x, y, x1, y1p, h);
                x = x1;
                y = y1c;
            }
            return y;
        }
    }
}

 POWER BY TRUFFER.CN

4 源代码

using System;
using System.Collections;
using System.Collections.Generic;namespace Legalsoft.Truffer.Algorithm
{/// <summary>/// 给定微分方程的一阶偏导方程/// </summary>/// <param name="x"></param>/// <param name="y"></param>/// <returns></returns>public delegate double SDE_Equation(double x, double y);/// <summary>/// 求解微分方程的预测校正法或修正欧拉法/// </summary>public static partial class Algorithm_Gallery{public static SDE_Equation dydx = null;private static double Predict(double x, double y, double h){double y1p = y + h * dydx(x, y);return y1p;}private static double Correct(double x, double y, double x1, double y1, double h){double e = 0.00001;double y1c = y1;do{y1 = y1c;y1c = y + 0.5 * h * (dydx(x, y) + dydx(x1, y1));}while (Math.Abs(y1c - y1) > e);return y1c;}public static double DE_Modified_Euler_Solve(double x, double xn, double y, double h){while (x < xn){double x1 = x + h;double y1p = Predict(x, y, h);double y1c = Correct(x, y, x1, y1p, h);x = x1;y = y1c;}return y;}}
}

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

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

相关文章

To 有缘看到的朋友,To myself

To 有缘看到的朋友&#xff0c;To myself 零、00时光宝盒 我们生而为人&#xff0c;而不是什么神仙妖怪&#xff0c;自然逃不脱凡尘种种不易。 世界并不完美&#xff0c;面对很多事情我们都很无奈甚至悲哀&#xff0c;但生活总要继续下去&#xff0c;当困难悄悄地来临&#xff…

【vue3之组合式API】

组合式API 一、setup1.写法2.如何访问3.语法糖4.同步返回对象 二、reactive()和ref()1.reactive()2.ref() 三、computed四、watch函数1侦听单个数据2.侦听多个数据3. immediate4. deep5.精确侦听对象的某个属性 五、生命周期函数六、组件通信1.父传子2. 子传父 七、模版引用1. …

shell脚本一键部署docker

Docker介绍 Docker 是一个开源的平台&#xff0c;用于开发、交付和运行应用程序。它利用容器化技术&#xff0c;可以帮助开发人员更轻松地打包应用程序及其依赖项&#xff0c;并将其部署到任何环境中&#xff0c;无论是开发工作站、数据中心还是云中。以下是 Docker 的一些关键…

declare module ‘global‘

在 TypeScript 中&#xff0c;declare module global 是一个特殊的模块声明&#xff0c;用于扩展全局作用域中的对象或接口。这通常用于为浏览器环境中的全局对象&#xff08;如 window&#xff09;或 Node.js 环境中的全局对象&#xff08;如 global&#xff09;添加类型声明。…

vue模拟数据mock放哪

在Vue中进行模拟数据的mock可以放在多个地方&#xff0c;具体取决于需求和个人偏好。以下是几种常见的方式&#xff1a; 在组件内部&#xff1a;可以在Vue组件的methods中编写模拟数据的逻辑。这种方式适用于只在当前组件中需要使用模拟数据的情况。 使用Mock.js库&#xff1a…

【Linux】软件包管理器yum

目录 一、yum是什么&#xff1f; 二、查看软件包 三、安装与卸载软件 1、如何安装软件 2、如何卸载软件 四、yum源的配置 一、yum是什么&#xff1f; 在Linux下安装软件, 一个通常的办法是下载到程序的源代码, 并进行编译, 得到可执行程序. 但是这样太麻烦了, 于是有些人…

Webserve(4): HTTP解析

// 写HTTP响应 bool http_conn::write() {int temp 0;if ( bytes_to_send 0 ) {// 将要发送的字节为0&#xff0c;这一次响应结束。modfd( m_epollfd, m_sockfd, EPOLLIN ); init();return true;}while(1) {// 分散写temp writev(m_sockfd, m_iv, m_iv_count);if ( temp <…

如何在华为云服务器部署安防监控EasyCVR平台?

随着视频技术的快速发展&#xff0c;安防视频汇聚平台EasyCVR可支持的协议也在不断拓展&#xff0c;平台兼容多类型的协议接入&#xff0c;包括&#xff1a;国标GB28181、RTSP/Onvif、RTMP&#xff0c;以及厂家的私有协议与SDK&#xff0c;如&#xff1a;海康ehome、海康sdk、大…

YOLOv7创新改进:SPPF创新涨点篇 | SPPELAN:SPP创新结合ELAN ,效果优于SPP、SPPF| YOLOv9

💡💡💡本文独家改进:新颖SPPF创新涨点改进,SPP创新结合ELAN,来自于YOLOv9,助力YOLOv7,将SPPELAN代替原始的SPPF 💡💡💡在多个私有数据集和公开数据集VisDrone2019、PASCAL VOC实现涨点 收录 YOLOv7原创自研 https://blog.csdn.net/m0_63774211/category…

C语言常见关键字:一文打尽

关键字 1. 前言2. 什么是关键字3. extern-声明外部符号4. auto-自动5. typedef-类型重定义&#xff08;类型重命名&#xff09;6. register-寄存器6.1 存储器6.2 register关键字的作用 7. static-静态7.1 static修饰局部变量7.1.1 代码对比7.1.2 原理分析 7.2 static修饰全局变…

Java中常见延时队列的实现方案总结

&#x1f3f7;️个人主页&#xff1a;牵着猫散步的鼠鼠 &#x1f3f7;️系列专栏&#xff1a;Java全栈-专栏 &#x1f3f7;️个人学习笔记&#xff0c;若有缺误&#xff0c;欢迎评论区指正 前些天发现了一个巨牛的人工智能学习网站&#xff0c;通俗易懂&#xff0c;风趣幽默&…

41. 【Linux教程】日志文件系统

前小节介绍了早期的 Linux 文件系统&#xff0c;本小节介绍 Linux 日志文件系统&#xff0c;日志文件系统是给 Linux 操作系统增加了一些安全相关能力&#xff0c;日志文件系统和早期的文件系统相比&#xff0c;它并不是先把数据写到中间存储设备&#xff0c;然后再写到存储设备…

2024/3/6打卡最短编辑距离---线性DP

题目&#xff1a; 给定两个字符串 A 和 B&#xff0c;现在要将 A 经过若干操作变为 B&#xff0c;可进行的操作有&#xff1a; 删除–将字符串 A 中的某个字符删除。插入–在字符串 A 的某个位置插入某个字符。替换–将字符串 A 中的某个字符替换为另一个字符。 现在请你求出&a…

springboot基于java的中医院门诊挂号诊断就诊系统ssm+jsp

主要研究内容&#xff1a; 医院门诊挂号系统分为护士&#xff0c;医生&#xff0c;药房&#xff0c;收费&#xff0c;管理员等权限。 护士&#xff1a;挂号、退号、查询病人。挂号——就诊科室(发热门诊、骨科、妇科等等)&#xff0c;就诊医生数据库获取&#xff0c;挂号类型—…

测试遍历1e5,1e8数组耗时

1e8大概0.38秒&#xff0c;即380ms 1e5耗时1ms左右&#xff1a; 代码使用方式来自&#xff1a;clock - C Reference (cplusplus.com)

Java 与 JavaScript 区别

Java 与 JavaScript 区别 制造商不同 Java 是 SUN Microsystems 公司推出的新一代面向对象的程序设计语言。javascript 是 Netscape 公司的产品&#xff0c;其目的是为了扩展 Netscape Navigator 功能&#xff0c;而开发的一种可以嵌入 Web 页面中的基于对象和事件驱动的解释…

python_读取txt文件绘制多条曲线III

先把文件中指定列&#xff0c;去重提取出来&#xff0c;然后根据指定列去匹配数据&#xff0c;最后完成多条数据的绘图&#xff1b; import matplotlib.pyplot as plt import re from datetime import datetime from pylab import mplmpl.rcParams["font.sans-serif"…

【python 1】----Pytest基础知识

定义 用于编写和执行Python测试全功能测试框架&#xff08;工具&#xff09;&#xff0c;是一个第三方库 安装 pip insatll pytest 安装pytest --version 校验 pytest的组成构成 不写调用语句也可以执行函数内容 在用例运行语句里面&#xff1a; -s:指的是开启与终端的…

vue3 setup函数与setup语法糖之间的区别

普通setup函数构建的组件 import {ref} from vueexport default {setup(){const countref(0)const handleUpdate()>{count.value}return{count,handleUpdate}}} </script>使用setup语法糖构建的组件 <script setup>import {ref} from vueconst countref(0)con…

devc++小游戏3.8.5

导航&#xff1a; Dev-c跑酷小游戏 1.0.0 devc跑酷小游戏1.2.5 devc跑酷游戏1.2.6 devc跑酷游戏2.0.0 devc跑酷游戏2.0.1 devc跑酷游戏2.4.0 devc跑酷小游戏3.5.0 更新内容 重磅回归&#xff0c;存档搞定&#xff01;&#xff01;&#xff01; 每一关需要前一关已…