当前位置： 首页 > wzjs >正文

开发技术网站开发技术路线俄罗斯搜索引擎入口 yandex

wzjs 2025/8/3 5:54:37
开发技术网站开发技术路线,俄罗斯搜索引擎入口 yandex,课程视频网站建设的必要性,什么视频直播网站做挣钱微分方程在自然科学、工程技术、社会科学等多个领域都有着广泛而重要的应用。而求解微分方程是数学与应用数据领域一大难题，对于一些复杂的微分方程无法通过计算推导计算其精确的方程表达式与 结果，因此，我们通过数学理论。迭代，微…

        微分方程在自然科学、工程技术、社会科学等多个领域都有着广泛而重要的应用。而求解微分方程是数学与应用数据领域一大难题，对于一些复杂的微分方程无法通过计算推导计算其精确的方程表达式与 结果，因此，我们通过数学理论。迭代，微分约等于思想寻求其精确解，随着计算机的快速发展，利用计算机程序求解微分方程非常重要，下面我们基于python代码分别利用显式欧拉方法、梯形法、改进欧拉方法、二阶龙格库塔方法、四阶龙格库塔方法求解微分方程。

1.问题描述

 已知微分方程如下

\begin{cases} y'=-2y-4x, & \ 0 \leq x \leq 1 \\ y(0)=2 \end{cases}

其方程精确解如下：

$ y(x)=e^{-2x}-2x+1 $

利用python编程技术，取步长h=0.1分别1.用显示欧拉方法，2.梯形法，3.改进欧拉方法，4.二阶龙格库塔方法，5.四阶龙格库塔方法计算数值解，并于精确解比较。

2.显式欧拉方法

2.1显式欧拉方法原理

        对于一阶常微分方程初值问题\frac{dy}{dt}=f(t,y)，y(t_0)=y_0其基本思想是在每个小区间上，用差商来近似代替导数。具体来说，在t_n时刻，已知y(t_n)=y_n，通过向前差分来近似t_{n + 1}=t_n + h时刻的y值，其中h为时间步长。根据导数的定义\frac{dy}{dt}\approx\frac{y_{n + 1}-y_n}{h}，将\frac{dy}{dt}=f(t,y)代入可得y_{n + 1}=y_n + hf(t_n,y_n)这就是显式欧拉方法的迭代公式，它通过已知的t_n和y_n，以及函数f(t,y)来直接计算出y_{n + 1}的值，从而逐步推进求解常微分方程在不同时刻的数值解。

2.2 python代码实现

import numpy as np
import matplotlib.pyplot as plt# 定义微分方程
def f(x, y):return -2 * y - 4 * x# 定义精确解
def exact_solution(x):return np.exp(-2 * x) - 2 * x + 1# 显式欧拉方法
def euler_method(x0, y0, h, n):x = np.zeros(n + 1)y = np.zeros(n + 1)x[0] = x0y[0] = y0for i in range(n):y[i + 1] = y[i] + h * f(x[i], y[i])x[i + 1] = x[i] + hreturn x, y
# 参数设置
x0 = 0
y0 = 2
h = 0.1
n = int(1 / h)# 计算数值解
x_euler, y_euler = euler_method(x0, y0, h, n)
# 计算精确解
x_exact = np.linspace(x0, x0 + 1, n + 1)
y_exact = exact_solution(x_exact)# 打印结果并绘图比较
print("显式欧拉方法结果：", y_euler)
plt.plot(x_exact, y_exact, label='Exact Solution')
plt.plot(x_euler, y_euler, 'o-', label='Euler Method')
plt.xlabel('x')
plt.ylabel('y')
plt.title('Numerical Solutions Comparison')
plt.legend()
plt.show()

运行结果

3.梯形法

3.1梯形法原理

考虑一阶常微分方程的初值问题：\frac{dy}{dt}=f(t,y)，y(t_0)=y_0 其中f(t,y) 是已知的函数，t_0和 y_0 分别是初始时刻和初始条件。在数值求解过程中，我们将时间区间进行离散化，设时间步长为 h，t_{n + 1}=t_n + h 。

显式欧拉方法是基于向前差分近似导数\frac{dy}{dt}\approx\frac{y_{n + 1}-y_n}{h}得到y_{n + 1}=y_n+hf(t_n,y_n)，它仅使用了当前时刻 t_n的斜率 f(t_n,y_n)。

