[MATLAB] 진동비에 따른 강제 조화 진동 (Harmonic Excitation Vibration) Graph

이 글은 진동학의 강제 조화 진동 (Harmonic Excitation of Damped system)의 그래프를 그린것입니다.

clc; clear;
r = linspace(0,2,500); % r의 범위 0~2
ze = [0.00; 0.20; 0.50; 0.70]; % zeta값 설정
G = 1./sqrt((1-r.^2).^2 + (2.*ze.*r).^2); 
figure(1) % 그래프 
plot(r,G,'Linestyle','-','Linestyle','--','Linestyle',':','Linestyle','-.') % 4개의 그래프 Plotting
xlim([0 2]) % x 그래프 범위 설정
ylim([0 10]) % y 그래프 범위 설정
hold on
grid on
title('Steady-state Response')
xline(1,'LineStyle','--','Label','\omega_n = \omega','LabelOrientation','horizontal') % wn = w 즉 r이 1인 라인 그리기
legend('\zeta = 0.00','\zeta = 0.20','\zeta = 0.50','\zeta = 0.70')
xlabel('r')
ylabel('X\omega_n^2/f_0')