【Matlab】某机器的振动规律为 (1)x=a sin⁡ω t+3a sin⁡3 ωt; (2)y=a sin⁡ω t+3a sin⁡( 1.01ωt) 编写matlab代码 (参数a=2)，分别作位移、速度和加速度随时间变化的曲线图。

在MATLAB中，你可以使用`sin`函数生成这两个振动方程，并通过微分计算出速度和加速度。以下是创建位移(x), 速度(dx/dt)和加速度(d²x/dt²)图形的步骤：

首先，我们需要定义一些变量和常数：

% 定义参数
a = 2; % 振幅
omega = 2*pi; % 角频率

% 时间范围
t = linspace(0, 2*pi/omega, 1000); % 从0到半个周期，采样点足够密集

% 方程1和2的位移
x1 = a*sin(omega*t) + 3*a*sin(3*omega*t);
y1 = a*sin(omega*t) + 3*a*sin(1.01*omega*t);

% 计算速度 (dx/dt)
dx1 = omega*a*cos(omega*t) + 9*omega*a*cos(3*omega*t);
dy1 = omega*a*cos(omega*t) + 3.03*omega*a*cos(1.01*omega*t);

% 计算加速度 (d²x/dt²)
da1 = -omega^2*a*sin(omega*t) - 27*omega^2*a*sin(3*omega*t);
db1 = -omega^2*a*sin(omega*t) - 3.0601*omega^2*a*sin(1.01*omega*t);

然后，绘制三个图形：

% 创建子图
figure;
subplot(3,1,1);
plot(t, x1, 'b', t, y1, 'r'); % 位移图，蓝色是第一个方程，红色是第二个方程
xlabel('Time');
ylabel('Displacement');
title('Position vs Time');

subplot(3,1,2);
plot(t, dx1, 'b', t, dy1, 'r');
xlabel('Time');
ylabel('Velocity');
title('Velocity vs Time');

subplot(3,1,3);
plot(t, da1, 'b', t, db1, 'r');
xlabel('Time');
ylabel('Acceleration');
title('Acceleration vs Time');

最后，运行这段代码，你就得到了位移、速度和加速度随时间变化的曲线图。