而梯形法结合了当前时刻 t_n 和下一时刻t_{n + 1} 的斜率信息。根据积分中值定理，y(t_{n + 1})-y(t_n)=\int_{t_n}^{t_{n + 1}}f(t,y(t))dt

梯形法使用梯形公式来近似这个积分，即\int_{t_n}^{t_{n + 1}}f(t,y(t))dt\approx\frac{h}{2}[f(t_n,y(t_n)) + f(t_{n + 1},y(t_{n + 1}))]。

由此得到梯形法的迭代公式：y_{n + 1}=y_n+\frac{h}{2}[f(t_n,y_n)+f(t_{n + 1},y_{n + 1})]

        与显式欧拉方法不同，梯形法的迭代公式是隐式的，因为y_{n + 1} 同时出现在等式的两边。这意味着在每一步求解 y_{n + 1}时，通常需要求解一个非线性方程。

3.2 python代码实现

import numpy as np
import matplotlib.pyplot as plt# 定义微分方程
def f(x, y):return -2 * y - 4 * x# 定义精确解
def exact_solution(x):return np.exp(-2 * x) - 2 * x + 1# 梯形法
def trapezoid_method(f, x0, y0, h, n):x_values = [x0]y_values = [y0]x = x0y = y0for _ in range(n):def g(y_next):return y_next - y - 0.5 * h * (f(x, y) + f(x + h, y_next))from scipy.optimize import rootsol = root(g, y)y = sol.x[0]x = x + hx_values.append(x)y_values.append(y)return x_values, y_values
# 参数设置
x0 = 0
y0 = 2
h = 0.1
n = int(1 / h)# 计算数值解
x_trapezoid, y_trapezoid = trapezoid_method(f, x0, y0, h, n)
# 计算精确解
x_exact = np.linspace(x0, x0 + 1, n + 1)
y_exact = exact_solution(x_exact)# 打印结果并绘图比较
print("梯形方法结果：", y_trapezoid)
plt.plot(x_exact, y_exact, label='Exact Solution')
plt.plot(x_trapezoid, y_trapezoid, 'o-', label='Trapezoid Method')
plt.xlabel('x')
plt.ylabel('y')
plt.title('Numerical Solutions Comparison')
plt.legend()
plt.show()

运行结果

 4.改进欧拉方法

4.1改进欧拉方法原理

改进欧拉方法的核心思想是综合考虑当前点(t_n,y_n)和预测的下一个点\(t_{n + 1},\overline{y}_{n + 1})处的斜率信息，通过取二者的平均值作为这一步的平均斜率，以此来得到更准确的y_{n + 1}的估计值。

具体原理和步骤

改进欧拉方法采用了预测 - 校正的策略，具体步骤如下：

  1. 预测阶段 利用显式欧拉方法对t_{n + 1}时刻的y值进行初步预测。设时间步长为h，t_{n + 1}=t_n + h，预测公式为：\overline{y}_{n + 1}=y_n+hf(t_n,y_n)这里\overline{y}_{n + 1}是对y(t_{n + 1})的一个初步估计，它基于当前点(t_n,y_n)处的斜率f(t_n,y_n)得到。
  2. 校正阶段 计算t_n时刻的斜率f(t_n,y_n))和预测的t_{n + 1}时刻的斜率f(t_{n + 1},\overline{y}_{n + 1})，然后取这两个斜率的平均值作为这一步的平均斜率，再用这个平均斜率来校正y_{n + 1}的值。校正公式为： y_{n + 1}=y_n+\frac{h}{2}[f(t_n,y_n)+f(t_{n + 1},\overline{y}_{n + 1})]从几何角度看，改进欧拉方法不再仅仅依赖于当前点的斜率，而是综合考虑了当前点和预测点的斜率，用这两个斜率的平均值来确定在区间[t_n,t_{n + 1}]上的近似直线的斜率，从而使近似曲线更接近真实的解曲线。

4.2 python代码实现

