%% prog20230524 arreglos y graficos

clear all;
clearvars;
close all;
clc;

%% graficos simples basados en datos

n1 = 0:1:25;              %arreglo fila de numeros del 0 al 25 incremento igual a 1
d1 = randi([0 10],26,1);  %arreglo columna numeros enteros aleatoreos entre 1 y 10

%grafico tipo plot - linea
figure(1)                                           %abre ventana de figura
plot(n1,d1);                                        %grafico de linea
title('este es mi primer grafico');                 %titulo del grafico
xlabel('n1 numeros del 0 al 25 (m)');               %titulo del eje x con sus unidades
ylabel('d1 numneros aleatoreos del 0 al 10 (m/s)'); %titulo del eje y con sus unidades
legend('d1(n1)');                                   %legenda
axis([min(n1),max(n1),min(d1),max(d1)]);            %ejes
grid on;                                            %grilla
box on;                                             %caja

%grafico tipo plot - linea modificada color y marcadores
figure(2)                                           %abre ventana de figura
plot(n1,d1,'r -o');                                 %r = red -o circulo marcador con raya
title('este es mi segundo grafico');                %titulo del grafico
xlabel('n1 numeros del 0 al 25 (m)');               %titulo del eje x con sus unidades
ylabel('d1 numneros aleatoreos del 0 al 10 (m/s)'); %titulo del eje y con sus unidades
legend('d1(n1)');                                   %legenda
axis([min(n1),max(n1),min(d1)-1,max(d1)+1]);            %ejes
grid on;                                            %grilla
box on;                                             %caja

%grafico tipo bar - barra
figure(3)                                           %abre ventana de figura
bar(n1,d1,'b');                                     %grafico de barra b = blue
title('este es mi segundo grafico');                %titulo del grafico
xlabel('n1 numeros del 0 al 25 (m)');               %titulo del eje x con sus unidades
ylabel('d1 numneros aleatoreos del 0 al 10 (m/s)'); %titulo del eje y con sus unidades
legend('d1(n1)');                                   %legenda
axis([min(n1)-1,max(n1)+1,min(d1)-1,max(d1)+1]);    %ejes
grid on;                                            %grilla
box on;                                             %caja

%% graficos comparativos
m1 = 1;
m2 = 4;
k1 = (2*pi)^2*2500;
k2 = (2*pi)^2*4900;
c1 = 5;
c2 = 10;
F1 = 1000;
F2 = 2000;

for n = 1:1001
    f(n) = (n - 1)/10;
    w = 2*pi*f(n);
    X1(n) = F1/sqrt( (k1 - w^2*m1)^2 + (w*c1)^2 );
    X2(n) = F2/sqrt( (k2 - w^2*m2)^2 + (w*c2)^2 );
end;

%plot dibuja graficos en escala lineal
figure(4)
plot(f,X1,'linewidth',1,'r',f,X2,'linewidth',1,'b')
title('respuesta de frecuencia de desplazamiento de dos sistemas vibratorios');
xlabel('frecuencia f (Hz)');
ylabel('desplazamiento X(f) (m)');
legend('X_1(f)','X_2(f)');
axis([min(f),max(f),min(min(X1),min(X2)),1.1*max(max(X1),max(X2))]);
grid on;
box on;

%plot dibuja graficos en escala lineal en eje x y logaritmica en eje y
figure(5)
semilogy(f,X1,'linewidth',1,'k',f,X2,'linewidth',1,'g')
title('respuesta de frecuencia de desplazamiento de dos sistemas vibratorios');
xlabel('frecuencia f (Hz)');
ylabel('desplazamiento X(f) (m)');
legend('X_1(f)','X_2(f)');
axis([min(f),max(f),min(min(X1),min(X2)),1.1*max(max(X1),max(X2))]);
grid on;
box on;

%plot dibuja graficos en escala logaritmica en eje x y lineal en eje y
figure(6)
semilogx(f,X1,'linewidth',1,'c',f,X2,'linewidth',1,'y')
title('respuesta de frecuencia de desplazamiento de dos sistemas vibratorios');
xlabel('frecuencia f (Hz)');
ylabel('desplazamiento X(f) (m)');
legend('X_1(f)','X_2(f)');
axis([min(f),max(f),min(min(X1),min(X2)),1.1*max(max(X1),max(X2))]);
grid on;
box on;

%plot dibuja graficos en escala logaritmica en eje x y logaritmica en eje y
figure(7)
loglog(f,X1,'linewidth',1,'b',f,X2,'linewidth',1,'r')
title('respuesta de frecuencia de desplazamiento de dos sistemas vibratorios');
xlabel('frecuencia f (Hz)');
ylabel('desplazamiento X(f) (m)');
legend('X_1(f)','X_2(f)');
axis([min(f),max(f),min(min(X1),min(X2)),1.1*max(max(X1),max(X2))]);
grid on;
box on;