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

以下是一个使用龙格-库塔方法（RK4）求解转子系统动力学问题的 MATLAB 代码示例：

% 定义系统参数
J = 2.5;   % 转动惯量
b = 0.1;   % 阻尼系数
K = 0.5;   % 刚度系数
t0 = 0;    % 初始时间
tf = 10;   % 最终时间
h = 0.01;  % 步长
w0 = 0;    % 初始角速度
theta0 = 0;% 初始角度

% 定义函数
f = @(t, w, theta) [-b/J*w - K/J*theta; w];

% 初始化
t = t0:h:tf;
w = zeros(size(t));
theta = zeros(size(t));
w(1) = w0;
theta(1) = theta0;

% RK4数值解
for i = 1:length(t)-1
    k1 = f(t(i), w(i), theta(i));
    k2 = f(t(i)+h/2, w(i)+h/2*k1(1), theta(i)+h/2*k1(2));
    k3 = f(t(i)+h/2, w(i)+h/2*k2(1), theta(i)+h/2*k2(2));
    k4 = f(t(i)+h, w(i)+h*k3(1), theta(i)+h*k3(2));
    w(i+1) = w(i) + h/6*(k1(1) + 2*k2(1) + 2*k3(1) + k4(1));
    theta(i+1) = theta(i) + h/6*(k1(2) + 2*k2(2) + 2*k3(2) + k4(2));
end

% 画图
figure();
subplot(2,1,1);
plot(t, w);
title('角速度 vs. 时间');
xlabel('时间');
ylabel('角速度');
subplot(2,1,2);
plot(t, theta);
title('角度 vs. 时间');
xlabel('时间');
ylabel('角度');

在这个例子中，我们使用了龙格-库塔方法（RK4）来求解转子系统的动力学问题。我们首先定义了系统的参数和一个函数 `f`，用来计算下一步的角速度和角度。接下来，我们初始化时间、角速度和角度，然后使用 RK4 方法进行数值解。最后，我们用 MATLAB 的 `plot` 函数来绘制角速度和角度随时间变化的图形。