Решение дифференциального уравнения можно записать в виде:
— корни характеристического уравнения. Они совпадают с собственными числами матрицы А, следовательно могут быть найдены из уравнения:
, где Е — единичная матрица.
Найдём корни полученного квадратного уравнения.
Подставим полученные корни уравнения, значения тока на катушке и напряжение на конденсаторе , найденные ранее, в решение дифференциального уравнения.
Найдем производную от полученной системы.
Используя начальные условия запишем:
Перепишем эту систему:
;
Решим систему этих уравнений:
Из последнего уравнения системы выразим , подставим его во второе уравнение системы, получим:
Из первого уравнения системы выразим , подставим это выражение в третье уравнение системы, получим:
Запишем окончательно решение дифференциального уравнения:
Определим:
Найдём время переходного процесса:
Решение уравнений состояния с помощью матричной экспоненты.
Запишем систему дифференциальных уравнений.
Запишем систему в виде:
, где ; ; , т.е.:
Запишем переменные состояния в виде:
, где
— матричная экспонента.
Следовательно:
Числа , входящие в матричную экспоненту — корни характеристического уравнения. Они совпадают с собственными числами матрицы А, следовательно могут быть найдены из уравнения:
, где Е — единичная матрица.
Эти корни были найдены в решении уравнений состояния классическим методом, и они равны: .
Подставим значения корней характеристического уравнения и элементы матрицы А в матричную экспоненту, получим:
, тогда
Откуда находим:
Преобразуем цепь после коммутации согласно правилам для операторного метода.
Запишем уравнения Кирхгофа для данной цепи:
Определение тока :
Подставим третье уравнение системы во второе, и выразим ток :
Подставим найденный ток в третье уравнение системы и выразим из него ток :
Ток подставим в первое уравнение системы и выразим из него ток :
Перепишем это уравнение относительно s:
Найдём корни уравнения
Подставим в уравнение численные значения параметров цепи:
; ;
Найдем производную :
Запишем выражение для числителя:
Подставим в уравнение численные значения параметров цепи:
k |
||||
1 |
-250 |
-25 |
-300 |
0.08 |
2 |
-625 |
-100 |
750 |
-0.13 |
3 |
0 |
150 |
500 |
0.3 |
Запишем окончательно по теореме разложения ток :
Определение напряжения :
Для определения напряжения определим вначале ток , а потом воспользуемся компонентным уравнением .
Из первого уравнения системы выразим ток : .
Из второго уравнения системы выразим ток : .
Из третьего уравнения системы выразим ток :. В полученное выражение подставим токи и :
, отсюда можем выразить ток через параметры цепи:
Перепишем числитель этого выражения:
Подставим в уравнение численные значения параметров цепи:
Знаменатель H(s) полученного выражения совпадает со знаменателем в выражении для тока при его определении в предыдущем пункте, поэтому, определённые ранее, корни H(s) и значения его производных можно сразу записать в таблицу значений.
k |
||||
1 |
-250 |
25 |
-300 |
-0.08 |
2 |
-625 |
25 |
750 |
0.03 |
3 |
0 |
150 |
500 |
0.3 |
Запишем по теореме разложения ток :
Численное интегрирование уравнений состояния.
Запишем уравнения состояния:
,
Проведём численное интегрирование уравнений состояния в среде MATLAB 6.1. Для этого создадим m-файл для вычисления значений правых частей (назовём его fun.m). Запишем в этот файл правую часть системы дифференциальных уравнений:
function dydt=fun(t,y)
dydt=[-750*y(1)-2.5*y(2)+375;25000*y(1)-125*y(2)];
где ,.
Сохраним этот файл в каталоге к которому имеет доступ MATLAB. В рабочем окне MATLAB обратимся к процедуре численного интегрирования, сохраняющую результаты интегрирования в рабочем пространстве в виде переменных:
[t,y]=ode45(@fun,[0 0.02],[0.25 50]);
При этом значения t будут выведены в виде вектора-столбца, а y-в виде матрицы.
Теперь для построения графиков переменных состояния можно воспользоваться командой plot:
для построения графика — plot(t,y(:,1))
для построения графика — plot(t,y(:,2)).
|
t
0
0.0002
0.0004
0.0006
0.0008
0.0012
0.0016
0.0019
0.0023
0.0027
0.0031
0.0035
0.0039
0.0044
0.0049
0.0054
0.0059
0.0064
0.0069
0.0074
0.0079
0.0084
0.0089
0.0094
0.0099
0.0104
0.0109
0.0114
0.0119
0.0124
0.0129
0.0134
0.0139
0.0144
0.0149
0.0154
0.0159
0.0164
0.0169
0.0174
0.0179
0.0184
0.0190
0.0195
0.0200
y(1) y(2)
0.2500 50.0000
0.2617 50.0297
0.2717 50.1122
0.2802 50.2388
0.2875 50.4014
0.2983 50.7765
0.3061 51.2195
0.3115 51.7055
0.3152 52.2110
0.3178 52.7643
0.3193 53.3100
0.3198 53.8402
0.3197 54.3462
0.3192 54.9152
0.3182 55.4412
0.3170 55.9248
0.3157 56.3645
0.3143 56.7661
0.3130 57.1274
0.3118 57.4518
0.3106 57.7417
0.3095 57.9997
0.3085 58.2294
0.3075 58.4337
0.3067 58.6150
0.3060 58.7756
0.3053 58.9179
0.3047 59.0440
0.3042 59.1555
0.3037 59.2541
0.3033 59.3413
0.3029 59.4184
0.3026 59.4865
0.3023 59.5467
0.3020 59.5998
0.3018 59.6467
0.3016 59.6882
0.3014 59.7248
0.3012 59.7571
0.3011 59.7856
0.3009 59.8108
0.3008 59.8339
0.3007 59.8542
0.3006 59.8721
0.3006 59.8877
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.