x = 2048*exp(-(1/2)*(((1:300).-150)./40).^2);
d = 20*sin((2*pi*(1:300)/3000)+pi/2);
n = randn(1,300)*sqrt(10);
y = x + d + n;

for i=1:299; z(i) = y(i+1) - y(i); end;
for i=1:298; w(i) = z(i+1) - z(i); end;

for i=1:299; q(i) = x(i+1) - x(i); end;
for i=1:298; r(i) = q(i+1) - q(i); end;

figure
subplot(3,2,1)
plot(x)
title("Signal")
subplot(3,2,2)
plot(y)
title("Signal mit Offset und Rauschen")
subplot(3,2,3)
plot(q)
title("Erste Ableitung von Signal")
subplot(3,2,4)
plot(r)
title("Zweite Ableitung von Signal")
subplot(3,2,5)
plot(z)
title("Erste Ableitung von Signal mit Offset und Rauschen")
subplot(3,2,6)
plot(w)
title("Zweite Ableitung von Signal mit Offset und Rauschen")

format long g

disp("Fläche unter dem Signal: ")
disp(sum(x))
disp("Fläche unter dem Signal mit Offset und Rauschen: ")
disp(sum(y))
