(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 15203, 503] NotebookOptionsPosition[ 14122, 461] NotebookOutlinePosition[ 14465, 476] CellTagsIndexPosition[ 14422, 473] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{"H", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}]}], " ", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"HT", " ", "=", " ", RowBox[{"Transpose", "[", "H", "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.4814591112804203`*^9, 3.4814591340935044`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"P0", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", SuperscriptBox["E", FractionBox[ RowBox[{"-", "r12"}], "d"]], ",", SuperscriptBox["E", FractionBox[ RowBox[{"-", "r13"}], "d"]]}], "}"}], ",", RowBox[{"{", RowBox[{ SuperscriptBox["E", FractionBox[ RowBox[{"-", "r12"}], "d"]], ",", "1", ",", SuperscriptBox["E", FractionBox[ RowBox[{"-", "r23"}], "d"]]}], "}"}], ",", RowBox[{"{", RowBox[{ SuperscriptBox["E", FractionBox[ RowBox[{"-", "r13"}], "d"]], ",", SuperscriptBox["E", FractionBox[ RowBox[{"-", "r23"}], "d"]], ",", "1"}], "}"}]}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.4814591369060764`*^9, 3.481459158031617*^9}, { 3.4814595200897274`*^9, 3.481459557075286*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", "P0", "]"}]], "Input", CellChangeTimes->{{3.482051782914525*^9, 3.482051786883301*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r12", "d"]}]], SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r13", "d"]}]]}, { SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r12", "d"]}]], "1", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r23", "d"]}]]}, { SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r13", "d"]}]], SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r23", "d"]}]], "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.4820517873051786`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"H", " ", ".", " ", "P0", " ", ".", " ", "HT"}], "//", "MatrixForm"}]], "Input", CellChangeTimes->{{3.481459578779105*^9, 3.481459591685768*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r23", "d"]}]]}, { SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r23", "d"]}]], "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.4814595832323728`*^9, 3.4814595919826527`*^9}, 3.482051814774104*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Inverse", "[", " ", RowBox[{"H", ".", "P0", ".", "HT"}], " ", "]"}], "//", "MatrixForm"}]], "Input", CellChangeTimes->{{3.4814595934983263`*^9, 3.48145960581122*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox["1", RowBox[{"1", "-", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]}]]}]], RowBox[{"-", FractionBox[ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r23", "d"]}]], RowBox[{"1", "-", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]}]]}]]}]}, { RowBox[{"-", FractionBox[ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r23", "d"]}]], RowBox[{"1", "-", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]}]]}]]}], FractionBox["1", RowBox[{"1", "-", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]}]]}]]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.481459601545459*^9, 3.4814596061237307`*^9}, 3.4820518832845263`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"K", "=", RowBox[{"Simplify", "[", RowBox[{"P0", " ", ".", " ", "HT", " ", ".", " ", RowBox[{"Inverse", "[", " ", RowBox[{"H", ".", "P0", ".", "HT"}], " ", "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"K", "//", "MatrixForm"}]}], "Input", CellChangeTimes->{{3.481468912833393*^9, 3.481468922786645*^9}, { 3.482053384410884*^9, 3.4820534122235613`*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"r12", "+", "r13", "-", "r23"}], "d"]}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", SuperscriptBox["\[ExponentialE]", FractionBox["r12", "d"]]}], "+", SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"r13", "+", "r23"}], "d"]]}], ")"}]}], RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]]}]], FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"r12", "+", "r13", "-", "r23"}], "d"]}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", SuperscriptBox["\[ExponentialE]", FractionBox["r13", "d"]]}], "+", SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"r12", "+", "r23"}], "d"]]}], ")"}]}], RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]]}]]}, {"1", "0"}, {"0", "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{ 3.481468925599181*^9, 3.48205195023808*^9, {3.4820533899265437`*^9, 3.4820534128798156`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "multiply", " ", "by", " ", "a", " ", "constant", " ", "to", " ", "look", " ", "at", " ", "the", " ", "simplified", " ", "structure"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1", "-", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]}]]}], ")"}], " ", "K"}], "//", "Simplify"}], "//", "MatrixForm"}]}]], "Input", CellChangeTimes->{{3.4820523677251267`*^9, 3.48205240808476*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r12", "d"]}]], "-", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"r13", "+", "r23"}], "d"]}]]}], RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["r13", "d"]}]], "-", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"r12", "+", "r23"}], "d"]}]]}]}, { RowBox[{"1", "-", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]}]]}], "0"}, {"0", RowBox[{"1", "-", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]}]]}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.482052403490981*^9, 3.4820524083035116`*^9}, 3.4820534371299706`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "compute", " ", "the", " ", "posteriori", " ", "state", " ", "estimate"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"xhatPlus", "=", RowBox[{ RowBox[{"K", " ", ".", " ", RowBox[{"{", RowBox[{ RowBox[{"{", "p2", "}"}], ",", RowBox[{"{", "p3", "}"}]}], "}"}]}], "//", "Simplify"}]}]}]], "Input", CellChangeTimes->{{3.481468934114915*^9, 3.4814689458650656`*^9}, { 3.48205322914789*^9, 3.4820532635699854`*^9}, {3.4820534637082663`*^9, 3.48205346506765*^9}, {3.4820537233193026`*^9, 3.482053724975563*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"r12", "+", "r13", "-", "r23"}], "d"]}]], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", SuperscriptBox["\[ExponentialE]", FractionBox["r12", "d"]]}], " ", "p2"}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"r13", "+", "r23"}], "d"]], " ", "p2"}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", FractionBox["r13", "d"]], " ", "p3"}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"r12", "+", "r23"}], "d"]], " ", "p3"}]}], ")"}]}], RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]]}]], "}"}], ",", RowBox[{"{", "p2", "}"}], ",", RowBox[{"{", "p3", "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.4814689474275856`*^9, 3.482053246351125*^9, {3.482053457801978*^9, 3.4820534653801517`*^9}, 3.4820537258193183`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Simplify", "[", RowBox[{ RowBox[{"(", RowBox[{"1", "-", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"2", " ", "r23"}], "d"]}]]}], ")"}], RowBox[{"xhatPlus", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4820537268193245`*^9, 3.4820537545695024`*^9}, 3.482053803882318*^9}], Cell[BoxData[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"r12", "+", "r13", "+", "r23"}], "d"]}]], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", SuperscriptBox["\[ExponentialE]", FractionBox["r12", "d"]]}], " ", "p2"}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"r13", "+", "r23"}], "d"]], " ", "p2"}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", FractionBox["r13", "d"]], " ", "p3"}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"r12", "+", "r23"}], "d"]], " ", "p3"}]}], ")"}]}]], "Output", CellChangeTimes->{{3.4820537291630898`*^9, 3.482053755147631*^9}, 3.482053804788574*^9}] }, Open ]] }, WindowSize->{616, 750}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (February 18, 2009)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[545, 20, 448, 13, 52, "Input"], Cell[996, 35, 946, 31, 44, "Input"], Cell[CellGroupData[{ Cell[1967, 70, 130, 2, 31, "Input"], Cell[2100, 74, 1225, 37, 104, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3362, 116, 183, 4, 31, "Input"], Cell[3548, 122, 869, 25, 78, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4454, 152, 209, 5, 31, "Input"], Cell[4666, 159, 1677, 51, 122, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6380, 215, 428, 10, 52, "Input"], Cell[6811, 227, 1990, 59, 118, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8838, 291, 600, 17, 68, "Input"], Cell[9441, 310, 1516, 46, 104, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10994, 361, 616, 15, 52, "Input"], Cell[11613, 378, 1230, 36, 73, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12880, 419, 431, 13, 47, "Input"], Cell[13314, 434, 792, 24, 47, "Output"] }, Open ]] } ] *) (* End of internal cache information *)