Re: Прямоугольный импульс в крокозябру.
От: Rostislav_Pro  
Дата: 13.04.10 18:21
Оценка:
Ради интереса набросал примерчик на матлабе, который восстанавливает импульсную характеристику из кракозябры:


%передаточная ф-я, которую надо найти (НЧ фильтр баттерворта например)
[b,a] = butter(4,.1);

%реальная импульсная характеристика
imp_real = filter(b,a,[1 zeros(1,100)])

%прямоугольный импульс, кот. подаем на вход
square_imp = [1 1 1 1 1 zeros(1,100)];

%кракозябра
craco = filter(b,a,square_imp);

%инверсный фильтр, который из прямоугольного импульса делает почти
%единичный
b_inv = [1 -1];
a_inv = [1 0 0 0 0 -.99 0 0 0 0 -.01];

%восстановленная импульсная характеристика
imp_restored = filter(b_inv,a_inv,craco);

%рисуем рядом восстановленную хар-ку и реальную, и видим что очень похожи
plot(1:100,imp_real(1:100),1:100,imp_restored(1:100)+.01)
 
Подождите ...
Wait...
Пока на собственное сообщение не было ответов, его можно удалить.