LU-DSP-b: Difference between revisions

From DiLab
Jump to navigation Jump to search
← Older editNewer edit →
Content deleted Content added
Line 403: Line 403:
Apgūt un praktizēties Octave vides un valodas lietošanā.
Apgūt un praktizēties Octave vides un valodas lietošanā.


====Uzdevums====
Pieejams e-studijās
Pieejams e-studijās


Line 409: Line 408:
Izmantojot Octave, konstruēt algoritmu, kas izgriež pauzes runas signālam.
Izmantojot Octave, konstruēt algoritmu, kas izgriež pauzes runas signālam.


Pieejams e-studijās
Skaņas faili, piemēram:
* [https://www.dropbox.com/s/34ait9wo4b1j1ld/test1.ogg?dl=1 test1.ogg] [https://www.dropbox.com/s/v8wyoidysq94nsh/test1.wav?dl=1 .wav]
* [https://www.dropbox.com/s/n7k971ppe1zfq9o/test2.ogg?dl=1 test2.ogg] [https://www.dropbox.com/s/vpr9dbd7famcjjk/test2.wav?dl=1 .wav]

Darbības plāns:
# Ielasīt audio failu
# Izveidot filtra kodolu lai nofiltrētu augstfrekvenču troksni, "pīķus" signālā, it īpaši paužu laikā
# pielietot šo filtru izmantojot konvolūciju
# uzlikt slieksni filtrētajam audio signālam: pauzes uz 0, pārējais uz 1.
# kombinējot oriģinālo un filtrēto signālus, konstruēt jaunu signālu bez pauzēm.

====Ieteikumi====
* Risinājumus saglabāt failā '''md_audio.m'''
* Rakstiet kodā arī komentārus
* Audio failu var ielasīt no faila un rakstīt failā ar komandām <code>audioread</code> un <code>audiowrite</code>
* Vairāk informāciju par komandām var uzzināt ar <code>help nezināmā-komanda</code>
* Reizēm ir Noderīgi apskatīties uz signāla grafiku ar komandu <code>plot(matrica)</code>. Ja vēlaties apskatīt vairākus grafikus vienlaicīgi, tad pirms plot() komandas dodiet komandu <code>figure();</code>

====Iesniegt====
Risinājumu ('''md_audio.m''' failu) iesniegt estudijās MD_Audio sadaļā.


==MD_Synth==
==MD_Synth==
Izmantojot Octave, uzrakstīt funkcijas vairākiem skaņas efektiem.
Izmantojot Octave, uzrakstīt funkcijas vairākiem skaņas efektiem.


Pieejams e-studijās
Turpmākajos uzdevumos:
* 'sample' ir skaņas matrica (viena kolonna ar mērījumiem - "samples"). Līdzīgi kā lasīts ar audioread, bet viens kanāls.
* 'fs' ir ''sample-rate'', piemēram 44100.
* 'len' ir skaņas ilgums, sekundēs

====Uzdevumi====
Visas funkcijas atgriež matricu ar skaņu (mono). Pašām funkcijām nav jāatskaņo, tās tikai rēķina.
# Izgriež 1 sekundi no skaņas
#* Funkcija <code>onesec(sample, fs)</code> - kas atgriež vektora matricu kas veidota no 'sample'
# Nospēlēt oktāvu: 8 skaņas no Do līdz Do, katru noti 0.25 sekundes ilgi.
#* Funkcija <code>oneoctave(fs)</code> - kas atgriež matricu ar skaņu.
# Atskaņot skaņu kas sastāv no 3 sinusoīdām ar maināmu frekvenci un amplitūdu
#* Funkcija <code>sine3(f1, f2, f3, a1, a2, a3, fs, len)</code> - kur <code>f*</code> ir frekvences un <code>a*</code> ir amplitūdas.
# Piemērot ADSR skaņas avotam, piemēram, sinusoīdai
#* Funkcija <code>adsr(sample, fs, attack, decay, sustain, release)</code> - kas pielieto ADSR "envelope" effektu signālam. "a,d,s,r" parametru summai jābūt 1.
# Izveidot sintētiskās skaņas instrumentu, pēc savas radošās izvēles.
#* Funkcija <code>mysynth(freq, fs)</code> - kas ņem vērā pamata frekvenci freq efektam un atgriež skaņas matricu. Matricas garums atkarīgs no jūsu efekta garuma.

====Resursi====
* [https://mybudgetstudio.com/what-is-sound-envelope-adsr-explained-with-example Envelope, ADSR]
* [http://digitalsoundandmusic.com/chapters/ch6/ Digital sound and music] - sadaļa 6.7.1 - par skaņas sintēzi.
* Ģenerēto skaņu iespējams rakstīt audio failā ar audiowrite funkciju.
* Skaņu iespējams atspēlēt arī no Octave. Sīkāk skat <code>help audioplayer</code>

====Iesniegt====
* '''md_synth.m''' failu ar visām funkcijām
* Estudijās, MD_Synth sadaļā.

====Piemēri====

Piemērs, kurā tiek ģenerēts 1 sekundi garš mazās oktāvas La signāls un saglabāts kā skaņas fails
// Example, La=220Hz, 1 second
fs = 44100;
freq = 220;
samples_per_period = fs/freq;
step = 2*pi / samples_per_period;
limit = fs/step;
t = 0: step: limit;
s = sin(t);
stereo = [s', s']
audiowrite ("la.ogg", stereo, fs);

Alternatīvs risinājums šim pašam uzdevumam
fs = 44100;
freq = 220;
t = linspace(0, 1, fs);
s = sin( t * 2*pi * freq );
stereo = [s', s'];

Veids, kā spēlēt signālu no Octave:
player = audioplayer (stereo, 44100, 8);
play (player);


==MD_Image==
==MD_Image==
Izmantojot Octave, uzrakstīt funkcijas attēlu apstrādei.
Izmantojot Octave, uzrakstīt funkcijas attēlu apstrādei.


Pieejams e-studijās
Turpmākajos uzdevumos:
* 'path' ir taciņa uz skaņas failu, iekļaujot faila vārdu.
* 'img' ir 2d matrica ar (bezkrāsainu) attēlu. Katra vērtība ir vesels skaitlis no 0 līdz 255.
* 'colorimg' ir matrica ar krāsainu RGB attēlu, līdzīgi ka IMG daļa no imread() funkcijas

====Uzdevumi====
# Ielasīt attēlu no faila (piemēram picture.jpg) un atgriezt 2D matricu ar melnbaltu attēlu.
#* Funkcija <code>file2gray(path)</code> - kas atgriež matricu ar attēlu, ko var rādīt ar imshow().
# Pielietot 2D konvolūciju attēlam. Funkcijai jāatgriež apstrādātais attēls.
#* Funkcija <code>filter1(img, kernel)</code> - kas atgriež matricu ar attēlu, ko var rādīt ar imshow().
# Pārbaudīt filtra kodolu. Ko dara sekojošais filtra kodols? Uzrakstīt funkciju kas to pārbauda, un komentāros aprakstīt rezultātu.
#* Funkcija <code>test3(img)</code> - kas atgriež matricu ar attēlu pēc tā apstrādes ar sekojošu filtra kodolu h:
#* <code> h = [4 3 2; 1 0 -1; -2 -3 -4] </code>
# Atrast filtru, kas veic sekojošu uzdevumu: iezīmē attēla kontūras.
#* Funkcija <code>edgifyer(img)</code> - kur <code>img</code> ir frekvences un <code>a*</code> ir amplitūdas.
# Atrast filtru, kas veic sekojošu uzdevumu: izveido attēla "negatīvu", ņemot vērā krāsas. Piemērs [https://www.dropbox.com/s/0yki43k1i03dwpx/im_positive.jpg?raw=1 pozitīvam] un [https://www.dropbox.com/s/jfx0j9aw1g3brwl/im_negative.jpg?raw=1 negatīvam] attēlam. Piemērs ņemts no [https://www.lunapic.com/editor/premade/negative.gif Lunapic redaktora vietnes].
#* Funkcija <code>negimage(colorimg)</code> - kur <code>img</code> ir frekvences un <code>a*</code> ir amplitūdas.

====Resursi====
* Funkcijas attēlu lasīšanai un rakstīšanai failos: imread(), imwrite()
* Funkcija attēlu rādīšanai uz ekrāna: imshow()
* 2D konvolūciju var veikt ar conv2() funkciju. Ņemiet vērā, ka tai ir trešais parametrs, ar kuru var norādīt kā rīkoties robežgadījumos. "same" ir variants kā paturēt rezultējošās matricas izmērus tādus pat kā sākotnējai.
* Ņemiet vērā, ka ielasot attēlu no faila tam var būt '''triju krāsu komponentes''' RGB. Apstrādājot ar conv2() nepieciešams vai no to darīt katrai krāsai atsevišķi, vai melnbaltā režīmā. Te var noderēt funkcija rgb2gray()
* Ņemiet vērā, ka '''imshow()''' funkcija sagaida attēla matricas elementus
** robežās no '''0 līdz 255''' ja tiem ir '''uint8''' tips, vai
** robežās no '''0.0 līdz 1.0''', ja tie ir '''double'''.

====Iesniegt====
* '''md_image.m''' failu ar visām funkcijām
* Estudijās, MD_Image sadaļā.

====Piemēri====

Piemērs, kurā tiek lasīts krāsains un attēlots attēls kā melnbalts
[im, map, alpha] = imread ("lowres1.jpg");
r = im(:,:,1);
g = im(:,:,2);
b = im(:,:,3);
y = 0.299 * r + 0.587 * g + 0.114 * b;
imshow (y)

Piemērs, kā konvertēt 2D matricu b uz tipu, ko sapratis imshow() - pie nosacījuma ka elementu vērtības ir starp 0..255
a = uint8(b)
imshow(a)


=Resursi=
=Resursi=

Revision as of 10:05, 23 February 2023

Īssaites: Kalendārs | Uzdevumi | Resursi | Šodiena: 11.11.26 : Rītdiena: 12.11.26 (ja ir lekcija)

Digitālā (ciparu) signālu apstrāde

LU DF bakalauru kurss.

  • Pasniedzējs Leo Seļāvo [vards.uzvards @gmail.com].

Mērķi

  • Iepazīties ar Digitālās Signālu Apstrādes pielietojumiem un principiem. Praktiskas darbības ar signālu apstrādi. Sistēmu Matlab un Octave iepazīšana un programmēšana tajās.

Motivācija

"Think how electronics has changed the world in the last 50 years. DSP will have the same role over the next 50 years. Learn it or be left behind!" - Steve Smith, autors "DSP Guide"

Mājas darbi un citi iesniedzamie darbi

  • Iesniegšanas termiņa laiks nozīmē ka uz šo datumu darbam jābūt jau gatavam un iesniegtam. Tātad jāiesniedz ir līdz IEPRIEKŠĒJĀS dienas beigām, ja netiek noteikts cits termiņš. Piemēram, ja termiņš ir ceturtdiena, tad darbs jāiesniedz pirms tās, tātad līdz trešdienas beigām.
  • Ja darbs tiek iesniegts ar novēlošanos, rezultāts tiek samazināts par 50%. Ja darbs iesniegts vairāk kā nedēļu pēc termiņa, pasniedzējs darbu var nepieņemt.
  • Darbi iesniedzami elektroniski, E-studijās, izņemot ja prasīts cits iesniegšanas veids.
  • Teksts noformējams PDF faila formātā. Piemēram MD2 fails būtu sekojošs: DSP_MD2_Vards_Uzvards.pdf
  • Ja iesniedzami vairāk par vienu failu, tie arhivējami ZIP formātā, sekojoši: DSP_MD2_Vards_Uzvards.zip Atspiežot šo failu tam jārada direktorija ar tādu pašu nosaukumu, kur atrodas visi faili.

Vērtējums kursā

Kursa galējais vērtējums sastāv no sekojošām komponentēm:

  • 10% - dalība lekcijās un diskusijās
  • 20% - mājas darbi, praktiskie darbi un testi
  • 20+20% - divi kontroldarbi semestra vidū
  • 30% - rakstisks gala eksāmens
  • Obligāta kursa aptaujas anketas aizpildīšana LUIS sistēmā


Kalendārs

Datums Tēma, saturs Uzdevumi

09.02.23
Digitālā signālu apstrāde, ievads.


DSP (Digital Signal Processing) ievadlekcija. DSP definīcija un pielietojumi.

Lasāmviela un resursi:

16.02.23
Statistika, varbūtība, trokšņi un gadījumsignāli


Analogais-ciparu pārveidotājs (ADC jeb ACP modelis) jeb signālu mērīšana (ciparu signāls = signālu diskretizācija + kvantēšana). Signālu klasifikācijas un galvenie parametri (vidējā vērtība, vidējā kvadrātiskā vērtība, SNR, enerģija, jauda, min un max vērtība, dinamiskais diapazons u.c.).

Lasāmviela un resursi:

23.02.23
Matlab un Octave


Svarīgāko Octave komandu apskats (ar kurām var konstruēt signālapstrādes algoritmus). Praktiskais darbs ar Octave.

  • OCTAVE instalējama Ubuntu sistēmā sekojoši:
sudo apt install octave
  • Octave darbināšana
    • grafiskā redaktora režīmā: octave --gui
    • terminālī: octave vai octave-cli

02.03.23
Lineāras sistēmas un signālu klasifikācija


Signālu klasifikācija un galvenie parametri (vidējā vērtība, vidējā kvadrātiskā vērtība, SNR, enerģija, jauda, min un max vērtība, dinamiskais diapazons u.c.).

Lasāmviela un resursi:

09.03.23
Konvolūcija


Konvolūcija un tās īpašības

Lasāmviela un resursi:

Termiņš MD_Sys (nav jāiesniedz, tiks pārrunāts lekcijā)

16.03.23
KD1

Vidus semestra kontroldarbs KD-1.

Lasāmviela un resursi:

  • KD1

23.03.23
Ievads digitālajos filtros, Moving average filtrs

Digitālie filtri, kustīgā vidējā (moving average) filtrs.

Lasāmviela un resursi:

30.03.23
Prakse ar Octave, skaņa

Praktiskā nodarbība. Izmantojot Octave, konstruēt algoritmu, kas apstrādā audio failu.

  • Uzdevumi, Octave vidē
    • Izgriež 1 sekundi no skaņas
    • Nospēlēt oktāvu
    • Atskaņot skaņu kas sastāv no 3 sinusoīdām ar maināmu frekvenci un amplitūdu
    • Piemērot ADSR skaņas avotam, piemēram, sinusoīdai
    • Izveidot savu sintētiskās skaņas instrumentu

Resursi:

06.04.23
Audio signāli

Audio uztvere in apstrāde

Lasāmviela un resursi:


13.04.23
Attēlu signāli

Attēlu uztvere un apstrāde

Lasāmviela un resursi:

20.04.23
Prakse ar Octave, attēli

Speciāla attēlu apstrāde.

Praktiskā nodarbība. Attēlu apstrāde.

27.04.23
KD2

Vidus semestra kontroldarbs KD2

  • Praktisku uzdevumu risināšana ar Octave.
  • Audio un attēlu apstrāde.
  • KD2

11.05.23
Furjē transformācija. Diskrētā Furjē transformācija

Video materiāls

Lasāmviela un resursi:

18.05.23
FIR un IIR


Rekursīvie filtri. FIR un IIR signālu filtri. Signālu filtrācija frekvenču apgabalā.

Pratiskais darbs: Izmantojot Octave, izveidot Short-time DFT jeb spektrogrammu runas signālam.


Lasāmviela un resursi:

25.05.23
Filtru salīdzinājums

Lasāmviela un resursi:

01.06.23
Kopsavilkums

Lasāmviela un resursi:

xx.06.23
Eksāmens
Laiks: xx:xx
Vieta: xx

Uzdevumi

MD_Rnd

Analizēt nejaušu skaitļu (RND) un to summas histogrammu sekojošiem gadījumiem:

  1. RND
  2. RND + RND
  3. RND + RND + RND + RND

Pētījumu veikt vismaz 100000 nejaušiem skaitļiem.

Ar RND saprotam funkciju, kas katrā izsaukumā atgriež nejaušu skaitli no 0 līdz 1.0.

Katram no iepriekš minētajiem gadījumiem uzrakstīt programmu, kas izrēķina un uzzīmē histogrammu.

Programmu ieteicams rakstīt Octave vidē, bet šoreiz atļauts arī kādā citā jums pazīstamā valodā.

E-studijās kā MD1 iesniegt PDF failu ar risinājumu, kurā katram no gadījumiem doti:

  • Programma, kas ģenerē histogrammu
  • Histogramma
  • Secinājumi, kas sanāca, un kāpēc tā.

MD_Sys

Mērķis: Treniņš lineāru un laika invariantu sistēmu atpazīšanā un pierādīšanā.

Uzdevums

  1. Pierādīt, ka sekojošas sistēmas ir (vai nav) lineāras.
  2. Pierādīt, ka sekojošas sistēmas ir (vai nav) laika invariantas.
  • y[n] = x[n] + 3
  • y[n] = x[n] * x[n]
  • y[n] = x[n] + x[n + 1]

Šis uzdevums paredzēts treniņam, un nav obligāti jāiesniedz.

MD_Octave

Mērķis: Apgūt un praktizēties Octave vides un valodas lietošanā.

Pieejams e-studijās

MD_Audio

Izmantojot Octave, konstruēt algoritmu, kas izgriež pauzes runas signālam.

Pieejams e-studijās

MD_Synth

Izmantojot Octave, uzrakstīt funkcijas vairākiem skaņas efektiem.

Pieejams e-studijās

MD_Image

Izmantojot Octave, uzrakstīt funkcijas attēlu apstrādei.

Pieejams e-studijās

Resursi

Kursa materiāli

"The Scientist and Engineer's Guide to Digital Signal Processing, copyright ©1997-1998 by Steven W. Smith. For more information visit the book's website at: www.DSPguide.com"



Furjē transformācija

FIR & IIR

Citi resursi

Retrieved from "https://andromeda.df.lu.lv/wiki/index.php?title=LU-DSP-b&oldid=10502"