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

 
 
 
Reply to this topicStart new topic
> фильтр бесселя в матлаб, цифровой
shf_05
сообщение Sep 19 2008, 05:05
Сообщение #1


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

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



можно ли рассчитать цифровой фильтр бесселя в матлабе в версии 7,0 нашел только аналоговый?
Go to the top of the page
 
+Quote Post
uriy
сообщение Sep 19 2008, 08:48
Сообщение #2


Гуру
******

Группа: Свой
Сообщений: 2 429
Регистрация: 30-11-05
Из: Ижевск
Пользователь №: 11 606



А как же fdatool целый пакет в матлабе для рассчета фильтров.
Go to the top of the page
 
+Quote Post
shf_05
сообщение Sep 19 2008, 09:18
Сообщение #3


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

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



Цитата(uriy @ Sep 19 2008, 14:48) *
А как же fdatool целый пакет в матлабе для рассчета фильтров.

список доступного в FDA БИХ фильтры:
buterwort
chebishev1
chebishev2
elliptic
max flat
least P-norm
const least P-norm
Go to the top of the page
 
+Quote Post
AsJohnAs
сообщение Sep 23 2008, 08:39
Сообщение #4


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

Группа: Свой
Сообщений: 125
Регистрация: 14-07-05
Из: Санкт-Петербург
Пользователь №: 6 793



Ну в Matlab есть неплохой хелп. Вот в нем указано что есть функция:
besselj

Bessel function of first kind
Syntax

J = besselj(nu,Z)
J = besselj(nu,Z,1)
[J,ierr] = besselj(nu,Z)
Go to the top of the page
 
+Quote Post
shf_05
сообщение Sep 26 2008, 07:18
Сообщение #5


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

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



Цитата(AsJohnAs @ Sep 23 2008, 14:39) *
Ну в Matlab есть неплохой хелп. Вот в нем указано что есть функция:
besselj

Bessel function of first kind
Syntax

J = besselj(nu,Z)
J = besselj(nu,Z,1)
[J,ierr] = besselj(nu,Z)


я знаю про неплохость хелпа в матлабе)). фильтр бесселя и функция бесселя к-го порядка от аргумента нужны для разных задач (если я не ошибаюсь конечно), см. все тот же хелп.
Go to the top of the page
 
+Quote Post
AsJohnAs
сообщение Sep 26 2008, 20:17
Сообщение #6


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

Группа: Свой
Сообщений: 125
Регистрация: 14-07-05
Из: Санкт-Петербург
Пользователь №: 6 793



Но ведь передаточная функция фильтров Бесселя записывается в виде отношения коэффициента к функци Бесселя n-го порядка.

Так что вроде все есть.
Go to the top of the page
 
+Quote Post
shf_05
сообщение Sep 30 2008, 11:35
Сообщение #7


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

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



Цитата(AsJohnAs @ Sep 27 2008, 02:17) *
Но ведь передаточная функция фильтров Бесселя записывается в виде отношения коэффициента к функци Бесселя n-го порядка.

Так что вроде все есть.


если можно опишите поподробнее, как из функции бесселя получить фильтр бесселя скажем 4-го порядка, частота Fs=8000, Fpass= 500?
Go to the top of the page
 
+Quote Post
AsJohnAs
сообщение Sep 30 2008, 18:49
Сообщение #8


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

Группа: Свой
Сообщений: 125
Регистрация: 14-07-05
Из: Санкт-Петербург
Пользователь №: 6 793



Цитата(shf_05 @ Sep 30 2008, 15:35) *
если можно опишите поподробнее, как из функции бесселя получить фильтр бесселя скажем 4-го порядка, частота Fs=8000, Fpass= 500?

Ну можно например так:
[b,a]=besself(5,500);
f=0:1:8000;
freqs(b,a,f);
Go to the top of the page
 
+Quote Post
shf_05
сообщение Oct 2 2008, 04:36
Сообщение #9


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

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



Цитата(AsJohnAs @ Oct 1 2008, 00:49) *
Ну можно например так:
[b,a]=besself(5,500);
f=0:1:8000;
freqs(b,a,f);