import numpy as np
import matplotlib.pyplot as plt# 定义微分方程
def f(x, y):return -2 * y - 4 * x# 定义精确解
def exact_solution(x):return np.exp(-2 * x) - 2 * x + 1# 改进的欧拉方法
def improved_euler_method(x0, y0, h, n):x = np.zeros(n + 1)y = np.zeros(n + 1)x[0] = x0y[0] = y0for i in range(n):k1 = f(x[i], y[i])k2 = f(x[i] + h, y[i] + h * k1)y[i + 1] = y[i] + h * (k1 + k2) / 2x[i + 1] = x[i] + hreturn x, y
# 参数设置
x0 = 0
y0 = 2
h = 0.1
n = int(1 / h)# 计算数值解
x_euler, y_euler =improved_euler_method(x0, y0, h, n)
# 计算精确解
x_exact = np.linspace(x0, x0 + 1, n + 1)
y_exact = exact_solution(x_exact)# 打印结果并绘图比较
print("改进的欧拉方法：", y_euler)
plt.plot(x_exact, y_exact, label='Exact Solution')
plt.plot(x_euler, y_euler, 'o-', label='Improved Euler Method')
plt.xlabel('x')
plt.ylabel('y')
plt.title('Numerical Solutions Comparison')
plt.legend()
plt.show()

运行结果

5.二阶龙格库塔方法

5.1二阶龙格库塔方法原理

二阶龙格 - 库塔方法的基本思想是在每个时间步长内，通过计算两个不同点的斜率估计值，然后对这些斜率进行加权平均，以此来得到更准确的\(y_{n + 1}\)的估计值。

二阶龙格 - 库塔方法的一般形式可以表示为：

k_1=hf(t_n,y_n)

k_2=hf(t_n+\alpha h,y_n+\beta k_1)

y_{n + 1}=y_n+(c_1k_1 + c_2k_2)

其中\alpha、\beta、c_1和c_2是待确定的参数，需要满足一定的条件以保证方法具有二阶精度。

5.2 python代码实现

import numpy as np
import matplotlib.pyplot as plt# 定义微分方程
def f(x, y):return -2 * y - 4 * x# 定义精确解
def exact_solution(x):return np.exp(-2 * x) - 2 * x + 1# 二阶龙格-库塔方法
def second_order_rk_method(x0, y0, h, n):x = np.zeros(n + 1)y = np.zeros(n + 1)x[0] = x0y[0] = y0for i in range(n):k1 = h * f(x[i], y[i])k2 = h * f(x[i] + h / 2, y[i] + k1 / 2)y[i + 1] = y[i] + k2x[i + 1] = x[i] + hreturn x, y
# 计算精确解
x_exact = np.linspace(x0, x0 + 1, n + 1)
y_exact = exact_solution(x_exact)
x_second_order_rk, y_second_order_rk = second_order_rk_method(x0, y0, h, n)
# 参数设置
x0 = 0
y0 = 2
h = 0.1
n = int(1 / h)
plt.plot(x_exact, y_exact, label='Exact Solution')
plt.plot(x_second_order_rk, y_second_order_rk, '^-', label='Second Order RK Method')
plt.xlabel('x')
plt.ylabel('y')
plt.title('Numerical Solutions Comparison')
plt.legend()
plt.show()

运行结果

6.四阶龙格库塔方法

6.1四阶龙格库塔方法原理

前面提到的显式欧拉方法仅使用当前点的斜率来估计下一个点的值，精度较低；二阶龙格 - 库塔方法通过计算两个不同点的斜率估计值并加权平均来提高精度。四阶龙格 - 库塔方法进一步拓展了这一思想，在每个时间步长内计算四个不同点的斜率估计值，然后对这些斜率进行加权平均，从而得到更准确的y_{n + 1}的估计值。

四阶龙格 - 库塔方法公式 

