Постоянно обновляющая передаточная функция

Question

У меня есть функция передачи, которую я пытаюсь использовать для фильтрации данных ускорения.

До сих пор я мог использовать lsim с волной sin около 10 секунд, и я получаю ожидаемый результат. Однако я не могу работать, как получить данные в функцию в режиме реального времени.

Чтобы уточнить, каждые 0,1 секунды я получаю значение ускорения от внешней программы. Мне нужно отфильтровать значения для удаления высокочастотных вариаций в данных. Мне нужно, чтобы это произошло для каждой точки данных, которую я получаю, когда я использую текущее значение ускоренного фильтра в дополнительных этапах обработки.

Как использовать функцию передачи в непрерывном режиме и обновлять выходное значение каждый раз, когда будут получены новые данные?

Answer 1

Это пример того, как сделать это с filter:

filter_state = []; % start with empty filter history 
while not_stopped() 
    x = get_one_input_sample(); 
    [y, filter_state] = filter(B, A, x, filter_state); 
    process_one_output_sample(y); 
end; 

Обратите внимание, что вам нужно использовать расширенную форму фильтра, в котором вы передаете историю фильтра от одной итерации к следующей используя переменную filter_state. B и A являются коэффициентами вашего дискретного фильтра. Если ваш фильтр непрерывный, вам необходимо сначала преобразовать его в дискретный фильтр, используя c2d или около того. Обратите внимание, что если ваш фильтр очень сложный, вам может потребоваться разделить фильтрацию в несколько этапов (например, один этап второго этапа на каждый фильтр), чтобы избежать числовых проблем.