Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Предоплата всего
Подписываем
МІНІСТЕРСТВО ОСВІТИ І НАУКИ, МОЛОДІ ТА СПОРТУ УКРАЇНИ
НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ «ЛЬВІВСЬКА ПОЛІТЕХНІКА»
ІКТА
кафедра ЗІ
З В І Т
до лабораторної роботи №1
з курсу: «Системи запису та відтворення інформації»
Виконав: ст. гр. ЗІ-31
Cтельмахович Н.П.
Прийняв: Совин Я.Р.
Львів 2013
МЕТА РОБОТИ
Отримати базові навики роботи в середовищі MATLAB, вивчити основні можливості MATLAB по обробці аудіоінформації.
ЗАВДАННЯ
|
Варіант |
Назва файлу |
Частота завади |
Амплітуда завади, В |
Порядок не рекурсивного фільтра N |
Тип рекурсивного фільтру |
Порядок рекурсивного фільтру |
|
13 |
Lab_1_13.wav |
3200 |
5 |
330 |
Еліптичний |
4 |
ТЕКСТ ПРОГРАМИ
Для рекурсивного фільтра
[y, Fs, bits] = wavread ('E:\Lab_1_13.wav');
Fs;
bits;
N=length(y);
t=(0:N-1)*(1/Fs);
Zavada=5*sin(2*pi*3200*t);
S1=y+Zavada';
N=4;
Fn = Fs/2;
F=(3100:3300);
[B,A]=ellip(N, 0.1, 100, [3150/Fn 3250/Fn], 'stop');
h=freqz(B, A, F, Fs);
y=filter(B, A, S1);
figure(1);
subplot(4,1,1);
plot(t,y); grid on; xlim([0 t(end)]); title('Оригінал');
subplot(4,1,2);
plot(t,Zavada); grid on; xlim([0 t(end)]); title('Завада');
subplot(4,1,3);
plot(t,S1); grid on; xlim([0 t(end)]); title('Шум');
subplot(4,1,4);
plot(t, y); grid on; title('Відфільтрований рекурсивний');
figure(2);
plot(F, abs(h)); grid on; title('Рекурсивний АЧХ');
Для не рекурсивного фільтра
[y, Fs, bits] = wavread ('E:\Lab_1_13.wav');
Fs;
bits;
N=length(y);
t=(0:N-1)*(1/Fs);
Zavada=5*sin(2*pi*3200*t);
S1=y+Zavada';
N=330;
Fn = Fs/2;
F=(3100:3300);
b=fir1( N, [3150/Fn 3250/Fn], 'stop');
h=freqz(b, 1, F, Fs);
y=filter(b, 1, S1);
figure(1);
subplot(4,1,1);
plot(t,y); grid on; xlim([0 t(end)]); title('Оригінал');
subplot(4,1,2);
plot(t,Zavada); grid on; xlim([0 t(end)]); title('Завада');
subplot(4,1,3);
plot(t,S1); grid on; xlim([0 t(end)]); title('Шум');
subplot(4,1,4);
plot(t, y); grid on; title('Відфільтрований не рекурсивний');
figure(2);
plot(F, abs(h)); grid on; title('Не рекурсивний АЧХ');
РЕЗУЛЬТАТ ВИКОНАННЯ ПРОГРАМИ
АЧХ рекурсивного фільтра:
АЧХ не рекурсивного фільтра:
ВИСНОВКИ
Завдання даної лабораторної роботи полягало у ознайомленні з середовищем MATLAB, вивченні основних можливостей MATLAB по обробці аудіоінформації.
Завантаживши і прослухавши wav-файл, за допомогою MATLAB вивів основну інформації про файл. Потім на корисний вхідний сигнал накладалася завада у вигляді синусоїдального сигналу з вказаними параметрами, а саме частотою 3200 Гц та амплітудою 5В. Для усунення завади розраховувалися коефіцієнти не рекурсивного фільтра та рекурсивного вказаного порядку N.
Для відповідних фільтрів було побудовано АЧХ.