用matlab写一段曲柄摇杆机构上连杆的一点的运动轨迹代码

在MATLAB中，你可以使用`animate`函数结合图形绘制和动画功能来模拟曲柄摇杆机构上连杆的一点的运动轨迹。这里是一个简单的例子，假设我们已经有了连杆长度、关节角度以及初始条件。请注意，你需要安装`Graphics`和`Animation`工具箱才能运行此代码。

% 定义常数和变量
L1 = 5; % 曲柄长度
L2 = 3; % 摇杆长度
theta = linspace(0, 2*pi, 100); % 关节角度数组
x = L1*sin(theta); % 连杆1端点x坐标
y = -L1*cos(theta) + L2*sin(theta); % 连杆1端点y坐标
initial_pos = [x(1), y(1)]; % 初始位置

% 创建figure和axis
fig = figure;
ax = axes('Position', [0 0 1 1], 'Box', 'off');

% 动画创建
for i = 1:length(theta)
    plot(ax, x(1:i), y(1:i), 'LineWidth', 2);
    hold on;
    pause(0.01); % 控制动画帧率
end

hold off;
xlabel('X');
ylabel('Y');
title('连杆端点运动轨迹');
annotation('line', initial_pos, [x(end), y(end)], 'Color', 'r', 'LineWidth', 2); % 标记初始和最终位置

% 清理多余信息
delete(findobj(fig, 'Type', 'line'));

% 添加动画显示
anim = VideoWriter('crank_rod_motion.mp4', 'MPEG-4');
open(anim);
for i = 1:length(theta)
    frame = getframe(fig);
    writeVideo(anim, frame);
end
close(anim);