绘制万向节的相对角速度

Question

我试图绘制Universal Joint相对于其输入天使的输出角速度，我称之为phi（在phi = landa_1上方的链接中）。 我正在使用MATLAB进行此操作

  

请注意链接中符号的翻译以及我的代码：

     

beta = Beta，phi = landa_1，Omega_2 = Omega_B，Omega_1 = Omega_A

这是我的代码：

  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

但是！情节不是it should look like（在上面的链接中可用） My current plot

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

输出：