kursovik (Исследование распределения температуры в тонком цилиндрическом стержне), страница 2
Описание файла
Документ из архива "Исследование распределения температуры в тонком цилиндрическом стержне", который расположен в категории "". Всё это находится в предмете "математика" из , которые можно найти в файловом архиве . Не смотря на прямую связь этого архива с , его также можно найти и в других разделах. Архив можно найти в разделе "рефераты, доклады и презентации", в предмете "математика" в общих файлах.
Онлайн просмотр документа "kursovik"
Текст 2 страницы из документа "kursovik"
Т.е. проверить, значимо ли уменьшение дисперсии при увеличении степени многочлена.
В
качестве статического критерия рассмотрим случайную величину, равную:
(2.9)
имеющую распределение Фишера с(r ; r1) степенями свободы. Выбираем уровень распределения Фишера, находим критическое значение F*, удовлетворяющее равенству: p(F>F*=
В нашем случае F=349.02, а F*=10,13.
Е
сли бы выполнилось практически невозможное соотношение F>F, имевшее вероятность 0,01, то гипотезу Н0 пришлось бы отклонить. Но в нашем случае можно ограничиться многочленом
, коэффициенты в котором неодинаковы.
3. Нахождение коэффициента теплопроводности .
Коэффициент вычислим по формуле (1.5), обозначим:
О
пределим допустимую абсолютную погрешность величины интеграла I, исходя из требования, чтобы относительная погрешность вычисления не превосходила 0,1%, т.е.:
(3.2)
Т
.к. из (3.1) очевидно, что , то условие (3.2) заведомо будет выполнено, если:
(3.3)
Т.е. в качестве предельно допустимой абсолютной погрешности вычисления интеграла I возьмём 0,001Т (3.4)
Т=218 оС, следовательно, 0,218 оС.
3.1 Вычисление интеграла I методом трапеции
Использование теоретической оценки погрешности
Д
ля обозначения требуемой точности количества частей n, на которые нужно разбить отрезок интегрирования [0;T] определяется по формуле:
, где M[f”(t)], t e [0;T], f(t)=e-bt3
У
читывая формулу (3.4) получаем:
(3.5)
А
необходимое условие экстремума: f”(t)-f’’’(t)=0, откуда получаем:
Далее вычисляем значения f’’(t) при t=t1, t=t2, t=0 и t=T, получаем:
f’’(t1)=1.5886 10-4
f’’(t2)=-1.6627 10-4
f’’(0)=0
f’’(T)=7.4782 10-6
Итак: M1,5886 10-4, откуда n=25.66; принимаем N=26.
Погрешность вычисления :
3.2 Вычисление интеграла I методом парабол
П
ри расчётах будем использовать теоретическую оценку погрешности с помощью правила Рунге. Для обеспечения заданной точности количество частей n, на которое следует разделить интервал интегрирования можно определить по формуле:
Нахождение М4 можно провести аналогично нахождению М2 в предыдущем пункте, но выражение для fIV(t) имеет довольно громоздкий вид. Поэтому правило Рунге – наиболее простой способ.
Обозначим через In и I2n значение интеграла I, полученное при разбиении промежутка интегрирования соответственно на n и 2n интервалов. Если выполнено равенство: |I2n-In| = 15 , то |I-I2n|=
Б
удем , начиная с n=2, удваивать n до тех пор, пока не начнёт выполняться неравенство (*1), тогда:
(3.6)
С
огласно формуле парабол (3.7):
Результаты вычислений сведём в таблицу:
n | In | I2n |
4 | 102.11 | |
8 | 101.61 | 0.5017 |
По формуле (3.7) I = 101,61 что в пределах погрешности совпадает со значением, полученным по методу трапеций
n=8 | n=4 | ||
ti (8) | y8 | ti (4) | y4 |
0 | 1 | 0 | 1 |
27.25 | 0.9864 | ||
54.5 | 0.8959 | 54.5 | 0.8959 |
81.75 | 0.6901 | ||
109 | 0.4151 | 109 | 0.4151 |
136.25 | 0.1796 | ||
163.5 | 0.0514 | 163.5 | 0.0514 |
190.75 | 0.0089874 | ||
218 | 0.00088179 | 218 | 0.00088179 |
4. Вычисление времени Т0 установления режима
4.1 Решение уравнения комбинированным методом
Время установления режима определяется по формулам (1.6) и (1.7).
Проведём сначала отделение корней. Имеем y = ctg(x) и y = Ax. Приведём уравнение к виду: A x sin(x)-cos(x) = 0. Проведём процесс отделения корня.
F(x) | -1 | -0.6285 | 0.4843 |
x | 0.01 | 0.05 | 0.1 |
т.е. с [0.01;0.05]
Убедимся, что корень действительно существует и является единственным на выбранном интервале изоляции.
f(a) f(b)<0 – условие существования корня выполняется
f’(x) на [a;b] – знакопостоянна: f’(x)>0 – условие единственности также выполняется. Проведём уточнение с погрешностью не превышающей
Строим касательные с того конца, где f(x) f”(x)>0
f
”(x)=(2A+1)cos(x) – A x sin(x). f”(x)>0 на (a;b), следовательно касательные строим справа, а хорды слева. Приближение корня по методу касательных:
В
ычисление ведём до того момента, пока не выполнится условие:
Результаты вычислений заносим в таблицу:
n | an | bn | f(an) | f(bn) |
0 | 0.05 | 0.1 | -0.6285 | 0.4843 |
1 | 0.07824 | 0.08366 | -0.0908 | 0.0394 |
2 | 0.08202 | 0.08207 | -9.1515 10-4 | 3.7121 10-4 |
3 | 0.08206 | 0.08206 | -8.4666 10-8 | 3.4321 10-8 |
Т0 = 72,7176 секунд.
4.2 Решение уравнения комбинированным методом
Приведём f(x) = 0 к виду x = (x). Для этого умножим обе части на произвольное число , неравное нулю, и добавим к обеим частям х:
X = x - f(x)
xx - A x sin(x) + cosx)
В качестве возьмём:
где М = max [f’(x)] на [a;b], а m = min [f’(x)] на [a’b]
В силу монотонности f’(x) на [a;b] имеем m = f’(а), М = f’(b). Тогда 0,045.
П
риближение к корню ищем по следующей схеме:
В
ычисление ведём до тех пор, пока не выполнится условие:
(q = max |’(x)| на [a’b])
’(x) на [a’b] монотонно убывает, поэтому максимум его модуля достигается на одном из концов.
’(0,05) = 0,3322 ’(0,1) = -0,3322, следовательно, q = 0.3322 < 1. В этом случае выполняется условие сходимости и получается последовательность:
i | xi | ( xi) | xi |
0 | 0.075 | 0.082392 | 0.00739 |
1 | 0.082392 | 0.082025 | 0.000367 |
2 | 0.082025 | 0.08206 | 3.54 10-5 |
3 | 0.08206 | 0.082057 | 3.33 10-6 |
4 | 0.082057 | 0.082057 | 3.15 10-7 |
Итак, с погрешностью, меньшей 10-4, имеем:
Т0 = 72,7176 с. , 0.03142
5. Решение краевой задачи
И
спользуем метод малого параметра. Краевую задачу запишем в виде:
(5.1)
В
ведя новую переменную y = (U - , запишем (5.1) в виде:
(5.2)
0.18L/2 =0.0193. В качестве малого параметра возьмём .
Тогда, подставив y(x) в уравнение (5.2) и перегруппировав члены при одинаковых степенях , получим:
(5.3)
О
граничимся двумя первыми членами ряда:
И
з (5.2) и (5.3) находим общее решение уравнения для y0:
где y0 с тильдой – частное решение данного неоднородного уравнения; y(1) и y(2) – линейно независимые решения однородного уравнения.
y0общ = 1 + c1ch(px)+c2sh(px), где p = 0.01953
0>