Построение относительной угловой скорости универсального шарнира

Question

Я пытаюсь построить выходную угловую скорость универсального шарнира относительно входного ангела, который я называю фи (в ссылке выше phi = landa_1). Я использую MATLAB, чтобы сделать это

обратите внимание на перевод обозначений в ссылке и мой код:

бета = бета, фи = landa_1, омега_2 = омега_B, омега_1 = омега_А

вот мой код:

  clear all, close all, clc
  phi=0:360; % one rotation of the input shaft
  Beta=60;
  Omega_A=1;
  Omega_B=(Omega_A*cos(Beta))./(1-((sin(Beta))^2)*((cos(phi)).^2))
  plot(phi,Omega_B,'LineWidth',2), grid on

НО! сюжет не тот, что должен выглядеть (который доступен по ссылке выше) Мой текущий сюжет

Answer 1

Вам необходимо преобразовать градусы в радианы, чтобы использовать sin и cos. Итак

clear all, close all, clc
phi=0:360; % one rotation of the input shaft
Beta=60;
Omega_A=1;
Omega_B=(Omega_A*cos(Beta/180*pi))./(1-((sin(Beta/180*pi))^2)*((cos(phi/180*pi)).^2))
plot(phi,Omega_B,'LineWidth',2), grid on

Выход: