Mitä rakennetta noudattavalla koodilla saisit tässä harjoituksessa käsitellyt ajan funktiona olevan PPG-signaalin ja tämän spektriestimaatin samaan kuvaikkunaan erillisiin alikuviin ("subplotteihin") alla olevan kuvan mukaisesti? Valitse kaikki oikeat vaihtoehdot.

Palauta tähän palautusboksiin kaikki työssä käyttämäsi MATLAB-koodit .m-päätteisinä tiedostoina. Tarvittaessa

voit palauttaa myös useamman .m-päätteisen tiedoston, jos kirjoitit harjoituksen eri tehtävien koodit eri tiedostoihin.

Mitä rakennetta noudattavalla koodilla saisit tässä harjoituksessa käsitellyt ajan funktiona olevan PPG-signaalin ja tämän spektriestimaatin samaan kuvaikkunaan erillisiin alikuviin ("subplotteihin") alla olevan kuvan mukaisesti? Valitse kaikki oikeat vaihtoehdot.

Piirrä jo aiemmin MATLABiin ladatun PPG-mittausdatan infrapunavalon aallonpituudella mitatun signaalin spektriestimaatti (data-matriisin 3. sarakkeen spektriestimaatti) pwelch-funktiota käyttäen. Vähennä signaalista keskiarvo mean(data(:, 3)) ennen spektrilaskentaa. Muistutuksena: signaalin näytteistystaajuus fs eli näytteenottotaajuus fs on 250 Hz.

Zoomaile tarvittaessa tuloksena saatavaa kuvaajaa ja valitse kaikki ne väittämät, jotka ovat totta:

Kerrataan kuvaajanpiirtoa yksinkertaisen tehtävän avulla

Koska tarkoitus on piirtää signaali ajan funktiona, luodaan aluksi aikavektori. Jotta aikavektorin luonti onnistuu alla olevan ohjeen mukaisesti, pitää data-matriisin olla luotuna eli mittausdatan olla ladattuna MATLABiin aiemmin esillä olleen ohjeistuksen mukaisesti. Tässä tapauksessa aikavektori aika voidaan luoda esimerkiksi seuraavasti:

 fs = 250; % Hz, signaalien näytteistystaajuus
% Varsinainen aikavektorin luonti:
aika = ((1:size(data,1))-1)/fs;  % sekuntia. Tässä siis luodaan ensin 
vektorillinen kokonaislukuja 1:stä
N:ään (N= signaalin pituus näytteissä eli matriisin data korkeus, jonka komento size(data,1)
 palauttaa ), 
vähennetään jokaisesta vektorin alkiosta 1
 (jotta saadaan aikavektori 
alkamaan nollasta) ja lopuksi jaetaan jokainen vektorin alkio 
näytteistystaajuudella fs
. (Yhdessä sekunnissa saadaan yhtä monta
näytettä kuin on näytteistystaajuuden suuruus hertseinä). 

Seuraavaksi piirrä MATLABiin seuraavien vaatimusten mukainen kuvaaja. Tarvittaessa hyödynnä 1. tietokoneharjoituksen osan 2 kuvaajanpiirtotehtävän tuloksia.

• Piirrä data-matriisin 3. sarake ajan funktiona eli kuvaaja, jossa on

  • y-koordinaatteina data-matriisin 3. sarake
  • x-koordinaatteina edellä mainittu aika-vektori, jossa aika-arvot ovat sekunteina.

• Rajaa kuvaa siten, että kuvassa näytetään signaali ainoastaan aikavälillä 20 s ... 26 s . Käytä esimerkiksi xlim-käskyä, ks. doc xlim tai kuvaikkunan zoomaustyökalua.
• X-akselin otsikko on "Aika (s)" ja y-akselin otsikko on "Sormen läpäisseen valon intensiteetti (a.u.)

  "

  • Edellä a.u. = arbitrary unit, 'mielivaltainen yksikkö'.  Tämä siis tarkoittaa, että signaalin voimakkuus on verrannollinen valoanturilla havaitun valon intensiteettiin, mutta kytkentää ei ole kalibroitu eli ei tiedetä, kuinka suurta intensiteettiä kukin signaaliarvo vastaa.  

 Tuloksena saatavassa kuvaajassa on/pitää olla lataamasi data-matriisin 3. sarakkeen arvot y-koordinaatteina ja edellä luomasi aikavektori aika x-koordinaatteina. Lisäksi rajaa kuvaaja siten, että näet siinä aikavälin 20 s ... 26 s  . Rajauksen voi esimerkiksi xlim-komennolla tai kuvaikkunan zoomaus-työkalulla.

Tutki tuloskuvaa ja valitse ne väittämät, jotka ovat totta: