function signalsampling
%Hauptprogramm zum Downsampling von Signalen und anschlie§ender Rekonstruktion
%mit der Funktion interpolation
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]; 
elseif input==2
  f=sin(0:0.2:3*pi);
end; 

% Downsampling
M=3;                       %Sampling-Faktor
[N1,N2] = size(f)

k=floor((N2-1)/M)
for j=0:k g(j+1)=f(1+M*j); end;
 


%===Algorithmen===% 


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


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