Курсовая Проектирование оптимальных структур активных RC-фильтров
Работа добавлена на сайт bukvasha.net: 2015-10-25Поможем написать учебную работу
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
от 25%
договор
1. Базовая структура нестационарных устройств
Важным классом современных устройств автоматики, технической диагностики и техники специальных измерений являются нестационарные блоки и подсистемы, обеспечивающие обработку сигналов в реальном масштабе времени. Такие наблюдатели, оцениватели и фильтры достаточно часто строятся на базе сигнальных процессоров и воспроизводят непосредственно систему дифференциальных уравнений, вытекающую из основных процессов. Внедрение в инженерную практику объектов микросистемной техники, создание новых версий систем управления, обеспечивающих работоспособность объектов в критических ситуациях, связано с разработкой нового поколения высокоточных и экономичных нестационарных устройств.
Результаты показывают, что перестраиваемые ARC-фильтры с собственной компенсацией обеспечивают существенное повышение точности преобразования сигнала при невысоких требованиях к частоте единичного усиления активных элементов и, следовательно, низкой потребляемой мощности. Таким образом, объединение в единую систему экономичного микроконтроллера с развитой архитектурой портов ввода/вывода и ARC-фильтра с цифроуправляемыми параметрами теоретически позволяет решить сформулированную задачу (рис. 1). В приведенной структуре микроконтроллер вырабатывает управляющие воздействия на ARC-схему с цифроуправляемыми параметрами и контролирует процесс оценки координат объекта. В частном случае, когда управляющие слова могут быть определены заранее, многоканальный АЦП может отсутствовать, и микроконтроллер работает в режиме логического управления под действием внутреннего таймера.
С точки зрения принципа обработки входных аналоговых сигналов такую систему уместно назвать гибридной.
Методики, разработанные для исследований и анализа процессов, протекающих в нестационарных цепях, как известно, достаточно сложные. В общем случае для анализа применяют двухмерные преобразования Лапласа, теорию уравнений Хилла, а также различные спектральные методы. Достаточно важными результатами в области анализа линейных нестационарных систем, получившими распространение в радиотехнических цепях, являются работы Л.А. Заде. Необходимо отметить, что в радиотехнике большое развитие получили в основном методы анализа нестационарных цепей с периодически изменяющимися коэффициентами. В случае анализа нестационарных систем с непериодическими параметрами обычно пользуются приближенными методами и оценками. По аналогии с линейными стационарными системами в нестационарных определяющими являются понятия сопряженной импульсной переходной характеристики и параметрической передаточной функции. Первая показывает реакцию предварительно невозбужденной («пустой») нестационарной системы в момент времени приложения единичного импульса и является функцией двух переменных – . Вторая является изображением по Лапласу отношения выходной реакции системы к ее входу в момент времени t и также является функцией двух переменных – .
Рис. 1. Структура нестационарных ARC-устройств
Ограничимся линейной версией системы, когда любое воздействие может быть пересчитано к одному из ее входов. Следовательно, система линейных уравнений n-го порядка, описывающая нестационарное устройство, может быть представлена в виде следующего дифференциального уравнения с нестационарными коэффициентами:
,(1)
где m £ n; – одна из выходных координат рассматриваемой системы; – приведенное эквивалентное входное воздействие, учитывающее скалярные сигналы аналоговой части, которое в символической операторной форме можно представить следующим образом:
, (2)
где ;
где ;
;
.
Как показано в [2], с учетом (2) параметрическую передаточную функцию можно определить из решения следующего дифференциального уравнения:
(3)
где , .
Следуя [2], приведем методику приближенного определения параметрической передаточной функции (3), идея которой принадлежит Л.А. Заде. Перепишем (3) в следующем виде:
, (4)
где .
Решение отыскивается в виде ряда
, (5)
где – «замороженная» передаточная функция ФКБ, .
В работе [5] отмечается, что ряд (5) хорошо сходится только в случае медленно изменяющихся коэффициентов и . Однако быструю сходимость ряда можно обеспечить, выполняя построение функции на ограниченных интервалах времени, где функциональные зависимости и можно аппроксимировать полиномами низкого порядка.
Создание нестационарных устройств в рамках экономичных и быстродействующих структур предполагает обработку сигнала аналоговым способом с применением для указанной цели ARC-схемы с цифроуправляемыми параметрами. В такой постановке задачи необходимо говорить о дискретно-непрерывной фильтрации. Очевидно, что максимального приближения характера поведения дискретно-непрерывной и непрерывной систем можно добиться, уменьшая интервал дискретизации, верхняя граница которого может быть определена из теоремы Котельникова, однако при этом возрастают требования предъявляемые к производительности цифровой части устройства.
Необходимо отметить, что параметры цифроуправляемой ARC-схе-мы на интервале времени, определяемом частотой дискретизации, остаются постоянными, то есть параметрическая передаточная функция (5) совпадает с передаточной функцией стационарной схемы при замораживании в ней на -м шаге коэффициентов
, (6)
где ;
(7)
– обозначает номер интервала дискретизации.
Как видно из приведенного анализа, задача синтеза аналоговой части нестационарной системы сводится к построению такой ARC-схемы, которая бы обеспечила на каждом -м шаге максимальное приближение к идеализированной замороженной передаточной функции (6). Следовательно, проектирование аналоговой части устройства возможно выполнить в рамках известных методов синтеза стационарных ARC-цепей, включая и частотные методы, рассмотренные ранее.
2. Обобщенный алгоритм решения задачи синтеза структур нестационарных ARC-схем
Полученный результат показывает, что задача синтеза структур нестационарных устройств сводится к аналогичной стационарной задаче в точке «наихудшего случая», когда совокупность управляющих параметров из множества допустимых параметрических воздействий приводит к максимальному отклонению частотных характеристик от желаемых. Таким образом, согласно предложенной в настоящей работе методологии синтеза структур рассматриваемую задачу можно разделить на три относительно самостоятельных этапа.
Первый этап заключается в синтезе исходной принципиальной схемы, получении набора локальных передаточных функций, определяющих функции активной составляющей чувствительности, и принятии решения о направлении проектных процедур. Настоящий этап состоит из ряда составляющих. Прежде всего, по модели нестационарного устройства синтезируется стартовая конфигурация принципиальной схемы. По стартовой конфигурации путем коммутации базисных структур строится принципиальная схема аналоговой части проектируемого устройства, воспроизводящая заданный набор «замороженных» передаточных функций. Выбор числа разрядов умножающих ЦАП, входящих в состав управляемых усилителей и интеграторов, может осуществляться по следующей оценочной формуле:
,
где d1 и d2 – верхняя и нижняя границы диапазона измеряемой величины; D – шаг квантования по уровню, который выбирается из соображений точности реализации требуемых коэффициентов.
Для определения набора локальных передаточных функций Fsi(p), Fkj(p), Hi(p), Hj(p), Fii(p), Fjj(p) по синтезированной схеме достаточно вычислить обратную матрицу. Получение последней в символьном виде позволяет не только повысить наглядность представляемой информации, но и обеспечивает на последнем этапе синтеза согласованных с Fii(p), Fjj(p), Hi(p), Hj(p) законов изменения дополнительных компенсирующих цепей обратных связей. На этом же этапе становится возможным вычисление коэффициентов , определяющих верхний уровень динамического диапазона во всех стационарных точках x. Этап завершается определением функций чувствительности к площади усиления всех активных элементов.
На втором этапе синтеза, с целью выбора предпочтительного варианта реализации компенсирующих контуров обратных связей, необходимо определить доминирующие активные элементы, параметры которых наибольшим образом оказывают влияние на достижимый частотный и динамический диапазон схемы. Для ранжирования степеней влияния каждого ОУ наиболее целесообразно, с точки зрения рассматриваемой концепции синтеза, произвести исследование наборов модулей функций чувствительности с целью определения их максимума. Для этого прежде всего необходимо определить область изменения параметров схемы, соответствующую «наихудшему случаю», когда отклонение реализуемых функций амплитудно-частотных характеристик (АЧХ) и фазочастотных характеристик (ФЧХ) от идеальных в полосе рабочих частот схемы окажется максимальным:
, (8)
(9)
где – максимальная граничная частота работы схемы.
Значение оценки верхней границы частотного диапазона схемы для можно определить по формуле
, (10)
где – свободный член полинома знаменателя идеализированной замороженной передаточной функции W0.
Отметим, что важно не только определить значение экстремума функций (8) и (9), но и найти координаты указанных глобальных экстремумов. На эти экстремальные задачи накладывается система ограничений в виде неравенств, следующая из максимально и минимально возможного коэффициента передачи ЦАП, масштабных усилителей, значения постоянной времени интеграторов и рабочего диапазона частот нестационарной схемы
(11)
В частном случае из (11) могут быть исключены ограничения, соответствующие неизменяемым параметрам коэффициентов передачи масштабных усилителей или постоянным времени интеграторов.
Таким образом, в результате решения экстремальных задач (8) и (9) с ограничениями (11) становится возможным определение следующего вектора оптимальных координат:
, (12)
соответствующего наихудшему случаю.
С учетом результатов (12) для ранжирования ОУ по степени их влияния находится решение следующей экстремальной задачи с соответствующей системой ограничений:
, (13)
(14)
где – относительный доверительный интервал решения экстремальной задачи.
Указанный доверительный интервал необходим вследствие того, что максимумы модулей функций чувствительности в общем случае не совпадают с определенной оптимальной точкой (12), а лишь находятся в ее окрестности [9].
Выражение (14) является интегральной оценкой, позволяющей произвести качественный анализ влияния площади усиления i-го ОУ на частотные свойства передаточных функций. Успешное решение экстремальных задач (8), (9) и (13) во многом зависит от специфики работы нестационарного устройства, диапазона изменения управляющих параметров и требований, предъявляемых к точности реализации.
Для визуальной оценки степени влияния параметров каждого ОУ по результатам проведенных исследований (13) строится набор диаграмм по каждому из выходов, из которых можно определить доминирующий активный элемент
(15)
Следует отметить, что функциональные зависимости коэффициента передачи ЦАП , масштабного усилителя и постоянной времени интегратора определяются дискретным моментом (интервалом) времени , поэтому функции (8), (9), (13) и соответствующие им системы ограничений (11) и (14), а также вектор оптимальных координат (12) яв-ляются значениями только двух параметров и . Однако для примене-ния оптимизационных методов решения этих экстремальных задач, не свя-занных в общем случае с прямым перебором возможных комбинаций зна-чений параметров оптимизируемых функций, целесообразной оказывается приведенная выше форма представления указанных выражений.
Задача третьего этапа синтеза связана с введением в схему дополни-тельных компенсирующих контуров обратной связи. В качестве нулевого приближения при проверке результатов синтезированных схемных решений целесообразно взять допустимое значение отклонения АЧХ реализуемой реальной схемы с замороженными коэффициентами от идеализированной. Допустимый коридор отклонений АЧХ определяется допустимыми значениями отклонений коэффициентов полиномов числителя и знаменателя замороженной пере-даточной функции (в точке (12)) от идеальных.
Чувствительность модуля передаточной функции Ф к изменению параметра можно представить следующим образом:
.
Представим замороженную передаточную функцию идеализированного устройства в следующем виде:
. (16)
Тогда допустимое значение отклонения АЧХ можно определить по формуле
, (17)
где , ; и – допустимые отклонения значений коэффициентов полинома числителя и знаменателя передаточной функции (16); – приведенная максимальная статическая погрешность умножающих ЦАП.
Значения и определяются допустимыми интервалами изменения элементов вектора отклонений, следующего из решения системной задачи.
Процедура синтеза низкочувствительной схемы заключается во введении в последнюю дополнительных компенсирующих контуров обратной связи и носит итерационный характер [8]:
- по оценкам, полученным на втором этапе синтеза, выбирается доминирующий активный элемент;
- для выбора предпочтительного варианта реализации компенсирующих контуров обратных связей производится поиск необходимого набора локальных передач, поиск производится по строке матрицы ; если необходимых передач нет, то последние ищутся по всей матрице, исключая элементы главной диагонали;
- для уменьшения влияния частотных свойств доминирующего ОУ на достижимый схемой частотный и динамический диапазоны в схему вводятся дополнительные компенсирующие контуры обратной связи и осуществляется ее параметрическая оптимизация;
- с целью проверки качества принятого в предыдущем пункте решения производится численное моделирование синтезируемой схемы в частотной области, например с помощью одного из современных пакетов программ;
- выход из алгоритма производится либо по достижении требуемых качеств проектируемого устройства (если они не удовлетворены, выбирается следующий по установленному ранжиру активный элемент), либо при исчерпывании всех степеней свободы схемы, количество которых определяется числом заземленных входов ОУ, при этом необходимо учитывать, что для обеспечения правильного режима работы схемы по постоянному току хотя бы один вход ОУ должен быть заземлен;
- с целью проверки качества синтезированных схемных решений производится численное моделирование синтезируемой схемы во временной области.
При получении неудовлетворительных результатов (невозможности достижения заданных требований к качеству проектируемого устройства) в исходной схеме, полученной на первом этапе синтеза, с целью перераспределения значений функций чувствительности необходимо выполнить иной выбор параметров базисных структур, после чего повторить приведенный выше алгоритм синтеза компенсирующих контуров обратных связей. Отмеченная ситуация, например, может возникнуть при синтезе компенсирующих контуров обратных связей, когда для достижения требуемого (достаточного) уровня компенсации влияния инерционных свойств i-го активного элемента на параметры схемы в дополнительном контуре обратной связи необходимо обеспечить большое усиление. Указанного можно достичь несколькими способами. В первом случае в схему вводится дополнительный усилитель, во втором в схеме выполняется иной выбор параметров базисных структур, который обеспечивает получение требуемых значений усиления в компенсирующих контурах схемы путем перераспределения усиления между ее функциональными узлами. В отличие от первого способа, второй не требует дополнительных аппаратных затрат.
Задача третьего этапа синтеза в части синтеза схемных решений не может быть полностью формализована – выбор предпочтительного варианта реализации компенсирующих контуров остается за проектировщиком.
Рассматриваемая задача может быть алгоритмизирована в виде некоторой экспертной системы, исходными данными для которой служат полный набор передаточных функций в символьном виде, полученных на первом этапе, и наборы оценок из второго этапа синтеза. Таким образом, в результате решения задачи последнего этапа проектирования возможно получить схемные решения, позволяющие создать схему с собственной компенсацией влияния частотных свойств активных элементов на ее параметры.
3. Пример синтеза структуры аналоговой части циклического фильтра Калмана – Бьюси
Исходными данными для синтеза схемы циклического фильтра Калмана – Бьюси (ФКБ) являются стартовая конфигурация его структурной схемы, значения коэффициентов усиления и времени функционирования на цикле .
Пусть необходимо производить измерения на интервале и значение скорости изменения входного сигнала на входе ФКБ не превышает , а интенсивность сигнала типа белого шума определяется значением , причем . Следуя методике, изложенной в работе [4], зададим начальное значение ковариационной матрицы следующим образом:
, (18)
где – некоторый коэффициент пропорциональности. Тогда значения коэффициентов усиления [4], определяемые решениями матричного уравнения Риккати, будут иметь следующий вид:
, (19)
. (20)
Дополнительные исследования показывают, что оптимальная точность фильтра достигается в случае, если функциональная зависимость этих коэффициентов для безразмерного времени q имеет вид, представленный на рис. 2.
Из уравнения Риккати [2] легко синтезируется исходная принципиальная схема фильтра рис. 3. На приведенной принципиальной схеме в начальном (некомпенсированном) варианте отсутствуют операционный усилитель А9 и резисторы R13–R17, номинал резистора , а неинвертирующие входы ОУ А2, А3 и А7 соединены с землей. При указанных на схеме номиналах резисторов и конденсаторов максимальный коэффициент передачи умножающих ЦАП (ОУ А2 и А3) не превышает единицы. По формуле (7) определяем, что для в обоих каналах ФКБ необходимо использовать 10-разрядные ЦАП, что позволяет воспроизводить характеристики (19) и (20) в каждый момент времени с высокой точностью, то есть фактически непрерывно.
Результаты численного моделирования схемы ФКБ (рис. 3) показывают, что в рассматриваемом случае достаточным является разбиение интервала времени цикла на 100 отсчетов. Таким образом, частота работы ЦАП составляет .
Рис. 3. Принципиальная схема гибридного циклического ФКБ 2-го порядка
Отметим, что в рассматриваемом случае каких-либо формальных строгих процедур определения допустимого интервала отклонения значений коэффициентов усиления нет. Знаменатель замороженной передаточной функции идеализированного ФКБ на -м фиксированном интервале времени следует из (6) и может быть представлен следующим образом:
. (21)
Из представленных на рис. 2 временных зависимостей коэффициентов усиления следует, что на каждом интервале времени () указанный знаменатель является гурвицевым, так как выполняется условие . С учетом частотных свойств ОУ, входящих в состав реального ФКБ, и их идентичности полином (21) можно записать следующим образом:
, (22)
где .
В выражении (22) не учтены все члены, обратно пропорциональные произведениям площадей усиления ОУ, влияние которых на свойства реализуемого ФКБ пренебрежимо мало. Используя результат [9], условие гурвицевости полинома (22) можно представить следующим образом:
. (23)
Учитывая, что , и пренебрегая членами второго порядка малости, неравенство можно записать в виде
. (24)
Как видно из (24), требования к минимально возможному значению площади усиления ОУ в основном определяются максимально возможным значением отношения коэффициентов усиления фильтра и могут быть снижены при компенсации (уменьшении) величины . Анализ неравенства (24) показывает, что в рассматриваемом случае условие гурвицевости полинома (22) при выполнении не зависит от вариаций приращений , и . Поэтому дальнейший синтез схемы будем производить таким образом, чтобы обеспечить минимальное отклонение АЧХ и переходных характеристик реального фильтра от идеализированного. В этом случае допустимое (минимально возможное) значение площади усиления ОУ может быть определено из анализа отклонений временных характеристик реального фильтра от идеализированного.
В соответствии с предложенной методикой определим необходимые для анализа схемы наборы локальных передач , , и . Для этого по синтезированной принципиальной схеме путем сопоставления локальных передач ветвей схемы с ветвями графа обобщенной структуры определим компоненты матриц и векторов, входящих в систему, и составим блочную матрицу , а также определим обратную :
, (25)
. (26)
нестационарный схема фильтр циклический
Из этой же матрицы легко определяется набор локальных передач , , и [4]. Результаты вычислений полиномов числителей указанных функций локальных передач схемы сведены в табл. 1. Знаменатель рассматриваемых передаточных функций определяется выражением (26).
Таблица 1. Наборы локальных передач схемы ФКБ второго порядка (рис. 3)
Вид локальной передачи | Числитель локальной передаточной функции |
1 | 2 |
| , |
| , , |
| , , , , , |
| , , , ,, , |
| , , ,, , , |