Мощность, выделяемая на резисторе, можно записать как:\[P=\fbox{$\frac{U_R^2}R$}=U_RI=U_R\dot Q_C=\frac1CU_R\dot U_C=\frac1CU_R\frac{\mathrm d}{\mathrm dt}\left[U-U_R\right]=-\frac1CU_R\dot U_R=\fbox{$-\frac1{2C}\frac{\mathrm d}{\mathrm dt}\left[U_R^2\right]$}.\]При этом мощность, рассеиваемая резистором в окружающее пространство, равна $\kappa\tau$, откуда:
Введём обозначения $c=a_1U^2$ и $\kappa=a_2 U^2$. Для поиска $a_1$ проведём касательную к графику $\tau(t)$ в точке $t=0$. Её угловой коэффициент будет равен $\fbox{$\frac{U^2}{cR_0}=\frac1{a_1R_0}$}$. Для поиска $a_2$ снимем несколько точек с "хвоста" графика и перестроим их в координатах $\ln\tau(t)$. Угловой коэффициент полученной прямой будет равен $\fbox{$-\frac\kappa c=-\frac{a_2}{a_1}$}$. Окончательный ответ:
Тепло, которое рассеется на термисторе, равно интегралу мощности теплообмена от $t=0$ до $t=+\infty$:\[Q_R=\fbox{$\kappa\int_0^{+\infty}\tau(\xi)~\mathrm d\xi$}.\]Работа источника за это время $A=QU=CU^2$, а изменение энергии конденсатора $W_C=\frac C2U^2$ $\implies$ на термисторе выделится $Q_R=A-W_C=\fbox{$\frac C2U^2$}$ $\implies$
Подставляя выражение для $\kappa$ из A2, имеем:
Из пунктов A1 и A2 получаем:\[P=\frac{U_R^2}R=c\dot\tau+\kappa\tau=\fbox{$\left[a_1\dot\tau+a_2\tau\right]U^2$}.\]
Вспомним, что мощность на термисторе равна $P=\frac C2\frac{\mathrm d}{\mathrm dt}\left[U_R^2\right]$, откуда в результате интегрирования получаем:\[\int_0^tP(\xi)~\mathrm d\xi=\frac C2\left[U^2-U_R^2(t)\right]=c\tau(t)+\kappa\int_0^{+\infty}\tau(\xi)~\mathrm d\xi\implies U_R=\sqrt{U^2-\frac2C\left[c\tau+\kappa\int_0^t\tau(\xi)~\mathrm d\xi\right]}.\]Подставляя выражение для $\kappa$ из A2, имеем:
Подставляя результаты A4 и A5 в выражение $P=\frac{U_R^2}R$, имеем:
Оригинальный график в максимальном разрешении можно найти в материалах к задаче.
Подставляем в выражения из A2:
По МНК нам нужно минимизировать следующую сумму:\[\mathcal S=\sum_i\left[A+B\ln R+C\ln^3R-\frac1{T_0+\tau}\right]^2.\]Для этого необходимо решить систему уравнений:
$\begin{cases}\frac{\partial\mathcal S}{\partial A}=0\\\frac{\partial\mathcal S}{\partial B}=0\\\frac{\partial\mathcal S}{\partial C}=0\end{cases}$
Беря производные и переписывая уравнения через средние по выборке, получим:\[\]
Выбор желаемого алгоритма решения оставим за читателем.
В стационарном режиме $\dot\tau=0$, поэтому:\[\frac{U^2}R=\kappa\tau=\kappa(T-T_0).\]Исключая отсюда $T$ с помощью уравнения Стейнхарта—Харта, имеем:
С помощью итеративных уравнений:
Нужно преобразовать уравнение к форме с $R$ в левой части так, чтобы обеспечить сходимость.
Параметризацией по $R$:
Куда проще можно поступить, если заметить, найти зависимость $U(R)$ очень легко. Вычислив пары точек, мы получим зависимость $R(U)$.
Оригинальные графики в максимальном разрешении можно найти в материалах к задаче.
С помощью итеративных уравнений:
При увеличении напряжения можно заметить, что сходимость при решении итеративным методом происходит всё медленнее и медленнее, пока по достижении некоторого напряжения сопротивление не "улетает" скачком в область очень малых значений (которые, как несложно проверить, соответствуют температурам порядка $\sim10^6~К$). Это – пробой термистора. Напряжение и ток, при котором он происходит:
Параметризацией по $R$:
При уменьшении $R$ при некотором значении $R_\mathrm {cr}=11.241~(256.5)~Ом$ ВАХ “разворачивается”, что, как известно, соответствует неустойчивому участку, а потому – пробою термистора. Отсюда сразу находим напряжение и ток, при котором он происходит:
С помощью итеративных уравнений:
Нужно преобразовать уравнение к форме с $R$ в левой части так, чтобы обеспечить сходимость. Температуру, которой соответствует конкретная точка на ВАХ, можно найти из уравнения Стейнхарта—Харта. При $T=T_{\max}$ ВАХ обрывается.
Параметризацией по $R$:
Сопротивление термистора $R_\min$ находим из уравнения Стейнхарта—Харта. Строим ВАХ от этого значения до найденного в предыдущем пункте.
Оригинальные графики в максимальном разрешении можно найти в материалах к задаче.
С помощью итеративных уравнений:
Температуру, которой соответствует конкретная точка на ВАХ, можно найти из уравнения Стейнхарта—Харта.
Параметризацией по $\tau$:
Ответ получается непосредственной подстановкой в уравнение.