标签：plt 画出 初值问题 step Delta 微分方程 欧拉

微分方程初值问题

初值问题\(\begin{cases}y^{\prime}=f(x, y)\\ y(x_{0})=y_{0}\end{cases}\)的解\(y=y(x)\)代表通过点\((x_0, y_0)\)的一条称为微分方程的积分曲线。积分曲线上的每一个点\((x, y)\)的切线斜率等于函数\(y^{\prime}\)在这点的值.

欧拉方法画出函数图像

在最一开始的\((x_0, y_0)\)位置开始，每次我们都给\(x\)一个步长\(\Delta x\)，如果这个步长比较小我们就可以认为\((x, g(x))\)和\((x+\Delta x, g(x+\Delta x))\)两点构成的直线的斜率近似为函数在\(x\)处的导数。

这样每次都画出这样一条“线段”，然后\(x\)跳到\(x+\Delta x\)这个位置，重复上面的过程就可以近似的画出这条直线的图像。

代码实现

import matplotlib.pyplot as plt

x, y = 0.0, 1.0
step = 0.1
X, Y = [x], [y]

def f(x:float, y:float):
    return y - 2 * x / y

while x < 10.0:
    k = f(x, y)
    y = y + k * step
    x = x + step
    X.append(x)
    Y.append(y)

print(X, Y, sep='\n')
plt.plot(X, Y)
plt.show()

标签：plt,画出,初值问题,step,Delta,微分方程,欧拉
来源： https://www.cnblogs.com/chantmee/p/15359882.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。