用matlab编写忽略重力的Jeffcott转子的振幅与时间的关系

以下是用MATLAB编写的忽略重力的Jeffcott转子的振幅与时间的关系的代码：

% Jeffcott转子的振幅与时间的关系
% 忽略重力

clear all;  % 清除所有变量
clc;        % 清空命令窗口

% 转子参数
m = 1;      % 转子质量
k = 1000;   % 转子刚度
c = 0.1;    % 转子阻尼系数
r = 0.5;    % 转子半径
e = 0.2;    % 不平衡质量离心率

% 初始条件
x0 = 0;     % 初始位移
v0 = 1;     % 初始速度

% 时间范围
t = 0:0.01:10;

% 计算振幅与时间的关系
omega_n = sqrt(k/m);                    % 自然频率
zeta = c/(2*sqrt(k*m));                 % 阻尼比
omega_d = omega_n*sqrt(1-zeta^2);        % 阻尼比下的阻尼振荡频率
A = e*r/(m*omega_n^2);                  % 不平衡质量的振幅
phi = atan(2*zeta*omega_n/omega_d);      % 相位角
x = A*exp(-zeta*omega_n*t).*sin(omega_d*t+phi);   % 振幅与时间的关系

% 绘制图形
plot(t,x,'LineWidth',2);                % 绘制振幅与时间的关系图像
title('Jeffcott转子的振幅与时间的关系');  % 添加标题
xlabel('时间');                           % 添加横坐标名称
ylabel('振幅');                           % 添加纵坐标名称
grid on;                                 % 添加网格线

运行以上代码，可以得到Jeffcott转子的振幅与时间的关系图像