%%%%% Práctica 2


close all,
% Ejercicio 1
load ecg_red.mat
figure, plot (ecg)
Fm=1000; %Frecuencia de muestreo
fe=50; %Eliminar frecuencia de 50 Hz
N=Fm/fe;
k=1:N-1;
f_eliminadas=Fm/N*k;
%espectro(ecg, Fm)
B=1/N*ones (1,N);
A=1;
y=filter (B,A,ecg);
hold on, plot(y, 'r'), legend ('ECG', 'Filtrada')

%Diagrama de ceros y polos
figure, zplane(B,A)

%Respuesta en frecuencia en modulo y fase
figure, freqz (B,A)

espectro(y, Fm)

%%%%%%% Ejercicio 2 %%%%%%%%%

close all, clear
x=load('milk.dat');
figure, plot(x)

%Promediador anual
N=12;
B=1/N*ones(1,N);