给定时间步长h，从t_n到t_{n + 1}=t_n + h，四阶龙格 - 库塔方法的更新公式如下：


 

  1. 计算k_1：k_1 = hf(t_n,y_n)
    • 这里k_1是在当前时间点t_n)和当前函数值y_n处的斜率乘以步长h，它反映了从当前点出发的初始斜率信息。
  2. 计算k_2：k_2 = hf(t_n+\frac{h}{2},y_n+\frac{k_1}{2})
    • k_2是在时间点t_n+\frac{h}{2}和预测的函数值y_n+\frac{k_1}{2}处的斜率乘以步长h。通过在时间区间的中点进行斜率估计，考虑了函数在该区间内的变化趋势。
  3. 计算k_3：k_3 = hf(t_n+\frac{h}{2},y_n+\frac{k_2}{2})
    • k_3同样是在时间点t_n+\frac{h}{2}处计算斜率，但使用了基于k_2更新后的预测函数值y_n+\frac{k_2}{2}，进一步细化了对区间内斜率的估计。
  4. 计算k_4：k_4 = hf(t_n + h,y_n + k_3)
    • k_4是在时间点t_{n + 1}=t_n + h和基于\(k_3\)更新后的预测函数值y_n + k_3处的斜率乘以步长h，反映了在该时间步结束时的斜率信息。
  5. 更新y_{n + 1}：y_{n+1}=y_n+\frac{1}{6}(k_1 + 2k_2+2k_3 + k_4)
    • 最终的y_{n + 1}是通过对k_1、k_2、k_3和k_4进行加权平均得到的，权重分别为\frac{1}{6}、\frac{2}{6}、\frac{2}{6}和\\frac{1}{6}。这种加权方式使得方法具有四阶精度。

6.2 python代码实现

import numpy as np
import matplotlib.pyplot as plt# 定义微分方程
def f(x, y):return -2 * y - 4 * x# 定义精确解
def exact_solution(x):return np.exp(-2 * x) - 2 * x + 1# 四阶龙格-库塔方法
def fourth_order_rk_method(x0, y0, h, n):x = np.zeros(n + 1)y = np.zeros(n + 1)x[0] = x0y[0] = y0for i in range(n):k1 = h * f(x[i], y[i])k2 = h * f(x[i] + h / 2, y[i] + k1 / 2)k3 = h * f(x[i] + h / 2, y[i] + k2 / 2)k4 = h * f(x[i] + h, y[i] + k3)y[i + 1] = y[i] + (k1 + 2 * k2 + 2 * k3 + k4) / 6x[i + 1] = x[i] + hreturn x, y# 计算精确解
x_exact = np.linspace(x0, x0 + 1, n + 1)
y_exact = exact_solution(x_exact)
x_fourth_order_rk, y_fourth_order_rk = fourth_order_rk_method(x0, y0, h, n)
# 参数设置
x0 = 0
y0 = 2
h = 0.1
n = int(1 / h)
plt.plot(x_exact, y_exact, label='Exact Solution')
plt.plot(x_fourth_order_rk, y_fourth_order_rk,'s-', label='Fourth Order RK Method')
plt.xlabel('x')
plt.ylabel('y')
plt.title('Numerical Solutions Comparison')
plt.legend()
plt.show()

运行结果

7.总结

为了直观的比较以上5种方法计算的准确性，将5种计算结果统一绘图比较分析。

