Разработать КИХ фильтр и реализовать это все дело в matlab-е с использованием окна Хэмминга

@sendrom · Регистрация: 13.12.2009

Author24 — интернет-сервис помощи студентам

Уважаемые форумчане помогите пожалуйста разобраться вот с этой лабораторной работой:
https://docs.google.com/docume... W1DodhTpV0

Здесь нужно разработать КИХ фильтр и реализовать это все дело в matlab-е с использованием окна Хэмминга.

Теоретически знаю вот что:
В матлабе открываем аудио файл, далее командой sound(y,32000) можно его воспроизвести. В этом файле искуственно создан шум который и нужно отфильровать КИХ фильтром. Препод сказал так: нужно взять 128 точек - это число самого фильтра далее там что-то вроде обратного преобразования фурье умножения на окно хэмминга в итоге прослушиваем этот сигнал и там играет музыка без шума. Вот как это сделать подскажите пожалуйста. Файл с аудио в формате матлаб прикрепляю.

@vital792 · 17.12.2012, 17:56

Matlab M

clc; clear all; close all;
 
load('noised_audio.mat');
 
% sp = abs(fft(y));
% plot(linspace(0,1,length(sp)/2), sp(1:end/2));
% посмотрели на спектр, увидели куда препод спрятал шум
Wn = 0.6; % и соответственно выбираем полосу фильтра
 
sound(y(1:8192)); % 1 секунда (зажать уши)
 
n = 128;
% b = fir1(n, Wn, 'low');
% как выглядит обратное фурье преобразование единичного импульса и так
% известно, можно его записать сразу:
b = Wn * sin(pi*(-(n/2+0.00001):(n/2-0.00001))*Wn) ./ (pi*(-(n/2+0.00001):(n/2-0.00001))*Wn);
b = b .* hamming(n+1)'; % домножаем на окно
b = b ./ sum(b); % нормируем
y1 = filtfilt(b, 1, y); % и фильтруем
 
% sp1 = abs(fft(y1));
% plot(linspace(0,1,length(sp1)/2), sp1(1:end/2));
 
sound(y1);

@sendrom · 17.12.2012, 20:43 **[ТС]**

Большое спасибо очень выручили, буду вникать в суть.

@Алексей888 · 09.06.2015, 20:22

подскажите а можно сюда запустит дискретный сигнал в этот фильтр (A*sin(2*pi*f0*t) нужно чтобы фильтр убирал помехи на частоте 0,2-2 Гц

@sendrom 1 / 1 / 0 Регистрация: 13.12.2009 Сообщений: 10
	17.12.2012, 20:43 [ТС]	3
	Большое спасибо очень выручили, буду вникать в суть. 0

@Алексей888 0 / 0 / 0 Регистрация: 09.03.2015 Сообщений: 8
	09.06.2015, 20:22	4
	подскажите а можно сюда запустит дискретный сигнал в этот фильтр (Asin(2pif0t) нужно чтобы фильтр убирал помехи на частоте 0,2-2 Гц 0