用matlab2022b编写在机械原理中的杆组运动，忽略杆的质量

以下是一个简单的Matlab程序，用于计算机械原理中的杆组运动，忽略杆的质量。注意，这个程序是基于理想情况下的假设，实际情况可能会因为杆的质量、空气阻力等因素而产生误差。

% 定义常数
L1 = 1; % 杆1的长度
L2 = 2; % 杆2的长度
theta1 = pi/4; % 杆1的初始角度
theta2 = pi/4; % 杆2的初始角度
omega1 = 0; % 杆1的初始角速度
omega2 = 0; % 杆2的初始角速度
tspan = [0 10]; % 时间范围

% 定义函数
f = @(t, y) [y(3); y(4); -9.81/L1*sin(y(1)-y(2)); 9.81/L2*sin(y(1)-y(2))];
% y(1) = theta1, y(2) = theta2, y(3) = omega1, y(4) = omega2

% 解微分方程
[t, y] = ode45(f, tspan, [theta1, theta2, omega1, omega2]);

% 绘制图形
figure;
plot(t, y(:,1), 'r-', t, y(:,2), 'b-');
legend('theta1', 'theta2');
xlabel('Time (s)');
ylabel('Angle (rad)');

这个程序使用ode45函数来解决微分方程。f是一个函数句柄，它返回微分方程的右手边。在这个例子中，我们有四个未知函数：theta1、theta2、omega1和omega2。它们的初始值被传递给ode45函数。在解决微分方程之后，我们使用plot函数来绘制杆的角度随时间的变化。