用matlab2022b写一份展现机械原理中不同杆组连接运动的代码

以下是一个简单的 MATLAB 代码，用于模拟不同杆组连接运动的机构。该机构由两个旋转连接杆和一个滑动杆组成。

%定义常数
L1 = 1; %旋转杆1长度
L2 = 2; %旋转杆2长度
L3 = 3; %滑动杆长度
theta1 = pi/4; %旋转杆1初始角度
theta2 = pi/3; %旋转杆2初始角度
t = linspace(0, 10, 1000); %时间范围

%计算各点坐标
x1 = L1*cos(theta1);
y1 = L1*sin(theta1);
x2 = x1 + L2*cos(theta2);
y2 = y1 + L2*sin(theta2);
d = sqrt((x2-x1)^2 + (y2-y1)^2);
alpha = acos((L1^2 + L2^2 - d^2)/(2*L1*L2));
theta3 = pi - alpha;
x3 = x1 + L3*cos(theta3);
y3 = y1 + L3*sin(theta3);

%绘图
figure;
for i = 1:length(t)
    %计算各点坐标
    theta1 = pi/4*cos(t(i));
    theta2 = pi/3*sin(2*t(i));
    x1 = L1*cos(theta1);
    y1 = L1*sin(theta1);
    x2 = x1 + L2*cos(theta2);
    y2 = y1 + L2*sin(theta2);
    d = sqrt((x2-x1)^2 + (y2-y1)^2);
    alpha = acos((L1^2 + L2^2 - d^2)/(2*L1*L2));
    theta3 = pi - alpha;
    x3 = x1 + L3*cos(theta3);
    y3 = y1 + L3*sin(theta3);
    
    %绘制机构
    plot([0 x1],[0 y1],'r-o');
    hold on;
    plot([x1 x2],[y1 y2],'b-o');
    plot([x1 x3],[y1 y3],'g-o');
    axis([-5 5 -5 5]);
    hold off;
    pause(0.01);
end

该代码使用了 MATLAB 中的绘图函数 `plot` 来绘制机构，并通过循环计算各点坐标来模拟机构运动。运行代码后，可以看到机构会按照预设的运动规律进行运动。

需要注意的是，该代码并没有考虑机构的动力学特性，仅仅是模拟了机构的几何运动。如果需要更加精确的模拟，需要考虑机构的质量、惯性、摩擦等因素，可以通过 MATLAB 中的动力学模拟工具箱来实现。