Difference between revisions of "LU-DSP-b"

From DiLab
Jump to: navigation, search
(18.03.2019.)
(13.03.24)
 
(232 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
__NOTOC__
 
__NOTOC__
 
<big>
 
<big>
'''Īsceļi:'''
+
'''Īssaites:'''
 
[[#Kalendārs | Kalendārs]] |  
 
[[#Kalendārs | Kalendārs]] |  
[[#Prakstiskie uzdevumi | Uzdevumi]] |  
+
[[#Uzdevumi | Uzdevumi]] |  
[[#Kursa materiāli | Resursi]]
+
[[#Resursi | Resursi]] |
 
</big>
 
</big>
 +
[[#{{LOCALDAY2}}.{{LOCALMONTH}}.{{#sub:{{LOCALYEAR}}|2}} | Šodiena: {{LOCALDAY2}}.{{LOCALMONTH}}.{{#sub:{{LOCALYEAR}}|2|3}}]] :
 +
[[#{{#expr:{{LOCALDAY2}}+1}}.{{LOCALMONTH}}.{{#sub:{{LOCALYEAR}}|2|3}} | Rītdiena: {{#expr:{{LOCALDAY2}}+1}}.{{LOCALMONTH}}.{{#sub:{{LOCALYEAR}}|2|3}}]]
 +
(ja ir lekcija)
 +
 
=Digitālā (ciparu) signālu apstrāde=
 
=Digitālā (ciparu) signālu apstrāde=
 
LU DF bakalauru kurss.
 
LU DF bakalauru kurss.
Line 17: Line 21:
 
"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"
 
"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 un citi iesniedzamie darbi====
+
====Mājas darbi un citi iesniedzamie darbi====
* Iesniegšanas termiņa laiks ir 30 minūtes pirms nākamās lekcijas sākumam, ja netiek noteikts cits termiņš.  
+
* 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'''.  
 
* 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, pa epastu pasniedzējam.
+
* Darbi iesniedzami elektroniski, E-studijās, izņemot ja prasīts cits iesniegšanas veids.
* E-pasta Subj. ir sekojošs: "'''DSP MD1 Vards Uzvards'''" - piemērs pirmajam mājasdarbam. Vārds uzvārds rakstāmi latīņu burtiem - bez garumzīmēm un tml. (failu sistēmu saderībai).
 
 
* Teksts noformējams PDF faila formātā. Piemēram MD2 fails būtu sekojošs: '''DSP_MD2_Vards_Uzvards.pdf'''
 
* Teksts noformējams PDF faila formātā. Piemēram MD2 fails būtu sekojošs: '''DSP_MD2_Vards_Uzvards.pdf'''
* Iesniedzamie projektu faili 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.
+
* 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.
 +
<!--
 +
* E-pasta Subj. ir sekojošs: "'''DSP MD1 Vards Uzvards'''" - piemērs pirmajam mājasdarbam. Vārds uzvārds rakstāmi latīņu burtiem - bez garumzīmēm un tml.
 +
-->
  
 
====Vērtējums kursā====
 
====Vērtējums kursā====
Line 34: Line 40:
  
  
== Kalendārs ==
+
=Kalendārs=
  
{| border=1 cellspacing=0 cellpadding=4
+
{| border=1 cellspacing=0 cellpadding=4 style = 'border-color: #ddd;'
 
|-
 
|-
! Datums, nedēļa
+
! Datums
! Kursa saturs
+
! Tēma, saturs
 
! Uzdevumi
 
! Uzdevumi
|-
+
|- style='vertical-align: top;'
 
|
 
|
==== 04.02.2019.====
+
==== 07.02.24 ====
 
|
 
|
[https://failiem.lv/u/g8watnvs Lekcijas slaidi ]
+
{| width='100%' style='background-color:#ddd;'
DSP (Digital Signal Processing) ievadlekcija. DSP definīcija un pielietojumi.  
+
|<big>'''Digitālā signālu apstrāde, ievads.'''</big>
 +
|}
 +
 
 +
 
 +
DSP (Digital Signal Processing) ievadlekcija. DSP definīcija un pielietojumi.
 +
 
 +
'''Lasāmviela un resursi''':
 +
* [https://youtu.be/d3gj0t4ddPE Ievadlekcijas video ieraksts]
 +
* [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatā] 1. nodaļa.
  
[http://www.dspguide.com/pdfbook.htm DSPGuide grāmatā] 1. nodaļa.
 
 
|
 
|
|-
+
 
 +
|- style='vertical-align: top;'
 
|
 
|
  
==== 11.02.2019.====
+
==== 14.02.24 ====
 
|
 
|
Statistika, varbūtība, trokšņi un gadījumsignāli.
+
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''Statistika, varbūtība, trokšņi un gadījumsignāli'''</big>
 +
|}
 +
 
  
 
Analogais-ciparu pārveidotājs (ADC jeb ACP modelis) jeb signālu mērīšana (ciparu signāls = signālu diskretizācija + kvantēšana).
 
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.).
 
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.).
  
[http://www.dspguide.com/pdfbook.htm DSPGuide grāmatā] 2-3. nodaļas.
+
'''Lasāmviela un resursi''':
 +
* [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatā] 2-3. nodaļas.
  
| Uzdots [[#MD1]]
 
|-
 
 
|
 
|
 +
* '''Uzdots [[#MD_Rnd | MD_Rnd]]'''
 +
 +
|- style='vertical-align: top;'
 +
|
 +
 +
==== 21.02.24 ====
 +
|
 +
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''Matlab un Octave'''</big>
 +
|}
 +
 +
 +
Svarīgāko Octave komandu apskats (ar kurām var konstruēt signālapstrādes algoritmus). Praktiskais darbs ar Octave.
 +
 +
* [http://www.octave.org OCTAVE] instalējama Ubuntu sistēmā sekojoši:
 +
sudo apt install octave
 +
* Octave darbināšana
 +
** grafiskā redaktora režīmā:  <code>octave --gui</code>
 +
** terminālī: <code>octave</code>  vai  <code>octave-cli</code>
 +
 +
* [https://www.dropbox.com/s/1f1q8m02gvaoox6/Octave_intro.pdf?raw=1 Ievads Octave], slaidi.
 +
 +
|
 +
* ''' '''Termiņš''' [[#MD_Rnd | MD_Rnd]]'''
 +
 +
* '''Uzdots [[#MD_Octave | MD_Octave]]'''
 +
 +
|- style='vertical-align: top;'
 +
|
 +
 +
==== 28.02.24 ====
 +
|
 +
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''Lineāras sistēmas un signālu klasifikācija'''</big>
 +
|}
 +
 +
 +
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''':
 +
* [http://www.dspguide.com/CH5.PDF DSPGuide 5. nodaļa]
  
==== 18.02.2019.====
 
 
|
 
|
Svarīgāko Octave komandu apskats (komandu pilna bāze, ar kurām var konstruēt signālapstrādes algoritmus). Praktiskais darbs ar Octave.
 
  
* M-valoda darbam ar OCTAVE datu apstrādes sistēmām: sudo apt install octave
+
* '''Termiņš''' [[#MD_Octave | MD_Octave]]
* Octave grafiskā redaktora palaišana Ubuntu terminālī: <code>octave</code>  vai  <code>octave-cli</code>
 
  
<!--
+
* Uzdots [[#MD_Sys | MD_Sys]] (nav jāiedniedz)
ADC darbības modelis Octave vidē. [http://failiem.lv/g/fjtyjrs Slaidi].
+
 
-->
+
|- style='vertical-align: top;'
| [[#MD1 | MD1]] termiņš.
 
|-
 
 
|
 
|
  
==== 25.02.2019.====
+
==== 06.03.24 ====
 
|
 
|
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.) [http://failiem.lv/g/fjtyjrs Slaidi].
+
{| width='100%' style='background-color:#ddd;'
Lineāras sistēmas ([http://www.dspguide.com/CH5.PDF DSPGuide 5. nodaļa])
+
|<big>'''Konvolūcija '''</big>
|
+
|}
|-
+
 
 +
 
 +
Konvolūcija un tās īpašības
 +
 
 +
'''Lasāmviela un resursi''':
 +
* [http://www.dspguide.com/CH6.PDF 6.nodaļa]
 +
* [http://www.dspguide.com/CH7.PDF 7.nodaļa]
 +
 
 
|
 
|
 +
'''Termiņš''' [[#MD_Sys | MD_Sys]] (nav jāiesniedz, tiks pārrunāts lekcijā)
  
==== 04.03.2019.====
+
|- style='vertical-align: top;'
 
|
 
|
Konvolūcija ([http://www.dspguide.com/CH6.PDF 6.nodaļa]). Konvolūcijas īpašības ([http://www.dspguide.com/CH7.PDF 7. nodaļa])
+
 
| '''MD-1 termiņš '''
+
==== 13.03.24 ====
|-
 
 
|
 
|
 +
{| width='100%' style='background-color:#ddd;'
 +
||''Attālināti'':  <big>'''KD1'''</big>
 +
|}
 +
 +
Vidus semestra '''kontroldarbs KD-1'''. 
 +
 +
'''Lasāmviela un resursi''':
 +
* 1. - 7. nodaļas no [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatas] un lekcijām.
  
==== 11.03.2019.====
 
 
|
 
|
Vidus semestra '''kontroldarbs KD-1'''. 1. - 7. nodaļas no [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatas] un lekcijām.
+
* '''KD1'''
<br>[http://failiem.lv/u/jdshznk Gadījuma skaitļu masīvs 1x2000]
+
 
| '''  '''
+
|- style='vertical-align: top;'
|-
 
 
|
 
|
  
==== 18.03.2019.====
+
==== 20.03.24 ====
 
|
 
|
Ievads digitālajos filtros.
+
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''Ievads digitālajos filtros, ''Moving average'' filtrs'''</big>
 +
|}
  
 
Digitālie filtri, kustīgā vidējā (moving average) filtrs.
 
Digitālie filtri, kustīgā vidējā (moving average) filtrs.
| Pasludināts [[#MD2 | MD2]]
+
 
|-
+
'''Lasāmviela un resursi''':
 +
* 14.,15. nodaļas no [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatas] un lekcijām.
 +
* [https://qr.ae/pNnblQ What does 50dB sound like?]
 +
 
 +
|
 +
* '''Uzdots [[#MD_Audio | MD_Audio]]'''
 +
 
 +
|- style='vertical-align: top;'
 +
|
 +
 
 +
==== 03.04.24 ====
 +
|
 +
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''Prakse ar Octave, skaņa'''</big>
 +
|}
 +
 
 +
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:
 +
* [https://en.wikipedia.org/wiki/Envelope_%28music%29 Envelope, ADSR]
 +
* [http://digitalsoundandmusic.com/chapters/ch6/ Digital Sound and Music] - Skaņas sintēze - 6.1.7
 +
<!-- * [https://mybudgetstudio.com/what-is-sound-envelope-adsr-explained-with-example Envelope, ADSR] (depreciated) -->
 +
 
 +
* Skaņas faili:
 +
** [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]
 +
 
 +
|
 +
* Uzdots [[#MD_Synth | MD_Synth]]
 +
 
 +
|- style='vertical-align: top;'
 +
|
 +
 
 +
==== 10.04.24 ====
 +
|
 +
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''Audio signāli'''</big>
 +
|}
 +
 
 +
Audio uztvere in apstrāde
 +
 
 +
'''Lasāmviela un resursi''':
 +
* [http://www.dspguide.com/pdfbook.htm DSPGuide 23., 24., 25. nodaļas]
 +
 
 +
|
 +
* ''' '''Termiņš''' [[#MD_Audio | MD_Audio]]'''
 +
 
 +
|- style='vertical-align: top;'
 +
|
 +
 
 +
==== 17.04.24 ====
 
|
 
|
 +
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''Attēlu signāli'''</big>
 +
|}
 +
Attēlu uztvere un apstrāde
 +
 +
'''Lasāmviela un resursi''':
 +
* [http://www.dspguide.com/pdfbook.htm DSPGuide 23., 24., 25. nodaļas]
  
==== 25.03.2019.====
 
 
|
 
|
Praktiskā nodarbība. Izmantojot Octave, konstruēt algoritmu, kas izgriež pauzes runas signālam.
+
* ''' '''Termiņš''' [[#MD_Synth | MD_Synth]]'''
 +
 
 +
|- style='vertical-align: top;'
 
|
 
|
|-
+
==== 24.04.24 ====
 
|
 
|
 +
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''Prakse ar Octave, attēli'''</big>
 +
|}
 +
 +
Speciāla attēlu apstrāde.
 +
 +
Praktiskā nodarbība.
 +
Attēlu apstrāde.
  
==== 01.04.2019.====
 
 
|
 
|
FIR un IIR signālu filtri. Šenona diskretizācijas teorēma. Saskaņotā filtrācija.
+
* Uzdots [[#MD_Image | MD_Image]]
|  '''Termiņš''' [[#MD2 | MD2]]
+
 
|-
+
|- style='vertical-align: top;'
 
|
 
|
  
==== 08.04.2019.====
+
==== 08.05.24 ====
 
|
 
|
Diskrētā Furjē transformācija [http://failiem.lv/u/ovrjrhg DFT-I]. Praktiskā implementācija Octave.
+
{| width='100%' style='background-color:#ddd;'
| '''   '''
+
|<big>'''KD2'''</big>
|-
+
|}
 +
 
 +
Vidus semestra kontroldarbs KD2
 +
* Praktisku uzdevumu risināšana ar Octave.
 +
* Audio un attēlu apstrāde.
 +
 
 
|
 
|
==== 15.04.2020.<br>-22.04.2020.====
+
* '''KD2'''
| ''Lieldienu brīvdienas''
+
 
 +
|- style='vertical-align: top;'
 
|
 
|
|-
+
 
 +
==== 15.05.24 ====
 
|
 
|
==== 29.04.2019.====
+
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''Furjē transformācija. Diskrētā Furjē transformācija '''</big>
 +
|}
 +
 
 +
[[#Furjē_transformācija | Video materiāls]]
 +
 
 +
'''Lasāmviela un resursi''':
 +
* [http://www.dspguide.com/CH8.PDF DSPGuide 8. nodaļa]
 +
* [https://qr.ae/pGpaXf Kā saprast Futjē transformāciju] (Quora)
 +
* [https://youtu.be/h7apO7q16V0 FFT] un idejas uz ko FFT balstās.
 +
 
 
|
 
|
DFT-II
+
* ''' '''Termiņš''' [[#MD_Image | MD_Image]]'''
| '''   '''
+
 
|-
+
|- style='vertical-align: top;'
 
|
 
|
==== 06.05.2019.====
+
 
 +
==== 22.05.24 ====
 
|
 
|
 +
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''FIR un IIR'''</big>
 +
|}
 +
 +
 +
Rekursīvie filtri. FIR un IIR signālu filtri.
 
Signālu filtrācija frekvenču apgabalā.
 
Signālu filtrācija frekvenču apgabalā.
| Pasludināts [[#MD3 | MD3]]
+
 
|-
+
Pratiskais darbs: Izmantojot Octave, izveidot Short-time DFT jeb spektrogrammu runas signālam.
 +
 
 +
 
 +
'''Lasāmviela un resursi''':
 +
* [http://www.dspguide.com/CH17.PDF DSPGuide 17. nodaļa] - Custom filters
 +
* [http://www.dspguide.com/CH19.PDF DSPGuide 19. nodaļa] - Recursive filters
 +
 
 +
* [https://www.youtube.com/watch?v=NvRKtdrssFA Intro to FIR]
 +
 
 
|
 
|
  
==== 13.05.2019.====
+
|- style='vertical-align: top;'
 
|
 
|
Izmantojot Octave, izveidot Short-time DFT jeb spektrogrammu runas signālam.  
+
 
| '''Termiņš''' [[#MD3 | MD3]]
+
==== 29.05.24 ====
|-
 
 
|
 
|
 +
{| width='100%' style='background-color:#ddd;'
 +
|<big>'''Filtru salīdzinājums'''</big>
 +
|}
  
==== 20.05.2019.====
+
'''Lasāmviela un resursi''':
 +
* [http://www.dspguide.com/CH21.PDF DSPGuide 21. nodaļa - filtru salīdzinājums]
 
|
 
|
Vidus semestra '''kontroldarbs KD-2'''. Visas nodaļas no [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatas] un lekcijām.
+
|- style='vertical-align: top;'
<br> [http://failiem.lv/u/pidnxjt Testa signāls]
 
| '''  '''
 
|
 
|-
 
 
|
 
|
==== 27.05.2019.====
+
 
 +
==== 29.06.24 ====
 
|
 
|
Vielas atkārtojums.
+
{| width='100%' style='background-color:#ddd;'
| ''' '''  
+
|<big>'''Kopsavilkums'''</big>
|-
+
|}
 +
'''Lasāmviela un resursi''':
 +
* Visas līdz šim kursā apskatītās nodaļas no [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatas] un lekcijām.'''
 +
 
 +
|
 +
|- style='vertical-align: top;'
 
|
 
|
==== 05.06.2019.====
+
 
 +
==== xx.06.24 ====
 
|
 
|
'''Eksāmens'''
+
{| width='100%' style='background-color:#ddd;'
:Laiks: 10:30
+
|<big>'''Eksāmens'''</big>
:Vieta: 312. auditorija
+
|}
  
* Svarīgi: pēc eksāmena pasniedzējs '''nevarēs izlikt kursā atzīmi''', ja nebūs elektroniski aizpildīta kursa vērtējuma aptauja (kā jau visos LU DF kursos). Tāpēc visiem studentiem lūdzu to savlaicīgi izdarīt.
+
:Laiks: xx:xx
 +
:Vieta: xx <!--E-studijas un tiešsaiste.-->
  
| '''  '''
+
|
|-
+
|- style='vertical-align: top;'
 
|}
 
|}
  
== Signāli ==
+
=Uzdevumi=
* '''S-1''': [http://failiem.lv/u/flhglcy 1. audio wav fails]
+
==MD_Rnd==
* '''S-2''': [http://failiem.lv/u/ikzoaqi 2. audio wav fails]
+
Analizēt nejaušu skaitļu (RND) un to summas histogrammu sekojošiem gadījumiem:
* '''S-3''': [http://failiem.lv/u/pidnxjt 3. audio wav fails]
+
# RND
* '''S-4''': [http://failiem.lv/g/lhnxysu bilde]
+
# RND + RND
* '''S-5''': [https://failiem.lv/u/jzr664nu bilde 2]
+
# 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 [https://www.gnu.org/software/octave/index 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====
 +
# Pierādīt, ka sekojošas sistēmas ir (vai nav) lineāras.  
 +
# Pierādīt, ka sekojošas sistēmas ir (vai nav) laika invariantas.
  
== Prakstiskie uzdevumi ==
+
* y[n] = x[n] + 3
====MD1====
+
* y[n] = x[n] * x[n]
Ģenerēt nejaušus gadījuma skaitļus (vismaz 80000) un izveidot tiem histogrammu (gadījuma procesa empīrisko varbūtību sadalījumu) grafikus gadījumos, kad
+
* y[n] = x[n] + x[n + 1]
# pastāv viens gadījuma skaitļu avots (to Octave simulē rand() funkcija),
 
# gadījuma process ir summa no 2 gadījuma skaitļu avotiem,
 
# gadījuma process ir summa no 4 gadījuma skaitļu avotiem.
 
  
Iesūtīt prasītos histogrammu grafikus gadījuma skaitļu avotiem un secinājumus pdf dokumentā.  
+
Šis uzdevums paredzēts treniņam, un nav obligāti jāiesniedz.
 +
<!--
 +
====Iesniegšana====
 +
Risinājums jāiesniedz e-Studijās kā PDF fails sadaļā MD_Sys.
 +
-->
  
Vēlamā programmēšanas vide: Octave.
+
==MD_Octave==
Iepazīties ar Octave funkciju rand().
+
'''Mērķis''':
 +
Apgūt un praktizēties Octave vides un valodas lietošanā.
  
====MD2====
+
Pieejams e-studijās
Veikt neliela audio faila (dažas sekundes) konvolūciju ar
 
* gausa funkciju,
 
* moving average filtru MAF (kad MAF koeficientu skaits N = 5 un N = 30).
 
Uzrakstīt secinājumus un  iesniegt kā PDF failu.
 
  
====MD3====
+
==MD_Audio==
 +
Izmantojot Octave, konstruēt algoritmu, kas izgriež pauzes runas signālam.
  
====Uzdevumi treniņam====
+
Pieejams e-studijās
* '''Uzd-2''': Izveidot simulatorus: (A) Monētas krišanai pa sazarotu celiņu (tas sakrīt ar šādu [http://www.cut-the-knot.org/Curriculum/Combinatorics/PascalTriangle2.gif varbūtību trīsstūri] ar 5 pakāpēm), kur varbūtība monētai aiziet pa kreiso zaru un pa labo zaru: 1)0.5, 0.5 un 2) 0.3 un 0.7. Secinājumi. (B) Objekta X diskrētai random walk, kur varbūtība iet vienu diskrēto lielumu uz priekšu vai atpakaļ ir 0.5 un 0.5, varbūtība iet vienu diskrēto lielumu pa kreisi, pa labi ir 0.5 un 0.5, attēlot grafiski šādas random walk realizācijas, kad gājienu skaits N=10, 100, 1000. Secinājumi.
+
 
* '''Uzd-3''': Audio faila konvolūcija ar gausa funkciju, moving average filtru MAF (kad MAF koeficientu skaits N = 5 un N = 30), secinājumi.
+
==MD_Synth==
* '''Uzd-4''': DFT matricu formā, kods, ilustratīvs piemērs, secinājumi.
+
Izmantojot Octave, uzrakstīt funkcijas vairākiem skaņas efektiem.
* '''Uzd-5''': Implementēt Kalmana filtru Octave vidē. Ģenerēt testa signālu ar troksni. Filtrēt to ar Kalmana filtru. Mainīt Kalman Gain un veikt secinājumus.
+
 
* '''Uzd-6''': Ģenerēt testa signālu y:
+
Pieejams e-studijās
** N=1000;
+
 
** sigmau=1.8;
+
==MD_Image==
** x=linspace(0,1,N)*.5;
+
Izmantojot Octave, uzrakstīt funkcijas attēlu apstrādei.
** y=x+sigmau*randn(1,N);
+
 
** Izmantojot signālu y, atdalīt no tā troksni, izmantojot (a) Moving Average, (b) Gausa un (c) Kalmana filtrus. Salīdzināt visu filtru rezultātus, izmantojot vidējās kvadrātiskās kļūdas kritēriju (ja tiek izmantots cits kritērijs, tad pamatot tā izvēli). Izdarīt secinājumus.
+
Pieejams e-studijās
* '''Uzd-7''': Izveidot ADC (analog-to-digital signālu pārveidotājs) Octave modeli. Ieejas analogais testa signāls [http://failiem.lv/g/lhnxysu šāds]. Grafiski attēlot ADC ieejas un izejas signālus vienā attēlā  pie parametriem: (a) diskretizācijas frekvence (sampling frequency) Fs = 5 Hz un 4 bitu kvantēšana (ADC kvantēšanas diapazons ir -2 līdz 2 ), (b) Fs = 25 Hz, 10 bitu kvantēšana. Iesūtīt grafikus, aprakstu un secinājumus pdf dokumentā.
+
 
* '''Uzd-8''': 1) Doti punkti x=[1  2  3  4  5  6  7] un y=[0.341627  0.662821  0.534425  0.270544  0.026933  0.937675  0.183618], izrēķināt/interpolēt punktu vērtības starp šiem punktiem, 2) veikt filtrāciju ar dažādām impulsa reakcijām h[n], kad tā ir 1, [-1 1], Gausa funkcija.
+
=Resursi=
* '''Uzd-9''': Izmantojot Octave, konstruēt algoritmu, kas izgriež pauzes runas signālam.
 
* '''Uzd-10''': Izmantojot Octave, konstruēt algoritmu, kas atrod cilvēku sejas attēlā.
 
* '''Uzd-11''': Izmantojot Octave: 1) ielasīt audio failu S-1, 2) pārrēķināt signāla nolases tā, lai tas būtu kvantēts visā signāla diapazonā ar 4 bitiem, 3) saglabāt un atskaņot jauniegūto signālu.
 
* '''Uzd-12''': Izmantojot Octave: 1) izveidot attēlu, kurā redzams jūsu vārda pirmais burts, 2) no iepriekšējā attēla izveidot jaunu, kam fonā būtu cits attēls S-4.
 
* '''Uzd-13''': Izveidot programmu, kas sintezē sinusoīdas signālu ar parametriem: Amplitūda 1.2, Frekvence 5 Hz, Fāze pi/4. Attēlot to uz ekrāna laikā periodā [0 1] (sek.)
 
* '''Uzd-14''': Izmantojot Octave, izveidot Short-time DFT jeb spektrogrammu runas signālam.
 
* '''Uzd-15''': Signālu S-1 pārdiskretizēt ar 0.8 no esošās diskretizācijas frekvences. Abus signālus attēlot grafiski.
 
* '''Uzd-16''': Veikt signālu filtrāciju frekvenču apgabalā. Ar joslas filtru no S-1 izfiltrēt frekvenču apgabalu 0.1 Fs līdz 0.4 Fs (kur Fs ir S-1 diskretizācijas frekvence). Orģinālo signālu un filtrēto signālu attēlot: 1) frekvenču apgabalā, 2) laika apgabalā.
 
* '''Uzd-17''': Implementēt konvolūciju Octave vidē.
 
* '''Uzd-18''': Mainīt iepriekš ģenerētā signāla amplitūdu tā, lai signāla enerģija laika periodā [0 1] (sek.) būtu 1 un tad 1.2 (Normēt signālu pēc tā enerģijas). Attēlot normētos signālus uz ekrāna.
 
* '''Uzd-19''': OCTAVE piemērs:
 
** Fs = 11025; % Definējam diskretizācijas frekvenci
 
** y  = wavrecord(5*Fs, Fs, 'double'); % Ierakstam 5 sekundes skaņas
 
** Fa = 11025; % Definējam ieraksta atskaņošanas frekvenci
 
** sound(y/max(y),Fa); % Atskaņojam ierakstīto
 
** plot(y); % Attēlojam signālu grafiski
 
** Klausīties ierakstu pie vairākām Fa vērtībām. Izdarīt secinājumus!
 
* '''Uzd-20''': Izrēķināt sinusa funkcijas jaudu pie dažādām signāla diskretizācijas frekvencēm fs laika intervālā [0 1 sek.]. Sinusa amplitūda 1.45, frekvence 4 Hz un fāze pi/2. Izdarīt secinājumus.
 
* '''Uzd-21''': Vai pēc diskretizācijas iespējams precīzi noteikt nepārtraukta signāla maksimālo vērtību, izmantojot tikai signāla diskrētās vērtības un papildus tās neapstrādājot? Pamatot atbildi.
 
* '''Uzd-22''': Ģenerēt attēlu (100x120), kas sastāv tikai no melniem vai baltiem pikseļiem. Melnā pikseļa parādīšanās varbūtība 0.7, bet baltā 0.3 Izvadīt attēlu uz datora ekrāna.
 
* '''Uzd-23''': Ielādēt attēlu (jebkādu) OCTAVE vidē. Izvadīt to uz datora ekrāna. Veikt tam 2D konvolūciju ar sekojošām impulsa reakcijām: h1=[ 0 0 0; 0 1 0; 0 0 0 ] un h2=[ 1 1 1; 1 0 1; 1 1 1 ]. Rezultātus izvadīt uz datora ekrāna
 
* '''Uzd-24''': Doti dati: X ass [0  1  2  3  4  5  6  7], Y ass [0    0.6700    3.9237  11.0340  22.9788  40.5927 64.6190  95.7359]. Izrēķināt vērtības pie sekojošām X vērtībām, pieņemot, ka punkti pieder funkcijai, kas ir nepārtraukta: X ass [0    1.3000    2.6000    3.9000    5.2000    6.5000], Y ass = ???
 
* '''Uzd-25''': Izpildīt Octave kodu: C( : , : , 1 )=[1 0]; C( : , : , 2 )=[1 0]; C( : , : , 3 )=[1 0]; image( C ) un veikt attēla C filtrāciju ar Gausa filtru.
 
  
 
== Kursa materiāli ==
 
== Kursa materiāli ==
 
* [http://www.dspguide.com/ DSP Guide] gramata tiešsaistē
 
* [http://www.dspguide.com/ DSP Guide] gramata tiešsaistē
 +
:: "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"
 +
* [https://www.dropbox.com/s/vvw92wtspjqg7ki/refcard-a4.pdf?dl=0 Octave 3.0.0 reference card]
 +
 +
 
* [http://failiem.lv/u/suipzqt Kompleksie skaitļi DSP]
 
* [http://failiem.lv/u/suipzqt Kompleksie skaitļi DSP]
 
* [http://failiem.lv/u/njjfvyx Kalmana filtrs]
 
* [http://failiem.lv/u/njjfvyx Kalmana filtrs]
<!-- * [http://failiem.lv/down.php?i=xigtwaj&n=dft.m DFT] -->
+
* [http://failiem.lv/down.php?i=xigtwaj&n=dft.m DFT]
 +
 
 +
 
 +
==== Furjē transformācija ====
 +
* [https://youtu.be/ds0cmAV-Yek Furjē sērija, jebkurš signāls no sinusoīdām]
 +
* [https://youtu.be/ykNtIbtCR-8 Ievada video]
 +
* [https://youtu.be/spUNpyF58BY Furjē transformācijas vizualizācija]
 +
* [https://youtu.be/1JnayXHhjlg lekcija par Inverso FT (1.daļa)]
 +
* [https://youtu.be/kKu6JDqNma8 lekcija par FT (2.daļa)]
 +
 
 +
* [https://youtu.be/XtypWS8HZco Ievads FFT, The Cooley-Tukey Algorithm]
 +
* [https://www.algorithm-archive.org/contents/cooley_tukey/cooley_tukey.html Cooley-Tukey apraksts] (Algoritmu arhīvs)
 +
 
 +
* [https://www.mathworks.com/help/matlab/ref/fft.html FFT pipemērs Matlab sistēmā]
 +
* [https://upload.wikimedia.org/wikiversity/en/d/dd/Octave.DFT.1.B.FFT.20170706.pdf FFT piemērs Octave sistēmā]
 +
 
 +
* [https://youtu.be/r7GdEWDHtbQ Furje transformācija, vizualizācija ar apļiem] (Smarter every day)
 +
* [https://youtu.be/r6sGWTCMz2k But what is a Fourier series? From heat flow to drawing with circles] (3Blue1Brown)
 +
 
 +
==== FIR & IIR ====
 +
* [https://www.youtube.com/watch?v=NvRKtdrssFA Intro to FIR]
 +
 
 +
==== Citi resursi ====
 +
* [https://ptolemy.berkeley.edu/eecs20/weekly.html Signals and Systems] - kurss Berkeley universitātē
 +
* [https://dspillustrations.com/pages/index.html Learning DSP illustrated]
 +
* [https://brilliant.org/wiki/linear-time-invariant-systems/ Lineāras laika invariantas sistēmas] @Brilliant wiki
 +
* [https://youtu.be/zMkXxI63_Og Laika invariantas sistēmas] - videolekcija ar pierādījumu piemēriem.
 +
 
 +
* [https://setosa.io/ev/image-kernels/ Attēlu filtru kodoli] - vizuāls skaidrojums
 +
 
 +
* [https://mybudgetstudio.com/what-is-sound-envelope-adsr-explained-with-example Envelope, ADSR]
 +
* [http://digitalsoundandmusic.com/chapters/ch1/ Digital sound and music]

Latest revision as of 15:07, 2 March 2024

Īssaites: Kalendārs | Uzdevumi | Resursi | Šodiena: 28.03.24 : Rītdiena: 29.03.24 (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

07.02.24

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


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

Lasāmviela un resursi:

14.02.24

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:

21.02.24

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

28.02.24

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:

06.03.24

Konvolūcija


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

Lasāmviela un resursi:

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

13.03.24

Attālināti: KD1

Vidus semestra kontroldarbs KD-1.

Lasāmviela un resursi:

  • KD1

20.03.24

Ievads digitālajos filtros, Moving average filtrs

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

Lasāmviela un resursi:

03.04.24

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:

10.04.24

Audio signāli

Audio uztvere in apstrāde

Lasāmviela un resursi:

17.04.24

Attēlu signāli

Attēlu uztvere un apstrāde

Lasāmviela un resursi:

24.04.24

Prakse ar Octave, attēli

Speciāla attēlu apstrāde.

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

08.05.24

KD2

Vidus semestra kontroldarbs KD2

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

15.05.24

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

Video materiāls

Lasāmviela un resursi:

22.05.24

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:

29.05.24

Filtru salīdzinājums

Lasāmviela un resursi:

29.06.24

Kopsavilkums

Lasāmviela un resursi:

xx.06.24

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