Labor7 fertig

%Marvin Meier
%Aufgabe 7 : Rechteck programmieren und transformieren in den
%Frequentbereich

%Abtastrate
fa = 11025;

%Länge des Signals
duration = 1.0;

%50Hz im Frequenzbereich
hzAngabe = 50;

%Definition des start und stop
start = 0.4;
stop = 0.6;

%Komplettes Signal wird auf 0 gesetzt
rectSignal = zeros(1,fa*duration);

%Wird der Abtastrate angepasst
stopFa = floor(stop*fa);
startFa = floor(start*fa);

%Zeitvektor
tVec = 0:1/fa:duration-1/fa;

%Stop und Start
%einsenVektor = ones(1,stopFa-startFa);

rectSignal(startFa:stopFa) = 1;

%In den Frequenzbereich umgewandelt
[specVec, phaVec] = funcRdft(rectSignal, length(rectSignal));

xFVec = 0:1:hzAngabe-1;
%Abbildung 1 Zeitbereich
figure(1);
subplot(3,1,1);
plot(tVec,rectSignal);
title('Zeitbereich');
xlabel('Zeit in Sekunden');
ylabel('Amplitude');

%Abbildung 1 Frequenzbereich
subplot(3,1,2);
plot(specVec);
title('Frequenzbereich');
xlabel('Hz bis fa/2');
ylabel('Amplitude');

%Offset nach oben schieben
%Abbildung 1 Frequenzbereich 25 Hz
subplot(3,1,3);
plot(xFVec,specVec(1:hzAngabe));
title('Frequenzbereich');
xlabel('Hz bis 50 obere Frequenz');
ylabel('Amplitude');