реклама на сайте
подробности

 
 
 
Reply to this topicStart new topic
> Реализация цифрового фильтра, на основе заданной передаточной функции
Impartial
сообщение Apr 12 2017, 22:58
Сообщение #1


Участник
*

Группа: Участник
Сообщений: 58
Регистрация: 11-06-10
Пользователь №: 57 871



Задана передаточная функция фильтра
W(p) = (0.067*p+1)^4/((0.5*p+1)^2 * (0.008*p+1)^2)
Необходимо сделать синтез для ПЛИС Циклон 2 16 битными числами.
Думаю, что необходимо как то ввести функцию в матлаб и получить коэффициенты.
Подскажите с чего начать, чтобы уметь это делать.
Go to the top of the page
 
+Quote Post
XVR
сообщение Apr 13 2017, 10:03
Сообщение #2


Гуру
******

Группа: Свой
Сообщений: 3 022
Регистрация: 7-04-07
Из: Химки
Пользователь №: 26 847



Для начала поставить Simulink (в Matlab)


И еще вам может пригодится FDATool (из того же Matlab'а)
http://www.chipinfo.ru/literature/chipnews/200109/1.html
Go to the top of the page
 
+Quote Post
Gruffly
сообщение Apr 13 2017, 16:43
Сообщение #3


Частый гость
**

Группа: Участник
Сообщений: 103
Регистрация: 6-04-17
Пользователь №: 96 386



Цитата(Impartial @ Apr 12 2017, 22:58) *
Подскажите с чего начать, чтобы уметь это делать.

Для столь простых случаев достаточно ручками использовать билинейное преобразование.
W(p) = W1*W1*W1*W1/(W2*W2*W3*W3) - исходная Ваша ПФ.
Go to the top of the page
 
+Quote Post
Impartial
сообщение Apr 15 2017, 21:49
Сообщение #4


Участник
*

Группа: Участник
Сообщений: 58
Регистрация: 11-06-10
Пользователь №: 57 871



Почитал про билинейное преобразование. Букв W1,W2,W3 не нашел. Объясните подробнее.
Go to the top of the page
 
+Quote Post
Gruffly
сообщение Apr 16 2017, 00:23
Сообщение #5


Частый гость
**

Группа: Участник
Сообщений: 103
Регистрация: 6-04-17
Пользователь №: 96 386



Цитата(Impartial @ Apr 15 2017, 22:49) *
Почитал про билинейное преобразование. Букв W1,W2,W3 не нашел. Объясните подробнее.

"Букв" - очень много надо учить.

W(p) = (0.067*p+1)^4/((0.5*p+1)^2 * (0.008*p+1)^2)

W1(p) = (0.067*p+1);
W2(p) = (0.5*p+1);
W3(p) = (0.008*p+1).

В итоге:
W(p) = (W1*W1*W1*W1)/(W2*W2*W3*W3).

Сообщение отредактировал Gruffly - Apr 16 2017, 00:24
Go to the top of the page
 
+Quote Post
Impartial
сообщение Apr 16 2017, 22:31
Сообщение #6


Участник
*

Группа: Участник
Сообщений: 58
Регистрация: 11-06-10
Пользователь №: 57 871



Я догадывался, что такое W.
А как перейти от оператора "p" к оператору "z" ?
Go to the top of the page
 
+Quote Post
Gruffly
сообщение Apr 16 2017, 23:08
Сообщение #7


Частый гость
**

Группа: Участник
Сообщений: 103
Регистрация: 6-04-17
Пользователь №: 96 386



Цитата(Impartial @ Apr 16 2017, 22:31) *
Я догадывался, что такое W.
А как перейти от оператора "p" к оператору "z" ?

Билинейное преобразование, к примеру - уже советовал.
Там еще будет "искажение" частот, но это известно - как исправлять.

Сообщение отредактировал Gruffly - Apr 16 2017, 23:35
Go to the top of the page
 
+Quote Post
HardEgor
сообщение Apr 17 2017, 01:17
Сообщение #8


Профессионал
*****

Группа: Свой
Сообщений: 1 720
Регистрация: 3-03-06
Из: Tomsk
Пользователь №: 14 925



Цитата(Impartial @ Apr 17 2017, 05:31) *
А как перейти от оператора "p" к оператору "z" ?

Найти книжку "Цифровая обработка сигналов. Практический подход Айфичер Э., Джервис Б. 2004" там на пальцах объясняется переход.
Go to the top of the page
 
+Quote Post
Impartial
сообщение Apr 17 2017, 09:55
Сообщение #9


Участник
*

Группа: Участник
Сообщений: 58
Регистрация: 11-06-10
Пользователь №: 57 871



Я цифровик и вместо математики учил ее дискретный аналог.
Мне проще перейти к более понятным вещам, например ПИД регулятору, хотя это тоже фактически фильтр.
Приведенная выше передаточная функция предназначена для компенсация реакции механической системы.
Возможно ли из этой функции вытащить момент инерции механической системы и максимальные ускорение и скорость?
Go to the top of the page
 
+Quote Post
AlexandrY
сообщение Apr 17 2017, 12:04
Сообщение #10


Mentor
******

Группа: Модераторы
Сообщений: 5 332
Регистрация: 19-01-05
Пользователь №: 2 050



Цитата(Impartial @ Apr 17 2017, 12:55) *
Приведенная выше передаточная функция предназначена для компенсация реакции механической системы.
Возможно ли из этой функции вытащить момент инерции механической системы и максимальные ускорение и скорость?

Преобразование делается в Mathcad-е элементарно:
Прикрепленное изображение

T - это период сэмплирования. Для простоты приравнял единице.
Но если его не укажете, то Mathcad вам даст общую формулу для T переменной.
Момент инерции, а тем более максимальные ускорения и скорость из этого регулятора узнать невозможно.
Пытаться можно узнать инерцию, но надо знать критерии компенсатора и природу реакции(скорость, позиция... или что? )
Go to the top of the page
 
+Quote Post
Impartial
сообщение Apr 17 2017, 15:24
Сообщение #11


Участник
*

Группа: Участник
Сообщений: 58
Регистрация: 11-06-10
Пользователь №: 57 871



Контролируется позиция.
С этим фильтром наблюдается болтанка при позиционировании вокруг заданной точки.
Хотелось перейти к RAMP генератору скоростной траектории разгона-торможения.
Но для этого расчета нужны характеристики макс. ускорений и скоростей механической системы.
Можно, конечно, на глазок выбрать, но хочется точных цифр.
Go to the top of the page
 
+Quote Post
Gruffly
сообщение Apr 17 2017, 21:06
Сообщение #12


Частый гость
**

Группа: Участник
Сообщений: 103
Регистрация: 6-04-17
Пользователь №: 96 386




На форуме есть раздел "Электрические машины", "Обратные связи" - Вам, туда.
Go to the top of the page
 
+Quote Post

Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 19th October 2017 - 23:58
Рейтинг@Mail.ru


Страница сгенерированна за 0.01339 секунд с 7
ELECTRONIX ©2004-2016