function signalsampling
%Hauptprogramm zum Downsampling von Signalen und anschlie§ender Rekonstruktion
%mit der Funktion interpolation2
clf reset
clc
clear all;
close all;
echo off

%input signal
input=2;

if input==1
  f=[1,1,2,3,2,0,-1,2,1,2,4,2,1];  % Beispielsignal f mit 13 Werten aus der Vorlesung
elseif input==2
  f=sin(0:0.1:3*pi);               % Sinussignal f
elseif input==3  
  x=0:pi/10:99.1*pi/10;
  f=(sin(x)).^2.*exp(-x);          % Chirp-Signal f
end; 

% Downsampling
M=3;                               % Sampling-Faktor (kann hier veraendert werden)
N2 = length(f);                    % N2 ist Laenge des Signals f

k=floor((N2-1)/M)
for j=0:k g(j+1)=f(1+M*j); end;    % g ist das reduzierte Signal nach Downsampling 
 


%===Algorithmen===% 

f1=interpolation1(g,f,M,N2,2);       % f wird hier nur zu Illustrationszwecken genutzt!

MSE=(norm(f1-f))^2/N2