Численное моделирование и анализ переходных процессов в электрической цепи

Автор работы: Пользователь скрыл имя, 01 Апреля 2012 в 23:13, курсовая работа

Описание

В ходе данной работы будут проанализированы переходные процессы в электрической цепи. Работа подразделена на 3 части:

I часть работы включает в себя решение системы дифференциальных уравнений, описывающих переходный процесс в электрической цепи.

II часть работы представляет из себя процесс решения задачи аппроксимации.

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

Содержание

Введение.
Постановка задачи, исходные данные.
Реализация поставленной задачи:

3.1. первый этап решения задачи: численная реализация системы дифференциальных уравнений:

3.1.1. алгоритмически, построив блок-схему и программу,

написанную на языке программирования Free Pascal,

используя алгоритм метода Рунге-Кутта ;

3.2. второй этап решения задачи:

3.2.1. решение задачи аппроксимации используя язык программирования Free Pascal;

3.3. третий этап решения задачи: вычисление интеграла, определяющего количество теплоты, выделяемое на резисторе R4 за период времени t1<t<t2:

3.3.1. алгоритмически, построив блок-схему и программу,

написанную на языке программирования Free Pascal,

используя метод трапеций;
Выводы по проделанной работе.

Скачать (181.26 Кб) Сколько стоит заказать работу?

Работа состоит из 1 файл

информатика.docx

— 196.75 Кб (Скачать документ)

Содержание:

Введение.
Постановка задачи, исходные данные.
Реализация поставленной задачи:

3.1. первый этап решения задачи: численная реализация системы дифференциальных уравнений:

3.1.1. алгоритмически, построив блок-схему и программу,

написанную на языке программирования Free Pascal,

используя алгоритм метода Рунге-Кутта ;

3.2. второй этап решения задачи:

3.2.1. решение задачи аппроксимации используя язык программирования Free Pascal;

3.3. третий этап решения задачи: вычисление интеграла, определяющего количество теплоты, выделяемое на резисторе R4 за период времени t1<t<t2:

3.3.1. алгоритмически, построив блок-схему и программу,

написанную на языке программирования Free Pascal,

используя метод трапеций;

Выводы по проделанной работе.

ВВЕДЕНИЕ

В ходе данной работы будут проанализированы переходные процессы в электрической цепи. Работа подразделена на 3 части:

I часть работы включает в себя решение системы дифференциальных уравнений, описывающих переходный процесс в электрической цепи.

II часть работы представляет из себя процесс решения задачи аппроксимации.

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

Постановка задачи, исходные данные

Дана схема электрической цепи:

Параметры элементов цепи:

E=E0*sin(w*t+φ) – гармонический источник тока

E0=15 В – амплитуда колебаний

w=2П*f – циклическая частота

t– текущий момент времени

f=30 Гц – линейная частота, φ=9П/5 – (параметры заданные по вариантам)

R1=30 Ом, R2=25 Ом , R3=50 Ом , R4=1.88 Ом

R5=15 Ом, R6=50 Ом - резисторы

L=5*57 мГн - катушка индуктивности

С=20 мкФ - конденсатор

В начальный момент t=t0=0 ключ находится в положении 1. При этом цепь разомкнута, напряжение на конденсаторе и ток на катушке равны нулю (U=0, I=0). Происходит первое переключение ключа (ключ мгновенно переводится в положение 2). При этом происходит заряд конденсатора, меняются значения U и I. В момент времени t=t1=0.01c ключ мгновенно переводится в положение 1. Конденсатор разряжается, вновь меняются значения U и I. Анализ схемы заканчивается в момент t=t1=0.02c

Реализация поставленной задачи

Вывод системы дифференциальных уравнений.

В соответствии со схемой запишем соотношения из законов Ома и Кирхгофа

Ключ в положении “2”:

i=i1+i2

ключ в положении “1”:

i=i1+i2

Исключив токи i1 и i2 получили ключ в положении “2”:

Для положения ключа “1”:

т.е. E(t)=0

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

Описание методов применяемых в курсовой работе.

Метод Рунге-Кутта:

Метод Рунге—Кутта 4 порядка столь широко распространён, что его часто называют просто методом Рунге—Кутта.

Рассмотрим задачу Коши

Тогда значение в следующей точке вычисляется по формуле:

где h — величина шага сетки по x и вычисление нового значение проходит в четыре этапа:

