Построение АЧХ по функции H(z)

Форум, посвященный общим вопросам, связанным с MATLAB

Модератор: Admin

=Anatoliy=
Пользователь
Сообщения: 160
Зарегистрирован: Вт ноя 18, 2008 1:53 pm

Построение АЧХ по функции H(z)

Сообщение =Anatoliy= » Вс фев 10, 2019 4:19 pm

Например есть передаточная функция фильтра H(z) = (1 - z^-D)/(1 - z^-1). Как по этому выражению в матлабе построить АЧХ этого фильтра?

Grizzly
Пользователь
Сообщения: 841
Зарегистрирован: Сб май 28, 2011 2:00 am

Re: Построение АЧХ по функции H(z)

Сообщение Grizzly » Вс фев 10, 2019 4:39 pm

Заменить z на exp(1j*w).

=Anatoliy=
Пользователь
Сообщения: 160
Зарегистрирован: Вт ноя 18, 2008 1:53 pm

Re: Построение АЧХ по функции H(z)

Сообщение =Anatoliy= » Вс фев 10, 2019 7:48 pm

Спасибо!
Ну тогда уже на exp(-1j*w/wa). Это единственный способ? Специальных функций нет?

barkot
Пользователь
Сообщения: 41
Зарегистрирован: Вс янв 14, 2018 2:28 pm

Re: Построение АЧХ по функции H(z)

Сообщение barkot » Вс фев 10, 2019 9:20 pm

Специальных функций нет?
bode

sandy
Эксперт
Сообщения: 5499
Зарегистрирован: Ср сен 22, 2004 4:49 pm

Re: Построение АЧХ по функции H(z)

Сообщение sandy » Вс фев 10, 2019 10:06 pm

=Anatoliy= писал(а):Специальных функций нет?

freqz, fvtool
С уважением

Александр Сергиенко

=Anatoliy=
Пользователь
Сообщения: 160
Зарегистрирован: Вт ноя 18, 2008 1:53 pm

Re: Построение АЧХ по функции H(z)

Сообщение =Anatoliy= » Пн фев 11, 2019 9:16 am

barkot писал(а):
Специальных функций нет?
bode

Замечательно! Очень похоже на то что я искал. Спасибо!

=Anatoliy=
Пользователь
Сообщения: 160
Зарегистрирован: Вт ноя 18, 2008 1:53 pm

Re: Построение АЧХ по функции H(z)

Сообщение =Anatoliy= » Пн фев 11, 2019 9:18 am

sandy писал(а):
=Anatoliy= писал(а):Специальных функций нет?

freqz, fvtool

Честно говоря не очень представляю как freqz можно подсунуть уравнение H(z) = (1 - z^-D)/(1 - z^-1). Нет ли какого либо примера?

=Anatoliy=
Пользователь
Сообщения: 160
Зарегистрирован: Вт ноя 18, 2008 1:53 pm

Re: Построение АЧХ по функции H(z)

Сообщение =Anatoliy= » Пн фев 11, 2019 11:18 am

barkot писал(а):bode

Сразу не спросил, но здесь же задаётся аналоговая функция, мне же нужна дискретная передаточная функция. c2d здесь не помогает так как коэффициенты у дискретной не совпадают с коэффициентами аналоговой. Можете что-то посоветовать? Я преобразовываю H() в Hd(), а мне нужно сразу задать Hd().Как его задать?

=Anatoliy=
Пользователь
Сообщения: 160
Зарегистрирован: Вт ноя 18, 2008 1:53 pm

Re: Построение АЧХ по функции H(z)

Сообщение =Anatoliy= » Пн фев 11, 2019 11:37 am

sandy писал(а):
=Anatoliy= писал(а):Специальных функций нет?

freqz, fvtool

Большое Вам спасибо! Наконец-то я всё понял.

barkot
Пользователь
Сообщения: 41
Зарегистрирован: Вс янв 14, 2018 2:28 pm

Re: Построение АЧХ по функции H(z)

Сообщение barkot » Пн фев 11, 2019 7:55 pm

=Anatoliy= писал(а):
barkot писал(а):bode

но здесь же задаётся аналоговая функция, мне же нужна дискретная передаточная функция.

bode работает с моделями непрерывного и дискретного времени. help bode:
For discrete-time models with sample time Ts, bode uses the transformation Z = exp(j*W*Ts) to map the unit circle to the real
frequency axis. The frequency response is only plotted for frequencies smaller than the Nyquist frequency pi/Ts, and the default
value 1 (time unit) is assumed when Ts is unspecified.

Пример сравнения АЧХ непрерывной H и дискретной Hd (Ts=0.1s) моделей

Код: Выделить всё

H = tf([1 0.1 7.5],[1 0.12 9 0 0]);
Hd = c2d(H,0.1,'zoh');
bode(H,'r',Hd,'b--')

=Anatoliy=
Пользователь
Сообщения: 160
Зарегистрирован: Вт ноя 18, 2008 1:53 pm

Re: Построение АЧХ по функции H(z)

Сообщение =Anatoliy= » Вт фев 12, 2019 9:30 am

barkot писал(а):

Извиняюсь, но Вы не поняли. В вашем примере передаточная функция задаётся для непрерывной модели, а потом с помощью c2d пересчитывается в дискретную модель. У меня же ПФ сразу дискретная, вопрос был как её подсунуть в bode. Александр подсказал ответ:
Изображение

barkot
Пользователь
Сообщения: 41
Зарегистрирован: Вс янв 14, 2018 2:28 pm

Re: Построение АЧХ по функции H(z)

Сообщение barkot » Вт фев 12, 2019 12:56 pm

=Anatoliy= писал(а):
У меня же ПФ сразу дискретная, вопрос был как её подсунуть в bode.


Пример:

Код: Выделить всё

Ts=0.001;
Dz=tf([1 1.1 1.2],[0.2 1.0 2.1],Ts)
bode(Dz)

=Anatoliy=
Пользователь
Сообщения: 160
Зарегистрирован: Вт ноя 18, 2008 1:53 pm

Re: Построение АЧХ по функции H(z)

Сообщение =Anatoliy= » Вт фев 12, 2019 3:04 pm

barkot писал(а):Пример:

Код: Выделить всё

Ts=0.001;
Dz=tf([1 1.1 1.2],[0.2 1.0 2.1],Ts)
bode(Dz)

Спасибо! Тоже хороший вариант.