что-то мы др. др. не понимаем:

besself- это аналоговый фильтр, а как из него сделать цифровой?

Вы писали про besselj, вот я спрашивал "если можно опишите поподробнее, как из функции бесселя получить фильтр бесселя скажем 4-го порядка, частота Fs=8000, Fpass= 500?"
Go to the top of the page
 
+Quote Post
AsJohnAs
сообщение Oct 3 2008, 09:03
Сообщение #10


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

Группа: Свой
Сообщений: 125
Регистрация: 14-07-05
Из: Санкт-Петербург
Пользователь №: 6 793



Мндя... да я не прав.
Да и ответ у меня тогда не такой уж и простой будет...
Правельно это сделать через расчет полюсов т.д. а там ведь функция бесселя у которой нули ваще какие-то не такие...

Вообщем есть такая идея:
[b,a]=besself(5,0.2); % надо нормированную частоту, а то весь матлаб плющит
randn('state',0);
noise = randn(50000,1);
x=filter(b,a,noise);
h = lpc(x,15);

Вообщем применям линейное предсказание. Вообщем задачу можно решить быстро и с нужной точностью, но так сказать "в лоб".
А и еще: h - это импульсная которую можно и в fdatool загнать
Go to the top of the page
 
+Quote Post
shf_05
сообщение Oct 3 2008, 11:59
Сообщение #11


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

Группа: Свой
Сообщений: 1 143
Регистрация: 22-04-08
Из: г. Екатеринбург
Пользователь №: 36 992



дак все это ясно, не хочется велосипед изобретать,....
PS: а есть еще ф-я impinvar.
Go to the top of the page
 
+Quote Post
PlainUser
сообщение Apr 19 2016, 06:39
Сообщение #12


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

Группа: Участник
Сообщений: 84
Регистрация: 23-06-05
Пользователь №: 6 244



Нужен фильтр бесселя в fdatools , не нашел.
Может он там есть но я его не вижу?

Ну скажем "maximaly flat" это не он случайно?

Сообщение отредактировал PlainUser - Apr 19 2016, 06:42
Go to the top of the page
 
+Quote Post
Fat Robot
сообщение Apr 19 2016, 06:58
Сообщение #13


ʕʘ̅͜ʘ̅ʔ
*****

Группа: Свой
Сообщений: 1 008
Регистрация: 3-05-05
Пользователь №: 4 691



doc besselap
doc bilinear
doc zp2tf
Go to the top of the page
 
+Quote Post
PlainUser
сообщение Apr 21 2016, 05:41
Сообщение #14


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

Группа: Участник
Сообщений: 84
Регистрация: 23-06-05
Пользователь №: 6 244



Цитата(PlainUser @ Apr 19 2016, 09:39) *
Нужен фильтр бесселя в fdatools , не нашел.
Может он там есть но я его не вижу?

Ну скажем "maximaly flat" это не он случайно?


Расследование выявило что в матлабе цифровой фильтр бесселя отсутствует принципиально.
Ввиду ....

" Аналоговые фильтры Бесселя характеризуются практически постоянной групповой задержкой в полосе пропускания, что позволяет сохранить форму пропускаемых через фильтр сигналов (если их спектр сосредоточен в полосе пропускания фильтра). Дискретные фильтры Бесселя не сохраняют это свойство, поэтому функция besself не реализует синтез дискретных фильтров Бесселя."

Его можно родить но это будет не он.
Go to the top of the page
 
+Quote Post
Fat Robot
сообщение Apr 21 2016, 08:50
Сообщение #15


ʕʘ̅͜ʘ̅ʔ
*****

Группа: Свой
Сообщений: 1 008
Регистрация: 3-05-05
Пользователь №: 4 691



Вполне сохраняют, если Fc < Fs/4

Цитата(PlainUser @ Apr 21 2016, 09:41) *
Дискретные фильтры Бесселя не сохраняют это свойство, поэтому функция besself не реализует синтез дискретных фильтров Бесселя.

Go to the top of the page
 
+Quote Post

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

 


RSS Текстовая версия Сейчас: 18th June 2025 - 08:27
Рейтинг@Mail.ru


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