Главная страница

Реферат Схемотехническое моделирование аналоговых устройств. Титульный лист Содержание


Скачать 0.51 Mb.
НазваниеТитульный лист Содержание
Дата15.02.2020
Размер0.51 Mb.
Формат файлаdoc
Имя файлаРеферат Схемотехническое моделирование аналоговых устройств.doc
ТипРеферат
#108570
страница4 из 10
1   2   3   4   5   6   7   8   9   10
Численные методы решения системы обыкновенных дифференциальных уравнений при моделировании переходных процессов

Задача расчета переходного процесса в схеме математически формулируется как задача Коши, т.е. ищется решение x(t) на интервале от 0 до t при заданных начальных условиях x(t0)=x0. Искомая функция x(t) определяется в отдельных точках (узлах) t1,t2,…,tn, на которые разбивается заданный интервал времени. Расстояние между узлами ∆t=tn+1–tn называется шагом интегрирования h, который определяется в процессе вычислений или может быть задан [4].

Таким образом, процесс решения задачи сводится к нахождению для очередного момента времени tn+1 значения xn+1=x(tn+1) по уже известным значениям узловых напряжений (или токов) и их производных в предыдущих временных узлах.

Методами решения уравнения является явный метод Эйлера: xn+1 = xn + hf(xn, tn) и неявный метод Эйлера: x n+1 = x n + hf(xn+1, tn+1).

Пример. Рассчитать переходной процесс в цепи, изображённой на рисунке 5, с помощью прямой формулы Эйлера.

В такой цепи происходит разряд конденсатора с ёмкостью C = 1 мкФ, заряженного до напряжения 1 В, через резистор с сопротивлением R = 1 кОм.



Рисунок 5

Если принять uc(t)=x(t), x(0)=1В и постоянную времени τ = RC = 10-3 сек, то процессы в этой цепи будут описываться дифференциальным уравнением



Если решать это уравнение с помощью прямой формулы Эйлера и выбрать h = 200 мкс с начальным условием x(0)=1, то

x(200 мкс)=x1=1-200*10-6*1000*1=0,8

x(400 мкс)=x2=0,8-200*10-6*1000*0,8=0,64 и т.д.

На практике модель схемы чаще представляется в неявной форме

. (*)

При этом модель формируется на основе замены дифференциальных и интегральных уравнений на соответствующие им конечно-разностные уравнения с помощью дискретизации компонентных уравнений емкости и индуктивности:

.

Таким образом, от системы (*) переходим к системе конечно-разностных алгебраических уравнений

Fi(xn+1,x,…,xn-k)=0, i=1,…,p.

Составленная модель решается относительно xn+1 каким-либо численным методом, обычно методом Ньютона.
  1. 1   2   3   4   5   6   7   8   9   10


написать администратору сайта