МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ «ЛЬВІВСЬКА ПОЛІТЕХНІКА» ІКТА Кафедра ЗІ / ЗВІТ До лабораторної роботи №1 Знайомство з середовищем Matlab. Дослідження методів обробки аудіо інформації в Matlab. / // Лістинг програми [y, Fs, bits]= wavread('H:\навчання\Лабораторні роботи СЗіВІ\Аудіосигнали\Lab_1_15.wav'); Fs bits Nr=320; a=6; f=1300; N = length(y); t = ( 0 : N - 1) / Fs; Noise = a * sin(2 * pi * f * t); S = Noise' + y; figure(2); subplot(2, 1, 1); plot( t, y); grid on; subplot(2, 1, 2); plot( t,S); grid on; fn=Fs/2; F=[1200 1400]; b=fir1(Nr,F/fn,'stop'); y1=filter(b,1,S); figure(3); subplot(3, 1, 1); plot( t, y1); grid on; F1=[800:1700]; subplot(3, 1, 2); plot( t,S); grid on; h=freqz(b,1,F1,Fs); figure(4); plot(F1,abs(h));grid on; wavwrite(y1,Fs,bits,'lab11'); Rp=0.1; Rs=60; Wp=[10 2000]/fn; Ws=[1295 1305]/fn; [n,Wn]=cheb2ord(Wp,Ws,Rp,Rs); n=5; [B,A]=cheby2(n,Rs,Wn,'stop'); F1=[800:1400]; h=freqz(B,A,F1,Fs); figure(5); plot(F1,abs(h));grid on;; y2=filter(B,A,S); figure(6); subplot(3,1,2) plot(t,y2);grid on; subplot(3,1,1) plot( t,S); grid on; wavplay(y2,Fs,'async') wavwrite(y2,Fs,bits,'lab12'); Графіки АЧХ Нерекурсивного фільтра Рекурсивного фільтра / / Сигнал до , після накладання завад та після фільтрації фільтрами Рекурсивним Нерекурсивним / / Висновок Фільтри призначені для послаблення сигналів або завад певної частоти. Фільтри бувають аналогові і цифрові. Аналогові поділяються на активні та пасивні. Цифрові - на рекурсивні та не рекурсивні. Перевага рекурсивних над нерекурсивними в тому, що для отримання певної АЧХ порядок рекурсивного має бути значно нижчим ніж не рекурсивного.