计算浮子和振子在波浪激励力 $f cos \omega t$（$f$ 为波浪激励力振幅，$\omega$ 为波浪频率）作用下前 $40$ 个波浪周期内时间间隔为 $0.2 ~ s$ 的垂荡位移和速度：
(1) 直线阻尼器的阻尼系数为 $10000 ~ N·s/m$；
(2) 直线阻尼器的阻尼系数与浮子和振子的相对速度的绝对值的幂成正比，其中比例系数取 $10000$，幂指数取 $0.5$。
将结果存放在 result1-1.xlsx 和result1-2.xlsx 中。在论文中给出 $10 s、20 s、40 s、60 s、100 s$ 时，浮子与振子的垂荡位移和速度。
# TODO import
import re
import os
import sys
import hmz
import pathlib
import mitosheet
import numpy as np
import pandas as pd
import matlab.engine
import scipy
from scipy.integrate import odeint
import time
import copy
import random
import sympy
from sympy import limit
from sympy import diff
from sympy import integrals
import sklearn
import graphviz
from sklearn import tree
from sklearn.model_selection import cross_val_score
from sklearn.model_selection import train_test_split
from sklearn.metrics import r2_score
from sklearn.metrics import mean_squared_error as MSE
from sklearn.metrics import mean_absolute_error as MAE
from sklearn.metrics import classification_report, roc_auc_score
import sko
from sko.GA import GA
import plotly
import plotly.express as px
import plotly.graph_objects as go
import plotly.figure_factory as ff
plotly.offline.init_notebook_mode()
import cufflinks as cf
cf.set_config_file(
offline=True,
world_readable=True,
theme='pearl', # cf.getThemes()
)
import matplotlib.pyplot as plt
plt.rcParams['font.sans-serif'] = ['SimHei'] # KaiTi
plt.rcParams['axes.unicode_minus'] = False
from IPython.core.interactiveshell import InteractiveShell
InteractiveShell.ast_node_interactivity = 'all'
# InteractiveShell.ast_node_interactivity = 'last'
import cv2 as cv
# import torch
# import torchvision
# import torch.nn as nn
# import torch.nn.functional as F
# import torch.utils.data as Data
# from torch.utils.data import DataLoader
# from torch.utils.data.dataset import Dataset
import pylatex
import latexify
import warnings
warnings.filterwarnings("ignore")