Exercise L3:
% 1 Bedenk 3 systemen en geef de input en output van deze systemen
Input: Systeem: Output:
Luchttrillingen Microfoon Geluid
Kracht in N Weegschaal Gewicht
Positie lichtknop Licht schakelaar Status van licht
aan/uit
% Opdracht 2a Plot een grafiek y tegen t, zonder de DFT te
berekenen. Schat welke
% frequentie vooral aanwezig is
plot (t,y)
xlabel ('Tijd [s]')
ylabel ('Positie van de pen [cm]')
title ('Pen beweging over tijd')
Er zijn ongeveer 2.7 perioden in 2 seconden tijd. Dus 2.7/2 = 1.4 Hz
% 2 b plot de PSD van het y signaal met pwelch kies een lege matrix
voor de 2e 3e en 4e input. Is de resolutie van de frequentie as hoog
genoeg?
fs = 86;
dt = 1/fs;
figure(2)
[Pxx, F] = pwelch (y-(mean(y)), [], [], [], fs);
plot (F, Pxx);
ylabel ('Positie van de pen [cm^2/Hz]')
xlabel ('Tijd [s]')
title ('PSD')
Nee de plot laat niet duidelijk zien welke frequentie aanwezig is,
dit komt omdat de frequentie te laag is.
% 2 c plot de PSD opnieuw maar deze keer met een hogere resolutie
fs = 171/2;
[Pxx, F] = pwelch (y-(mean(y)),130, [], [], fs);
plot (F, Pxx)
ylabel ('Positie van de pen [cm^2/Hz]')
xlabel ('Tijd [s]')
title ('PSD')
% 2 d Welke frequentie is het meest aanwezig in het plot bij c, komt
het overeen met je schatting bij a
Ja de piek vormt zich rond 1.4 Hz
% 2 e plot de PSD opnieuw maar deze keer zonder het gemiddelde van
het
% signaal af te trekken. Zorg dat de resolutie hoog genoeg is. Wat
is het
% verschil
fs = 171/2;
[Pxx, F] = pwelch (y,130, [], [], fs);
plot (F, Pxx)
ylabel ('Positie van de pen [cm^2/Hz]')
xlabel ('Tijd [s]')
title ('PSD')
Er vormen zich 2 pieken aan het begin, eentje bij 0 Hz
% 1 Bedenk 3 systemen en geef de input en output van deze systemen
Input: Systeem: Output:
Luchttrillingen Microfoon Geluid
Kracht in N Weegschaal Gewicht
Positie lichtknop Licht schakelaar Status van licht
aan/uit
% Opdracht 2a Plot een grafiek y tegen t, zonder de DFT te
berekenen. Schat welke
% frequentie vooral aanwezig is
plot (t,y)
xlabel ('Tijd [s]')
ylabel ('Positie van de pen [cm]')
title ('Pen beweging over tijd')
Er zijn ongeveer 2.7 perioden in 2 seconden tijd. Dus 2.7/2 = 1.4 Hz
% 2 b plot de PSD van het y signaal met pwelch kies een lege matrix
voor de 2e 3e en 4e input. Is de resolutie van de frequentie as hoog
genoeg?
fs = 86;
dt = 1/fs;
figure(2)
[Pxx, F] = pwelch (y-(mean(y)), [], [], [], fs);
plot (F, Pxx);
ylabel ('Positie van de pen [cm^2/Hz]')
xlabel ('Tijd [s]')
title ('PSD')
Nee de plot laat niet duidelijk zien welke frequentie aanwezig is,
dit komt omdat de frequentie te laag is.
% 2 c plot de PSD opnieuw maar deze keer met een hogere resolutie
fs = 171/2;
[Pxx, F] = pwelch (y-(mean(y)),130, [], [], fs);
plot (F, Pxx)
ylabel ('Positie van de pen [cm^2/Hz]')
xlabel ('Tijd [s]')
title ('PSD')
% 2 d Welke frequentie is het meest aanwezig in het plot bij c, komt
het overeen met je schatting bij a
Ja de piek vormt zich rond 1.4 Hz
% 2 e plot de PSD opnieuw maar deze keer zonder het gemiddelde van
het
% signaal af te trekken. Zorg dat de resolutie hoog genoeg is. Wat
is het
% verschil
fs = 171/2;
[Pxx, F] = pwelch (y,130, [], [], fs);
plot (F, Pxx)
ylabel ('Positie van de pen [cm^2/Hz]')
xlabel ('Tijd [s]')
title ('PSD')
Er vormen zich 2 pieken aan het begin, eentje bij 0 Hz