function K5_ma
% K5_ma 
%   ei input/output muuttujia

% Jukka Parviainen, 2003-2005
t = [2001 : 1/365 : 2004];
x = 5*cos(2*pi*1*t) + 0.2*randn(1,length(t)) + 40;

figure(1); clf;    % avaa/aktivoi ikkuna nro 1
plot(t, x); grid on; 
xlabel('aika (vuosia)'), ylabel('GWh'); title('Energian kulutus');

for k = 1:length(x)
  y(k) = 0.5 * (x(k) + x(k-1));
end;
hold on; plot(t, y, 'r');   % piirretään suodatettu punaisella
legend({'Alkuper.','Suodatettu'});

figure(2); clf;
w = [0 : pi/64 : pi];  
H = 0.5 * (1 + exp(-j*w));
plot(w, abs(H));            % magnitudivaste
grid on; title('Amplitudivaste');


