
- •1 Цель работы
- •2 Краткие теоретические сведения
- •3 Составление блок-схем алгоритмов программ (БСА)
- •3.1 Описание символов, применяемых в БСА программ
- •3.2 Правила применения символов и выполнения схем
- •3.3 Примеры выполнения схем
- •4 Программа работы
- •Лабораторная работа №2 «Решение систем линейных уравнений»
- •1 Цель работы
- •2 Краткие теоретические сведения
- •2.1 Метод Крамера
- •2.2 Метод исключений Гаусса
- •2.3 Метод LU–разложения (LU–факторизации)
- •3 Порядок работы
- •4 Методические указания к выполнению работы
- •5 Содержание отчета о лабораторной работе
- •1 Цель работы
- •2 Краткие теоретические сведения
- •Листинг программы интегрирования СДУ (2.16)
- •3 Порядок работы
- •4 Методические указания к выполнению работы
- •5 Содержание отчета о лабораторной работе
- •Литература
3 Порядок работы
3.1Ознакомиться с методами численного интегрирования систем дифференциальных уравнений первого порядка, а также с методами формирования уравнения цепи в виде системы дифференциальных уравнений первого порядка.
3.2Для указанной преподавателем схемы составить математическую модель в виде системы дифференциальных уравнений первого порядка. Рассчитать начальные условия (напряжение конденсаторов и токи катушек до коммутации). Схемы для анализа приведены на рисунке 2.
3.3Составить блок-схему алгоритма программы и программу для численного интегрирования уравнения состояния методом РунгеКутта четвёртого порядка. Значения параметров элементов схемы должны вводиться с клавиатуры.
3.4Представить программу для контроля преподавателю.
3.5Произвести расчет переходного процесса для произвольных значений параметров элементов и записать результаты (для внесения в отчёт о лабораторной работе).
3.6Оформить отчет о лабораторной работе (содержание отчета приводится в п.5).
4 Методические указания к выполнению работы
Схемы для расчета приведены на следующем рисунке. Необходимо выбрать схему, указанную преподавателем.
50.

L L
|
|
R2 |
R1 |
E |
|
|
|
S |
|
S |
|
|
|
|
C |
E |
C |
R2 |
|
|
R1
|
1 |
|
|
2 |
|
|
|
|
|
L |
|
R1 |
R2 |
|
R1 |
S |
|
S |
|
|
|
|
|
E |
L |
C |
E |
R2 |
C |
|
|
3 |
4 |
Рис. 2. Варианты схем для расчета
В качестве примера покажем вывод основных расчетных формул для анализа переходного процесса в схеме, представленной на рисунке 3.
1 |
2 |
|
|
E |
S |
|
|
|
|
R1 |
L |
|
C |
|
|
|
|
|
0 |
3 |
|
R2 |
Рис. 3. К примеру расчета
Начальные условия для анализа переходного процесса определяем при разомкнутом состоянии ключа S. Очевидно,
51.

UC 0 |
= E |
(4.1) |
|
IL0 |
= 0 |
||
|
Далее, в предположении, что ключ замкнут, составим направленный граф цепи и выделим на нём нормальное дерево (E, C, R2).
1 |
2 |
3 |
|
E |
L |
|
|
C |
|
R1 |
R2 |
0
Рис. 4. Направленный граф для схемы примера
Топологическая матрица F будет иметь вид:
|
|
|
E |
C |
R2 |
F = |
R |
−1 |
1 |
0 |
|
1 |
|
0 |
−1 |
|
|
|
L |
|
1 |
Топологическое уравнение (2.1) для данного примера будет иметь вид:
IE |
|
|
0 0 0 −1 0 |
|
UE |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
IC |
|
|
0 0 0 1 |
−1 |
UC |
|
||||||||||||
I2 |
|
= 0 |
0 |
|
|
0 |
0 |
1 |
|
U2 |
|
(4.2) |
||||||
|
|
|
|
−1 0 0 0 |
|
|
I1 |
|
|
|||||||||
U1 |
|
|
1 |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
−1 0 0 |
|
|
IL |
|
|
||||||
UL |
|
0 1 |
|
|
|
|
|
|||||||||||
В правую часть этого уравнения введём компонентные зависимости ре- |
||||||||||||||||||
зистивных элементов и источников: |
|
|
|
|
|
|
|
|
|
|
|
|||||||
U |
E |
= E; |
U |
2 |
= I |
R ; |
I |
1 |
= |
U1 |
, |
|
||||||
|
|
|||||||||||||||||
|
|
|
|
|
|
|
2 2 |
|
|
|
R1 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
а в левую часть — компонентные зависимости реактивных элементов: |
|
|||||||||||||||||
|
IC |
= C |
dUC |
|
=CUC ; |
UL = LIL . |
|
|||||||||||
|
|
|
||||||||||||||||
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
В результате будем иметь уравнение:
52.

|
IE |
|
0 0 0 |
−1 0 |
|
|
E |
|
|
||||
|
|
|
|
|
|
|
1 |
|
|
UC |
|
|
|
CUC |
0 0 0 |
−1 |
|
(4.3) |
|||||||||
|
I |
|
|
= 0 |
0 |
0 |
0 |
1 |
|
I |
R |
|
|
|
|
2 |
|
|
−1 0 |
0 0 |
|
|
2 2 |
|
|
||
|
U1 |
1 |
U1 R1 |
|
|
||||||||
|
|
|
|
|
|
−1 0 0 |
|
|
IL |
|
|
||
|
LIL |
0 1 |
|
|
|
|
Далее при ручных расчетах удобнее всего расписать систему (4.3) по строкам в виде отдельных уравнений:
|
IE |
= −U1 R1 |
|
|
|||||||
CU |
C |
=U |
1 |
R |
− I |
L |
|||||
|
|
|
1 |
|
|
||||||
|
I2 |
= IL |
|
|
|
|
(4.4) |
||||
|
|
|
|
|
|||||||
|
U |
|
= E −U |
|
|
|
|
||||
|
LI |
|
1 |
=U |
|
|
C |
R |
|
||
|
L |
C |
− I |
|
|||||||
|
|
|
|
|
2 |
|
2 |
За основу уравнения состояния выберем второе и пятое уравнения системы (4.4), поскольку они содержат производные переменных состояния —
UC и IL . Из правых частей этих уравнений исключим «лишние» переменные
(должны остаться только переменные состояния UC, IL и независимые источники E). Для этого воспользуемся оставшимися уравнениями системы (4.4):
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
CUC = (E −UC ) R1 − IL |
|
|||||||||||||||
|
|
|
|
|
=U |
|
− I |
|
|
R |
|
|
|
||||
|
LI |
L |
C |
L |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|||||
Оставляя в левой части только производные, запишем |
|||||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
1 |
|
|
||
UC = − |
|
|
|
UC |
− |
|
|
|
IL + |
|
|
E |
|||||
R1C |
|
C |
R1C |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
(4.5) |
|||||||
|
|
1 |
|
|
|
|
|
R2 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
IL = |
|
UC |
− |
|
|
IL |
|
|
|
|
|
||||||
L |
L |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Это и есть уравнение состояния исследуемой цепи. Для его интегрирования можно использовать программу, листинг которой приведён выше, дополнив её модулем ввода с клавиатуры численных значений параметров элементов.
53.