O
omidi_sbu
Guest
Ця програма MATLAB для прийняття відповідного ядра вейвлет
Воно випливає фільтр низьких шляхом відбору проб з аналогового косінусоідальним функції в частотній області і беручи ОБПФ цієї функції поступитися імпульсного відгуку.(Іншими фільтрами, виводяться за допомогою ортогональних відносини).
Перш за все я хочу знати, якщо метод проектування фільтр має яких-небудь проблем?
Фільтр має дві параметрів в частотній області Ф і Ф B. це число зразків, які будуть довжини фільтра і B, який контролює чіткість кривої в частотній області.
Ці фільтри мають досить велику погрішність реконструкції.навіть тоді, коли довжина фільтра великий (Fs = 32 і B =. 15)
Коли я приймаю БПФ цей фільтр з великою кількістю точки ШПФ (FFT (TZ, 8000)) результат показує, дрібна брижі в смузі.Я перевірив ж в БД фільтрів і ці фільтри не було цих пульсацій.
Я заплутався!======================================= М-файл
% Програма для створення відповідного ядра для кодера
очистити всі;
закрити всі;
Fs = 16;% вибірки frequncy від аналогового фільтру T (F)
% B =. 1505, вхідний параметр%
B =. 15;
р = 1 / 2;% г (другий параметр)
А = Fs / 2;
для / = 1: K 1% вибірки аналогового фільтра в частотній області
якщо (F-1) <= (РБ) * А
T (F) = 1;
ELSEIF (F-1) <= (г B) * А
T (F) = ство (р / (4 * B) * ((F-1) / А-Г B));
ELSEIF (F-1) <= А
T (F) = 0;
цілі;
цілі;
х = (0: А) / Ф;
ділянка (х, T);
для / = (А 2): Fs% рішень другій половині зразків
T (F) = T (Fs 2- е);
цілі;
T1 = реал (ОБПФ (T));% оригінальний КИХ-фільтра
для я = 1: K% TL випадково фільтр
TL (я K) = T1 (I);
цілі;
для я = 1: K
TL (я) = T1 (K I);
цілі;
TX = TL;
%%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$%%
[Y, FS] = wavread ('taha.wav');% отримувати вхідний сигнал (16KHz WAV-файл)
ділянка (у);
Назва ('Original');
вхід = у ';
%------------------------------------------------- -------------------------
для я = 1: довжина (TX) звітності% зразків фільтра (ядро)
(Г) = fprintf ('% 0,16 F', штат Техас (я));
цілі;
[Лод, HID, Лор, Hir] = orthfilt (TX);% рішень іншими фільтрами
QL = filterx1 (введення, Лод);% аналізу вхідного сигналу
QH = filterx1 (введення, HID);
QLD = d_convrt (ЯЖ);% до перетворення
QHD = d_convrt (Qh);
qldu = u_convrt (QLD);% до перетворення
qhdu = u_convrt (QHD);
outl = filterx1 (qldu, Лор);% синтезу фільтрації
ЖНИЙ = filterx1 (qhdu, Hir);
final_out = outl ЖНИЙ;% підсумовування вихідних синтезу фільтрів
тримати;
% Дільниці (final_out, 'G');
ділянки (final_out (Fs: Ф-1 довжина (вхід)), 'м'); резюме доповіді%
Назва ('реконструйованих та оригінальний сигнал');
% Пауза;
фігури;
ділянки (final_out (Fs: Ф-1 довжина (вхід))-вхід, 'м');
заголовок ('Помилка');
у.о. = final_out (Fs: Ф-1 довжина (вхід))-вхід;
L = довжина (Е);
E = SQRT (1 / L * сума (Є. ^ 2))
Воно випливає фільтр низьких шляхом відбору проб з аналогового косінусоідальним функції в частотній області і беручи ОБПФ цієї функції поступитися імпульсного відгуку.(Іншими фільтрами, виводяться за допомогою ортогональних відносини).
Перш за все я хочу знати, якщо метод проектування фільтр має яких-небудь проблем?
Фільтр має дві параметрів в частотній області Ф і Ф B. це число зразків, які будуть довжини фільтра і B, який контролює чіткість кривої в частотній області.
Ці фільтри мають досить велику погрішність реконструкції.навіть тоді, коли довжина фільтра великий (Fs = 32 і B =. 15)
Коли я приймаю БПФ цей фільтр з великою кількістю точки ШПФ (FFT (TZ, 8000)) результат показує, дрібна брижі в смузі.Я перевірив ж в БД фільтрів і ці фільтри не було цих пульсацій.
Я заплутався!======================================= М-файл
% Програма для створення відповідного ядра для кодера
очистити всі;
закрити всі;
Fs = 16;% вибірки frequncy від аналогового фільтру T (F)
% B =. 1505, вхідний параметр%
B =. 15;
р = 1 / 2;% г (другий параметр)
А = Fs / 2;
для / = 1: K 1% вибірки аналогового фільтра в частотній області
якщо (F-1) <= (РБ) * А
T (F) = 1;
ELSEIF (F-1) <= (г B) * А
T (F) = ство (р / (4 * B) * ((F-1) / А-Г B));
ELSEIF (F-1) <= А
T (F) = 0;
цілі;
цілі;
х = (0: А) / Ф;
ділянка (х, T);
для / = (А 2): Fs% рішень другій половині зразків
T (F) = T (Fs 2- е);
цілі;
T1 = реал (ОБПФ (T));% оригінальний КИХ-фільтра
для я = 1: K% TL випадково фільтр
TL (я K) = T1 (I);
цілі;
для я = 1: K
TL (я) = T1 (K I);
цілі;
TX = TL;
%%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$%%
[Y, FS] = wavread ('taha.wav');% отримувати вхідний сигнал (16KHz WAV-файл)
ділянка (у);
Назва ('Original');
вхід = у ';
%------------------------------------------------- -------------------------
для я = 1: довжина (TX) звітності% зразків фільтра (ядро)
(Г) = fprintf ('% 0,16 F', штат Техас (я));
цілі;
[Лод, HID, Лор, Hir] = orthfilt (TX);% рішень іншими фільтрами
QL = filterx1 (введення, Лод);% аналізу вхідного сигналу
QH = filterx1 (введення, HID);
QLD = d_convrt (ЯЖ);% до перетворення
QHD = d_convrt (Qh);
qldu = u_convrt (QLD);% до перетворення
qhdu = u_convrt (QHD);
outl = filterx1 (qldu, Лор);% синтезу фільтрації
ЖНИЙ = filterx1 (qhdu, Hir);
final_out = outl ЖНИЙ;% підсумовування вихідних синтезу фільтрів
тримати;
% Дільниці (final_out, 'G');
ділянки (final_out (Fs: Ф-1 довжина (вхід)), 'м'); резюме доповіді%
Назва ('реконструйованих та оригінальний сигнал');
% Пауза;
фігури;
ділянки (final_out (Fs: Ф-1 довжина (вхід))-вхід, 'м');
заголовок ('Помилка');
у.о. = final_out (Fs: Ф-1 довжина (вхід))-вхід;
L = довжина (Е);
E = SQRT (1 / L * сума (Є. ^ 2))