Реферат: Температурный расчет с помощью вычислений информационной математики
![](/file/7d5aacc146456be8f3642.jpg)
👉🏻👉🏻👉🏻 ВСЯ ИНФОРМАЦИЯ ДОСТУПНА ЗДЕСЬ ЖМИТЕ 👈🏻👈🏻👈🏻
По
длинной квадратного
сечения трубе
течет горячая
жидкость. Труба
наполовину
погружена в
ледяную ванну,
так, что температура
нижней половины
поверхности
трубы равна
0 0
С. Верхняя плоскость
трубы имеет
постоянную
температуру
100 0
С. На участке
между ледяной
ванной и верхней
плоскостью
температура
наружной поверхности
трубы изменяется
линейно по
высоте от 0
0 С до
100 0
С. Жидкость
внутри трубы
имеет температуру
200 0
С.
Распределение
температуры
в
теле трубы
удовлетворяет
уравнению
С
погрешностью
не более 0,5
0 С вычислить
распределение
температуры
в теле трубы.
Метод
последовательной
верхней релаксации
Библиотечная
графическая
подпрограмма
Решить
диф.уравнение
в частных
производных:
с
задаными началиными
условиями на
границах области
дифференцирования.
При
решении уравнения
приблизительно
заменю производные
второго порядка
конечно-разностными
отношениями:
в
результате
чего диф.уравнение
преобразуется
в 5-ти диаганальную
систему алгеброических
уравнений n-го
порядка.
Систему
алгеброических
уравнений буду
решать методом
Зейделя.
Погрешность
решения задачи
найду по формуле:
где,
и
-решения,полученные
для одной и той
же точки с разными
шагами.
Так назван
метод решения
краевых задач,
основанный
на приближенной
замене производных,
входящих в
дифференциальные
уравнения
и краевые условия,
нонечно-разностными
отношениями.
Эта замена
позволяет
свести краевую
задачу к задаче
решения системы
алгебраических
уравнений.
Конечные
разности и
производные.Пусть
некоторая
функция y(x)
задана на отрезке
[a,b].
Будем считать,
что она непрерывна
и многократно
дифференцируема
на этом отрезке.
Разделим отрезок
на равные части
длиною h
и обозначим
точки деления
x0,x1,...,xi,...,xn.Значения
функции в этих
точках обозначим
соответственно
y0,y1,...,yi,...,yn.Первой
центральной
разностью в
i-й
точке
(i=1,2,...,n-1) называют
разность:
С помощью
этой разности
можно приближенно
вычислить
значение первой
производной
у`
в i-й
точке.
Разложим
функцию y(x)
в степенной
ряд. приняв за
центр разложения
точку
xi и ограничившись
четырьмя членами:
Аналогично
найдем значение
ф-ции и в точке ,отстоящей
от центра разложения
на шаг (-h):
Таким
образом,производная
y`
приближонно
заменяется
конечно-разностным
отношением
с ошибкой порядка
h*h:
Второй
центральной
разностью ф-ции
y(x)
в i-й
точке называют
величину:
С
помощью этой
разности можно
приближонно
вычислить
значение второй
производной
y``
в i-й
точке.Используем
теперь 5 членов
разложения
в ряд Тейлора:
Таким
образом,вторая
производная
y``
с ошибкой порядка
h*h
может быть
приближонно
заменена
конечно-разностным
отношением:
При определении
разностей в
i
-и точке использовались
значения функции
в точках, расположенных
симметрично
относительно
xi
. Поэтому
эти разности
называются
центральными.
Существуют
также левые
и правые разности,
использующие
точки, расположенные
соответственно
левее и правее
точки xi.
С помощью этих
разностей можно
также приближенно
вычислять
значения производных,
но погрешность
при этом будет
больше
-порядка
h.
Разностные
системы уравнений
составляются
в следующем
порядке.
1. Исходное
дифференциальное
уравнение
преобразуют
к такой форме,
чтобы затем
получить из
него наиболее
простую разностную
систему уравнений.
При этом учитывают,
что коэффициенты
при производных
войдут в разностную
схему одновременно
в несколько
ее членов и
затем будут
распространены
на всю систему
уравнений.
Поэтому желательно
иметь единичные
коэффициенты
при производных
в исходном
уравнении.
2. На
интервале
интегрирования
исходного
уравнения
устанавливают
равномерную
сетку с шагом
h и записывают
разностную
схему, приближенно
заменяя производные
соответствующими
центральными
конечно-разностными
отношениями.
3.Применяя
разностную
схему для узлов
сетки записывают
разностные
уравнения. При
этом можно
получить уравнения
содержащие
так называемые
внеконтурные
неизвестные,
то есть неизвестные
в точках, лежащих
за пределами
установленной
сетки.
4.В разностной
форме записывают
краевые условия
и составляют
полную систему
разностных
уравнений.
Оценка
погрешности
решения краевой
задачи
Решение
разностной
системы уравнений
дает приближенное
решение
краевой задачи.
Поэтому возникает
вопрос о точности
этого приближенного
решения.
Для линейных
краевых задач
доказана теорема
о том, что порядок
точности решения
краевой задачи
не ниже порядка
точности
аппроксимации
производных
конечно-разностными
отношениями.
Оценку погрешности
производят
приемом Рунге.
Краевую задачу
решают дважды:
с шагом сетки
h и
с шагом сетки
H=kh, погрешность
решения с малым
шагом h
оценивают
по формуле:
где y(h)
и y(H)
- решения,
полученные
для одной и той
же точки
-xi отрезка
интегрирования
с разными шагами.
Относительную
погрешность
E
оценивают по
формуле:
Если при
составлении
разностной
системы уравнений
используются
левые или правые
разности, то
погрешность
решения будет
выше, порядка
0(h),
и для ее оценки
в формулах
следует заменить
k*k
на k
.
Применение
метода конечных
разностей для
решения уравнений
в частных проиэводных
Для
применения
разностного
метода в области
изменения
независимых
переменных
вводят некоторую
сетку. Все
производные,
входящие в
уравнение и
краевые условия,
заменяют разностями
значений функции
в узлах сетки
и получают
таким образом
алгебраическуго
систему уравнений.
Решая эту систему,
находят приближенное
решение задачи
в узлах сетки.
c------------------------------------------------------------------
c
ПОДПРОГРАММА
СОСТАВЛЕНИЯ
СИСТЕМЫ УРАВНЕНИЙ
c
МЕТОДОМ
КОНЕЧНЫХ РАЗНОСТЕЙ
c
real N-количество
уравнений(примерное
число,желательно
N=M*P)
c
real y(6,N)-выходной
массив уравнений,содержащий
следующие поля:
c
y(3,N)-коэфициен
уравнения для
Q(y(1,N)-1,y(2,N))
c
y(3,N)=h^2/(2*(h^2+k^2))
c
y(4,N)-коэфициен
уравнения для
Q(y(1,N),y(2,N)-1)
c
y(4,N)=k^2/(2*(h^2+k^2))
c
y(5,N)-коэфициен
уравнения для
Q(y(1,N)+1,y(2,N))
c
y(5,N)=h^2/(2*(h^2+k^2))
c
y(6,N)-коэфициен
уравнения для
Q(y(1,N),y(2,N)+1)
c
y(6,N)=k^2/(2*(h^2+k^2))
c
integer N-выходное
количество
получившихся
уравнений
c------------------------------------------------------------------
integer M,P,IIX,IIY,NN,N,KR1,KR2,KR3
c-----------------------------------------------------------------
c-----------------------------------------------------------------
c-----------------------------------------------------------------
c
присваивание
начальных
значений
c
iix-номер текущего
узла по оси X
c
iiy-номер текущего
узла по оси Y
c-----------------------------------------------------------------
c-----------------------------------------------------------------
c
проверка
границы трубы
с жидкостью
c-----------------------------------------------------------------
if
((IIY.ge.KR1).and.(IIY.le.KR2).and.(IIX.ge.KR3)) then
c-----------------------------------------------------------------
c-----------------------------------------------------------------
elseif (((IIY.lt.KR1).or.(IIY.gt.KR2)).and.(IIX.eq.M))
then
c-----------------------------------------------------------------
c
составление
уравнений во
внутренних
точках фигуры
c-----------------------------------------------------------------
c-----------------------------------------------------------------
c
присваивание
начальных
значений на
границе фигуры
c------------------------------------------------------------------
Для
тестирования
составлю разностную
систему с шагом
вдоли оси X
и Y=0.05
Неизвестные
значения в
узлах матрицы
находящихся
внутри фигуры
высчитываются
по формуле:
Неизвестные
значения в
узлах матрицы
находящихся
на оси симметрии
высчитываются
по формуле:
Метод Зейделя
относится к
числу итерационных
методов, в которых
принципиально
отсутствует
фактор накопления
погрешностей.
Поэтому он
широко применяется
для решения
больших систем
уравнений.
Будем рассматривать
корни решаемой
системы как
компоненты
некоторого
вектора у
. Основная
идея всех
итерационных
методов заключается
в том, что берется
приближенное
значение вектора
у и по формулам,
составленным
на основании
решаемых уравнений,
вычисляется
новое приближенное
значение вектора
у .
Назовем эти
приближенные
значения
y(k) и
y(k+1) соответственно.
Поскольку
исходное
приближение
выбиралось
произвольно,
то у(k+1)в
свою очередь
может послужить
исходным для
получения по
тем же формулам
нового приближения
y(k+2)
. Очевидно,
этот процесс
можно продолжать
сколь угодно
долго. Говорят,
что процесс
итераций сходится,
если получаемая
при этом последовательность
векторов у(k)
(к=0,1,2,...}
имеет своим
пределом вектор
y,являющийся
точным решением
системы:
На практике
невозможно
достигнуть
этого предела,
но можно приблизиться
н нему с любой
наперед заданной
точностью. Так
и поступают
задаются некоторой
погрешностью,
вектором начального
приближения
и получают
последовательные
приближения
до тех пор, пока
действительная
погрешность
корней не станет
меньше заданной.
Различные
методы отличаются
друг от друга
способом вычисления
очередного
приближения,
но во всех методах
существуют
две главные
проблемы:
обеспечение
сходимости
процесса итераций;
Предполагая,
что диагональные
коэффициенты
разрешим
первое уравнение
относительно
y1 ,
второе
- относительно
y2
и т.д.
Такую систему
будем в
дальнейшем
называть приведенной.
Метод Зейделя
заключается
в следующем.
Выбрав вектор
начального
приближения
подставим
его компоненты
в правую часть
первого уравнения
системы
и вычислим
первую компоненту
y`1
нового вектора
y`(ср)
. В правую
часть второго
уравнения
подставим
компоненты
(y`1,y2,y3,...,yn)
и вычислим
вторую компоненту
y`2'нового
вектора. В третье
уравнение
подставим
(y`1,y`2,y3,...,yn) и т.д.
Очевидно,подстановкой
в каждое
уравнение мы,
дойдя до последнего
уравнения,
обновим все
компоненты
исходного
вектора и получим
первое приближение
к решению
Далее
, взяв в
качестве исходного
вектор у`(ср)
, выполним
вторую
итерацию
и получим y``(ср).
Этот процесс
будем продолжать
до достижения
заданной степени
точности.
Оценка
погрешности
приближений
процесса Зейделя
Для оценки
погрешности
прежде всего
вычисляют
показатель
скорости сходимости
То есть для
каждой строки
матрицы коэффициентов
системы
вычисляется
сумма модулей
коэффициентов,
лежащих правее
главной диагонали
:
и сумма модулей
коэффициентов,
лежащих левее
главной диагонали:
Для каждой
i-й
строки
(i =1,2,...,n
) вычисляется
отношение
и в качестве
берется
максимальное
из этих отношений.
Чем меньше
окажется
,
тем большей
будет скорость
сходимости.
Для процесса
Зейделя справедлива
следующая
оценка погрешности
К-го приближения:
то есть модуль
отклонения
любого
i -го корня
системы в К-м
приближении
от точного
значения того
же корня
не больше, чем
умноженное
на множитель
максимальное
из приращений
корней, полученных
в результате
перехода от
(K-1)
-го приближения
к К-му.
Если задаться
абсолютной
погрешностью
и
потребовать
выполнения
условия
то есть заданная
степень точности
на К-й итерации
будет достигнута.
На практике
это означает,
что после каждой
итерации необходимо
выделить тот
корень, изменение
которого по
сравнению с
предыдущим
значением
оказалось
наибольшим
по модулю. Модуль
приращения
этого корня
необходимо
умножить на
и сравнить
результат с
выбранной
абсолютной
погрешностью .
Достаточные
условия сходимости
процесса Зейделя
Если модули
коэффициентов
системы
удовлетворяют
хотя бы одному
из условий
то процесс
Зейделя для
соответствующей
приведенной
системы сходится
к её единственному
решению при
любом выборе
начального
вектсра y(ср)
Такие системы
называют системами
с диагональным
преобладанием.
Метод Зейдедя
имеет свойство,
позволяющее
обеспечить
сходимость
процесса для
любых систем
уравнений с
неособенной
матрицей
коэфициентов.
Если обе
части систем
с неособенной
матрицей коэфициентов
А=[aij]
умножить слева
на транспонировнную
матриц A*[aij]
, то будет
получена новая,
равносильная
исходной система,
которая называется
нормальной.
Процесс Зейделя
для приведенной
системы, полученной
из нормальной,
всегда сходится
независимо
от выбора нача
льного приближения.
c-----------------------------------------------------------------
cПОДПРОГРАММА
РЕШЕНИЯ СИСТЕМЫ
УРАВНЕНИЙ
МЕТОДОМ ЗЕЙДЕЛЯ
с
integer N-входное
количество
уравнений
c
real y(6,N)-входной
массив уравнений,содержащий
следующие поля:
c
y(3,N)-коэфициен
уравнения для
Q(y(1,N)-1,y(2,N))
c
y(3,N)=h^2/(2*(h^2+k^2))
c
y(4,N)-коэфициен
уравнения для
Q(y(1,N),y(2,N)-1)
c
y(4,N)=k^2/(2*(h^2+k^2))
c
y(5,N)-коэфициен
уравнения для
Q(y(1,N)+1,y(2,N))
c
y(5,N)=h^2/(2*(h^2+k^2))
c
y(6,N)-коэфициен
уравнения для
Q(y(1,N),y(2,N)+1)
c
y(6,N)=k^2/(2*(h^2+k^2))
c
real Q(M,P)-входной
массив начальных
значений Y
c
real Q(M,P)-выходной
массив вычисленых
значений Y
c------------------------------------------------------------------
c------------------------------------------------------------------
c
вычисление
коэфициента
сходимости
процесса
C
НО Т.К. MJ=0.5 ТО KM=1 И
СЛЕДОВАТЕЛЬНО
ЕГО МОЖНО ОПУСТИТЬ
c-----------------------------------------------------------------
c
KM=(y(5,1)+y(6,1))/(1-y(5,1)+y(6,1))
c------------------------------------------------------------------
c------------------------------------------------------------------
c------------------------------------------------------------------
c
вычисление
погрешности
на данной итерации
c------------------------------------------------------------------
if
(abs(NEXTQ-q(y(1,i),y(2,i))).gt.EI)
+
EI=abs(NEXTQ-q(y(1,i),y(2,i)))
c
print *,'x=',y(1,i),' y=',y(2,i)
c
print '(16h Итерация
номер ,i5,13h погрешность=,E15.7)',S,EI
c
print '(21e10.3)',(q(j,i),j=1,M)
В
качестве теста
выполним одну
итерацию для
системы , полученной
в предыдущем
пункте.
c------------------------------------------------------------------
c
ПОДПРОГРАММА
СОСТАВЛЕНИЯ
СИСТЕМЫ УРАВНЕНИЙ
c
МЕТОДОМ
КОНЕЧНЫХ РАЗНОСТЕЙ
c
real N-количество
уравнений(примерное
число,желательно
N=M*P)
c
real y(6,N)-выходной
массив уравнений,содержащий
следующие поля:
c
y(3,N)-коэфициен
уравнения для
Q(y(1,N)-1,y(2,N))
c
y(3,N)=h^2/(2*(h^2+k^2))
c
y(4,N)-коэфициен
уравнения для
Q(y(1,N),y(2,N)-1)
c
y(4,N)=k^2/(2*(h^2+k^2))
c
y(5,N)-коэфициен
уравнения для
Q(y(1,N)+1,y(2,N))
c
y(5,N)=h^2/(2*(h^2+k^2))
c
y(6,N)-коэфициен
уравнения для
Q(y(1,N),y(2,N)+1)
c
y(6,N)=k^2/(2*(h^2+k^2))
c
integer N-выходное
количество
получившихся
уравнений
c------------------------------------------------------------------
integer M,P,IIX,IIY,NN,N,KR1,KR2,KR3
c-----------------------------------------------------------------
c-----------------------------------------------------------------
c-----------------------------------------------------------------
c
присваивание
начальных
значений
c
iix-номер текущего
узла по оси X
c
iiy-номер текущего
узла по оси Y
c-----------------------------------------------------------------
c-----------------------------------------------------------------
c
проверка
границы трубы
с жидкостью
c-----------------------------------------------------------------
if
((IIY.ge.KR1).and.(IIY.le.KR2).and.(IIX.ge.KR3)) then
c-----------------------------------------------------------------
c-----------------------------------------------------------------
elseif (((IIY.lt.KR1).or.(IIY.gt.KR2)).and.(IIX.eq.M))
then
c-----------------------------------------------------------------
c
составление
уравнений во
внутренних
точках фигуры
c-----------------------------------------------------------------
c-----------------------------------------------------------------
c
присваивание
начальных
значений на
границе фигуры
c------------------------------------------------------------------
c-----------------------------------------------------------------
c
ПОДПРОГРАММА
РЕШЕНИЯ СИСТЕМЫ
УРАВНЕНИЙ
МЕТОДОМ ЗЕЙДЕЛЯ
с
integer N-входное
количество
уравнений
c
real y(6,N)-входной
массив уравнений,содержащий
следующие поля:
c
y(3,N)-коэфициен
уравнения для
Q(y(1,N)-1,y(2,N))
c
y(3,N)=h^2/(2*(h^2+k^2))
c
y(4,N)-коэфициен
уравнения для
Q(y(1,N),y(2,N)-1)
c
y(4,N)=k^2/(2*(h^2+k^2))
c
y(5,N)-коэфициен
уравнения для
Q(y(1,N)+1,y(2,N))
c
y(5,N)=h^2/(2*(h^2+k^2))
c
y(6,N)-коэфициен
уравнения для
Q(y(1,N),y(2,N)+1)
c
y(6,N)=k^2/(2*(h^2+k^2))
c
real Q(M,P)-входной
массив начальных
значений Y
c
real Q(M,P)-выходной
массив вычисленых
значений Y
c------------------------------------------------------------------
c------------------------------------------------------------------
c
вычисление
коэфициента
сходимости
процесса
C
НО Т.К. MJ=0.5 ТО KM=1 И
СЛЕДОВАТЕЛЬНО
ЕГО МОЖНО ОПУСТИТЬ
c-----------------------------------------------------------------
c
KM=(y(5,1)+y(6,1))/(1-y(5,1)+y(6,1))
c------------------------------------------------------------------
c------------------------------------------------------------------
c------------------------------------------------------------------
c
вычисление
погрешности
на данной итерации
c------------------------------------------------------------------
if
(abs(NEXTQ-q(y(1,i),y(2,i))).gt.EI)
+
EI=abs(NEXTQ-q(y(1,i),y(2,i)))
c
print *,'x=',y(1,i),' y=',y(2,i)
c
print '(16h Итерация
номер ,i5,13h погрешность=,E15.7)',S,EI
c
print '(21e10.3)',(q(j,i),j=1,M)
c------------------------------------------------------------------
c
ПОДПРОГРАММА
АЛФАВИТНО-ЦИФРОВОГО,МОЗАИЧНОГО
c
real Q(M,P)-входной
массив значений
Y
c------------------------------------------------------------------
character
a(11)/'.','+','*','','','-','-','-','' ,'-','-'/
print '(400A2)',(a(int(q(I,J)/21)+1),I=1,M),
+
(a(int(q(I,J)/21)+1),I=M-1,1,-1)
print *,'''',a(I),'''','---> от
',20*(I-1),', до ',
c------------------------------------------------------------------
c
ПОДПРОГРАММА
ВЫЧИСЛЕНИЯ
ОШИБКИ
c
real q-массив значений
Y с шагом =2*h
c
real qq-массив значений
Y с шагом =h
c------------------------------------------------------------------
real
qq(M,P),q(int(M/2)+1,int(P/2)+1),max,E,other
c------------------------------------------------------------------
c------------------------------------------------------------------
real
y(6,90000),q(300,300),H/.05/,K/.05/,E/.5/,qq(300,300)
c
print *,'Введите шаг
вдоль оси X '
c
print *,'Введите шаг
вдоль оси Y '
c
print *,'Введите
точность вычислений
'
1.
И.С.Березин,Н.П.Жидков
’Методы
вычислений’,том
1,М.,1966,632 стр.
2.’
Численные
методы решения
задач на ЭВМ
’ , Учебное
пособие , Г.Н.Рубальченко
, К. ,
1989
, 148 стр.
3.’Справочник
языка ФОРТРАН’
, М.,1996 ,106 стр.
Температурный
расчет с помощью
вычислений
информационной
математики.
Название: Температурный расчет с помощью вычислений информационной математики
Раздел: Рефераты по математике
Тип: реферат
Добавлен 06:08:14 15 августа 2005 Похожие работы
Просмотров: 96
Комментариев: 16
Оценило: 4 человек
Средний балл: 5
Оценка: неизвестно Скачать
Срочная помощь учащимся в написании различных работ. Бесплатные корректировки! Круглосуточная поддержка! Узнай стоимость твоей работы на сайте 64362.ru
Привет студентам) если возникают трудности с любой работой (от реферата и контрольных до диплома), можете обратиться на FAST-REFERAT.RU , я там обычно заказываю, все качественно и в срок) в любом случае попробуйте, за спрос денег не берут)
Да, но только в случае крайней необходимости.
Реферат: Температурный расчет с помощью вычислений информационной математики
Сочинение Герой Моего Времени 9 Класс
Учебное пособие: Методические указания по выбору и применению асинхронного частотно-регулируемого электропривода мощностью до 500 квт врд 39 10-052-2001
Контрольная работа по теме Маркетинговые коммуникации холдинга Яки’Да Restaurant Group
Реферат: И.С.Тургенев. Скачать бесплатно и без регистрации
Интернет Источники В Курсовой Работе
Курсовая работа по теме Образ Петербурга в произведениях Н.В. Гоголя
Отчет По Практике Бухгалтера 2022
Сочинение Рассуждение На Тему Сильный Характер
Контрольная Работа 1 По Математике 5 Класс
Реферат: Tolerance Happens Essay Research Paper Tolerance HappensEffective
Реферат: Тэги списков. Скачать бесплатно и без регистрации
Отчет По Учебной Практике Страховое Дело Заочное
Дипломная Работа На Тему Бизнес-План Производственного Предприятия
Курсовая работа по теме Особенности газетных заголовков в английском языке
Вступление Сочинения Конфликт Между Поколениями
Доклад по теме Березневі статті
Контрольная работа: Понятие и методы решения задач бухгалтерского учета
Реферат по теме Асаблівасці фарміравання беларускай нацыi
Реферат: Диагностика готовности ребенка к школе
Курсовая работа по теме Создание и редактирование таблиц Microsoft Word
Реферат: Внутренняя политика 60-90 гг.
Доклад: Спорт и окружающая среда: перспективы развития
Сочинение: Я люблю Пушкина Цветаевой (цикл М. Цветаевой «Стихи Пушкину»)