22.01.2014 Views

TFE4120 Elektromagnetisme Laboppgave 1

TFE4120 Elektromagnetisme Laboppgave 1

TFE4120 Elektromagnetisme Laboppgave 1

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Appendiks<br />

Her nnes det tilleggsinformasjon som kan være nyttig for å gjennomføre<br />

laboppgavene.<br />

A<br />

analyserlyd.m<br />

Denne Matlab-funksjonen håndterer plotting av tids- og frekvensdata fra et<br />

opptak. Fila ligger på fagets hjemmeside, men koden nnes her også i tilfelle<br />

det er noen problemer med nettet.<br />

function analyserlyd(y)<br />

%ANALYSERLYD(Y) Finner og plotter tids- og frekvensdataen til et<br />

%lydopptak Y.<br />

Fs = 44100;<br />

N = length(y);<br />

t = 0 : 1/Fs : (N-1)/Fs;<br />

f = linspace(0, Fs/2, round(N/2));<br />

% Anta samplingrate paa 44100 Hz<br />

% Antall samples<br />

% Tidsvektor<br />

% Frekvensvektor<br />

% Spill av lyden<br />

soundsc(y, Fs)<br />

% Plott tidsdata<br />

subplot(2, 1, 1)<br />

plot(t, y)<br />

xlabel('Tid (s)')<br />

ylabel('y(t)')<br />

zoom on<br />

% Plott frekvensdata med bruk av Hann-tidsvindu<br />

subplot(2, 1, 2);<br />

ywin = y .* ( sin( pi*(0:length(y)-1) / (length(y)-1) ) )'.^2;<br />

Y = abs(fft(ywin)) ./ N;<br />

semilogy(f, Y(1:round(N/2)));<br />

xlabel('Frekvens (Hz)')<br />

ylabel('Y(f)')<br />

xlim([0, 3000])<br />

zoom on<br />

10

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!