%% areglos parte 2

clearvars
close all; 
clc;

%% prealocar juan more time
N = 101;
%ceros
w1 = zeros(N,1); %vector columna de ceros
w2 = zeros(1,N); %vector fila de ceros

%unos
w3 = ones(N,1); %vector columna de unos
w4 = ones(1,N); %vector fila de unos

%valores elegidos
w5(1:N,1) = 36.8278; %rellenar con cualquier numero columna
w6(1,1:N) = 9999.1; %rellenar con cualquier numero fila

%numeros consecutivos
for n = 1:N % esto es mas lento
    w7a(n) = n;  %arreglo de numeros del 1 al 101 fila
    w7b(n,1) = n;%arreglo de numeros del 1 al 101 columna
end;
%esto es mas rapido
w8a = 1:N; %arreglo de numeros del 1 al 101 fila
w8b = transpose(w8a);%arreglo de numeros del 1 al 101 columna
w8c(:,1) = 1:N; %arreglo de numeros del 1 al 101 columna
w8d = (1:N)'; %el signo (') comilla simple trasnpone y conjuga

%numeros aleatoreos punto flotante
w9 = rand(N,1); %columna de valores al azar entre 0 y 1 punto flotante
w10 = rand(1,N);%fila de valores al azar entre 0 y 1  punto flotante

%numeros aleatoreos enteros
w11 = randi(10,N,1); %columna de valores al azar entre 1 y 10 entero
w12 = randi(10,1,N);%fila de valores al azar entre 1 y 10  entero

%numeros aleatoreso entre cero y 10
w13 = randi(11,N,1) - 1;
w14 = round(10*rand(N,1));

%datos ordenados
w15 =  1:N;             %de uno a 101;
w16 =  0:N-1;           %de cero a 100
w17 = -N:N;             %de -101 a 101
w18 =  0:(1/100):1;     %de 0 a 1 cada 0.01
w19 =  0:0.001:3;       %de cero a 3 cada 0.01 
w20 =  0:(1/44100):1;   %de cero a 1 cada 1/44100

%% graficar

%graficar una funcion senoidal con ciclo for
f = 2;%frecuencia
for n = 1:N
    t1(n) = (n-1)/(N-1); %de cero a 1 cada 0.01
    f1(n) = sin(2*pi*f*t1(n));
end;

%con prelocacion
t2 = 0:(1/(N-1)):1;
f2 = sin(2*pi*f*t2);

%grafico
figure(1)%numero de figura
plot(t2,f2); %grafica ambos ejes en escala lineal
title('onda sonora senoidal'); %t itulo del grafico
xlabel('tiempo t (s)')%titulo del eje x, simbología y unidades
ylabel('presion p (Pa)');%titulo del eje y, simbología y unidades
legend('sin(2\pift)');%leyenda
grid on;%grilla
box on;%caja y bordes
axis([0 1 -1.5 1.5]);%limites de los ejes [xmin xmax ymin ymax]

figure(2)
plot(t1,f1)%no hacer graficos sin informacion

%datos aleatoreos
lala = rand(N,1);
lolo = rand(N,1);

tata = rand(N,1);
tato = rand(N,1); 

%graficamos dos series de datos con simbolos
figure(3)
plot(lala,lolo,'+',tata,tato,'o');
title('datos aleatoreos');
xlabel('x (m)');
ylabel('y (m)');
legend('y(x)','u(v)')
grid on;
box on;
axis([0 1 0 1]);

%% graficos en escalas lineales y logaritmicas

Nfreq = 1000;
f = (0:Nfreq)/100;  %cuidado he actualizado el valor y 
                    %la estructura de esta variable
                    %es arreglo
w = 2*pi*f;         %es arreglo

%sistema vibratorio masa resorte y amortuguador
m = 1;          %masa
k = (2*pi)^2;   %resorte
c = 0.1;        %amortiguador
F = 10;         %fuerza

%respuesta de frecuencia
H1 = F./sqrt( (k - w.^2*m).^2 + (w*c).^2 );

%H1 se puede calcular usando for, es lo mismo que H2
for n = 1:Nfreq+1
    H2(n) = F/sqrt( (k - w(n)^2*m)^2 + (w(n)*c)^2 );
end;

%grafico escala lineal
figure(4)
plot(f,H1); 
title('respuesta de frecuencia escala lineal-lineal');
xlabel('f (Hz)');
ylabel('H(f)(N/m)');
legend('H(f)');
grid on;
box on;
axis([0 10 0 20]);