clear; close all; m = 1; k = 1; wn = sqrt(k/m); wd = wn; c = 0.1; F = 1; zeta = c/(2*sqrt(m*k)); beta = zeta*wn; t = linspace(0,100,500); ws = [1,0.2,5]; for idx=1:length(ws) w = ws(idx); A = 1/wd*(m*w^2-k+beta*c)*w/((-m*w^2+k^2)^2+c^2*w^2…

clear; close all; zetas = [1,0.5,0.2,0.1,0.01]; w_wn = linspace(0,2,100); Mds=[]; figure(); hold on; for idx=1:length(zetas) zeta = zetas(idx); Md = 1./sqrt((1-w_wn.^2).^2+4*zeta^2.*w_wn.^2); plot(w_wn, Md,'DisplayName',sprintf('zeta=%.2g'…