Разработка, производство и продажа радиоэлектронной аппаратуры
|
Карта сайта
|
Пишите нам
|
В избранное
Требуется программист в Зеленограде
- обработка данных с датчиков; ColdFire; 40 тыс.
e-mail:
jobsmp@pochta.ru
Телесистемы
|
Электроника
|
Конференция «Цифровые сигнальные процессоры (DSP) и их применение»
ничего странного, зная, что в нуле ноль. 6 герц при 27 мкс - это нормированная переходная полоса 0.000081! Аццкая крутизна и соотв. добротность.
Отправлено
SM
16 января 2008 г. 19:30
В ответ на:
7 порядков - это очень странно IMHO
отправлено Oldring 16 января 2008 г. 19:19
Составить ответ
|
Вернуться на конференцию
Ответы
Автор, во-первых, неправильно посчитал, во-вторых, 7 порядков у него за счет малых коэффициентов числителя и больше не из-за чего, корректируется выносом к-та усиления за скобки. Он там маленький должен быть, этот к-т усиления, порядка 1/g, g - к-т связи аналоговых и цифровых частот=ctg(pi*FS/FD)~=2000 при его цифрах, соответвественно к-т усиления ЦФ~=0.0005. Вынесет его за скобки и исчезнут 7 порядков.
—
-=ВН=-
(16.01.2008 21:28:26
193.125.71.140
,
пустое
)
Хочу покаяться. Я больший тормоз, чем матлаб2007 на моей машине. Частота ср...{+}
—
visual_wind
(17.01.2008 13:32:46
217.12.205.67
, 418 байт)
Про матлаб ничего не подскажу. На глаз устроит, если циклические частоты с обычными не путать, я так думаю :-).
—
-=ВН=-
(17.01.2008 15:05:7
193.125.71.140
, 1126 байт)
И чтоб Вы больше не мучались:-)
—
-=ВН=-
(17.01.2008 15:34:14
193.125.71.140
, 813 байт)
Благодарю, пересчитал (мне легче было поменять l=2*pi*ctg(pi/Qa)). Резонанс...{+}
—
visual_wind
(17.01.2008 15:31:45
217.12.205.67
, 389 байт)
Вынесите усиление за скобки, в который уже раз пишу, ниже все Вам посчитал:-)))
—
-=ВН=-
(17.01.2008 15:36:47
193.125.71.140
,
пустое
)
Если я вынесу коэффициент за скобки, то влезу в 32 разряда, однако, как мне...{+}
—
visual_wind
(17.01.2008 15:41:2
217.12.205.67
, 357 байт)
С какой радости мой фильтр и с какой радости первого порядка - это Ваш фильтр:-)) Это все зависит от представлений данных, считайте, что входные данные у Вас 32-х разрядные
—
-=ВН=-
(17.01.2008 15:54:51
193.125.71.140
, 541 байт)
Т.е. выше
—
-=ВН=-
(17.01.2008 15:37:16
193.125.71.140
,
пустое
)
Построил я фильтр, который вы предложили, только исправил знак у C2 на прот...{+}
—
visual_wind
(17.01.2008 15:59:59
217.12.205.67
, 422 байт)
Не надо знак у C2 исправлять:-)
—
-=ВН=-
(17.01.2008 16:01:47
193.125.71.140
,
пустое
)
Так у него тогда ФЧХ от 180 до 270 градусов поднимается. Где тут "инте...{+}
—
visual_wind
(17.01.2008 16:04:16
217.12.205.67
, 165 байт)
Ничего там не поднимается. Вам просто нужно было посчитать точные значения к-тов, я же Вам дал все формулы, а численные значения привел для ориентировки. Не ленитесь, посчитайте.:-)
—
-=ВН=-
(17.01.2008 16:27:2
193.125.71.140
,
пустое
)
Построил. Получился практически один к одному обычный интегратор, во всяком...{+}
—
visual_wind
(17.01.2008 16:59:3
217.12.205.67
, 373 байт)
Так у и аналоговый фильтр, который Вы привели именно такой, у него тоже будет постоянная составляющая на выходе. Но Вы ошибаетесь, отождествляя с интеграторм.
—
-=ВН=-
(17.01.2008 17:26:10
193.125.71.140
, 693 байт)
xbcnsq==чистый:-)
—
-=ВН=-
(17.01.2008 17:27:37
193.125.71.140
,
пустое
)
=>
—
Oldring
(17.01.2008 17:04:16
91.76.19.190
,
пустое
,
ссылка
)
частота среза для разных типов фильтров определена по-разному, поэтому общего способа рассчета этой частоты по коэффициентам нет.
—
SM
(17.01.2008 14:02:56
80.92.255.53
,
пустое
)
А матлаб умеет находить корни кравнений высоких степеней?
—
visual_wind
(17.01.2008 14:25:13
217.12.205.67
,
пустое
)
да. roots
—
SM
(17.01.2008 15:47:16
80.92.255.53
,
пустое
)
Наверное можно в качестве универсального критерия использовать -3dB если не указано противное.
—
Oldring
(17.01.2008 14:07:4
91.76.19.190
, 222 байт)
Забыл сказать, в моем ФВЧ k=0.4
—
visual_wind
(17.01.2008 13:34:55
217.12.205.67
,
пустое
)
Что это за зверь?
—
Oldring
(17.01.2008 14:12:3
91.76.19.190
,
пустое
)
То есть ФВЧ A(p) = (0.4p+1)/(0.16p+0.4р+1)
—
visual_wind
(17.01.2008 14:24:26
217.12.205.67
,
пустое
)
Что еще нужно для анализа?
—
Oldring
(17.01.2008 14:37:47
91.76.19.190
, 677 байт)
можно еще добавить bode( F )
—
Oldring
(17.01.2008 14:40:25
91.76.19.190
,
пустое
)
А частота среза какая по мнению матлаба?
—
visual_wind
(17.01.2008 15:07:50
217.12.205.67
,
пустое
)
Посмотрите на график.
—
Oldring
(17.01.2008 15:28:56
91.76.19.190
, 106 байт)
А как плясать? Из вашей дискуссии с SM вчера вечером я мало что понял. Вы п...{+}
—
visual_wind
(17.01.2008 15:35:8
217.12.205.67
, 313 байт)
Я предлагаю подумать
—
Oldring
(17.01.2008 16:06:24
91.76.19.190
, 1199 байт)
Правда, с d в знаменателе будут некоторые неприятности - но в случае первого порядка они были легко обходимы
—
Oldring
(17.01.2008 16:36:12
91.76.19.190
, 381 байт)
И, кстати, у прототипа нет в нуле нуля - так что, вероятно, Вы что-то не поняли.
—
Oldring
(17.01.2008 16:42:24
91.76.19.190
, 125 байт)
Использовали. А Вы не могли бы уточнить, что значит нет нуля в нуле?
—
visual_wind
(17.01.2008 17:04:52
217.12.205.67
,
пустое
)
Это значит что у него усиление по DC равно единице, а дальше спадает как у интегратора.
—
Oldring
(17.01.2008 17:05:38
91.76.19.190
,
пустое
)
Первое приближение к прототипу (хуже по ФЧХ) - простая RC-цепь c постоянной времени полсекунды.
—
Oldring
(17.01.2008 17:15:16
91.76.19.190
,
пустое
)
Правильно я понимаю, что усиление по постоянному току, это значение на АЧХ ...{+}
—
visual_wind
(17.01.2008 17:36:9
217.12.205.67
, 178 байт)
Правильно.
—
Oldring
(17.01.2008 17:37:13
91.76.19.190
,
пустое
)
ЧТо же тогда получается, у интегратора это усиление равно бесконечности с ...{+}
—
visual_wind
(17.01.2008 17:43:41
217.12.205.67
, 174 байт)
У чистого интегратора - да. У прототипа оно равно 1.
—
Oldring
(17.01.2008 17:44:45
91.76.19.190
,
пустое
)
Получается, что для нулевого усиления по постоянному току в аналоговой пере...{+}
—
visual_wind
(17.01.2008 18:01:44
217.12.205.67
, 353 байт)
Скажу за себя - вижу по графику ФЧХ. Она хуже :)
—
Oldring
(17.01.2008 18:03:20
91.76.19.190
,
пустое
)
Из этого следует, что Вы не учли механизм устранения смещения нуля в аналоговой системе управления.
—
Oldring
(17.01.2008 17:39:49
91.76.19.190
,
пустое
)
да есть чисто цифровой путь. Сразу сделать фильтр, с двумя компл. сопряженными полюсами около ед. окр. на частоте, соотв. этим +-5 герц. И убрать интегратор. Т.е. (+)
—
SM
(17.01.2008 15:49:51
80.92.255.53
, 110 байт)
Выполнил я указанный фильтр Баттерворта
—
visual_wind
(17.01.2008 18:54:52
217.12.205.67
, 352 байт)
Можете просто добавить ПОСЛЕ него обычный интегратор. Он сам убъет свой ноль.
—
SM
(18.01.2008 14:55:14
80.92.255.53
,
пустое
)
Спасибо за roots. А Вы не могли бы поподробнее про фразу насчет уборки нуля.
—
visual_wind
(17.01.2008 16:02:52
217.12.205.67
,
пустое
)
Он там в явном виде в числителе сидит. Берете и поганой метлой его оттуда физически ликвидируете :)
—
SM
(17.01.2008 17:26:2
80.92.255.53
,
пустое
)
Объясните, пожалуста, как Вы его там видите и почму метла должна быть именн...{+}
—
visual_wind
(17.01.2008 17:38:21
217.12.205.67
, 169 байт)
Как я его там вижу? Смотрю передаточную ф-цию, и вижу в числителе "1-z". А метла поганой просто так быть должна. она нули лучше выметает.
—
SM
(18.01.2008 14:18:48
80.92.255.53
,
пустое
)
т.е. конечно 1-z^-1
—
SM
(18.01.2008 14:19:15
80.92.255.53
,
пустое
)
Только плох этот фильт по ФЧХ
—
Oldring
(17.01.2008 17:32:44
91.76.19.190
, 54 байт)
Совет+
—
andyP
(16.01.2008 20:20:4
89.18.130.241
, 144 байт)
И давно мы начали измерять ширину полосы в линейном масштабе?
—
Oldring
(16.01.2008 20:00:0
91.76.19.190
, 93 байт)
А в каком, пардоньте, еще полосы измерять? Я саму добротность не называл.
—
SM
(16.01.2008 20:10:31
80.92.255.53
,
пустое
)
Какая здесь добротность?
—
Oldring
(16.01.2008 20:16:52
91.76.19.190
, 153 байт)
PS кстати, как в Матлабе из скрипта заменить у этого графика единицы измерения частоты на герцы?
—
Oldring
(16.01.2008 20:22:47
91.76.19.190
,
пустое
)
не в курсе, я всегда только plot использую
—
SM
(16.01.2008 20:28:42
80.92.255.53
,
пустое
)
добротность полюса смысле его близости к ед. окр.
—
SM
(16.01.2008 20:19:7
80.92.255.53
,
пустое
)
Нельзя ли числом?
—
Oldring
(16.01.2008 20:20:57
91.76.19.190
,
пустое
)
Нельзя, я в уме такое не умею сосчитать. А опыта синтеза ФВЧ с такими резкими спадами в ноль знаю, какие там разрядности и коэффициенты получаются.
—
SM
(16.01.2008 20:28:2
80.92.255.53
,
пустое
)
Так я Вам привел пример фильтра ;)
—
Oldring
(16.01.2008 20:47:42
91.76.19.190
, 300 байт)
Я точно об этом же и говорил, там, где сказал про полосу переходную в 0.000081 и что 7 порядков вовсе не странно, (и как следствие 16 бит маловато).
—
SM
(16.01.2008 20:51:29
80.92.255.53
,
пустое
)
Да нет, 16 бит маловато по другим причинам
—
Oldring
(16.01.2008 20:58:49
91.76.19.190
, 189 байт)
Этого да, не знаем, но я не вижу способа обеспечить ноль в нуле и околоравномерную в 6..доконца герц полосу с только действительными полюсами при биквадратном фильтре. А кстати насчет КИХ - почему бы и нет :) только дюже длинный он будет.
—
SM
(16.01.2008 21:07:29
80.92.255.53
,
пустое
)
Чем приведенный мною пример фильтра не устраивает?
—
Oldring
(16.01.2008 21:10:16
91.76.19.190
,
пустое
)
Тем, что на нулевой частоте у него коэффициент передачи не ноль. И тем, что он не в том домене.
—
SM
(16.01.2008 21:19:57
80.92.255.53
,
пустое
)
Не ноль? Не в том домене? Уверены? ;)
—
Oldring
(16.01.2008 21:22:33
91.76.19.190
,
пустое
)
:)
—
Oldring
(16.01.2008 21:24:30
91.76.19.190
, 272 байт)
Сорри, моя виноват, я не заметил T четвертым параметром zpk. Как этот гнусный объект преобразовать в привычные [Num, Den]? zp2tf не арбайтен.
—
SM
(16.01.2008 21:34:51
80.92.255.53
,
пустое
)
Так не надо его туда преобразовывать.
—
Oldring
(16.01.2008 21:47:23
91.76.19.190
, 394 байт)
не устраивает тем, что мой любимый fvtool резльтат этого zpk не ест.
—
SM
(16.01.2008 21:56:2
80.92.255.53
,
пустое
)
Можете сами приготовить ему корм.
—
Oldring
(16.01.2008 22:05:4
91.76.19.190
, 230 байт)
да уже вроде разобрался, только все равно хня какая-то (+)
—
SM
(16.01.2008 22:09:41
80.92.255.53
, 114 байт)
В нуле нуль. Один. :)
—
Oldring
(16.01.2008 22:22:11
91.76.19.190
, 95 байт)
Да глазами я вижу, что в нуле нуль. А в fvtool этого не вижу. И не пойму, где я нечестно перешел из объекта zpk в [num,den]. Вопрос уже чисто технический.
—
SM
(16.01.2008 22:24:49
80.92.255.53
,
пустое
)
Поставьте логарифмический масштаб, 65536 точек и укажите частоту выборок.
—
Oldring
(16.01.2008 22:30:18
91.76.19.190
,
пустое
)
У меня и так логарифмический масштаб. По Y. А по X не могу, тогда нулевую частоту не увижу, где все должно уходить в минус бесконечноть децибел.
—
SM
(16.01.2008 22:50:54
80.92.255.53
,
пустое
)
Видимо это глюки fvtool
—
Oldring
(16.01.2008 23:03:17
91.76.19.190
, 98 байт)
Я же говорю: не так частоту считаете
—
Oldring
(16.01.2008 21:33:40
91.76.19.190
, 108 байт)
все я так считаю и в том масштабе (+)
—
SM
(16.01.2008 21:39:40
80.92.255.53
, 215 байт)
Еще раз: чем мой фильтр не устраивает?
—
Oldring
(16.01.2008 21:49:14
91.76.19.190
, 104 байт)
не в коэффициентах, а в расстоянии от полюса до ед.окр., отсюда будет следовать большой коэффициент передачи рекурсивной части, и отсюда нехватка разрядности.
—
SM
(16.01.2008 21:53:15
80.92.255.53
,
пустое
)
Расстояние - десять в минус третьей.
—
Oldring
(16.01.2008 21:55:45
91.76.19.190
,
пустое
)
Вообще-то в минус пятой :) А точнее 0.1799E-4, тут я просто допускаю погрешность у автора, он частотную "на глаз гнул"
—
SM
(16.01.2008 21:57:31
80.92.255.53
,
пустое
)
В моем примере - в минус третьей. Фильтр с действительными полюсами, характеристикой интегратора в рабочей полосе и нулем в нуле.
—
Oldring
(16.01.2008 21:59:46
91.76.19.190
, 63 байт)
Это в моем примере. Получено как 1-abs(p). Для фильтра, который ставится ПЕРЕД интегратором, для того, чтобы интегратор (голый полюс в нуле) потом не улетал никуда.
—
SM
(16.01.2008 22:06:21
80.92.255.53
,
пустое
)
Это, конечно, откровение.
—
Oldring
(16.01.2008 22:26:40
91.76.19.190
, 196 байт)
Можно и так, спору нет. Но человек просил "что-то сделать", чтобы его уже имеющийся интегратор никуда не улетал.
—
SM
(16.01.2008 22:29:41
80.92.255.53
,
пустое
)
Это все равно плохое решение
—
Oldring
(16.01.2008 22:32:54
91.76.19.190
, 122 байт)
Должен там быть ноль. Правда это уже от реализации зависит. Если сначала рекурсивная часть фильтра, потом нерекурсивная, потом его интегратор - то полюс интегратора совершенно точно убъет ноль фильтра.
—
SM
(16.01.2008 22:56:40
80.92.255.53
,
пустое
)
Это да. Один дифференциатор убивает один интегратор.
—
Oldring
(16.01.2008 23:01:46
91.76.19.190
, 75 байт)
И если
—
Oldring
(16.01.2008 23:12:57
91.76.19.190
, 99 байт)
И полюса там будут отстоять от единичной окружности если не на минус седьмой порядок, то наверное на минус пятый гарантировано.
—
SM
(16.01.2008 21:43:8
80.92.255.53
,
пустое
)
Подскажите, где я все-таки ошибся в коэффициентах при переводе
—
visual_wind
(16.01.2008 20:03:47
217.12.205.67
, 989 байт)
Я же Вам написал ниже, нормированная лапласовская переменная. Т.е. если
—
-=ВН=-
(16.01.2008 21:55:18
193.125.71.140
, 1186 байт)
:-))Еще описки: C1=-2*R*COS(FI, D2=-DNUL
—
-=ВН=-
(17.01.2008 11:53:51
193.125.71.140
,
пустое
)
Еще описка. Вместо Qa конечно ctg(pi/Qa)
—
-=ВН=-
(17.01.2008 10:22:48
193.125.71.140
,
пустое
)
C0=1, само собой. И C0*y(n)=D0*v(n)+D1*v(n-1)+D2*v(n-2)-C1*y(n-1)-C2*y(n-2). v(n)=B*x(n). x(n) - вход.
—
-=ВН=-
(16.01.2008 22:02:40
193.125.71.140
,
пустое
)
Прежде всего, в выборе названия переменной l ;)
—
Oldring
(16.01.2008 20:18:29
91.76.19.190
,
пустое
)
Отправка ответа
Имя*:
Пароль:
E-mail:
Тема*:
Сообщение:
Ссылка на URL:
URL изображения:
если вы незарегистрированный на форуме пользователь, то
для успешного добавления сообщения заполните поле, как указано ниже:
введите число 47:
Перейти к списку ответов
|
Конференция
|
Раздел "Электроника"
|
Главная страница
|
Карта сайта
Web
telesys.ru