next up previous contents
Next: Delay and Dispersion Correction Up: Mathematical Analysis Previous: Double-Weighted Integration Method

Triple-Weighted Integration Method

In order to model the complete two-compartment system, we must be able to solve equation (4). In section 2.2, we solved the equation by multiplying by two different weights and then dividing, and we can take a similar approach with the full two-compartment equation. We can weight equation (4) with three different weights and then integrate:


 \begin{displaymath}\int_{0}^{T} w_{1}A^{*}(t)dt = K_{1} \int_{0}^{T} w_{1}
\frac...
...{0}^{T}w_{1} \frac{\int_{T_1}^{T_2}
C_{a}(u) du}{T_2 - T_1} dt
\end{displaymath} (7)


 \begin{displaymath}\int_{0}^{T} w_{2}A^{*}(t)dt = K_{1} \int_{0}^{T} w_{2}
\frac...
...{0}^{T}w_{2} \frac{\int_{T_1}^{T_2}
C_{a}(u) du}{T_2 - T_1} dt
\end{displaymath} (8)


 \begin{displaymath}\int_{0}^{T} w_{3}A^{*}(t)dt = K_{1} \int_{0}^{T} w_{3}
\frac...
...{0}^{T}w_{3} \frac{\int_{T_1}^{T_2}
C_{a}(u) du}{T_2 - T_1} dt
\end{displaymath} (9)

By multiplying equation (7) by $V_{0}
\int_{0}^{T} w_{3} {(\int_{T_1}^{T_2}C_{a}(u)du)}/{(T_2 - T_1)} dt$and equation (9) by $V_{0} \int_{0}^{T} w_{1}
{(\int_{T_1}^{T_2}C_{a}(u)du)}/{(T_2 - T_1)} dt$, and then subtracting the two, we may eliminate the V0 term. A similar operation can be performed on equation (8) and equation (9). This leaves two equations that do not contain V0. They may then be divided to produce:


 \begin{displaymath}\hspace*{-1.0cm}
\frac{
\begin{array}{lcr}
\multicolumn{2}{l}...
...(T_2 - T_1)}\, dt} \end{array} \end{array} \!\!\!\!
\right\} }
\end{displaymath} (10)

The K1 term cancels out of both the numerator and denominator of equation (10), leaving an equation that only involves k2. As with the equation in section 2.2, this is very difficult to solve for k2. Therefore, a look-up table was again used. Once the table matching values of k2 with values of the right hand side of equation (10) has been created, we may evaluate k2 through simple lookup. With the k2 data computed, finding K1 is simply a matter of evaluating either the numerator or denominator of equation (10) without cancelling K1. With both K1 and k2 known, we may find V0 by evaluating equation (4).


next up previous contents
Next: Delay and Dispersion Correction Up: Mathematical Analysis Previous: Double-Weighted Integration Method
Mark Wolforth <wolforth@bic.mni.mcgill.ca>
Greg Ward <greg@bic.mni.mcgill.ca>
Sean Marrett <sean@bic.mni.mcgill.ca>