Differences between revisions 5 and 9 (spanning 4 versions)
⇤ ← Revision 5 as of 2005-05-25 13:05:58 →
Size: 2013
Editor: RadekPodgorny
Comment:
← Revision 9 as of 2005-05-26 10:27:16 → ⇥
Size: 2335
Editor: RadekPodgorny
Comment:
Deletions are marked like this. Additions are marked like this.
Line 29: Line 29:
[signal,fs,b] = wavread('Sound.wav');
N = length(signal);
figure; plot(signal);		 [i,fs,b] = wavread('Sound.wav');
N = length(i);
figure; stem((0:length(i)-1)/fs,i,'.');
Line 38: Line 38:
attachment:in.png
Line 41: Line 43:
filtr = dlmread('filtr.dat', '\n');
figure; plot(filtr);		 f = dlmread('filtr.dat', '\n');
figure; stem((0:length(f)-1)/fs,f);
Line 49: Line 51:
attachment:imp.png
Line 52: Line 56:
vystup = conv(signal, filtr);
wavwrite(vystup, fs, b, 'vystup.wav');
figure; plot(vystup);		 o = conv(i, f);
wavwrite(o, fs, b, 'vystup.wav');
figure; stem((0:length(o)-1)/fs,o,'.');
Line 61: Line 65:
attachment:out.png
Line 64: Line 70:
spektrum = fft(signal);
figure; plot(abs(spektrum));		 is = fftshift(fft(i));
figure; stem((-length(is)/2:length(is)/2-1)*fs/length(i),abs(is),'.');
Line 71: Line 77:
attachment:in-s.png
Line 74: Line 82:
ospektrum = fft(vystup);
figure; plot(abs(ospektrum));		 os = fftshift(fft(o));
figure; stem((-length(os)/2:length(os)/2-1)*fs/length(o),abs(os),'.');
Line 81: Line 89:
attachment:out-s.png
Line 84: Line 94:
Ein = sumsqr(signal)
Eout = sumsqr(vystup)		 Ein = sumsqr(i)
Eout = sumsqr(o)

Ein = 5.643853576660156e+003
Eout = 1.793025481292911e+002
Line 87: Line 100:

=== Zaver ===

Je videt, ze se jedna o horni propust.

Semestralni prace z predmetu "Signaly asoustavy"

Zadani

  • Soubor Sound.wav obsahuje 8-mi bitové vzorky audiosignálu, jenž je zpracováván obvodem na úpravu zvuku, jehož impulsová odezva je v souboru filtr.dat.
  • Zobrazte vstupní a výstupní signál a jejich spektra.
  • Určete energii vstupního a výstupního signálu.
  • Z tvaru a průběhu vstupního, výstupního signálu a impulsové odezvy filtru (případně z jejich spekter) usudte na typ zvukové úpravy.

Prilohy

  • attachment:2_P010-3.pdf
  • attachment:filtr.dat

Reseni

Nejdrive provedeme inicializaci prostedi v Matlabu. 

clc;
clear all;
close all;
format long;

Nacteme vstupni signal a zobrazime ho. 

[i,fs,b] = wavread('Sound.wav');
N = length(i);
figure; stem((0:length(i)-1)/fs,i,'.');
title('Vstupni signal');
xlabel('t[s]');
ylabel('s(t)');
grid on;

attachment:in.png

Dale nacteme impulzovou odezvu filtru a zobrazime ji. 

f = dlmread('filtr.dat', '\n');
figure; stem((0:length(f)-1)/fs,f);
title('Impulzova odezva filtru');
xlabel('t[s]');
ylabel('s(t)');
grid on;

attachment:imp.png

Dalsim krokem je aplikace impulzove odezvy filtru na vstupni signal. To provedeme pomoci konvoluce. Vystup opet zobrazime. 

o = conv(i, f);
wavwrite(o, fs, b, 'vystup.wav');
figure; stem((0:length(o)-1)/fs,o,'.');
title('Vystupni signal');
xlabel('t[s]');
ylabel('s(t)');
grid on;

attachment:out.png

Spektrum vstupniho signalu ziskame rychlou fourierovou transformaci. 

is = fftshift(fft(i));
figure; stem((-length(is)/2:length(is)/2-1)*fs/length(i),abs(is),'.');
title('Spektrum vstupniho sugnalu');
xlabel('f[Hz]');
ylabel('S(f)');

attachment:in-s.png

Obdobne i pro vystupni signal. 

os = fftshift(fft(o));
figure; stem((-length(os)/2:length(os)/2-1)*fs/length(o),abs(os),'.');
title('Spektrum vystupniho signalu');
xlabel('f[Hz]');
ylabel('S(f)');

attachment:out-s.png

Nakonec spocitame pozadove hodnoty energii pro vstupni a vystupni signal. 

Ein = sumsqr(i)
Eout = sumsqr(o)

Ein = 5.643853576660156e+003
Eout = 1.793025481292911e+002

Zaver

Je videt, ze se jedna o horni propust.

Prilohy

  • attachment:sempr2.m
  • attachment:vystup.wav