MATLAB输出一个龙格库塔方法求解转子动力学建模代码

以下是一个使用龙格-库塔方法（RK4）求解转子动力学建模的MATLAB代码：

%% 转子动力学建模
clear; clc;

%% 参数设置
M = 1.5; % 质量
k = 0.8; % 刚度
c = 0.2; % 阻尼
F = 0.5; % 外力
tspan = [0 10]; % 时间范围
y0 = [0 0]; % 初始状态

%% 龙格-库塔方法求解
h = 0.01; % 步长
t = tspan(1):h:tspan(2);
n = length(t);
y = zeros(n, length(y0));
y(1,:) = y0;

for i = 1:n-1
    % RK4方法
    k1 = h * f(t(i), y(i,:));
    k2 = h * f(t(i) + h/2, y(i,:) + k1/2);
    k3 = h * f(t(i) + h/2, y(i,:) + k2/2);
    k4 = h * f(t(i) + h, y(i,:) + k3);
    y(i+1,:) = y(i,:) + 1/6 * (k1 + 2*k2 + 2*k3 + k4);
end

%% 绘图
plot(t, y(:,1), 'r', t, y(:,2), 'b');
title('转子动力学建模');
xlabel('时间');
ylabel('位移/速度');
legend('位移', '速度');

%% 动力学方程
function dydt = f(t, y)
    dydt = zeros(2,1);
    dydt(1) = y(2);
    dydt(2) = 1/M * (F - k*y(1) - c*y(2));
end

代码中，首先定义了转子动力学建模需要用到的参数，包括质量、刚度、阻尼、外力、时间范围和初始状态。然后使用龙格-库塔方法（RK4）求解微分方程，最后绘制了位移和速度随时间变化的图像。注意，在求解微分方程时，需要定义动力学方程f(t, y)。