t = linspace(0, 10, 1000); A0 = 1;%начальная амплитуда(cм) k = 9.9;%коэффициент упругости пружины(cН/cм) ro = 11300/(10^-6);%плотностьсвинца(кг/cм^3) m = 1;%масса свинца(кг) nu = 35;%вязкость глицерина(cН*с/cм^2) V = (m/ro);%объём свинцового шара(cм^3) r = (((3*V)/(4*pi))^(1/3));%радиус свинцового шара(м) b = 6*pi*r*nu;%коэффициент сопротивления глиценрина alpha = b/(2*m);%коэффициент затухания w = sqrt(k/m-alpha^2); A = A0*exp(-alpha*t); x = A.*cos(w*t); figure("Name", "График затухающих колебаний", "NumberTitle", "off"); hold on; xlabel("t, с"); ylabel("x(t),A(t), см"); grid minor; set(gca, 'xtick', t(1):0.5:t(end), 'xlim', [t(1) t(end)]); set(gca, 'ytick', linspace(-A0, A0, 21), 'ylim', [-A0 A0]); %set(get(gca,'ylabel'),'rotation', -1); title("График затухающих колебаний"); plot(t, x); plot(t, A); plot(t, -A); legend("x(t)", "A(t)", "-A(t)", -1);