ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 07.12.2021
Просмотров: 2280
Скачиваний: 1
81
7.
Для
формирования
файла
отчета
запускаем
двойным
щелчком
левой
кнопки
мыши
на
рабочем
столе
программу
Microsoft Word,
если
же
ярлык
отсутствует
,
то
открываем
программу
через
кнопку
«
Пуск
».
Открываем
новый
документ
.
В
начале
документа
необходимо
оформить
титульный
лист
,
описать
математическую
постановку
задачи
и
результаты
выполнения
подготовительных
расчетов
.
Затем
скопировать
основные
результаты
расчетов
из
программы
Parab.mcd
в
документ
и
оформить
итоговый
отчет
.
Копирование
– ‘Ctrl’+’Insert’,
вставка
– ‘Shift’+’Insert’.
Сохранить
документ
как
«
ФамилияСтудента
_
группа
_Parab.doc»
и
распечатать
.
Файл
отчета
оформить
аналогично
приложению
А
,
описывающему
выполнение
лабораторной
работы
№
1
.
3.5.
Порядок
выполнения
лабораторной
работы
1.
Повторить
главу
1.
Изучить
разделы
3.1–3.6
данной
главы
и
подготовить
ответы
на
контрольные
вопросы
из
раздела
3.9.
2.
Пройти
собеседование
с
преподавателем
,
получить
допуск
к
выполнению
работы
на
ЭВМ
,
номер
варианта
задания
,
значение
параметра
n
и
указания
по
выбору
пробных
и
поверочных
функций
.
3.
Выполнить
первый
пункт
задания
,
связанный
с
построением
ряда
Фурье
для
точного
решения
задачи
)
,
(
t
x
U
и
нахождением
длины
отрезка
этого
ряда
,
обеспечивающую
точность
решения
0,001.
4.
Выполнить
подготовительный
шаг
алгоритма
метода
Галеркина
и
,
если
)
(
0
x
u
не
является
точным
решением
задачи
,
подготовить
все
числовые
и
строчные
данные
для
расчетов
и
в
пункте
«
Постановка
задачи
»
программы
Parab.mcd
ввести
их
вместо
данных
примера
,
введенных
изначально
.
5.
В
пункте
«
Получение
точного
решения
»
программы
ввести
число
,
намного
превышающее
найденное
в
3-
м
пункте
число
слагаемых
в
разложении
точного
решения
в
тригонометрический
ряд
Фурье
(
чтобы
гарантировать
достаточную
точность
решения
и
в
дальнейшем
считать
его
точным
).
Скопировать
график
получившегося
точного
решения
)
,
(
T
x
U
в
файл
отчета
.
6.
В
пункте
«
Получение
приближенного
решения
»
рассмотрено
применение
трех
систем
пробных
и
поверочных
функций
.
По
заданию
преподавателя
ввести
(
вместо
уже
введенных
для
примера
)
системы
пробных
)
,
1(
x
k
V
и
поверочных
)
,
(
x
k
W
функций
,
указанных
во
2-
м
пункте
(
см
.
раздел
3.6).
Выполнить
построение
n-
го
пробного
решения
задачи
.
Следует
скопировать
в
файл
отчета
вектор
коэффициентов
)
(
T
v
k
(
элементы
вектора
k
Y
,
100
программы
)
пробных
решений
и
набрать
в
отчете
решение
с
этими
коэффициентами
.
Так
же
необходимо
скопировать
в
этот
файл
пункт
«
Выводы
».
7.
Оформить
и
распечатать
файл
отчета
по
лабораторной
работе
,
который
должен
содержать
титульный
лист
,
математическую
постановку
задачи
и
ее
физическую
интерпретацию
,
результаты
выполнения
подготовительных
82
расчетов
,
основные
результаты
расчетов
на
ЭВМ
,
выводы
о
возможностях
использованных
систем
пробных
и
поверочных
функций
и
наиболее
приближенное
к
точному
аналитическое
решение
.
8.
Защитить
отчет
.
3.6.
Программа
в
системе
MathCAD
и
тестирующий
пример
В
данном
пункте
приведен
текст
программы
Parab.mcd,
разработанной
для
решения
начально
-
краевой
задачи
для
одномерного
параболического
уравнения
методом
Галеркина
.
В
тексте
разбирается
получение
значений
пробного
решения
)
,
(
5
t
x
u
при
1
t
задачи
:
найти
функцию
)
,
(
t
x
u
,
удовлетворяющую
в
области
}
1
0
,
0
:
)
,
{(
2
t
x
R
t
x
D
уравнению
2
2
1
,
0
x
u
t
u
(3.28)
и
условиям
,
1
)
,
0
(
t
u
,
2
)
,
(
t
u
(3.29)
.
8233
.
2
1
1
1
)
(
)
0
,
(
2
2
x
x
x
x
x
f
x
u
(3.30)
Задача
(3.28)–(3.30)
является
частным
случаем
задачи
(3.25)–(3.27)
при
1
,
1
,
2
,
1
,
1
,
0
4
3
2
1
T
t
c
c
c
с
.
Использовать
три
системы
пробных
и
поверочных
функций
:
1.
Пробные
и
поверочные
функции
–
многочлены
(2.28);
2.
Пробные
функции
–
многочлены
(2.28),
поверочные
функции
–
многочлены
Лежандра
(2.31);
3.
Пробные
и
поверочные
функции
–
тригонометрические
функции
x
k
)
1
2
sin(
.
Лабораторная
работа
«
Решение
начально
-
краевой
задачи
для
одномерного
параболического
уравнения
методом
Галеркина
»
Задание
на
лабораторную
работу
1.
В
пункте
«
Постановка
задачи
»
ввести
вместо
данных
примера
непрерывные
функции
уравнения
K(x) (K>0),
(x), g(x), f(x)
и
числовые
параметры
задачи
a, b, a0, a1, a2, b0, b1, b2, c1, c2, c3, c4
своего
варианта
.
2.
В
пункте
«
Получение
точного
решения
»
программы
ввести
число
слагаемых
в
разложении
решения
в
ряд
,
намного
превышающее
найденное
аналитически
число
,
обеспечивающее
точность
решения
0.001.
Скопировать
график
полученной
интегральной
кривой
в
файл
отчета
.
3.
В
пункте
«
Получение
приближенного
решения
»
выполнить
построение
n-
го
пробного
решения
задачи
тремя
системами
пробных
и
поверочных
83
функций
.
Скопировать
в
файл
отчета
вектор
коэффициентов
Y
100,k
пробного
решения
и
набрать
в
отчете
решение
с
этими
коэффициентами
при
t=T.
4.
Скопировать
результаты
пункта
«
Выводы
»
в
файл
отчета
,
и
,
анализируя
их
,
сделать
в
файле
отчета
выводы
о
точности
построенных
решений
.
Постановка
задачи
Требуется
в
двумерной
области
D={(x, t) | a
x
b
t
0
}
найти
решение
U(x, t)
дифференциального
уравнения
t
U
d
d
K x t
(
)
2
x
U
d
d
2
x
K x t
(
)
x
U
d
d
d
d
x t
(
)
U
g x t
(
)
удовлетворяющее
двум
краевым
условиям
a0 U a t
(
)
a1
x
U a t
(
)
d
d
a2 t
( )
b0 U b t
(
)
b1
x
U b t
(
)
d
d
b2 t
( )
и
начальному
условию
U x
0
(
)
f x
( )
Рассмотрим
случай
,
когда
функции
K,
, g, a2, b2
не
зависят
от
t.
Введите
непрерывные
функции
уравнения
K(x) (K>0),
(x), g(x), f(x)
и
числовые
параметры
задачи
a, b, a0, a1, a2, b0, b1, b2, c1, c2, c3, c4
c1
0.1
c2
1
c3
2
c4
1
K x
( )
c1
x
( )
0
g x
( )
0
a
0
b
a0
1
a1
0
a2
c2
b0
1
b1
0
b2
c3
f x
( )
c4 x
2
c3
c2
c4 b
2
b
x
c2
то
есть
f x
( )
x
2
1
2
x
1
Проверим
соответствие
граничных
и
начальных
условий
if a0 f a
( )
a1
a
f a
( )
d
d
a2
"Yes"
"No"
"Yes"
if b0 f b
( )
b1
b
f b
( )
d
d
b2
"Yes"
"No"
"Yes"
Если
хотя
бы
одно
условие
не
выполняется
(="No"),
то
задача
поставлена
не
корректно
.
84
Введите
конечный
момент
времени
,
до
которого
необходимо
провести
исследование
для
вашего
варианта
T
1
Получение
точного
решения
Найдем
точное
решение
U(x, t),
используя
разложение
функции
в
ряд
Фурье
.
Если
const
c
x
K
x
x
1
)
(
,
0
)
(
,
0
)
(
,
то
решение
имеет
вид
U(x, t)=U0(x, t)+
1
M
k
A
k
e
c1
k
2
2
b a
(
)
2
t
a
b
x
k
sin
.
Введите
число
слагаемых
,
обеспечивающих
достаточно
большую
точность
решения
(
для
примера
M=3
обеспечивает
точность
0,001,
поэтому
возьмем
число
,
намного
превышающее
данное
,
например
, M=30)
M
30
Если
a1=0, b1=0,
то
функцию
U0(x, t)
можно
взять
в
виде
U0 x
( )
b0 a2
b
b2 a0
a
a0 b0
b
a
(
)
b2 a0
b0 a2
(
)
x
a0 b0
b
a
(
)
U0 x
( )
1
x
Вычислим
коэффициенты
A
k
i
1
M
A2
i 1
2
b
a
a
b
x
f x
( )
U0 x
( )
(
)
sin
i
x
b
a
d
f x
( )
U0 x
( )
x
2
1
2
x
x
Следовательно
,
точное
решение
U(x, t)
имеет
вид
UT x t
(
)
U0 x
( )
1
M
k
A2
k
1
e
c1
k
2
2
b a
(
)
2
t
sin
k
x
b
a
График
точного
решения
при
t=T
0
2
2
0
2
UT x T
(
)
x
Скопируйте
график
полученной
интегральной
кривой
в
файл
отчета
.
85
Получение
приближенного
решения
Введите
порядок
пробного
решения
U
n
=
V(0,x)
+
1
n
k
V k x
(
)
H k t
(
)
.
n
5
1.
Введите
систему
пробных
функций
:
k
1
n
V0 k x
(
)
x
a
(
)
k
x
b
(
)
Нормируем
их
.
Для
этого
вычислим
нормировочные
коэффициенты
i
1 6
VV
i
1
a
b
x
V0 i x
(
)
(
)
2
d
Получили
нормированные
пробные
функции
V k x
(
)
if k
0
V0 k x
(
)
VV
k
1
b0 a2
b
b2 a0
a
a0 b0
b
a
(
)
b2 a0
b0 a2
(
)
x
a0 b0
b
a
(
)
Введите
функции
V1(k,x)
и
V2(k,x),
равные
первой
и
второй
производной
от
функции
V(k,x)
k
1
n
V1 k x
(
)
if k
0
x
a
(
)
k
x
b
(
)
k
x
a
(
)
k
1
VV
k
1
b2 a0
b0 a2
a0 b0
b
a
(
)
V2 k x
(
)
if k
0
if k
1
2 k
x
a
(
)
k
1
x
b
(
)
k
k
1
(
)
x
a
(
)
k
2
VV
k
1
2
VV
k
1
0
Введите
систему
поверочных
функций
:
W k x
(
)
V k x
(
)
т
.
е
.
для
примера
в
качестве
поверочных
возьмем
пробные
функции
.
Найдем
коэффициенты
системы
дифференциальных
уравнений
A
t
H
d
d
C H
B
для
отыскания
функций
H
k
(t)
с
начальными
условиями
A H
0
( )
D1:
i
1
n
B
i
1
a
b
x
K x
( )
V2
0
x
(
)
x
K x
( )
V1
0
x
(
)
d
d
x
( )
V
0
x
(
)
g x
( )
W i x
(
)
d
i
1 n
j
1 n
A
i
1
j
1
a
b
x
V j x
(
)
W i x
(
)
d