function quantisierung1D
%Hauptprogramm zur skalaren Quantisierung
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]; 
  x=[0,1,2,3,4,5,6,7,8,9,10,11,12];
  [N1,N2] = size(f);
elseif input==2
  f=sin(0:0.1:3*pi);
  [N1,N2] = size(f)
  x=0:0.1:(N2-1)/10;
  f=double(f);
end; 

% skalare gleichmŠ§ige Quantisierung
m1=floor(min(f))
m2=ceil(max(f))

A=3 % Anzahl der Quantisierungsintervalle
h=(m2-m1)/A
for j=1:N2
  g(j)=(floor((f(j)-m1)/h)+0.5)*h+m1;
end;



plot(x,f)
hold on
plot(x,f,'ro')
plot(x,g,'bd')
plot(x,g,'b')
hold off
 


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