Перейти к содержанию
    

Узнать полюсы и нули фильтра по измеренной АЧХ (Матлаб)

Здравствуйте!

 

Есть снятая с реального устройства АЧХ. Хочется узнать полюсы и нули фильтра. То есть проанализировать весь тракт, представив его как фильтр.

Это вообще возможно? и как? Применительно к Матлабу, если можно.

Достаточно ли самой АЧХ для этого?

 

Чего сделали:

1. с помощью polyfit() вычислена кривая по точкам АЧХ.

2. через invfreqz() получаем коэффициенты

3. через fvtool() разглядываем что получилось.

Так правильно?

 

исходные данные (56 точек АЧХ) DataIn.mat

Результат fvtool() ниже.

(АЧХ правильная)

image.thumb.png.fd5be14aac861d1ba18553b320dbd5f2.png

 

Полюсы и нули. Это нормально? Так может быть?

image.thumb.png.1222d8de1896fd737984437627aa9b1e.png

 

image.png.02f26c195e098f60896e8a63145cfd92.png

 

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

1, 2 вроде норм.

функцией в п.3 не пользовался.

На 100% не уверен, но последняя картинка при вашей АФЧХ  сомнительно, что может быть так что все точки на действительной отрицательной оси.

 

по идее надо создать объект для системы (например функцией tf)

А дальше pole, zero...

 

 

Изменено пользователем aBoomest

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

10 minutes ago, Lmx2315 said:

имхо по одной АЧХ нельзя, нужно ещё ФЧХ.

Можно. Функция КЧХ является аналитической. Ваш, КО.

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

1 минуту назад, blackfin сказал:

Можно. Функция КЧХ является аналитической. Ваш, КО.

Где вы увидели в первом сообщении что у них КЧХ, я не вижу. А АЧХ получается как модуль КЧХ , фарш назад не провернуть.

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

32 minutes ago, Lmx2315 said:

Где вы увидели в первом сообщении что у них КЧХ, я не вижу. А АЧХ получается как модуль КЧХ , фарш назад не провернуть.

ФЧХ можно легко вычислить, зная АЧХ. Учите мат-часть: Соотноше́ния Кра́мерса — Кро́нига

 

PS. Нет, я был не прав..

 

Зная Re[K(w)] можно вычислить Im[K(w)] и наоборот. С модулем K(w) такой фокус, кажется, не пройдет.

 

PPS. Но в любом случае, судя по графику из первого поста, ФЧХ известна.

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

4 минуты назад, blackfin сказал:

ФЧХ можно легко вычислить, зная АЧХ. Учите мат-часть: Соотноше́ния Кра́мерса — Кро́нига

спасибо.

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

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

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

1 час назад, Lmx2315 сказал:

 

 

1 час назад, blackfin сказал:

Можно. Функция КЧХ является аналитической. Ваш, КО.

Нельзя. Модуль.

Изменено пользователем thermit

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

03.03.2022 в 20:45, Ruslan1 сказал:

исходные данные (56 точек АЧХ) DataIn.mat

Результат fvtool() ниже.

(АЧХ правильная)

Исходные данные содержат модуль АЧХ, это только у меня так? Или такие исходные данные? Тогда ловить нечего...

Да, результат неправильный, АЧХ может и правильная, но вот ФЧХ реализуемой системы не может возрастать, это собственно и на карте нулей/полюсов видно, у Вас один полюс вне единичной окружности, такая система физически невозможна...

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

6 минут назад, thermit сказал:

Возможна. К физической реализуемости это никакого отношения не имеет.

Да-да, Матлаб все стерпит...

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

22 hours ago, Lmx2315 said:

имхо по одной АЧХ нельзя, нужно ещё ФЧХ.

Вроде написал АФЧХ, подразумевая несомненно комплексы.

Изменено пользователем aBoomest

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Присоединяйтесь к обсуждению

Вы можете написать сейчас и зарегистрироваться позже. Если у вас есть аккаунт, авторизуйтесь, чтобы опубликовать от имени своего аккаунта.

Гость
К сожалению, ваш контент содержит запрещённые слова. Пожалуйста, отредактируйте контент, чтобы удалить выделенные ниже слова.
Ответить в этой теме...

×   Вставлено с форматированием.   Вставить как обычный текст

  Разрешено использовать не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отображать как обычную ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставлять изображения напрямую. Загружайте или вставляйте изображения по ссылке.

×
×
  • Создать...