tp/informe/plot/script.m (3457B)
1 pkg load control 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', 1.00); 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 set(0, 'DefaultAxesMinorGridAlpha', 0.25); 43 44 45 % num = [6/250 6]; 46 % den = conv([1/500 1], [1/power(1/4000, 2) 1/9000 1]); 47 48 p1=-500 49 complex_roots = roots([1 power(4000, 2)/9000 power(4000, 2)]); 50 p2=complex_roots(1) 51 p3=complex_roots(2) 52 53 k1=-250/(500^2 + ((4000^2)/9000)*-500 + (4000^2)) 54 k2=(p2+250) / ((p2+500)*(p2-(conj(p2)))) 55 k3=conj(k2) 56 57 H1=tf([1 250], [1 500]) 58 H2=power(4000, 2)*tf(1, [1 power(4000, 2)/9000 power(4000, 2)]) 59 60 H = 12*H1*H2; 61 62 figure(1); 63 set(gcf, "paperunits", "inches"); 64 set(gcf, "papersize", [10 4]); 65 set(gcf, "paperposition", [0 0 10 4]); 66 67 % bode(H); 68 % grid on; 69 70 % Cambiar tamaño de las fuentes 71 % set(gca, "fontsize", 12); % ejes actuales 72 % set(findall(gcf,"type","text"), "fontsize", 10); % todos los textos 73 % set(findall(gcf,"type","line"), "linewidth", 1.25); 74 75 % diagrama de bode transferencia original 76 % print("bode_plot.png", "-dpng", "-r1000"); 77 78 set (0, "defaultlinelinewidth", 3.00); 79 set (0, "defaultaxesfontsize", 16); 80 set (0, "defaulttextfontsize", 16); 81 set(0, "DefaultAxesLineWidth", 1.00); 82 83 set (0, "defaultaxesxgrid", "on"); 84 set (0, "defaultaxesygrid", "on"); 85 set (0, "defaultaxeszgrid", "on"); 86 set (0, "defaultaxesxminorgrid", "on"); 87 set (0, "defaultaxesyminorgrid", "on"); 88 set (0, "defaultaxesxminortick", "on"); 89 set (0, "defaultaxesyminortick", "on"); 90 91 set (0, 'DefaultAxesGridAlpha', 0.30); 92 set (0, 'DefaultAxesMinorGridAlpha', 0.25); 93 94 set (gca, "minorgridlinestyle", "--"); % grilla mayor con línea de guiones 95 96 t=0:0.00001:0.02; 97 98 plot(t*1e3, 12.*4000^2.*(-16.3e-6.*exp(-500.*t) 99 + 2.*abs(k2).*exp(real(p2).*t).*cos(imag(p2).*t+arg(k2)))) 100 axis tight; 101 xlim([0 10]) 102 xticks(0:1:10); 103 ylim([-35000 35000]) 104 105 ylabel("AMPLITUD h(t)"); 106 xlabel("TIEMPO t [ms]"); 107 legend("Respuesta al impulso h(t)"); 108 109 graphics_toolkit("qt"); 110 set(gcf, "renderer", "painters"); 111 print("impulse.png", "-dpng", "-r500"); 112 113 k1=(p1+250)/((p1^2 + (4000^2/9000)*p1 + 4000^2) *p1) 114 k2=(p2+250) / ((p2+500)*(p2-conj(p2))*p2) 115 k3=conj(k2) 116 k4=k4=(250)/(500*(4000^2)) 117 118 plot(t*1e3, (6.24*exp(-500.*t))+(12.41*exp(-889.*t).*cos(3900.*t + 2.98)) + (6.01)) 119 axis tight; 120 xlim([0 10]) 121 xticks(0:1:10); 122 ylim([0 20]) 123 124 ylabel("AMPLITUD y(t)"); 125 xlabel("TIEMPO t [ms]"); 126 legend("Respuesta al escalón y(t)"); 127 128 print("heaviside.png", "-dpng", "-r500"); 129 130 131 % pause() 132 % cla