请扮演一个船舶运动仿真模拟工程师，帮我写一段模拟船舶横摇运动时历曲线的代码

好的，以下是一个 MATLAB 代码，可以模拟船舶的横摇运动时历曲线：

% 定义船舶参数
m = 5000; % 质量，单位 kg
Iz = 100000; % 惯性矩，单位 kg*m^2
B = 10; % 横向稳定力臂，单位 m
L = 50; % 船长，单位 m

% 定义时间间隔和模拟时间
dt = 0.01; % 时间间隔，单位 s
t = 0:dt:60; % 模拟时间，单位 s

% 定义初值
phi0 = deg2rad(5); % 初始横摇角，单位 rad
p0 = 0; % 初始横摇角速度，单位 rad/s

% 定义控制输入
u = zeros(size(t)); % 没有控制输入

% 模拟船舶横摇运动
phi = zeros(size(t));
p = zeros(size(t));
phi(1) = phi0;
p(1) = p0;
for i = 2:length(t)
    phi(i) = phi(i-1) + p(i-1)*dt;
    p(i) = p(i-1) + (1/Iz)*(B*L*m*9.81*sin(phi(i-1)) - B*p(i-1)*abs(p(i-1)))*dt;
end

% 绘制模拟结果
figure;
plot(t, rad2deg(phi));
xlabel('Time (s)');
ylabel('Roll angle (deg)');
title('Ship roll simulation');

% 计算时历曲线的相关参数
max_phi = max(rad2deg(phi));
max_p = max(rad2deg(p));
period = t(find(phi==max(phi),1,'last'))*2;

% 显示时历曲线的相关参数
fprintf('The maximum roll angle is %0.2f degrees.\n', max_phi);
fprintf('The maximum roll angular velocity is %0.2f degrees/s.\n', max_p);
fprintf('The period of the roll motion is %0.2f s.\n', period);

这个代码基于前面的船舶横摇数学模型，模拟了船舶的横摇运动，并且绘制了横摇角随时间的变化曲线。在模拟完成后，代码还会计算时历曲线的相关参数，包括最大横摇角、最大横摇角速度和横摇周期等。

希望这个代码能够帮助你更好地理解船舶横摇运动的数学模型和仿真方法。