МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ "ЛЬВІВСЬКА ПОЛІТЕХНІКА" / Звіт до лабораторної роботи № 1 з курсу: "Системи запису та відтворення інформації" Львів – 2009р. Мета роботи:
Отримати базові навики роботи в середовищі MatLab, вивчити основні можливості MatLab по обробці аудіоінформації. Повний текст завдання: Ознайомитись з теоретичним матеріалом. Завантажити wav-файл вказаний в завданні відповідно до свого варіанту. Вивести основну інформацію про файл: кількість вибірок, кількість каналів, частота дискретизації, кількість біт на вибірку. Прослухати аудіофайл. Сформувати заваду у вигляді синусоїдального сигналу з вказаними параметрами і накласти її на аудіосигнал. Прослухати зашумлений аудіосигнал. Розрахувати коефіцієнти нерекурсивного фільтра вказаного порядку N для усунення завади. Побудувати АЧХ розрахованого фільтра. Відфільтрувати сигнал. Прослухати очищений аудіосигнал. Зберегти результат в wav-файлі. Розрахувати коефіцієнти рекурсивного фільтра вказаного порядку N та типу для усунення завади (пульсації в зоні пропускання прийняти рівними 0,1 дБ, мінімально допустиме загасання в смузі затримки 60 дБ). Побудувати АЧХ розрахованого фільтра. Відфільтрувати сигнал. Прослухати очищений аудіосигнал. Зберегти результат в wav-файлі. Варіант  Назва файлу  Частота завади, Гц  Амплітуда завади, В  Порядок нерекурсивного фільтру N  Тип рекурсивного фільтру  Порядок рекурсивного фільтру N  4 Lab_1_12.wav 4000 5 230 Чебишева 1 7   Лістинг програми: clear all; [y, Fs, bits] = wavread('G:\С-ми запису і відтворення інф\Lab_1_11.wav'); f = 4000; Amp = 5; porNerTupy = 230; porRecTupy = 7; wavplay(y,Fs); t = 0 : length(y)-1; t = t / Fs; Noise = Amp * sin( 2*pi*f*t ); Ywym=y+Noise'; wavplay(Ywym,Fs); figure(1); strips(y(1:length(y),:),10000); figure(2); strips(Ywym(1:length(Ywym),:),10000); %розрахунок нерекурсивного фильтра fn1 = 3800; fn2 = 4200; b = fir1(porNerTupy, [fn1*2/Fs fn2*2/Fs], 'stop', 'scale'); yFil = filter(b, 1, Ywym); figure(3); strips(yFil(1:length(yFil),:),10000); %розрахунок ачх нерекурсивного фильтра figure(4); wavplay(yFil, Fs); f = 0 : 1000; h = freqz(b, 1, f, Fs); semilogx(f, abs(h)); grid on; % wavwrite(Ywym, Fs, bits, 'd:\school\Системи запису та відтворення інформації\new11_wym.wav'); % wavwrite(yFil, Fs, bits, 'd:\school\Системи запису та відтворення інформації\new11_fil.wav'); %розрахунок рекурсивного фильтра figure(5); fn1 = 3800; fn2 = 4200; Wp = fn1 * 2/Fs; Ws = fn2 * 2/Fs; [n, Wn] = cheb1ord(Wp, Ws, 0.1, 60); [B, A] = cheby1(porRecTupy,60, [fn1*2/Fs fn2*2/Fs], 'stop'); yFil2 = filter(B, A, Ywym); strips(yFil2(1:length(yFil2),:),10000); %розрахунок ачх рекурсивного фильтра figure (6); wavplay(yFil, Fs); f = 10 : 0.1:14000; h = freqz(B, A, f, Fs); semilogx(f,abs(h));grid on; Графіки сигналу і АЧХ фільтрів: Початковий сигнал: / Сигнал після накладання завади: / АЧХ нерекурсивного фільтра: / Сигнал після фільтрації нерекурсивним фільтром: / АЧХ рекурсивного фільтра: / Сигнал після фільтрації рекурсивним фільтром: / Висновок: Після виконання цієї лабораторної роботи я ознайомився із середовищем MatLab та отримав базові навики роботи, вивчив основні можливості MatLab по обробці аудіоінформації, а також дослідив основні принципи і методи, які дозволяють реалізувати теоретичні положення про обробку аудіосигналів на практиці. Я побачив, як працюють цифрові фільтри, і процес очищення звуку від шуму. Після виконання цієї лабораторної роботи, я переконався, щодо переваг цифрових фільтрів. На основі отриманих АЧХ можна краще зрозуміти принцип роботи рекурсивних і нерекурсивних фільтрів, смугових і режекторних фільтрів, і причину їх використовують для виділення певного сигналу на робочих частотах.