Этот метод имеет четвёртый порядок точности, т.е. суммарная ошибка на конечном интервале интегрирования имеет порядок O(h4) (ошибка на каждом шаге порядка O(h5))

Метод трапеций (ЛП):

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

Площадь трапеции на каждом отрезке:

Погрешность аппроксимации на каждом отрезке:

Полная формула трапеций в случае деления всего промежутка интегрирования на отрезки одинаковой длины h:

Погрешность формулы трапеций:

Метод Гаусса для решения систем линейных уравнений:

Пусть исходная система выглядит следующим образом

Тогда согласно свойству элементарных преобразований над строками эту систему можно привести к трапециальному виду:

Переменные называются главными переменными. Все остальные называются свободными.

Если , то рассматриваемая система несовместна.

Предположим, что .

Перенесём свободные переменные за знаки равенств и поделим каждое из уравнений системы на свой коэффициент при самом левом (, где — номер строки):

где

Если свободным переменным системы (2) придавать все возможные значения и вычислить через них главные переменные, то мы получим все решения этой СЛАУ. Так как эта система получена путём элементарных преобразований над исходной системой (1), то по теореме об эквивалентности при элементарных преобразованиях полученное нами решение является решением системы (1).

Метод парабол (Симпсона):

Основан на интерполяции подынтегральной функции y=f(x) параболой на паре соседних частичных интервалов (i=1,3), т.е. интерполяционным полиномом второй степени

Тогда оказывается, что частичный интеграл для интерполяционного полинома вычисляется по формуле

Выбирая число частичных интервалов n четным, получим суммированием составную

формулу для приближенного вычисления всего интеграла методом парабол

Аппроксимация полиномом 1-ой степени:

меру близости искомой прямой к графику узловых точек выразим величиной:

Найдем а₁ и а₂ из условия минимума суммы квадратов отклонений:

-положительно определенная квадратичная форма, то

минимум достигается в единственной точке.

Из этого условия получаем систему линейных уравнений для нахождения коэффициентов а₁,а₂ аппроксимирующей прямой.

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

Аппроксимация полиномом 2-ой степени

В матричной форме решение будет иметь вид:

А=В^-1С

В А С

Составим систему нормальных уравнений в матричной форме:

Решение системы дифференциальных уравнений на языке Free Pascal:

2. Решение для системы дифференциальных уравнений (3):

program cursc3;

var a,b,m,d,f,t,i,u,h,k1,k2,k3,k4,m1,m2,m3,m4,r1,r2,r3,r4,r5,r6,l,c:real;

begin

writeln(' vvod R1,R2,R3,R4,R5,R6');

READLN(R1,R2,R3,R4,R5,R6);

l:=5.57*0.001;

c:=20*0.000001;

A:=(R5+R6)/(R3+R5+R6);

B:=R1*R2/(R1+R2);

M:=1/(R3+R5+R6);

D:=R4+B+A*r3;

F:=R2/(R1+R2);

h:=0.00025;

i:=0;

u:=0;

t:=0;

while (t<=0.01001) do begin

k1:=h*(15*sin(62.8*t+7*pi/5)*f/l-d*i/l-a*u/l);

m1:=h*(i*a/c-m*u/c);

k2:=h*(15*sin(62.8*(t+h/2)+7*pi/5)*f/l-d*(i+k1/2)/l-a*(u+m1/2)/l);

m2:=h*((i+k1/2)*a/c-m*(u+m1/2)/c);

k3:=h*(15*sin(62.8*(t+h/2)+7*pi/5)*f/l-d*(i+k2/2)/l-a*(u+m2/2)/l);

m3:=h*((i+k2/2)*a/c-m*(u+m2/2)/c);

k4:=h*(15*sin(62.8*(t+h)+7*pi/5)*f/l-d*(i+k3)/l-a*(u+m3)/l);

m4:=h*((i+k3)*a/c-m*(u+m3)/c);

i:=i+(k1+2*k2+2*k3+k4)/6;

u:=u+(m1+2*m2+2*m3+m4)/6;

t:=t*1000;

writeln(' pri t=',t:2:2,'*10E-3',' u=',u:2:2,' i=',i:4:4);

t:=t/1000;

t:=t+h;

end;

readln

end.

Результаты работы программы:

Возьмем начальные данные для решения системы дифференциальных уравнений (2) из результата решения системы дифференциальных уравнений (3).

1. Решение для системы дифференциальных уравнений (2):

program cursc;