Разработка алгоритма анализа массива ЯЭФП

Вычисление Y-матрицы ЯЭФП, в общем случае, производится по следующему алгоритму. Для каждой из частот, на которой необходимо вычислить параметры структуры, производят действия:

1.

2. Производится нумерация узлов структуры;

3. Заполняется глобальная Y-матрица описывающая весь элемент;

4. Стягиваются узлы, объединенные контактными площадками;

5. Производится перенумерация узлов,

6. Порядок Y-матрицы уменьшается до необходимых размеров.

Глобальная матрица представляется в виде двумерного массива. Количество строк матрицы равно количеству узлов, а количество столбцов задается динамически и зависит от количества связей данного узла С соседними узлами. Примененный алгоритм нумерации узлов позволяет исключить шаги 4 и 5.

Заполнение глобальной Y-матрицы производится по следующему алгоритму. Для каждого конечного элемента (имеющего номер /) выполняем следующие действия:

1. Вычисляем матрицу проводимости конечного элемента

2. Для каждого номера- порядок матрицы проводимости

КЭ, выполняем шаги 3, 4.

Рис. 3.3 Фрагмент схемы замещения двумерного ЯЭФП.

Понижение порядка матрицы производится по следующему алгоритму (рисунок

3.4):

Рис. 3.4 Алгоритм понижения порядка матрицы.

1. Выбирается последний столбец в матрице, обозначим его В. Допустим, что число элементов в столбце п.

2. Элементы столбца, кроме последнего нумеруются: j = l..n-l.

3. Последний элемент столбца обозначим D

4. Если D не равно 0, то для всех j от 1 до п-1 выполняем пункты 5, 6, 7.

5. j-ю строку обозначим А.

6. Для все:выполняем пункт 7.

7.

8. Последний столбец отбрасывается, и итерация повторяется с пункта 1, пока порядок матрицы не уменьшиться до гп, где m - количество контактных площадок.

Рис. 3.5 Схема оаботы поогоаммы анализа ЭЯФП.

Схема работы программы, включающая формирование и вычисление матрицы проводимости двумерного ЯЭФП и вычисление передаточной функции устройства, изображена на рисунке 3.5.

3.3.