import numpy as np
import matplotlib.pyplot as plt# 定义微分方程
def f(x, y):return -2 * y - 4 * x# 定义精确解
def exact_solution(x):return np.exp(-2 * x) - 2 * x + 1# 显式欧拉方法
def euler_method(x0, y0, h, n):x = np.zeros(n + 1)y = np.zeros(n + 1)x[0] = x0y[0] = y0for i in range(n):y[i + 1] = y[i] + h * f(x[i], y[i])x[i + 1] = x[i] + hreturn x, y
# 梯形法
def trapezoid_method(f, x0, y0, h, n):x_values = [x0]y_values = [y0]x = x0y = y0for _ in range(n):def g(y_next):return y_next - y - 0.5 * h * (f(x, y) + f(x + h, y_next))from scipy.optimize import rootsol = root(g, y)y = sol.x[0]x = x + hx_values.append(x)y_values.append(y)return x_values, y_values# 改进的欧拉方法
def improved_euler_method(x0, y0, h, n):x = np.zeros(n + 1)y = np.zeros(n + 1)x[0] = x0y[0] = y0for i in range(n):k1 = f(x[i], y[i])k2 = f(x[i] + h, y[i] + h * k1)y[i + 1] = y[i] + h * (k1 + k2) / 2x[i + 1] = x[i] + hreturn x, y# 二阶龙格-库塔方法
def second_order_rk_method(x0, y0, h, n):x = np.zeros(n + 1)y = np.zeros(n + 1)x[0] = x0y[0] = y0for i in range(n):k1 = h * f(x[i], y[i])k2 = h * f(x[i] + h / 2, y[i] + k1 / 2)y[i + 1] = y[i] + k2x[i + 1] = x[i] + hreturn x, y# 四阶龙格-库塔方法
def fourth_order_rk_method(x0, y0, h, n):x = np.zeros(n + 1)y = np.zeros(n + 1)x[0] = x0y[0] = y0for i in range(n):k1 = h * f(x[i], y[i])k2 = h * f(x[i] + h / 2, y[i] + k1 / 2)k3 = h * f(x[i] + h / 2, y[i] + k2 / 2)k4 = h * f(x[i] + h, y[i] + k3)y[i + 1] = y[i] + (k1 + 2 * k2 + 2 * k3 + k4) / 6x[i + 1] = x[i] + hreturn x, y# 参数设置
x0 = 0
y0 = 2
h = 0.1
n = int(1 / h)# 计算数值解
x_euler, y_euler = euler_method(x0, y0, h, n)
x_improved_euler, y_improved_euler = improved_euler_method(x0, y0, h, n)
x_second_order_rk, y_second_order_rk = second_order_rk_method(x0, y0, h, n)
x_fourth_order_rk, y_fourth_order_rk = fourth_order_rk_method(x0, y0, h,n)
x_trapezoid, y_trapezoid = trapezoid_method(f, x0, y0, h, n)# 计算精确解
x_exact = np.linspace(x0, x0 + 1, n + 1)
y_exact = exact_solution(x_exact)# 打印结果并绘图比较
print("显式欧拉方法结果：", y_euler)
print("改进的欧拉方法结果：", y_improved_euler)
print("二阶龙格-库塔方法结果：", y_second_order_rk)
print("四阶龙格-库塔方法结果：", y_fourth_order_rk)
print("精确解结果：", y_exact)plt.plot(x_exact, y_exact, label='Exact Solution')
plt.plot(x_euler, y_euler, 'o-', label='Euler Method')
plt.plot(x_improved_euler, y_improved_euler, 'x-', label='Improved Euler Method')
plt.plot(x_second_order_rk, y_second_order_rk, '^-', label='Second Order RK Method')
plt.plot(x_fourth_order_rk, y_fourth_order_rk,'s-', label='Fourth Order RK Method')
plt.plot(x_trapezoid, y_trapezoid, '*-', label='Trapezoid Method')
plt.xlabel('x')
plt.ylabel('y')
plt.title('Numerical Solutions Comparison')
plt.legend()
plt.show()

运行结果

对比五种计算方法结果相差较小， 二阶龙格库塔方法计算结果最精确，

http://www.dtcms.com/wzjs/199209.html

相关文章：

  • 网站在哪里搜索做一个公司网站大概要多少钱
  • 网站开发流程简述自媒体发布软件app
  • 91手表网重庆seo优化公司
  • 南县做网站完整html网页代码案例
  • 最好用的网站每日财经要闻
  • 网站内容的建设媒体网站
  • 国内 上市网站建设公司免费的模板网站
  • 做百度网站要多少钱郑州seo团队
  • 小说网站开发中遇到的问题房地产网站模板
  • 网站开发什么语言比较好3小时百度收录新站方法
  • 长沙网站优化效果快速推广
  • 重庆市政府渝快办电商seo搜索优化
  • 威海专业做网站公司点石关键词排名优化软件
  • 网站内容怎么做备份成功营销十大经典案例
  • 嘉兴seo优化优化绿松石什么意思
  • 互联网网站建设哪里好seo排名方案
  • 免费企业网站系统源码宁德市人民政府
  • 网站的登陆注册页面怎么做有品质的网站推广公司
  • 做网站与做网页的区别房地产销售
  • 百度收录网站技巧输入关键词搜索
  • 怎样建外贸公司网站太原seo网络优化招聘网
  • 高端定制手机网站怎么把产品推广到各大平台
  • 建设一个网站要多少钱如何创建网站
  • 杰诚网站建设搜狗搜索引擎优化
  • 后台管理系统界面培训如何优化网站
  • 定制做网站费用20个排版漂亮的网页设计
  • 怀化主要网站磁力屋torrentkitty
  • 做一手机网站需要多少钱软文写作经验
  • 优质网站建设是哪家搜索引擎优化方法有哪几种
  • 域名买好了怎么做网站百度手机seo