Difference between revisions of "LU-DSP-b"

From DiLab
Jump to: navigation, search
(Redirecting to LU-DSP-b13)
 
(Removed redirect to LU-DSP-b13)
(Tag: Removed redirect)
Line 1: Line 1:
=Digitālā (ciparu) signālu apstrāde=
#REDIRECT [[LU-DSP-b13]]
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

====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 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ņš.
* 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.
* 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'''
* 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.

====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 ==

{| border=1 cellspacing=0 cellpadding=4
|-
! Nedēļas datumi
! Kursa saturs
|-
| 30.01.2018. - 04.02.2018. || Reģistrācijas nedēļa
|-
| 05.02.2018. - 11.02.2018.
C: 08.02.2018.
| [https://failiem.lv/u/g8watnvs Lekcijas slaidi ]
DSP (Digital Signal Processing) ievadlekcija. DSP definīcija un pielietojumi. [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatā] 1. nodaļa.
<br>M-valoda darbam ar OCTAVE datu apstrādes sistēmām: sudo apt install octave
<br>Octave grafiskā redaktora palaišana Ubuntu terminālī: octave --force-gui
|-
| 12.02.2018. - 18.02.2018.
C: 15.02.2018.
| 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.).
<br>Ievads MATLAB un OCTAVE datu apstrādes sistēmās: [http://failiem.lv/down.php?i=levllia&n=matlab-slides.pdf M-valodas Slaidi 1] [http://failiem.lv/u/desmkgz M-valodas Slaidi 2]
<br> [http://failiem.lv/u/flhglcy 1. wav fails]
<br> [http://failiem.lv/u/ikzoaqi 2. wav fails]
<br> DSP uzdevumi
| ''' '''
|-
| 19.02.2018. - 25.02.2018.
C: 22.02.2018.
| Svarīgāko Octave komandu apskats (komandu pilna bāze, ar kurām var konstruēt signālapstrādes algoritmus). Praktiskais darbs ar Octave.
ADC darbības modelis Octave vidē. [http://failiem.lv/g/fjtyjrs Slaidi].
<br>Statistika, varbūtība, trokšņi un gadījumsignāli. [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatā] 1-3. nodaļa.
| ''' '''
|-
| 26.02.2018. - 04.03.2018.
C: 01.03.2018.
| 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].
Lineāras sistēmas [http://www.dspguide.com/CH5.PDF DSPGuide 5. nodaļa !!!]
<br> Konvolūcijas Octave implementācija.
<br> [http://selavo.lv/~kursi/dsp/2011/Lekcija5_dekompozicija-un-sinteze.pdf Sintēze un dekompozīcija. ]
<br> Pietiekama diskretizācijas frekvence, uzklāšanās.
| '''MD-1 (Uzd-1) '''
|-
| 05.03.2018. - 11.03.2018.
C: 08.03.2018.
| Konvolūcijas Octave implementācija.
| '''MD-1 termiņš '''
|-
| 12.03.2018. - 18.03.2018.
C: 15.03.2018.
| Konvolūcijas Octave implementācija. Šenona diskretizācijas teorēma. Saskaņotā filtrācija.
| ''' '''
|-
| 19.03.2018. - 25.03.2018.
C: 22.03.2017.
| FIR un IIR signālu filtri. Šenona diskretizācijas teorēma. Saskaņotā filtrācija.
| '''MD-2 (Uzd-17 un Uzd-3) '''
|-
| 26.03.2018. - 02.04.2018.
| ''Lieldienu brīvdienas''
| ''' '''
|-
| 03.04.2018. - 08.04.2018.
C: 05.04.2018.
| Praktiskā nodarbība. Izmantojot Octave, konstruēt algoritmu, kas izgriež pauzes runas signālam.
| '''MD-2 termiņš '''
|-
| 09.04.2018. - 15.04.2018.
C: 12.04.2018.
| Vidus semestra '''kontroldarbs KD-1'''. 1. - 7. nodaļas no [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatas] un lekcijām.
<br>[http://failiem.lv/u/jdshznk Gadījuma skaitļu masīvs 1x2000]
| ''' '''
|-
| 16.04.2018. - 22.04.2018.
C: 19.04.2018.
| Diskrētā Furjē transformācija [http://failiem.lv/u/ovrjrhg DFT-I]. Praktiskā implementācija Octave.
| ''' '''
|-
| 23.04.2018. - 29.04.2018.
C: 26.04.2018.
| DFT-II
| ''' '''
|-
| 30.04.2018. - 06.05.2018.
C: 03.05.2018.
| Signālu filtrācija frekvenču apgabalā.
| ''' MD-3 (Uzd-4) '''
|-
| 07.05.2018. - 13.05.2018.
C: 10.05.2018.
| Izmantojot Octave, izveidot Short-time DFT jeb spektrogrammu runas signālam.
| ''' MD-3 termiņš '''
|-
| 14.05.2018. - 20.05.2018.
C: 17.05.2018.
| Vidus semestra '''kontroldarbs KD-2'''. Visas nodaļas no [http://www.dspguide.com/pdfbook.htm DSPGuide grāmatas] un lekcijām.
<br> [http://failiem.lv/u/pidnxjt Testa signāls]
| ''' '''
|
|-
| 21.05.2018. - 27.05.2018.
C: 24.05.2018.
| Vielas atkārtojums.
| ''' '''
|-
| 28.05.2018. - 03.06.2018.
C: 31.05.2018.
| '''Eksāmens'''
* 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.

| ''' '''
|-
|}

== Signāli ==
* '''S-1''': [http://failiem.lv/u/flhglcy 1. audio wav fails]
* '''S-2''': [http://failiem.lv/u/ikzoaqi 2. audio wav fails]
* '''S-3''': [http://failiem.lv/u/pidnxjt 3. audio wav fails]
* '''S-4''': [http://failiem.lv/g/lhnxysu bilde]
* '''S-5''': [https://failiem.lv/u/jzr664nu bilde 2]

== Prakstiskie uzdevumi ==
* '''Uzd-1''': Iepazīties ar Octave funkciju rand() un ģ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 1) pastāv viens gadījuma skaitļu avots (to Octave simulē rand funkcija), 2) gadījuma process ir summa no 2 gadījuma skaitļu avotiem, 3) 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ā. Programmēšanas vide Octave. (Vecā versija: Ģenerēt nejaušus skaitļus (vismaz 80000) un histogrammas sekojošām funkcijām: RND, RND+RND, RND+RND+RND+RND. Izdarīt secinājumus par histogrammām. Iesūtīt grafikus un secinājumus pdf dokumentā. Programmēšanas vide Octave.)
* '''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.
* '''Uzd-4''': DFT matricu formā, kods, ilustratīvs piemērs, secinājumi.
* '''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:
** N=1000;
** sigmau=1.8;
** x=linspace(0,1,N)*.5;
** 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.
* '''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.
* '''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 ==
* [http://www.dspguide.com/ DSP Guide] gramata tiešsaistē
* [http://failiem.lv/u/suipzqt Kompleksie skaitļi DSP]
* [http://failiem.lv/u/njjfvyx Kalmana filtrs]
<!-- * [http://failiem.lv/down.php?i=xigtwaj&n=dft.m DFT] -->

Revision as of 20:54, 26 January 2019

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

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 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ņš.
  • 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.
  • 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
  • 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.

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

Nedēļas datumi Kursa saturs
30.01.2018. - 04.02.2018. Reģistrācijas nedēļa
05.02.2018. - 11.02.2018.

C: 08.02.2018.

Lekcijas slaidi

DSP (Digital Signal Processing) ievadlekcija. DSP definīcija un pielietojumi. DSPGuide grāmatā 1. nodaļa.
M-valoda darbam ar OCTAVE datu apstrādes sistēmām: sudo apt install octave
Octave grafiskā redaktora palaišana Ubuntu terminālī: octave --force-gui

12.02.2018. - 18.02.2018.

C: 15.02.2018.

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.).
Ievads MATLAB un OCTAVE datu apstrādes sistēmās: M-valodas Slaidi 1 M-valodas Slaidi 2
1. wav fails
2. wav fails
DSP uzdevumi

19.02.2018. - 25.02.2018.

C: 22.02.2018.

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

ADC darbības modelis Octave vidē. Slaidi.
Statistika, varbūtība, trokšņi un gadījumsignāli. DSPGuide grāmatā 1-3. nodaļa.

26.02.2018. - 04.03.2018.

C: 01.03.2018.

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.) Slaidi.

Lineāras sistēmas DSPGuide 5. nodaļa !!!
Konvolūcijas Octave implementācija.
Sintēze un dekompozīcija.
Pietiekama diskretizācijas frekvence, uzklāšanās.

MD-1 (Uzd-1)
05.03.2018. - 11.03.2018.

C: 08.03.2018.

Konvolūcijas Octave implementācija. MD-1 termiņš
12.03.2018. - 18.03.2018.

C: 15.03.2018.

Konvolūcijas Octave implementācija. Šenona diskretizācijas teorēma. Saskaņotā filtrācija.
19.03.2018. - 25.03.2018.

C: 22.03.2017.

FIR un IIR signālu filtri. Šenona diskretizācijas teorēma. Saskaņotā filtrācija. MD-2 (Uzd-17 un Uzd-3)
26.03.2018. - 02.04.2018. Lieldienu brīvdienas
03.04.2018. - 08.04.2018.

C: 05.04.2018.

Praktiskā nodarbība. Izmantojot Octave, konstruēt algoritmu, kas izgriež pauzes runas signālam. MD-2 termiņš
09.04.2018. - 15.04.2018.

C: 12.04.2018.

Vidus semestra kontroldarbs KD-1. 1. - 7. nodaļas no DSPGuide grāmatas un lekcijām.


Gadījuma skaitļu masīvs 1x2000

16.04.2018. - 22.04.2018.

C: 19.04.2018.

Diskrētā Furjē transformācija DFT-I. Praktiskā implementācija Octave.
23.04.2018. - 29.04.2018.

C: 26.04.2018.

DFT-II
30.04.2018. - 06.05.2018.

C: 03.05.2018.

Signālu filtrācija frekvenču apgabalā. MD-3 (Uzd-4)
07.05.2018. - 13.05.2018.

C: 10.05.2018.

Izmantojot Octave, izveidot Short-time DFT jeb spektrogrammu runas signālam. MD-3 termiņš
14.05.2018. - 20.05.2018.

C: 17.05.2018.

Vidus semestra kontroldarbs KD-2. Visas nodaļas no DSPGuide grāmatas un lekcijām.


Testa signāls

21.05.2018. - 27.05.2018.

C: 24.05.2018.

Vielas atkārtojums.
28.05.2018. - 03.06.2018.

C: 31.05.2018.

Eksāmens
  • 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.

Signāli

Prakstiskie uzdevumi

  • Uzd-1: Iepazīties ar Octave funkciju rand() un ģ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 1) pastāv viens gadījuma skaitļu avots (to Octave simulē rand funkcija), 2) gadījuma process ir summa no 2 gadījuma skaitļu avotiem, 3) 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ā. Programmēšanas vide Octave. (Vecā versija: Ģenerēt nejaušus skaitļus (vismaz 80000) un histogrammas sekojošām funkcijām: RND, RND+RND, RND+RND+RND+RND. Izdarīt secinājumus par histogrammām. Iesūtīt grafikus un secinājumus pdf dokumentā. Programmēšanas vide Octave.)
  • Uzd-2: Izveidot simulatorus: (A) Monētas krišanai pa sazarotu celiņu (tas sakrīt ar šādu 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.
  • Uzd-4: DFT matricu formā, kods, ilustratīvs piemērs, secinājumi.
  • 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:
    • N=1000;
    • sigmau=1.8;
    • x=linspace(0,1,N)*.5;
    • 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.
  • Uzd-7: Izveidot ADC (analog-to-digital signālu pārveidotājs) Octave modeli. Ieejas analogais testa signāls šā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.
  • 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

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