matlab相位图

% 清空工作空间和命令窗口

clear;

clc;

% 模拟生成时间t，位移y(t)和角位移theta(t)

t = linspace(0, 100, 1000); % 时间从0到100，包含1000个点

y = 1e-5 * sin(2 * pi * 0.1 * t) .* exp(-0.01 * t); % 位移y(t) 振荡衰减

theta = 1e-6 * cos(2 * pi * 0.1 * t) .* exp(-0.01 * t); % 角位移theta(t) 振荡衰减

% 绘制时程图（线性位移与角位移随时间）

figure;

subplot(2,1,1);

plot(t, y, 'k', 'LineWidth', 1.5);

xlabel('time (s)');

ylabel('y (m)');

title('a) 线性位移随时间');

grid on;

subplot(2,1,2);

plot(t, theta, 'k', 'LineWidth', 1.5);

xlabel('time (s)');

ylabel('\theta (rad)');

title('b) 角位移随时间');

grid on;

% 绘制相位图（位移-角位移的相图）

figure;

plot(y, theta, 'k', 'LineWidth', 1.5);

xlabel('y (m)');

ylabel('\theta (rad)');

title('相位图：角位移 \theta 与线性位移 y 的关系');

grid on;