% script faltning.m

close all
clear all

%% Sätt upp en del parameterar

f = 10000;                   % Signalens frekvens
fs = 44100;                 % Sampelfrekvens
Ts = 1/fs;                  % Sampelintervall
k = 0:1:(0.01*fs);             % Diskret tidsvariabel över 0.01 sek
s1 = cos(2*pi*f*k*Ts);       % Signalen skapas

frek=0:22050;
h=[0.5 0.5];                % filtrets impulssvar
H=h(1)+h(2)*exp(i*2*pi*frek*Ts);    % (4.61) i kompendiet
figure(12);plot(frek,abs(H));

figure(11);
subplot(2,1,1);plot(s1,'-');    % Signalen ritas upp
s2= conv(s1,h,'same');
subplot(2,1,2);plot(s2,'-');


%% Längre filter

f = 2*22050/5;                   % Signalens frekvens
fs = 44100;                 % Sampelfrekvens
Ts = 1/fs;                  % Sampelintervall
k = 0:1:(0.01*fs);             % Diskret tidsvariabel över 0.01 sek
s1 = cos(2*pi*f*k*Ts);       % Signalen skapas

h = [1 1 1 1]/4;

H=abs(h(1)*exp(-i*4*pi*frek*Ts)+...
      h(2)*exp(-i*2*pi*frek*Ts)+...
      h(3)*exp(0*pi*frek*Ts)+...
      h(4)*exp(i*2*pi*frek*Ts));
figure(12);plot(frek,abs(H));

figure(11);
subplot(2,1,1);plot(s1,'o-');    % Signalen ritas upp
s2= conv(s1,h,'same');
subplot(2,1,2);plot(s2,'o-');

%%
