phase angle v/s frequency ratio of a damped Sdof system under harmonic force

9 次查看（过去 30 天）
KIRAN MUKUND
KIRAN MUKUND 2017-9-10
评论： william thomas 2022-10-10
clear all;clc
q=1;
for zeta=[0.05,0.15,.375,1];
p=1;
for r=0:0.001:5
phi(p)=atan((2*zeta*r)/(1-(r^2)));
R(p)=r;
p=p+1;
end
plot(R,phi,'color',rand(1,3),'linewidth',2)
pause(1) %just to see it being draw
hold on;
q=q+1;
end
grid on
% This is my code for plotting phase angle v/s frequency ratio of a damped Sdof system under harmonic force
% But the graph obtaining is wrong
% Error starts from r=1
% The graph i expect is attached with this mail.
% Pls help me in correcting the error...
% This is the graph am Expecting......
  4 个评论
显示 2更早的评论
David Goodmanson
David Goodmanson 2022-9-11
编辑：David Goodmanson 2022-9-11
If in KSSV's code you replace the atan line with
phi(i,:)=atan2((2*zeta(i)*r),(1-(r.^2)));
then the plot looks good. Compared to atan, atan2 keeps track of which quadrant to be in.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Shan Min Hii
Shan Min Hii 2021-5-19
clear all;close all;clc
for zeta=[0:0.2:1]
r=[0:0.1:3];
theta= atand((2*zeta*r)./(1-r.^2));
for j=1:31
if theta(j)<0
theta(j)=theta(j)+180;
end
end
plot(r,theta)
hold on
end
grid on
legend('\zeta=0','\zeta=0.2','\zeta=0.4','\zeta=0.6','\zeta=0.8','\zeta=1',"Location","southeast")
xlabel('w/w0')
ylabel('Phase Angle (\theta)')
title('Phase Angle')
This might help you.

类别

在 Help Center 和 File Exchange 中查找有关 Acoustics, Noise and Vibration 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by