Разработка программно-методического комплекса для анализа линейных эквивалентных схем в частотной области для числа узлов <=500

i

i U(t)

j -U(t) Этот вектор почти нулевой

Y

Принцип построения аналогичен матрице С.

Характеристики модели в ОКБ.

Достоинства:

- Метод построения прост, обладает низкой трудоемкостью.

- Матрицы, как правило, хорошо обусловлены, результатом чего является высокая точность решения.

Недостатки:

- Используется только один вид зависимых источников.

- Наличие интегральных уравнений.

2) Построение модели в РОКБ с помощью ММУП.

Цель - избавиться от интегральных уравнений и оставить только дифференциальные уравнения.

1. Записывается модель в ОКБ.

2. Избавляемся от интегральных членов уравнения ( вида 1/pL, т. к. 1/р - оператор интегрирования), преобразовывая их в новые неизвестные (например, токи).

3. Получим систему вида:

Это система линейных обыкновенных дифференциальных уравнений 1-го порядка с постоянными коэффициентами в неявной форме.

4. Решаем полученную систему.

{???????????????????????????????????????????????????????????????????}

Достоинства:

1. В модели могут быть любые типы источников.

2. Низкая трудоемкость (т. к. метод прост).

3. Отсутствуют интегральные уравнения.

Недостатки:

Выросла размерность решаемых задач.

3) Построение модели в СГКБ с помощью МПС

МПС сложен для осмысления и для реализации. МПС можно построить, если в схеме нет топологических выражений (это контуры из емкостей или звезды из индуктивностей).

Чтобы выйти из этой ситуации, в схему вводят дополнительные элементы, но снижается точность вычислений.

Вывод: модели СГКБ имеют смысл, когда <= 100, где и - собственные значения матрицы (А- Е).

Определение квазистатических (частотных) характеристик линейных эквивалентных схем.

{??????????????????????????????????????????????????????????????????????}

Для большинства линейных схем характерными являются такие показатели, как добротность, полоса пропускания, равномерность усиления в некотором частотном диапазоне и другие, определяемые по АЧХ и ФЧХ.

Основными широко применяемыми при “ручных” расчетах схем являются методы операционного исчисления, и в частности, спектральный (частотный) метод Фурье.

С помощью преобразований Лапласа решения системы линейных дифф. уравнений переводятся в область комплексной переменной p=Y+jw, показываемой комплексной частотой.

Функция от t, к которой применено преобразование Лапласа, называется оригиналом, а соответствующая функция от р - изображением. Связь между ними определяется формулами:

Основная цель этих преобразований - сведение дифференциальных уравнений к чисто алгебраическим относительно комплексной частоты р. Так, при нулевых начальных условиях операция дифференцирования соответствует умножению на р-изображение, следовательно, при х= уравнения системы:

х = Ах + f(t) х = х

х(t) - вектор переменных состояния,

А - матрица размерностью n x n,

х - вектор начальных значений

будут иметь вид: р Х(р) = А Х(р) - F(р),

а решение исходной системы вида:

х(t) = e x + e f(S) dS, где е = (матричная экспонента)

будет иметь вид:

Х(р) = (рЕ - А) * F(p) = K(p) F(p)

Так как выходные токи и напряжения линейным образом выражаются через переменные состояния и входные воздействия, то вектор выходных переменных z = Bx + Cf , где В, С - матрицы. Тогда матрица В(рЕ - А) + С соответствует матричной передаточной функции, обозначаемой обычно К(р). Отношения любых переменных вектора неизвестных называются схемными функциями. Численный расчет или формирование аналитических выражений для схемных функций составляют основу задачи анализа линейных эквив. схем в частотной области. Согласно правилам Крамера, эти функции описываются линейной комбинацией отношений алгебраических дополнений матрицы А. Таким образом, в общем случае схемные функции есть дробно-рациональные выражения относительно комплексной частоты. Форма их представления называется символьной (буквенной), если коэффициенты при различных степенях р определены через параметры элементов схемы. Если коэффициенты получены в численном виде, то такую форму представления принято называть символьно-численной или аналитической.

К достоинствам методов определения схемных функций на ЭВМ можно отнести: получение конечного результата анализа в аналитическом виде; возможность быстрого дальнейшего расчета значений схемных функций на заданных частотах; удобство при решении задачи оптимизации и определения устойчивости схемы.

К недостаткам при решении задачи на ЭВМ можно отнести: огромный порядок (до нескольких десятков) полиномов схемных функций, диапазон изменения коэффициентов полиномов может превышать возможности представления чисел в разрядной сетке ЭВМ, что требует проведения соответствующей нормировки и счета с удвоенной точностью. Это объясняется влиянием всех элементов схемы во всем частотном диапазоне.

