Исследование распределения температуры в тонком цилиндрическом стержне
Решая эту систему методом Гаусса, получим:
(2.8)
Чем лучше функция регрессии описывает эксперимент, тем меньше для неё должна быть оценка дисперсии отдельного измерения Ui, т.к. при плохом выборе функции в дисперсию войдут связанные с этим выбором дополнительные погрешности. Поэтому для того, чтобы сделать выбор между функциями U(x) и U(1)(x) нужно проверить значимость различия между соответствующими оценками дисперсии, т.е. проверить гипотезу:
Н0 – альтернативная гипотеза
Т.е. проверить, значимо ли уменьшение дисперсии при увеличении степени многочлена.
В качестве статического критерия рассмотрим случайную величину, равную:
(2.9)
имеющую распределение Фишера с(r ; r1) степенями свободы. Выбираем уровень распределения Фишера, находим критическое значение F*, удовлетворяющее равенству: p(F>F*=
В нашем случае F=349.02, а F*=10,13.
Если бы выполнилось практически невозможное соотношение F>F, имевшее вероятность 0,01, то гипотезу Н0 пришлось бы отклонить. Но в нашем случае можно ограничиться многочленом
, коэффициенты в котором неодинаковы.
3. Нахождение коэффициента теплопроводности .
Коэффициент вычислим по формуле (1.5), обозначим:
(3.1)
Определим допустимую абсолютную погрешность величины интеграла 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)-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.
Далее вычислим интеграл I:
Погрешность вычисления :
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) 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)