导数应用(一):差分计算(导数)

2024-06-16 22:08
文章标签 计算 应用 导数 差分

本文主要是介绍导数应用(一):差分计算(导数),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

导数应用(一):差分计算(导数)

  • 1.数学背景
  • 2.代码

1.数学背景

导数: d y d x = y ( x i ) − y ( x i − 1 ) x i − x x − i \frac{dy}{dx} = \frac{y(x_i) - y(x_{i-1})}{x_i - x_{x-i}} dxdy=xixxiy(xi)y(xi1)

差分: Δ Y Δ X = Y i − Y i − 1 X i − X i − 1 \frac{\Delta Y}{\Delta X} = \frac{Y_i- Y_{i-1}}{X_i - X_{i-1}} ΔXΔY=XiXi1YiYi1

然而由于数据特性,大多数情况 Δ X \Delta X ΔX 远远无法逼近 d x dx dx
由此衍生出各种差分形式,但主要思想万变不离其宗

2.代码

按照惯例,先放图:

  • 蓝色曲线为 y 1 = s i n ( x ) y_1 = sin(x) y1=sin(x)
  • 黄色曲线为 y 2 = ( s i n ( x ) ) ′ = c o s ( x ) y_2 = (sin(x))' = cos(x) y2=(sin(x))=cos(x)
  • 红色曲线为 y 3 = ( s i n ( x ) ) ′ ′ = − s i n ( x ) y_3 = (sin(x))'' = -sin(x) y3=(sin(x))=sin(x)

在这里插入图片描述

import matplotlib.pyplot as plt
import pandas as pd
import numpy as np# 一阶中心差分
def diff_1st_mid(series,dx):arr = []for i in range(0,len(series),1):if i ==0:arr.append((series[1] - series[0])/dx)elif i== len(series)-1:arr.append((series[len(series)-1] - series[len(series)-2])/dx)else:arr.append((series[i+1] - series[i-1])/ (2 * dx) )return arr# 二阶中心差分
def diff_2nd_mid(series,dx):arr = []dx2 = dx * dxfor i in range(0,len(series),1):if i ==0:arr.append((series[2] -2*series[1]  + series[0])/dx2)elif i== len(series)-1:arr.append((series[len(series)-1] - 2*series[len(series)-2] +series[len(series)-3])/dx2)else:arr.append((series[i+1] - 2 *series[i] +  series[i-1])/ (dx2) )return arrX=np.linspace(-np.pi,np.pi,1000,endpoint=True)Y = np.sin(X)
Diff1 = np.cos(X)DD1 = diff_1st_mid(Y,np.pi * 2 /(len(X)))
DD2 = diff_2nd_mid(Y,np.pi * 2 /(len(X)))plt.plot(X,Y)plt.plot(X,DD1, color = 'yellow', linewidth=1, linestyle="--")
plt.plot(X,DD2, color = 'red', linewidth=1, linestyle="--")print(DD1)plt.show()

这篇关于导数应用(一):差分计算(导数)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

PostgreSQL的扩展dict_int应用案例解析

《PostgreSQL的扩展dict_int应用案例解析》dict_int扩展为PostgreSQL提供了专业的整数文本处理能力,特别适合需要精确处理数字内容的搜索场景,本文给大家介绍PostgreS... 目录PostgreSQL的扩展dict_int一、扩展概述二、核心功能三、安装与启用四、字典配置方法

Python中re模块结合正则表达式的实际应用案例

《Python中re模块结合正则表达式的实际应用案例》Python中的re模块是用于处理正则表达式的强大工具,正则表达式是一种用来匹配字符串的模式,它可以在文本中搜索和匹配特定的字符串模式,这篇文章主... 目录前言re模块常用函数一、查看文本中是否包含 A 或 B 字符串二、替换多个关键词为统一格式三、提

Java MQTT实战应用

《JavaMQTT实战应用》本文详解MQTT协议,涵盖其发布/订阅机制、低功耗高效特性、三种服务质量等级(QoS0/1/2),以及客户端、代理、主题的核心概念,最后提供Linux部署教程、Sprin... 目录一、MQTT协议二、MQTT优点三、三种服务质量等级四、客户端、代理、主题1. 客户端(Clien

CSS中的Static、Relative、Absolute、Fixed、Sticky的应用与详细对比

《CSS中的Static、Relative、Absolute、Fixed、Sticky的应用与详细对比》CSS中的position属性用于控制元素的定位方式,不同的定位方式会影响元素在页面中的布... css 中的 position 属性用于控制元素的定位方式,不同的定位方式会影响元素在页面中的布局和层叠关

SpringBoot3应用中集成和使用Spring Retry的实践记录

《SpringBoot3应用中集成和使用SpringRetry的实践记录》SpringRetry为SpringBoot3提供重试机制,支持注解和编程式两种方式,可配置重试策略与监听器,适用于临时性故... 目录1. 简介2. 环境准备3. 使用方式3.1 注解方式 基础使用自定义重试策略失败恢复机制注意事项

Python并行处理实战之如何使用ProcessPoolExecutor加速计算

《Python并行处理实战之如何使用ProcessPoolExecutor加速计算》Python提供了多种并行处理的方式,其中concurrent.futures模块的ProcessPoolExecu... 目录简介完整代码示例代码解释1. 导入必要的模块2. 定义处理函数3. 主函数4. 生成数字列表5.

Python使用Tkinter打造一个完整的桌面应用

《Python使用Tkinter打造一个完整的桌面应用》在Python生态中,Tkinter就像一把瑞士军刀,它没有花哨的特效,却能快速搭建出实用的图形界面,作为Python自带的标准库,无需安装即可... 目录一、界面搭建:像搭积木一样组合控件二、菜单系统:给应用装上“控制中枢”三、事件驱动:让界面“活”

如何确定哪些软件是Mac系统自带的? Mac系统内置应用查看技巧

《如何确定哪些软件是Mac系统自带的?Mac系统内置应用查看技巧》如何确定哪些软件是Mac系统自带的?mac系统中有很多自带的应用,想要看看哪些是系统自带,该怎么查看呢?下面我们就来看看Mac系统内... 在MAC电脑上,可以使用以下方法来确定哪些软件是系统自带的:1.应用程序文件夹打开应用程序文件夹

Python Flask 库及应用场景

《PythonFlask库及应用场景》Flask是Python生态中​轻量级且高度灵活的Web开发框架,基于WerkzeugWSGI工具库和Jinja2模板引擎构建,下面给大家介绍PythonFl... 目录一、Flask 库简介二、核心组件与架构三、常用函数与核心操作 ​1. 基础应用搭建​2. 路由与参

Spring Boot中的YML配置列表及应用小结

《SpringBoot中的YML配置列表及应用小结》在SpringBoot中使用YAML进行列表的配置不仅简洁明了,还能提高代码的可读性和可维护性,:本文主要介绍SpringBoot中的YML配... 目录YAML列表的基础语法在Spring Boot中的应用从YAML读取列表列表中的复杂对象其他注意事项总