Вывод: используя метод оределения схемных функций, можно достичь в приемлемое время результатов для схем небольших размерностей.

Наряду с методами символьного анализа существуют методы численных решений или расчета тех же схемных функций по точкам. Целью анализа в том случае является получение набора численных значений схемных функций на заданных частотах путем многократного решения системы линейных алгебраических уравнений с комплексными коэффициентами. В процессе расчета необходимо учитывать разреженность матрицы и оптимальный порядок исключения переменных. Алгоритмы численных методов расчета схемных функций, как правило, легче реализуются на ЭВМ и требуют меньших объемов машинной памяти и используются при этом для расчета достаточно больших схем , имея при этом удовлетворительную погрешность и приемлемое время.

Численный метод.

Идея: Выбирается диапазон частот, для каждого значения частоты решают комплексное уравнение.

Достоинства и недостатки метода:

1. Можно работать с переменным шагом частоты. Чем сильнее меняются характеристики, тем меньше шаг, это может привести к огромному количеству шагов.

2. Трудоемкость линейно зависит от количества шагов.

Линейно-аналитический метод.

Идея метода: Определить выходные характеристики в аналитическом виде (т. е. как функция от р, где р - буква). Далее вместо р подставлять конкретное значение частоты и получать характеристики.

А х = Y ; [A ... A ] х = Y

Будем считать, что в схеме имеется единственный источник входных сигналов.

[ Cp + G ] x = Y - исходная модель

где А - большая матрица, в которой вычитаем строку и столбец,

А - алгебраическое дополнение, оставшееся после вычитания строки и столбца,

- говорят о номерах вычеркнутых строк и столбцов, многочлен имеет ровно столько корней, какова его степень. корни могут быть вещественными и/или комплексно сопряженными.

{ - константы = к,

z ,... ,z - нули,

р ,... ,р - полюсы,

к уровень ??? }

- формула вычисления частотных характеристик

Достоинства и недостатки:

- Нули и полюсы заранее известны по виду функции (больше полезной информации).

- Точное решение многочлена высокой степени (>4) не может быть получено, а вычисление значений многочлена степени >30 приведет к погрешности >50%.

- Нули и полюсы вычисляются как собственные значения матриц (числителя и знаменателя).

- Трудоемкость этой задачи 2 * n (n - порядок матрицы), и 4/3 * n - для вычислений в одной точке по частоте.

Вывод: применяется для задач малой размерности.

1. Обзор методов

Цель метода:

1. Составляем (или уже имеем) эквив. схему.

Эквив. схема отображает: способ связи элементов друг с другом, физическая сущность отдельных элементов, граф же только - способ связи.

Введем правила построения эквив. схем:

1) Эквив. схема, как и граф, состоит из множества ветвей и узлов.

2. Каждая ветвь относится к одному из 5-ти возможных типов:

а. б. в. г. д. е. ж. з.

II IU UU

3) Каждой ветви соответствует компонентное уравнение:

а.

dU

I=C*

dt

I, U - фазовые переменные типа потока и разности потенциалов (напряжения) в рассматриваемой ветви, С - емкость.

б.

dI

U=L*

dt

L - индуктивность

в.

U=R*I

R - сопротивление

г.

U=f1(V,t)

U - вектор фазовых переменных,

t - время, в частном случае возможное U=const

д.

I=f2(V,t)

U - вектор фазовых переменых,

I - м.б. I=const

Зависимая ветвь - ветвь, параметр которой зависит от фазовых переменных.

4) Каждому узлу схемы соответствует определенное значение фазовой переменной типа потенциала, каждой ветви - значения переменных I и U, фигурирующих в компонентных уравнениях. Соединение ветвей друг с другом (т.е. образование узлов) должно отражать взаимодействие элементов в системе. Выполнение этого условия обеспечивает справедливость топологических уравнений для узлов и контуров.

В качестве фазовых переменных нужно выбирать такие величины, с помощью которых можно описывать состояния физических систем в виде топологических и компонентных уравнений.

В ЭВМ эта схема представляется в табличном виде на внутреннем языке.

Граф электрич. схем характеризуется некоторыми так называемыми топологическими мат-рицами, элементами которых являются (1, 0, -1). С помощью них можно написать независимую систему уравнений относительно токов и напряжений ветвей на основании законов Кирхгофа. Соединения ветвей с узлами описываются матрицей инциденции А . Число ее строк равно числу узлов L, а число столбцов - числу ветвей b. Каждый элемент матрицы a(i, j):

м -1 - i-я ветвь входит в j-й узел,

a(i, j) = н 1 - i-я ветвь выходит из j-го узла,