一个四杆曲柄摇杆机构，四个顶点分别为A,B,C,D,其中机架AD=150mm，曲柄AB=15mm，摇杆CD=25mm，BC=149mm，E点在CD延长线上，DE=100mm，当AB杆做匀速圆周运动时，求E点的运动轨迹，请问这个MATLAB代码该怎么写呢

根据机构分析，我们可以得到以下公式：

- 角速度 ω = const
- 角度 θ = ωt
- 曲柄角度 β = θ
- 摇杆角度 γ = arccos((AB^2 + BC^2 - CD^2) / (2AB * BC)) + atan((BC * sin(γ)) / (AB + BC * cos(γ)))
- E点坐标 x = AD + DE * cos(γ + β)
- E点坐标 y = DE * sin(γ + β)

将上述公式转化为MATLAB代码，可以写成以下形式：

% 机构参数
AD = 150; % 机架AD
AB = 15; % 曲柄AB
CD = 25; % 摇杆CD
BC = 149; % 摇杆BC
DE = 100; % E点到CD延长线的距离

% 角速度
omega = 2*pi/10; % 假设匀速圆周运动的周期为10s

% 时间范围
t = 0:0.01:10; % 假设时间步长为0.01s

% 初始化E点坐标数组
E_x = zeros(size(t));
E_y = zeros(size(t));

% 计算E点坐标
for i = 1:length(t)
    theta = omega * t(i);
    beta = theta;
    gamma = acos((AB^2 + BC^2 - CD^2) / (2*AB*BC)) + atan((BC * sin(gamma)) / (AB + BC * cos(gamma)));
    E_x(i) = AD + DE * cos(gamma + beta);
    E_y(i) = DE * sin(gamma + beta);
end

% 绘制E点运动轨迹
plot(E_x, E_y);
xlabel('x');
ylabel('y');
title('E点运动轨迹');

执行上述代码，即可得到E点的运动轨迹图像。