====== Zpracování biologických signálů ====== Pro úlohu použijeme data z [[https://www.physionet.org/physiobank/database/ptbdb/|PTB Databáze]]. ===== Zadání ===== Doporučený postup: Předzpracování: - Odstranění rušení ze sítě (50Hz) - Odstranění driftu (Nízké frekvence / DC složka) - Odstranění artefaktů Zpracování: - Detekce QRS komplexu - Detekce začátku P-vlny - Detekce konce T-vlny Požadovaný výstup: - Prvních 10 rozměřených záznamů - Zobrazení signálu s rozměřenými komplexy - Číselný výpis parametrů pro každý záznam ve tvaru mean(std) pro následující příznaky: * Srdeční frekvence * Počet anomálních (V, S atp) beatů * Šířka QRS komplexu * Šířka QT intervalu * Šířka PQ intervalu * Sklon a plocha ST-úseku * Plocha T vlny ===== Načtení dat ===== fid_record = fopen('RECORDS','r'); fileName = fgetl(fid_record); headerFile = fullfile(path,[fileName '.hea']); ecgFile = fullfile(path,[fileName '.dat']); %------ LOAD HEADER ------------------------------------------- fprintf(strcat('Loading new record #',num2str(counter),': ',fileName)); fprintf('\n'); fid_header = fopen(headerFile,'r'); z = fgetl(fid_header); fclose(fid_header); %------ LOAD DATA -------------------------------------------- fid_ecg = fopen(ecgFile,'r'); ECG = fread(fid_ecg, [12, inf], 'int16'); fclose(fid_ecg); ... ... ===== Odevzdání ===== * __Výsledky:__ viz Požadovaný výstup. * __Forma:__ Skript v Matlabu, který bude načítat data umístěná v ../Data. * Vypracujte do 21.4! Odevzdejte prostřednictvím [[http://cw.felk.cvut.cz/upload|UploadSystemu]].