四阶龙格库塔方程解二阶常微分方程组并计算船舶在迎浪下的纵摇埀荡耦合运动方程-附Python代码

最新推荐文章于 2024-03-26 04:54:26 发布

Mezikov

最新推荐文章于 2024-03-26 04:54:26 发布

阅读量2.8k

点赞数 8

分类专栏：四阶龙格库塔计算四阶龙格库塔方程文章标签： matlab

本文链接：https://blog.csdn.net/Mezikov/article/details/107896332

版权

该博客记录了使用四阶龙格库塔方法在Python中解决二阶常微分方程组，以模拟船舶在迎浪时的纵摇和埀荡耦合运动。作者将MATLAB代码移植到Python，并对比了两者2%左右的计算误差。详细代码和结果展示在博客中。

摘要由CSDN通过智能技术生成

0 写在前面

这篇博客是在将我上一篇博客的matlab代码移植到python中，应为后续要开展深度强化学习下的船舶减摇研究，总的来说还是在python上做这项工作比较合适。具体方程的解法和背景不在赘述，见https://blog.csdn.net/Mezikov/article/details/107461970

1 代码

import math
import numpy as np
import matplotlib.pyplot as plt
from math import e
from numpy import reshape

plt.rcParams['font.sans-serif'] = ['SimHei']  # 用来正常显示中文标签
plt.rcParams['axes.unicode_minus'] = False  # 用来正常显示负号

##################   写在前面：这个代码只适合用于已知浪向角的规则波   ########################

ts, ys, zs, ys_, zs_ = [], [], [], [], []
# all below is coefficients of wave
fi = math.pi / 3
B = 1
d = 0.341
delta = 1.76 * 10 ** 3
V = 1.76
g = 9.8
L = 6
lamda = 1.64
rho = 1.025 * 10 ** 3
A_w = 5.448
A_m = 0.346
C_w = A_w / (L * B)
C_p = V / (L * A_m)
C_b = delta / (1000 * L * B * d)
v = 4.44
T = 0.8 * math.sqrt(lamda)
W_0 = 2 * math.pi / T
W_e = W_0 * (1 + v * W_0 * np.cos(fi) / g)
x_b = -0.23
GM_l = (L ** 2 * (5.55 * C_w + 1) ** 3 / 3450) / (C_b * d)
M = 1.6 * 10 ** 3
H_0 = B / (2 * d)
a_zz = 0.8 * H_0 * C_w * M
b_zz = (5.4 * (C_w / C_p) * (H_0 ** 0.5) - 4.7) * delta / ((g * L) ** 0.5)
c_zz = rho * g * A_w
a_z0 =