tp/scripts/plot_responses.m (2593B)
1 clc; clear all; close all; 2 3 %% Figure light theme 4 5 set(0, "defaultAxesColorOrder", [ 6 1.0000 0.0000 0.0000; # red 7 0.0000 0.0000 1.0000; # blue 8 1.0000 0.4980 0.0549; # orange 9 0.1725 0.6275 0.1725; # green 10 0.5804 0.4039 0.7412; # purple 11 0.5490 0.3373 0.2941; # brown 12 0.8902 0.4667 0.7608; # pink 13 0.4980 0.4980 0.4980; # gray 14 ]); 15 16 set(0, "DefaultFigureColor", [1 1 1]); 17 set(0, "DefaultAxesColor", [1 1 1]); 18 set(0, "DefaultAxesXColor", [0 0 0]); 19 set(0, "DefaultAxesYColor", [0 0 0]); 20 set(0, "DefaultTextColor", [0 0 0]); 21 22 set(0, "DefaultAxesGridColor", [0 0 0]); 23 set(0, 'DefaultAxesGridAlpha', 1.00); 24 25 set(0, "DefaultAxesMinorGridColor", [0 0 0]); 26 set(0, 'DefaultAxesMinorGridAlpha', 0.20); 27 28 set(0, "DefaultLineLinewidth", 3.00); 29 set(0, "DefaultAxesFontSize", 16); 30 set(0, "DefaultTextFontSize", 16); 31 set(0, "DefaultAxesLineWidth", 1.00); 32 33 set(0, "DefaultAxesXGrid", "on"); 34 set(0, "DefaultAxesYGrid", "on"); 35 set(0, "DefaultAxesZGrid", "on"); 36 set(0, "DefaultAxesXMinorGrid", "on"); 37 set(0, "DefaultAxesYMinorGrid", "on"); 38 set(0, "DefaultAxesXMinorTick", "on"); 39 set(0, "DefaultAxesYMinorTick", "on"); 40 41 set(0, 'DefaultAxesGridAlpha', 0.50); 42 43 omega = 2*pi*250; 44 45 t = linspace(0.00, 0.10, 4000); 46 47 % y = (-0.49 .* exp(-500 .* t) + ... 48 % 12.4 .* exp(-889 .* t) .* cos(3900 .* t - 1.24) + ... 49 % 7.37 .* cos(6283 .* t + 2.05)) .* (t >= 0); 50 51 y = (12*4000^2*1e-9) .* ... 52 ( -9.41 .* exp(-500 .* t) + ... 53 29.39 .* exp(-889 .* t) .* cos(3900 .* t + 1.108) + ... 54 69.83 .* cos(omega .* t - 1.62)) .* (t >= 0); 55 56 y_sine = sin(omega * t); 57 58 figure(1); 59 grid on; 60 set(gcf, "paperunits", "inches"); 61 set(gcf, "papersize", [10 4]); 62 set(gcf, "paperposition", [0 0 10 4]); 63 64 % bode(H); 65 % grid on; 66 67 % Cambiar tamaño de las fuentes 68 % set(gca, "fontsize", 12); % ejes actuales 69 % set(findall(gcf,"type","text"), "fontsize", 10); % todos los textos 70 % set(findall(gcf,"type","line"), "linewidth", 1.25); 71 72 % diagrama de bode transferencia original 73 % print("bode_plot.png", "-dpng", "-r1000"); 74 75 set (gca, "minorgridlinestyle", "--"); % grilla mayor con línea de guiones 76 77 plot(t*1e3, y_sine, "b", "LineWidth", 3); 78 hold on; 79 plot(t*1e3, y, "r", "LineWidth", 3); 80 81 % title("Time Response"); 82 xlim([0 16]) 83 xticks(0:1:16); 84 ylim([-15 15]) 85 86 xlabel("TIEMPO t [ms]"); 87 ylabel("AMPLITUD y(t)"); 88 legend("Señal senoidal 250Hz", "Respuesta a señal senoidal ", 89 "Location", "northeast"); 90 91 graphics_toolkit("qt"); 92 % set(gcf, "renderer", "painters"); 93 set(gcf, "renderer", "opengl"); 94 print("sine_response_omega1.png", "-dpng", "-r500");