Будь умным!


У вас вопросы?
У нас ответы:) SamZan.ru

Лабораторная работа 16 Решение дифференциальных уравнений и систем Цель работы- научиться решать

Работа добавлена на сайт samzan.ru: 2015-07-05

Лабораторная работа № 16

Решение дифференциальных уравнений и систем

Цель работы: научиться решать  дифференциальные уравнения и системы дифференциальных уравнений.

Задание к работе:

К работе допущен:

Работу выполнил:

Работу защитил:

Теоретическая часть

Решение дифференциальных уравнений и систем

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

Odesolve(x,b,[step]) - возвращает функцию, которая является решением дифференциального уравнения. Используется в блоке с оператором Given.

x - переменная интегрирования, действительное число

b - конечная точка отрезка интегрирования

step - величина шага по переменной интегрирования (необязательный аргумент)

Rkadapt(y,x1,x2,n,F) - возвращает матрицу решений методом Рунге-Кутта с переменным шагом для системы обыкновенных дифференциальных уравнений с начальными условиями в векторе y, правые части которых записаны в символьном векторе F, на интервале от x1 до x2 при фиксированном числе шагов n;

rkfixed(y,x1,x2,n,F) - возвращает матрицу решений методом Рунге-Кутта системы обыкновенных дифференциальных уравнений с начальными условиями в векторе y, правые части которых записаны в символьном векторе F, на интервале от x1 до x2 при фиксированном числе шагов n.

Для численного решения одиночного дифференциального уравнения в MathCAD имеется функция Odesolve, с помощью которой может быть решена как задача Коши для обыкновенного дифференциального уравнения, так и граничная задача. Эта функция входит в состав блока решения и является его заключительным ключевым словом.

Системы линейных дифференциальных уравнений первого порядка решаются с помощью функции Rkfixed.

На рис.2 приведен пример применения функции rkfixed для решения дифференциального уравнения, описывающего процесс свободных затухающих колебаний величины электрического заряда q (К) на конденсаторе с емкостью С (Ф), включенного в замкнутый контур, содержащий также сопротивление R (Ом) и индуктивность L (Гн).

Этот процесс описывается дифференциальным уравнением второго порядка

где =d2q/dt2 – ускорение изменения заряда, К/с2;

=dq/dt   – скорость изменения заряда, К/с;

b – коэффициент затухания, 1/с, ;

wc– круговая частота собственных колебаний контура, 1/с,  

Исходные данные к решению задачи:

Начальное условие: t=0, Vq=0, q=q0.

Номер

варианта

R, Ом

L, Гн

C, Ф

q0, K

1

2

3

4

5

1

3

4

6

8

5

15

25

40

55

0,0050

0,0035

0,0040

0,0075

0,0070

1

2

3

4

5

Процесс затухания колебаний рассчитать до tk

Исходное дифференциальное уравнение второго порядка может быть преобразовано в систему дифференциальных уравнений первого порядка.

Для этого введем подстановки:

q0=q

q1=

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

 

Правые части системы дифференциальных уравнений записываются в вектор правых частей системы уравнений D(t,q).

Матрица Z размерности n строк по числу точек вывода результатов решения и m+1 столбцов, равным числу уравнений в системе. В столбцах матрицы содержатся значения переменных соответственно t, ,. На рис.2 представлен график изменения заряда от времени.

Практическая часть.

1 Найти частное решение y(x) дифференциального уравнения для своего варианта при произвольных начальных условиях и построить график решения.

2 Решите систему дифференциальных уравнений для своего варианта на отрезке [0,3]. Выведите значения искомых функций и их производных в точке с координатой х=1.5

Заключение

Я научился решать  дифференциальные уравнения и системы дифференциальных уравнений.


L

R

C


1. 1неразвитыеэто задачи крые характеризуются след чертами- а это естандартная задача для решения которой не
2. Реферат- Трудовой кодекс башкортостана
3. Каролина - Уголовно-Судебное уложение феодальной германи
4. Реферат- Экономическое развитие германских городов Берлина и Бранденбурга
5. реферат дисертації на здобуття наукового ступеня кандидата технічних наук КИЇВ ~ Дисе1
6. Лабораторна робота 2 Дослідження ефекту Холла у напівпровіднтках Виконав- Студент групи Ф
7. государственная политика занятости и ее эффективность
8. IIIIVV пальцев левой руки рубцовосгибательная контрактура I II III IV V пальцев левой руки
9. проистекают то из психологии то из социологии то из экономики но неизменно не видят систему управления1 ка
10. Лабораторная работа 5
11. Я живу в твоем подвале 1 Не делай так
12. Пояснительная записка Основная цель курса дальнейшее развитие коммуникативной компетентности
13. учебное пособие по основам международного права представляет собой переработанные лекции прочитанные на ф
14. Металлургиялы~ процестерді~ теориясы п~ні бойынша тест с~ра~тары 5В070900 Металлургия маманды~ы студе
15. . To show dvntges nd disdvntges of online shopping 3.
16. Факультету-інституту курс гр
17. Производство сорбентов на основе фрезерного торфа
18. Компьютерные технологии в преподавании математики
19. Методические рекомендации по написанию контрольных работ по направлению подготовки 030909 Юриспруденция
20. 1 2009 г ФЕДЕРАЛЬНЫЙ ГОСУДАРСТВЕННЫЙ ОБРАЗОВАТЕЛЬНЫЙ СТАНДАРТ НАЧАЛЬНОГО П