(* Content-type: application/vnd.wolfram.cdf.text *) (*** Wolfram CDF File ***) (* http://www.wolfram.com/cdf *) (* CreatedBy='Mathematica 9.0' *) (*************************************************************************) (* *) (* The Mathematica License under which this file was created prohibits *) (* restricting third parties in receipt of this file from republishing *) (* or redistributing it by any means, including but not limited to *) (* rights management or terms of use, without the express consent of *) (* Wolfram Research, Inc. For additional information concerning CDF *) (* licensing and redistribution see: *) (* *) (* www.wolfram.com/cdf/adopting-cdf/licensing-options.html *) (* *) (*************************************************************************) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 1063, 20] NotebookDataLength[ 76802, 2081] NotebookOptionsPosition[ 73590, 1959] NotebookOutlinePosition[ 74052, 1977] CellTagsIndexPosition[ 74009, 1974] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Reconstrucci\[OAcute]n de datos reales a partir de los observados con el m\ \[EAcute]todo de M\[AAcute]xima Entrop\[IAcute]a (MEM)\ \>", "Title"], Cell[CellGroupData[{ Cell["\<\ Trabajo de curso de doctorado realizado por Victor de Buen (Bayes \ Inference), 2001\ \>", "Author", CellChangeTimes->{{3.6007630446317153`*^9, 3.6007630456157312`*^9}, { 3.6007630803268604`*^9, 3.600763131447829*^9}, {3.6007632938352556`*^9, 3.600763383647336*^9}, {3.600764059686153*^9, 3.600764133153102*^9}, { 3.60076422444149*^9, 3.6007642287540026`*^9}}], Cell[TextData[{ "Realizado dentro del programa de doctorado de Modelizaci\[OAcute]n y \ Simulaci\[OAcute]n impartido por Guillermo S\[AAcute]nchez (", ButtonBox["http://diarium.usal.es/guillermo", BaseStyle->"Hyperlink", ButtonData->{ URL["http://diarium.usal.es/guillermo"], None}, ButtonNote->"http://diarium.usal.es/guillermo"], "). Universidad de Salamanca. \nEl trabajo original est\[AAcute] hecho con ", StyleBox["Mathematica", FontSlant->"Italic"], " 3. He hecho (G. S\[AAcute]nchez) los cambios imprescindibles para que se \ ejecute en ", StyleBox["Mathematica", FontSlant->"Italic"], " 9." }], "Text", CellChangeTimes->{{3.6007630446317153`*^9, 3.6007630456157312`*^9}, { 3.6007630803268604`*^9, 3.600763131447829*^9}, {3.6007632938352556`*^9, 3.600763383647336*^9}, {3.600764059686153*^9, 3.600764215536167*^9}}], Cell["\<\ Empezaremos por generar una figura con 2 picos,de la que obtendremos 64 \ puntos.Esto puede corresponder a las emisiones una muestra radiactiva.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"n", " ", "=", " ", "64"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f", "[", "x_", "]"}], ":=", " ", RowBox[{"N", "[", RowBox[{ RowBox[{"2", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", " ", SuperscriptBox["x", "2"]}]}], ")"}]}], "+", RowBox[{"2", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", "3"}], ")"}], "2"]}]}], ")"}]}], "+", "0.7"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"underlying", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"f", "[", "x", "]"}], ",", " ", RowBox[{"{", RowBox[{"x", ",", " ", RowBox[{"-", "5"}], ",", " ", "10", ",", " ", RowBox[{"15", "/", RowBox[{"(", RowBox[{"n", " ", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"ourhope", " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"underlying", ",", " ", RowBox[{"Joined", "\[Rule]", "True"}], ",", " ", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "3"}], "}"}]}], ",", " ", RowBox[{"PlotStyle", " ", "->", " ", RowBox[{"RGBColor", "[", RowBox[{"0", ",", " ", "0", ",", " ", "1"}], "]"}]}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.600763408069225*^9, 3.600763408616103*^9}}], Cell[BoxData[ GraphicsBox[{{}, {RGBColor[0, 0, 1], LineBox[CompressedData[" 1:eJw10nsslXEcx/EjRCWT0RSSy0xCnHLp4JyP437cWrKwYloWJYcdMZfOnsyW jYn+oLlMWG6zEWZpSxe306Y0ly3LZJkhwsk1l6T5/p49z569nvc5v/2ePV/T W9KrsYd4PB72r//3g2NZ1F/SXNNgOSuiB4hY8bR8Hcysifw9wcqzLGYdONfU uWe8ZNZDsWZ89MwiswEyQyd6xp3myEZo75xSrOYxn0XYtq1LrZLZDGZDsQFV CT/JFojUmJy4oTpPtsT0uGOfsovZCpz/mCytaYFsjfq1C96u33+RbWDUGBMm TVki2+GpZGq3eWCZbA/FgNb6kuQ32QF50rD1uMRVMh9dO5Z/owo3yBfx6F2T 7LrVDvkSXM6b8MOFKjiwI8rzQ1MXWtXITtCtkw/aJh8mO6NF0ruRssq6C9qD Kupf9bL/XwaXxBdEjO7S+gLI1QLizbb/kF1hoyVfGlKw/bihtUlP1vJjnewO x+NJdrqGrAvRof5Rlv5wkyxC29KH8WNB22Rg4LF66dvZvQNzwHKa9EnqSVXa jwcMRkXy6Gn1A3MeyK3Q/6YwZO8jhs1waZxHDP2eE8OPJ9Oe6Kb1eJ7Ils6N RBhv0fqecCup0v58ju3XC4ap4t62jBXqXij5YlKtUa2k7o3UroS2yGv0vThv JJ7q2T8XqfsgN8s++XkkfW/OB+1F8tN3fNg8+ELljcC0O5HmhfOF7ljNC0E3 my8/9FdqZ+S4kDk/NJxIz+roY/Poj/l88/cWcWTOH7v6d60bDViXoHA3XOgw SvPOSZAt6lSLLWPzH4DNnMTYzHjWA6A0v2esI2Q9EAXlW+1HDVgPBMTz03nr M9SDUOVbcGbyK5kLwsSRteIHXawHQxgVujhSy3owpp3FpWVFrIdA0cZvSJOz HoLB25/kFvdZv4LK4VxJ3M0Z0T9LGC70 "]]}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{{0, 64.}, {0, 3}}, PlotRangeClipping->True, PlotRangePadding->{{1.28, 1.28}, {0., 0.}}]], "Output", CellChangeTimes->{3.6007634926557436`*^9, 3.6007636460985813`*^9, 3.600763954252386*^9}] }, Open ]], Cell["\<\ Supongamos que la muestra anterior es medida con un detector multicanal con \ resoluci\[OAcute]n gausiana que simulamos por:\ \>", "Text"], Cell[BoxData[ RowBox[{ RowBox[{"m", " ", "=", " ", "32"}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"resolution", " ", "=", " ", RowBox[{"Fourier", "[", RowBox[{"RotateRight", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", " ", "5"}], SuperscriptBox["x", "2"]}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", " ", RowBox[{"-", "4"}], ",", " ", "4", ",", " ", RowBox[{"8", "/", RowBox[{"(", RowBox[{"n", " ", "-", "1"}], ")"}]}]}], "}"}]}], "]"}], ",", "m"}], "]"}], "]"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"R", "[", "x_List", "]"}], ":=", " ", RowBox[{"Chop", "[", RowBox[{"InverseFourier", "[", RowBox[{"resolution", " ", RowBox[{"Fourier", "[", "x", "]"}]}], " ", "]"}], "]"}]}], ";"}]], "Input"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Rt", "[", "x_List", "]"}], ":=", " ", RowBox[{"Chop", "[", RowBox[{"Fourier", "[", RowBox[{"resolution", " ", RowBox[{"InverseFourier", "[", "x", "]"}]}], " ", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Observed", " ", "=", RowBox[{"Rt", "[", RowBox[{"R", "[", "underlying", "]"}], "]"}]}], ";"}]}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Obs", " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"Observed", ",", " ", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", " ", "4"}], "}"}]}], ",", " ", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"PlotStyle", " ", "->", " ", RowBox[{"RGBColor", "[", RowBox[{"0", ",", " ", "0", ",", " ", "1"}], "]"}]}]}], "]"}]}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.6007634151317806`*^9, 3.600763415663047*^9}}], Cell[BoxData[ GraphicsBox[{{}, {RGBColor[0, 0, 1], LineBox[CompressedData[" 1:eJw102tIU2Ecx/GlBYEU89bUEvKSSxPzMi+Vup/XeVcmaL7Q1JLQtswSs5U5 TCM1rKVUXlpJIaYR01LKSxqBE7MSX4RpVjPLsrzNBWlOM3r+Bw6Hz/nynOfF +T826dniDAMOh4P1+9/z/zUvTH6juukuHReyF0ivNLp3o4m8GbWtbisFOjIX zbtqjqRGfmA2QxZylr1byBaIzTDetuz4kXkHzqVxAoy7yDsxsXJyWObzidkW cXYOyWtNZHtk1E3mPfTXMDug7NHgwew18m6MWVvaa/5MMDvhmPnjNr/gSWZn VIoMy+qnvjC7QKY1jJHOTTG7orHx1Am3pGlmN8yq9KebuDPM7jgcMZCU8HOO 2QOa+Gyj1sQFZgEaBnPWIku1zJ4wkfPD7mQuMnshmOe75D1B9saALn1BOkv2 wVJgichaSd6Hq/rVrlET8n6cKc9vzsqn7x+AzqMnt2GZ9vdFX0uKZ3kn2Q9W Ytug3pdkf8wXraRs5NN6IfrtNdPmbWSAm+ko0saw/eRA91Fpfoma9g/A5VYJ V/+LegDKdfxbLj3UA/FV1quQbKceiEtayW8ej74fBNvcYBfO63nWg9BRrXg3 2jfLejBqrDus4hd/sL7uF2PH9yR+Zz0E6u6t3iJT9r/kIRAPWQqdBPQ/Q1Fa 8Naxr/0z66EY71cPqytoHkRIU0xlej1g8yIX4VlH2VPxGs1XGM4fGqlQypjl YWiofaUsNKUejpoCnsShmM2rPByChoVriRtofiMwfFtrUXKdzbs8AnXTWS2r fjT/kRiJrnN4vsjOizwSoaq7CQZtdH6iECJQPpkppB6FYpv6sxfiqUfDbkup c9le6tGYLshbLTSjHgMXcdXFfEPqMfim8FeN6t+zHov7V9qHDDZRj0WRa2qA xpLWx0Fbxa/u9B0X/gWinzb+ "]]}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{{0, 64.}, {0, 4}}, PlotRangeClipping->True, PlotRangePadding->{{1.28, 1.28}, {0., 0.}}]], "Output", CellChangeTimes->{3.6007634932494736`*^9, 3.600763646223561*^9, 3.6007639544086227`*^9}] }, Open ]], Cell["\<\ Simularemos que la se\[NTilde]al anterior est\[AAcute] enmascarada con un \ ruido que simularemos por n\[UAcute]meros aleatorios generados con una funci\ \[OAcute]n normal\ \>", "Text"], Cell[BoxData[ RowBox[{ RowBox[{"s", "=", " ", "0.01"}], ";"}]], "Input"], Cell[BoxData[ RowBox[{" ", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Random", "[", RowBox[{"NormalDistribution", "[", RowBox[{"0", ",", "s"}], " ", "]"}], "]"}], ",", " ", RowBox[{"{", "n", "}"}]}], "]"}], ";"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"noise", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"Random", "[", RowBox[{"NormalDistribution", "[", RowBox[{"0", ",", " ", "s"}], "]"}], "]"}], ",", " ", RowBox[{"{", "n", "}"}]}], "]"}]}], ";"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"observed", " ", "=", " ", RowBox[{"Observed", "+", "noise"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"obs", " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"observed", ",", RowBox[{"Joined", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "4"}], "}"}]}], ",", " ", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"PlotStyle", " ", "->", " ", RowBox[{"RGBColor", "[", RowBox[{"1", ",", " ", "0", ",", " ", "0"}], "]"}]}]}], "]"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.6007634190068045`*^9, 3.600763426647458*^9}}], Cell[BoxData[ GraphicsBox[{{}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJw103lIFGEYx/FVIy3vMjLUMrMtRGstr1xtf66364WlQbrg2UFtihGFtDa4 i0lKuhhRpLW2CZl2ebBqiCKpm5RS/tGh5rWaqamZhSlKRe8zMAyf+TLvCzPP 7EzNjMsw5PF4+Hv+u/4/5kXn3c62J5UNiNgNWMQeW1GVDzGbILL596t4s8/M VvhU6ZrlVky2gVocXirtItti3WT2ePw1et4eeTb8xCVTsiOcLK4EyMfITmjU q1L1BsPMztCjRWqWOcLMh6DefsNs8CjzXrRrChSB1mPMLnhYqektaya7wndy 05rUeIJ5H3av+G9sm/7CLICDRaxa2zPF7A5bfcKAp+k35gNYydJJ7jrPMx+E yGWPpfTEd2YPuGn7ky/VLDB7wtQo1/JB6A9mLygM15tlbF5k9ka65y1l9hB1 H7zQn5oWdZEPIeVx8arLLrIvbpxzsPaPoPWFeHIyf3Wlg/b3Q3/VT0HfM7I/ 8F6S4/WBfBgGz5cieu3oeRG28AsX5AlkYG7/+JuXzmw/DjB0GDZImaD9A9DY ndCWqKUegLqYyw2OHdTFmGvoPL5mTl0MaV7P8qqQ1g9Eae47c8UMe39cIPoF 18tlmGU9CEdHbl+o002zHoReLyvrEulX1oMhn/FWoIp9Py4Y7YWvVR+jxlkP wR0Ttax2jn1vLgTtF++nS5Q0H6HQCQtqrt5j88OFQp5bkr9sSfMUhooj+VOy JjZvXBi0qh7XiiKav3Bo8rYvtdD8cuEQqYOMBW8HWY9ApvjpL+Eom3cuAmla Tbm8iboE5nZ9VkVbqUuQlLLDL9KGeiTiph4l26VRj0S1rldmnEM9CrJO95ba bvY/clEYuunVekZCPRqt5YqcSh9mLholRbUingf1GFSnabxPN1OPQZORvVKZ RT0Wyu56/rbFQdEfG94vVQ== "]]}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{{0, 64.}, {0, 4}}, PlotRangeClipping->True, PlotRangePadding->{{1.28, 1.28}, {0., 0.}}]], "Output", CellChangeTimes->{3.6007634934056864`*^9, 3.600763646333022*^9, 3.6007639544711075`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"Obs", ",", "obs"}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{{}, {RGBColor[0, 0, 1], LineBox[CompressedData[" 1:eJw102tIU2Ecx/GlBYEU89bUEvKSSxPzMi+Vup/XeVcmaL7Q1JLQtswSs5U5 TCM1rKVUXlpJIaYR01LKSxqBE7MSX4RpVjPLsrzNBWlOM3r+Bw6Hz/nynOfF +T826dniDAMOh4P1+9/z/zUvTH6juukuHReyF0ivNLp3o4m8GbWtbisFOjIX zbtqjqRGfmA2QxZylr1byBaIzTDetuz4kXkHzqVxAoy7yDsxsXJyWObzidkW cXYOyWtNZHtk1E3mPfTXMDug7NHgwew18m6MWVvaa/5MMDvhmPnjNr/gSWZn VIoMy+qnvjC7QKY1jJHOTTG7orHx1Am3pGlmN8yq9KebuDPM7jgcMZCU8HOO 2QOa+Gyj1sQFZgEaBnPWIku1zJ4wkfPD7mQuMnshmOe75D1B9saALn1BOkv2 wVJgichaSd6Hq/rVrlET8n6cKc9vzsqn7x+AzqMnt2GZ9vdFX0uKZ3kn2Q9W Ytug3pdkf8wXraRs5NN6IfrtNdPmbWSAm+ko0saw/eRA91Fpfoma9g/A5VYJ V/+LegDKdfxbLj3UA/FV1quQbKceiEtayW8ej74fBNvcYBfO63nWg9BRrXg3 2jfLejBqrDus4hd/sL7uF2PH9yR+Zz0E6u6t3iJT9r/kIRAPWQqdBPQ/Q1Fa 8Naxr/0z66EY71cPqytoHkRIU0xlej1g8yIX4VlH2VPxGs1XGM4fGqlQypjl YWiofaUsNKUejpoCnsShmM2rPByChoVriRtofiMwfFtrUXKdzbs8AnXTWS2r fjT/kRiJrnN4vsjOizwSoaq7CQZtdH6iECJQPpkppB6FYpv6sxfiqUfDbkup c9le6tGYLshbLTSjHgMXcdXFfEPqMfim8FeN6t+zHov7V9qHDDZRj0WRa2qA xpLWx0Fbxa/u9B0X/gWinzb+ "]]}, {}}, {{}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJw103lIFGEYx/FVIy3vMjLUMrMtRGstr1xtf66364WlQbrg2UFtihGFtDa4 i0lKuhhRpLW2CZl2ebBqiCKpm5RS/tGh5rWaqamZhSlKRe8zMAyf+TLvCzPP 7EzNjMsw5PF4+Hv+u/4/5kXn3c62J5UNiNgNWMQeW1GVDzGbILL596t4s8/M VvhU6ZrlVky2gVocXirtItti3WT2ePw1et4eeTb8xCVTsiOcLK4EyMfITmjU q1L1BsPMztCjRWqWOcLMh6DefsNs8CjzXrRrChSB1mPMLnhYqektaya7wndy 05rUeIJ5H3av+G9sm/7CLICDRaxa2zPF7A5bfcKAp+k35gNYydJJ7jrPMx+E yGWPpfTEd2YPuGn7ky/VLDB7wtQo1/JB6A9mLygM15tlbF5k9ka65y1l9hB1 H7zQn5oWdZEPIeVx8arLLrIvbpxzsPaPoPWFeHIyf3Wlg/b3Q3/VT0HfM7I/ 8F6S4/WBfBgGz5cieu3oeRG28AsX5AlkYG7/+JuXzmw/DjB0GDZImaD9A9DY ndCWqKUegLqYyw2OHdTFmGvoPL5mTl0MaV7P8qqQ1g9Eae47c8UMe39cIPoF 18tlmGU9CEdHbl+o002zHoReLyvrEulX1oMhn/FWoIp9Py4Y7YWvVR+jxlkP wR0Ttax2jn1vLgTtF++nS5Q0H6HQCQtqrt5j88OFQp5bkr9sSfMUhooj+VOy JjZvXBi0qh7XiiKav3Bo8rYvtdD8cuEQqYOMBW8HWY9ApvjpL+Eom3cuAmla Tbm8iboE5nZ9VkVbqUuQlLLDL9KGeiTiph4l26VRj0S1rldmnEM9CrJO95ba bvY/clEYuunVekZCPRqt5YqcSh9mLholRbUingf1GFSnabxPN1OPQZORvVKZ RT0Wyu56/rbFQdEfG94vVQ== "]]}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{{0, 64.}, {0, 4}}, PlotRangeClipping->True, PlotRangePadding->{{1.28, 1.28}, {0., 0.}}]], "Output", CellChangeTimes->{3.6007634934212985`*^9, 3.6007636463641877`*^9, 3.6007639545179834`*^9}] }, Open ]], Cell[TextData[{ "Para reconstruir los datos originales aparentemente ser\[IAcute]a \ suficiente utilizar la trasformada inversa de Fourier (", StyleBox["InverseFourier[observed/resolution]", FontFamily->"Courier"], ") sin embargo falla debido a que la funci\[OAcute]n de convoluci\[OAcute]n \ tiene valores muy peque\[NTilde]os para frecuencias altas. Si dividimos en \ intervalos peque\[NTilde]os y eliminamos las altas frecuencias el tema puede \ ir mejor" }], "Text"], Cell[BoxData[ RowBox[{ RowBox[{"crudeinverse", " ", "=", " ", RowBox[{"Map", "[", RowBox[{ RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Abs", "[", "#", "]"}], ">", "0.1"}], ",", RowBox[{"1", "/", "#"}], ",", "0"}], "]"}], "&"}], ",", " ", "resolution"}], "]"}]}], ";"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"recon", " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Chop", "[", RowBox[{"Fourier", "[", RowBox[{ RowBox[{"InverseFourier", "[", "observed", "]"}], " ", "crudeinverse"}], "]"}], "]"}], " ", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"PlotStyle", " ", "->", " ", RowBox[{"RGBColor", "[", RowBox[{"1", ",", " ", "0", ",", " ", "0"}], "]"}]}]}], "]"}]}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{3.600763447381892*^9}], Cell[BoxData[ GraphicsBox[{{}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJw10gtIU2EYxnEVCgPzRplNi2XDtFTUNe+2Z/M2N2+ZUaKJClLJVEQrllLH GgtNTTJvRKmZYakQImLeho0wMZMUisqZttRZqSuzTfOS0fcdOBx+58/hvPC9 B1KzYtNMjIyMsHX/e/6/dPwrDwVRnHINn7xA6WR2y/oktSnuLLd22yxQW6Lz o61K1ke9C02rVUvXL1Hbgmf3mHM3htoem2dKPgzJqNmIS1bEXN72hdgBN0tX k5VT1BzIDQnGuTPTxI4oXA9eM/yYIXaCLjFLfmi7lvgwyqxbnEUn54hdIFTd 6GrkfSN2Q5Ri5z3BxjyxO6Dd2xQyryP2QELdrXRW3BKxJ4ZMRxzP6peJuXBm id/55hmIjyJEILLwy1sl5sFYbGey0fiH2AvNDbrBzDFqb2RXP+IpLah9IH3h WfCqeoXYFxdH49M0I3piPxR36ZuTOL+J/VGSIlWz134RB6C2KzAs+D51IEqe trYzB+m8x9BeURvqmkW/52M5sX2cxaXzA6rkdL3pJvk/A7zJKOxVxdP5BEjx LOVWThAzAhSrpx6wTtEuxOuGNvgP0u+FqBxzt/PPpfMHgVNR5ObRS+ZjgqAd zc9JtPpJejBsaiXKjbJF0oMRNFz+6bzDd9JDoJhK9Zt8T86TCYGrep9y9gI9 71BYr2j71VWzpIdu7cO5t3Nf6X6EwWxeVi8MIGbCwBq2l3pdpfskwoDLaqhT P9k3RgQzg9xQs4fuXzhO11X2cG+TfWXCcS3/eFUnj+6vGG3ZHeNlG59JFyOd 7f58YJrYSIIiQU2PUkO7BB38HEgXaY9ABKdKVL9GewRazOMsfddpj4Rswruo Zob2SPgcsdrBeUZ7FDozM7olctqjoFjwDnyZSns0npjnxjZl0h6NJA/bEwVt tMegb7+yejdbw/8LV+847g== "]]}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0.5745583720849516}, Method->{}, PlotRange->{{0, 64.}, {0.5745583720849516, 1.7418424619565018`}}, PlotRangeClipping->True, PlotRangePadding->{{1.28, 1.28}, {0.023345681797431003`, 0.023345681797431003`}}]], "Output", CellChangeTimes->{3.600763493483892*^9, 3.600763646411065*^9, 3.6007639545492477`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"recon", ",", "ourhope"}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{{}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJw10gtIU2EYxnEVCgPzRplNi2XDtFTUNe+2Z/M2N2+ZUaKJClLJVEQrllLH GgtNTTJvRKmZYakQImLeho0wMZMUisqZttRZqSuzTfOS0fcdOBx+58/hvPC9 B1KzYtNMjIyMsHX/e/6/dPwrDwVRnHINn7xA6WR2y/oktSnuLLd22yxQW6Lz o61K1ke9C02rVUvXL1Hbgmf3mHM3htoem2dKPgzJqNmIS1bEXN72hdgBN0tX k5VT1BzIDQnGuTPTxI4oXA9eM/yYIXaCLjFLfmi7lvgwyqxbnEUn54hdIFTd 6GrkfSN2Q5Ri5z3BxjyxO6Dd2xQyryP2QELdrXRW3BKxJ4ZMRxzP6peJuXBm id/55hmIjyJEILLwy1sl5sFYbGey0fiH2AvNDbrBzDFqb2RXP+IpLah9IH3h WfCqeoXYFxdH49M0I3piPxR36ZuTOL+J/VGSIlWz134RB6C2KzAs+D51IEqe trYzB+m8x9BeURvqmkW/52M5sX2cxaXzA6rkdL3pJvk/A7zJKOxVxdP5BEjx LOVWThAzAhSrpx6wTtEuxOuGNvgP0u+FqBxzt/PPpfMHgVNR5ObRS+ZjgqAd zc9JtPpJejBsaiXKjbJF0oMRNFz+6bzDd9JDoJhK9Zt8T86TCYGrep9y9gI9 71BYr2j71VWzpIdu7cO5t3Nf6X6EwWxeVi8MIGbCwBq2l3pdpfskwoDLaqhT P9k3RgQzg9xQs4fuXzhO11X2cG+TfWXCcS3/eFUnj+6vGG3ZHeNlG59JFyOd 7f58YJrYSIIiQU2PUkO7BB38HEgXaY9ABKdKVL9GewRazOMsfddpj4Rswruo Zob2SPgcsdrBeUZ7FDozM7olctqjoFjwDnyZSns0npjnxjZl0h6NJA/bEwVt tMegb7+yejdbw/8LV+847g== "]]}, {}}, {{}, {RGBColor[0, 0, 1], LineBox[CompressedData[" 1:eJw10nsslXEcx/EjRCWT0RSSy0xCnHLp4JyP437cWrKwYloWJYcdMZfOnsyW jYn+oLlMWG6zEWZpSxe306Y0ly3LZJkhwsk1l6T5/p49z569nvc5v/2ePV/T W9KrsYd4PB72r//3g2NZ1F/SXNNgOSuiB4hY8bR8Hcysifw9wcqzLGYdONfU uWe8ZNZDsWZ89MwiswEyQyd6xp3myEZo75xSrOYxn0XYtq1LrZLZDGZDsQFV CT/JFojUmJy4oTpPtsT0uGOfsovZCpz/mCytaYFsjfq1C96u33+RbWDUGBMm TVki2+GpZGq3eWCZbA/FgNb6kuQ32QF50rD1uMRVMh9dO5Z/owo3yBfx6F2T 7LrVDvkSXM6b8MOFKjiwI8rzQ1MXWtXITtCtkw/aJh8mO6NF0ruRssq6C9qD Kupf9bL/XwaXxBdEjO7S+gLI1QLizbb/kF1hoyVfGlKw/bihtUlP1vJjnewO x+NJdrqGrAvRof5Rlv5wkyxC29KH8WNB22Rg4LF66dvZvQNzwHKa9EnqSVXa jwcMRkXy6Gn1A3MeyK3Q/6YwZO8jhs1waZxHDP2eE8OPJ9Oe6Kb1eJ7Ils6N RBhv0fqecCup0v58ju3XC4ap4t62jBXqXij5YlKtUa2k7o3UroS2yGv0vThv JJ7q2T8XqfsgN8s++XkkfW/OB+1F8tN3fNg8+ELljcC0O5HmhfOF7ljNC0E3 my8/9FdqZ+S4kDk/NJxIz+roY/Poj/l88/cWcWTOH7v6d60bDViXoHA3XOgw SvPOSZAt6lSLLWPzH4DNnMTYzHjWA6A0v2esI2Q9EAXlW+1HDVgPBMTz03nr M9SDUOVbcGbyK5kLwsSRteIHXawHQxgVujhSy3owpp3FpWVFrIdA0cZvSJOz HoLB25/kFvdZv4LK4VxJ3M0Z0T9LGC70 "]]}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0.5745583720849516}, Method->{}, PlotRange->{{0, 64.}, {0.5745583720849516, 1.7418424619565018`}}, PlotRangeClipping->True, PlotRangePadding->{{1.28, 1.28}, {0.023345681797431003`, 0.023345681797431003`}}]], "Output", CellChangeTimes->{3.600763493515049*^9, 3.600763646457944*^9, 3.6007639545961323`*^9}] }, Open ]], Cell["\<\ Se observa que hay dos picos claramente visibles, sin embargo aparecen adem\ \[AAcute]s otros que corresponden a espureos. Vamos ahora a aplicar el MME \ usando el algoritmo del gradiente conjugado. \ \>", "Text"], Cell[TextData[{ "La entrop\[IAcute]a de un vector ", Cell[BoxData["x"]], " se define como" }], "Text"], Cell[BoxData[ RowBox[{ RowBox[{"S", RowBox[{"(", StyleBox["x", FontWeight->"Bold"], ")"}]}], "=", RowBox[{"-", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"k", "=", "1"}], "n"], RowBox[{ SubscriptBox[ StyleBox["x", FontWeight->"Bold"], "k"], "log", RowBox[{"(", SubscriptBox[ StyleBox["x", FontWeight->"Bold"], "k"], ")"}]}]}]}]}]], "Text", CellChangeTimes->{3.6007640319204082`*^9}], Cell[TextData[{ "Evidentemente la entrop\[IAcute]a solo est\[AAcute] definida para valores \ de ", Cell[BoxData[ SubscriptBox["x", "k"]]], " positivos. Para evitar problemas num\[EAcute]ricos sustituiremos esta \ funci\[OAcute]n por otra definida en todo", Cell[BoxData[ FormBox[ RowBox[{" ", SuperscriptBox["\[DoubleStruckCapitalR]", "n"]}], TraditionalForm]]], " que tenga el m\[AAcute]ximo para valores de ", Cell[BoxData[ SubscriptBox["x", "k"]]], " positivos. Por ejemplo" }], "Text"], Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox[ StyleBox["s", FontWeight->"Bold"], "k"], RowBox[{"(", StyleBox["x", FontWeight->"Bold"], ")"}]}], "=", RowBox[{ StyleBox["{", FontSize->16], RowBox[{GridBox[{ { RowBox[{ RowBox[{"-", SubscriptBox[ StyleBox["x", FontWeight->"Bold"], "k"]}], "log", RowBox[{"(", StyleBox["x", FontWeight->"Bold"], ")"}]}], RowBox[{"\[ForAll]", " ", RowBox[{ SubscriptBox[ StyleBox["x", FontWeight->"Bold"], "k"], ">", "0"}]}]}, { SubscriptBox[ StyleBox[ RowBox[{"\[ScriptCapitalK]", StyleBox["x", FontWeight->"Bold"]}]], "k"], RowBox[{"\[ForAll]", " ", RowBox[{ SubscriptBox[ StyleBox["x", FontWeight->"Bold"], "k"], "\[LessEqual]", "0"}]}]} }], " ", "\[And]", RowBox[{"\[ScriptCapitalK]", ">", "0"}]}]}]}]], "Text", CellChangeTimes->{3.6007640388892183`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"S", RowBox[{"(", StyleBox["x", FontWeight->"Bold"], ")"}]}], "=", RowBox[{"-", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"k", "=", "1"}], "n"], RowBox[{ SubscriptBox[ StyleBox["s", FontWeight->"Bold"], "k"], RowBox[{"(", "x", ")"}]}]}]}]}]], "Text", CellChangeTimes->{3.6007640388892183`*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"\[ScriptCapitalK]", " ", "=", " ", "1000"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"ScalarEntropy", "[", "x_", "]"}], " ", ":=", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"x", "<", "0"}], ",", RowBox[{"\[ScriptCapitalK]", " ", "x"}], ",", " ", RowBox[{ RowBox[{"-", "x"}], " ", RowBox[{"Log", "[", "x", "]"}]}]}], "]"}]}], " ", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"ScalarEntropiaD1", "[", "x_", "]"}], " ", ":=", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"x", "<", "0"}], ",", "\[ScriptCapitalK]", " ", ",", " ", RowBox[{ RowBox[{"-", "1"}], "-", " ", RowBox[{"Log", "[", "x", "]"}]}]}], "]"}]}], " ", ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"nEntropy", "[", "x_List", "]"}], ":=", RowBox[{"Apply", "[", RowBox[{"Plus", ",", RowBox[{"Map", "[", RowBox[{"ScalarEntropy", ",", "x"}], "]"}]}], "]"}]}], ";", RowBox[{ RowBox[{"EntropyGradient", "[", "x_List", "]"}], ":=", RowBox[{"Map", "[", RowBox[{"ScalarEntropiaD1", ",", "x"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"nEntropy", "[", RowBox[{"{", RowBox[{"x", ",", RowBox[{"1", "-", "x"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", ".1"}], ",", "1.1"}], "}"}], ",", RowBox[{"PlotStyle", " ", "->", " ", RowBox[{"RGBColor", "[", RowBox[{"0", ",", " ", "0.8", ",", " ", "1"}], "]"}]}]}], "]"}]}], "Input", CellChangeTimes->{ 3.6007634526712646`*^9, {3.6007636118132067`*^9, 3.60076361887127*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[0, 0.8, 1], LineBox[CompressedData[" 1:eJw12Hk4VO/7B3BJQiotJCQiS6lEG+EdSaRF2gshSqV904Ks2RKJUhI+9r3s ZUtCSUmWsc2cGXNGCkOFLPF7vn/8/prrNXPmmvs85zn3/T6jaH/B0lFQQEBA f5qAwP9erWyD59Vk7i7T17UWsJrsKyvrSlpVOHd32cai2+lM+74yjR/L/yos NCxL2XPbq+8lVXay4oGYqr1C2SrP1te/ppeX/ZpXTMtUrQX3qmJm0Q0HBAna Jsczt0HhlUVsY+F1zJ0c6F9ybi9O75G890/KA4yxl/49Pw6By1SUXO/pi7iR y8p5TjZYs8Zy85bcQDx1GdNQeeWIi/lx2pOfHmJy+OjYrEEnPHke3hG+Khwn rr2pHlzjjO96VyTD7z6GlI7BS+bBS7jus/1h8fQoUOunzFMbbyBvl+iRjbvi UPS4vFzByROvPMLWuG1NhWWB9Lq2jV7g/44982hDGnqbLyaFzfQGh3ew7Ypq OuSlFIOFk3zA0wp60i2YCc9HHsd6aT+I6WlusP6ajZ2hRsMFJx5AsuY119cs Fxe2ccKT1kRAizUW6lfyGsnXtb6Hh0TgecFBoaTFb0Aleep6D0Ygd8VFoWPX 3sBSTIllm/MYlZpibgEri7Hhywk1mY2RMMyyn70urAQTh+nXQQZRuLk93cZF rRzbNvzhF6rEIs0q5aRUdQW+mY0sjF6fhO3qrt550dWY7602pXQmCaLjM7JS qqphWXqkJyU6CXJL88fd+qvRoFVckjczGQWx9VN1+jVokPU4WdeaDIu0dWFC 7cR9YoUTrqlgXpV/NFPiI+pD5Y8dq8qA05sr9o2On9DdG2cRtPIVQk5xhrsd 6rFJICM57nkh6vPOeTsYNWLXitGAhphSTFQaLZ640ILTleu29z2sQOO5RAk5 jzbMnpo9Z/L5e4xET81eJdGJYy8Y1m8ia3CnIHK2lCUL7CPrSx2aayFkor1Q tYLCh6U5XU/bP2Nc7HppRDobQ+7sIKbEV8htTc4LCeagfZl977b0Bggsib+Z t7cLl82vZP5e2YjbppLaioNdWO5R57LoUxP+pDh7evtyMV/DQ1G9oQm9LIul Pn5cCLSs+6jLaILnMtur3gFcdGhEydpwmzDdIWu7ezAXDxlnSuMnmjCR0fTY MYKLyTUiQlqrmsGkeqcPJXDRwjR6sDO4GZ9T56xVqOTCX68wwcOiBSrig4+z J7hIm3VfRfJQC1IPvwzynuSirs0uKdm6BTKBy8QPCdCYd3NWypczLXCz/Jg2 Mp1GZJ5N+hKfFiz2ELmiOotGyuoZOUVFLRAev77dcDGNDwqW5YPLGBB/+EPz +DoaP/kqhj7qDIielf4itIHG7LLxt9KaDFg81vmXvJHGXuuEd9Bn4OKi36J9 ujRaI/9W3T/EgPNPn54ThjR65r+oUw8in38v4YjvpiE6o7fNbogBXdGKtztP 0hjOb8htHmdAICG+oeQUDa5TUbC5YCv6rbV9Vp2mUfrJ12j93FbMGGv+IuJM 4+ojxVQR9VaEvX8zmXmJBlPp0M0sq1bIzchkFt2hkWtUsejfu1acTwt2/RhK I/ZP8uClj62ILNARkAijEZz4oJZX34odJlZq+x/ROC1m7V7f2QqjXS6nGRE0 5BtGeP+NtKLSfNGr+mc0/O1W5e9Y2YYM+c6rdxNo2Nx9sv/JwzboXlp0s7OA hv1VUdPoJ21IflA9zC2kcdLp1ub46Dakbck3+1lE47yF1bLs1Db4Gj+/8OcN DXcFhYHqijY4/cgo+VtGI6Y8KXDkVxu2S5+98qGaBlugoOLgvnY0Hx6uDGym wfujmm91pB1WcRfmnG+h8eP7kxT74+34LHd6zR4GjV/1t0LOn22HSor839lt NARjDWx8vdrxV9060bWTxjLDqtG8nHYclVX2Xs6l4eDRpCm5sAO7p83IfjxA Iyry1ESUTAcqOWK9DoM0ml6OVisrdsBi8m2d5i8aJpwlx9et7kC/o79n5W8a altP3t9n2oGAcquc9mEavUIjPaF3OjCamrnlwwQN5SX+eYu9OkC5ZbDu/qNh vV7WI9a/A3ZRjss3TtL47AjplxEd6PJFwospGtlV90zqszsgts7YxVaQh2t+ 0v/NoTsQtMrNJXYmD5kxqecjfnbgmvO8gK0iPPAK9XTlf3Xgrcf2rVziwz22 X1dNdSDzckO6ohgPujtSpnYu7oR5HsM6SJyHyVm6VoG7OmFM9RspzuPB98Ex SZH8TkR5nT5ntpiHpp4bYwrFnQg5tjarkFjZ+BFLp6ITFnclqlRkeHj391PK 2c+dcGtZd2WKWMDewOALrxN6vwdr4+R4uL1O8VTEIiYGaptdXy7l4WOw/s6s JUx8X5y6cI4CDzI9R9bWKDGRvm2Wx2niwucPx0fXMCHfFNomp8jDsLBQiLUp E1WHHEOvLePhYiuvUPkmE5a0/wK+Mg9l2tOj9d2Z4CzSENZZzsOc4KVeB32Y kE4v/uxBnGF0eJdfKBO0TCl7rgoPP9I+sH+mMBGh/XWuvCoPOsK8KqFsJtRW X39lQ+xnK5i+JJ+J6eFrV0YTq0puvr6nggmpOZPRi9V4cHBLF8tpY+J+/ukQ QXUechg1/FqKCcE8A64OsaA23cjlMaGXv3POReLY7iUvpH4zEX1eiWYQMy0e aN+axQIjVzk6agUPq9PSpMPmsRBwfGX8R2LXGTX/0haxoKLtd3eEWO61QE2n Egu3zPcW71rJwxGly1aGeiy8NrixpZf4qm3M081GLLhIVm2V0OAh+PlnxnpT FjRMJNW1iSsXrTywYj8L++1KI68RM/cfCVM+ykI2X255OPFo6L2v8rYsNC/1 eZhDvFqcu2uBMwuWX28t6CM2NZsfNPsyCxEHxVVFVvFwwnfLx5kuLHj6ZMkt I34i8NxkwosFC3GtTEviV/q13sP+LBjVLNxxmrju1mjFwAMW/g7NrXEjFhw6 APoZC7JhB84kEstpebuyYln4WBUXXEi88cKrN61JLNT+mR3+gdgynRr9lsHC gV1PXFuJnXvmbPqcw0J6r5Hpd+J7KvrXa4pYkGmfMzREHHfibG5FGQtj4v+8 BFfzUBwT+av4PQup9nOGZhO3dFZrFtSyEJW81VSa+JfM8PmXX1kwefzcVZFY /LByRloLC6u/yYWrE6uGW/5M6GSRufs2WJPYqOGuekwXC3Ps7p/ZQGw1N+vU 0x5y/k/c1TYT39jZmfCIz4J995MaA+KH/rO4wUMsiD5o3WFInF6ls8x/nNRb bJhlRNw1bV32WwEKl+99HfqfLyqYfymdTkFg2E/+f540sO9/I0yh+O8p9S3E gdY3ZxeJUshwPi+lTyx9J0QjX5zCd5koehNx4tMk85y5FO4XDkRoE2sXlZ7J nk/hE+/silXE5S1N/hmSFDhas2OWE+8e7k1OlaawQKpxRI64faFQTZIshbyh srULiJ20Zbvj5SmoK9XvEiEe3qslHKdIoaN9+s4Jst5eF82Wv1CmYGx2VINP LPHA1jhKlQL3bkM/RbzyU7BnxCoKCw2VZd4SF/5IiA3TpMAMGffMIjYRLSkP 0aZwQGegLorY1uTnZIAOBbXtGmJXifsdBOX99CgIN14ctya+7bVY3wcUHonU fTYhjijfftt9GwXpjveyksRKLJund0wpVAgeDx8j+zf737Wim+YUAiiRQSZx rW78yJW9FMSHfHfHEwvkCVx1sqYQFDgYq0ic7KjmnmhLof+Yrdw4ub/2SFkE cE9Q4MevN/lGHH3jRYztGQruZ4WT3Yn1dA3qDt+gsN+Y+V8NuV+5PxwZEbco 6EoujXxMHPjsflejK4XKZSn6jsRtEx2jFt5kvUz8Hk2QfuBSdlvFLJTCUPuI iyyxwsX/1t57RMHaQ+Q7i/SbaoVavfePKbitNxD+j1jKU2afYTSF1XFCFsuJ c42L3HXSKKTZnJGWIf1q4MMwQ/09hU128TcqSX98ckuee6qGQmbAheYrxFhp wk+opeDiJshVJA4OCp+xrIHCeIzDrFtKPGjsWaclwyL779WBABnSb50aLwWK jVLY3XDOSI70bwmfyAjTCQoi2BWQL89Dwfq3sb5TFCbVbt/dQyz8WKJQUJiN HJOYG65LeIg/ksUdm8+G9EaZkA+yPLCZvfo/Ndjotai0WybNw6Ps1qbBNWwY uad+yFpE9oNn1bm/WmzsnDkmpkecvjzm+Qwd4ricxXulyHo57/u3dBsbB/Nj T55bSOoZKyreb83GsqdJ+eclyDyqTdx/zJaNH+qiL3vm8nD9eViv3Qk2AjWN vU8Qd2w5J3fhNBvFGUvz9s3hIemewh3/a+T3fx4kmZEHA8l7m0vvs8F1bOgN I/N1kHflW2UIG383zb44IUzOp9D2bG0YGxsKb9eeIBaz0n3GiGRj/tZqsTUz yLyM6xv7lcAG5dCRUkDmtbPm/teqpWxcDmg8Hkrm/TNzRZ3QfjbUxOYdOEny g1OWfAljkI13ona9kSRfbFggZ6gwxMaTIkeVTyR/1LdJmWWNs6GpZOKpwach dEb8SJ0IB4mCZ1s6ftI45zfiIqrEgWC7bPc4yS94X1focYiDqdItC7VJPhJX r9X7cJQDjX8OnbpNNNqCasolbDi4k1NjYthI8uL+dzUvHDioPzEyatZAI7mr kFFyiYPOZ36sHZ9JPhaM/zsayIH45LswsSoaXQa3Nl0u56DtuMTE3VwaUo8E Tqe+42Dg8PioXQ4N0x7fSE4VB5qvBeMMX9HIDAsbs6zjYJ+l0MBEFo1b39OL tds4eC4n63EijYbEQ9aWod8crFO/YDcQS0OPNjZzUe1CuX2ZvUgwjfDAuUfd grvwtGBMd6U9jcZpxxa5h3ahynf3/C+2NOa7JDa6P+qCcvVo0aXjJL866O/x eEq+b/rzZo4VDV/9s8beiV3Qmr+tSO0QDRd+1eqA0i6IPfuypM2cxtF9btMj +rtw5OiZwkMkzyvI9qVn7OaCc8oxqWAaOV7J7CB/Lxf5m6JLfkxx8XVlwpTm AS66quoc5cjzhKeetWXuMS6+CZ84e3uMvG9dN/LaiYuJeU9N1H5zkRiTYVTj SZ5jLM2l9LlcaKicZ3DyubCpWdO29z0XOmv506WX0hiPKD+yxIuLFQE6QQls Gm7nNx6KGumCr175BpsYHvblap+QOdqFi9fO0nsvdONTjsVrTgQHWoM27500 v2Ngg7Z0Tg0bnr+vp7ye1oNA19Au3TrSFza1lK5534Opmln2ZkEsGLRcGf7i /QNzDSQ/Pt7dieMT67TS9v9E40m7jJskpx9xMHiRItULwcHzs3QrW8C2dJbo +t6LqAPGXxP3NuLb0+b/thf0waP306mSa/WYdOipMHvdh9tO3WZmJ+uxYs0E 27ykD5uOCUd+OliPu+8UFSze9UG0bfb+oo310Oh1jjr8pQ/SM9fL7//1hVyf 6eGnu/tg45vtV/LnMzZRa3wDF/Wjofm4vvvPT4he7u/05UY/iqVMX1k8rQHX I/3K/Nv92ObHDTe4UgN15he3A279qCu5IzF3Zw1yI6TC2737ker2oOXgv2p8 mBn/tju0H0ppOz/dsqnGn55SGYH0fvh6XndtlanCjsw/nzSpfmztCxb+dvMd QsSkGVe6+nH48r4Wn93v0HRyc1c+rx+Ghb8nJJTewXap56heXz+6v4+7ZtZW 4FrIXBWzsX6sT3trMCpTgZjLK9ztFvKhnJD69ujLcoxssF370JSPV4VF/VWX 3sA01/TDc3M+rlq9Gw+RfIMna9fapezmI2PmNht+4WvoagiGlu/nw5KRw7Cb LIKbYvxAvy0fIYXKbTyvQgiJf88yv8nH+JHw+Bs38nDAv377oTt8ON8XOeC1 IA+JwkUse3c+7lw2d5zIzMX2aQFzb/nwsd5B/qsaNwcBwysvJIeSepp/zRDY /goS7AurZ6TycS/j0qElTRmwszlcJZHBR/oW/68Fw+l41b7FRi6bj3ZOo6uX VDr2Nc8L1s7nQ638gvopi1RE1Ob02VXwwetfyliYlohu0yifc+/5kNbXMPzz JgGbqryX3KzhoyhZdLXJx3i0lh/YFfKZD20v/b+ZzDjI5Y+kl7bx8XH0gciv iuc4p00Zf+zkgxsWe7jw5TOUZtd0NFF81C6YyHj/OBLH0yLF+3h8vJnSVJE3 CUe2mmf83x4+ckuSH0iKh0Eg8YyeUB8fDVtVaffKEMTFbHaW/c2H+N+vqhrC fvi9RFlIdZgPWV/uuUvTvWD8TDxKa5TUa5/zr8T3Nh5JD2kbTJD6MIN+vPci 6PDOWrMpPnL26m7pvrELz0S7Xq6YNoBtHqHnq9BQttW0xTVbcABlP0xCN4nQ ZSX//jGapg/g//+/+z+8wgUc "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{{-0.1, 1.1}, {-0.937299729071874, 0.6931471635000277}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.600763494311929*^9, 3.600763646645522*^9, 3.6007639548461065`*^9}] }, Open ]], Cell["\<\ Imponemos como restricci\[OAcute]n que nuestra soluci\[OAcute]n cumpla\ \>", "Text"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["\[ScriptCapitalX]", "2"], RowBox[{"(", StyleBox["x", FontWeight->"Bold"], ")"}]}], "=", RowBox[{ RowBox[{"-", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"k", "=", "1"}], "n"], SuperscriptBox[ RowBox[{"(", FractionBox[ RowBox[{ SubscriptBox[ StyleBox["y", FontWeight->"Bold"], "k"], "-", SubscriptBox[ StyleBox["o", FontWeight->"Bold"], "k"]}], SubscriptBox[ StyleBox["\[Sigma]", FontWeight->"Bold"], "k"]], ")"}], "2"]}]}], "=", " ", "n"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", "donde", " ", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ StyleBox["y", FontWeight->"Bold"], " ", "=", " ", RowBox[{"R", "[", StyleBox["x", FontWeight->"Bold"], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"{", SubscriptBox[ StyleBox["o", FontWeight->"Bold"], "k"], "}"}], " ", "son", " ", "los", " ", "datos", " ", "observados"}]}], "\[IndentingNewLine]"}]], "NumberedEquation"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"ChiSqr", "[", "x_List", "]"}], " ", ":=", " ", RowBox[{"Apply", "[", RowBox[{"Plus", ",", SuperscriptBox[ RowBox[{"(", " ", RowBox[{ RowBox[{"R", "[", "x", "]"}], "-", "observed"}], ")"}], "2"]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"ChiSqrGradient", "[", "x_List", "]"}], " ", ":=", " ", RowBox[{"2", " ", RowBox[{"Rt", " ", "[", RowBox[{ RowBox[{"R", "[", "x", "]"}], "-", "observed"}], "]"}]}]}], ";"}]}], "Input"], Cell["\<\ Tomaremos como soluci\[OAcute]n inicial la distribuci\[OAcute]n uniforme del \ total observado, es decir, el vector cuyas componentes son todas la media de \ los valores observados\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"L0", " ", "=", " ", RowBox[{"Apply", "[", RowBox[{"Plus", ",", " ", "observed"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"f0", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"L0", "/", "n"}], ",", RowBox[{"{", "n", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x0old", "=", " ", RowBox[{"Chop", "[", RowBox[{"Fourier", "[", RowBox[{ RowBox[{"InverseFourier", "[", "observed", "]"}], " ", "crudeinverse"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"S0", "=", RowBox[{"nEntropy", "[", "f0", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"C0", "=", RowBox[{"ChiSqr", "[", "f0", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ SubscriptBox["C", "aim"], "=", RowBox[{"n", " ", SuperscriptBox["s", "2"]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"{", " ", RowBox[{"L0", ",", "S0", ",", "C0", ",", SubscriptBox["C", "aim"]}], " ", "}"}]}], "Input", CellChangeTimes->{3.6007636714102097`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ "48.585907542730936`", ",", "13.387825579634745`", ",", "7.785071268521833`", ",", "0.0064`"}], "}"}]], "Output", CellChangeTimes->{3.600763494436926*^9, 3.6007636861914997`*^9, 3.6007639548929844`*^9}] }, Open ]], Cell["\<\ Las direcciones conjugadas en la m\[EAcute]trica de la entrop\[IAcute]a se \ calculan por el m\[EAcute]todo explicado en el art\[IAcute]culo Maximum entropy image reconstruction: general algorithm J.Skilling and R.K. Bryan Departement of Applied Mathematics and Theoretical Phisics, Silver Street, Cambridge CB3 9EW Accepted 1984 May 14 \ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"\[Del]", "S"}], "=", RowBox[{"EntropyGradient", "[", "f0", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Del]", "C"}], "=", RowBox[{"ChiSqrGradient", "[", "f0", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ SubscriptBox[ RowBox[{"L", " "}], "S"], "=", " ", SqrtBox[ RowBox[{ RowBox[{"\[Del]", "S"}], ".", RowBox[{"(", RowBox[{ RowBox[{"\[Del]", "S"}], "/", "f0"}], ")"}]}]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ SubscriptBox[ RowBox[{"L", " "}], "C"], "=", " ", SqrtBox[ RowBox[{ RowBox[{"\[Del]", "C"}], ".", RowBox[{"(", RowBox[{ RowBox[{"\[Del]", "C"}], "/", "f0"}], ")"}]}]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"e", "=", " ", RowBox[{"{", " ", RowBox[{ RowBox[{"f0", " ", RowBox[{"\[Del]", "S"}]}], ",", RowBox[{"f0", " ", RowBox[{"\[Del]", "C"}]}], ",", RowBox[{"2", " ", "f0", " ", RowBox[{"Rt", "[", RowBox[{"R", "[", RowBox[{ FractionBox[ RowBox[{"f0", " ", RowBox[{"\[Del]", "S"}]}], SubscriptBox[ RowBox[{"L", " "}], "S"]], "-", FractionBox[ RowBox[{"f0", " ", RowBox[{"\[Del]", "C"}]}], SubscriptBox[ RowBox[{"L", " "}], "C"]]}], "]"}], "]"}]}]}], " ", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{"{", RowBox[{ RowBox[{"e", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"e", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"e", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], "]"}]}], "Input", CellChangeTimes->{{3.6007636954884*^9, 3.6007636961602783`*^9}, { 3.600763736159358*^9, 3.6007637373468337`*^9}, {3.600763768380354*^9, 3.600763779317831*^9}}], Cell[BoxData[ GraphicsBox[{{}, { {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData[" 1:eJxdyb1KwgEUxuGXpkYHh4aGFIkICTPTPvWffZhWJo2NQXMX0OJteCOtXkt3 4CVoFIHPgcPhOb/K28fL+0aSYrU/93cWvc+vxfx19j3/exTr3sQlXMZbeBvv 4Cqu4V28h/dxHR/gBj7ETXyEW/gYt3EHn+BTfIbP8QW+xF3cw8W6pzhXdJw+ HeeajnNDx7ml49zRcQZ0nHs6zpCOM6LjPNBxHuk4T3ScMR3nmY4z+fcSLeeF /A== "]]}, {RGBColor[0.6, 0.24, 0.4428931686004542], PointBox[CompressedData[" 1:eJw103ss1WEcx/GTtaZZy3XphDDkFuV+3M7HwXFcj5JbhUOYsYOYJRXHbFjU NBu1WI0yRI5GF7FDR4lQOiutVnPJ6SJC/7CaVfN9fttvv72e937P88/3MUvJ OZqmweFw8O/9/916Vvgd9YEu3ptP+LQARXbiUHEJsybc1aJG17UhsjbmfdPW xbHM+lgarFs53jJINsTIaNeXC7UKshFSdW2+C30HyKY4Z7fc07T+iGwOj4hq QUl2L9kCOyebS+O775KtYCa0/Vj5rZFsDUuFp5dWVpliy7bYiN2zWBjVSbaH QUtBpbxogOyAX0JJUWLBCPkQBN1dYpv2V+TDWJBInqVtviE7of3FQfOM0x/I ztAqq1pXxcyQXXBvd2VbYccc2RUVoquc5zWfyW5wajDLd1QvkN0hXNWrU9xU kz3Qf+SyPH2MmQfdUe6sbxmzJ0TJa+q4P+x/L+Sr33ONpMzeiA5wmzffweyD tccNwdZqdr4vesSj6SH6rPPRbVwQryphBiIidap2bafzZIBznMUEP4ed74eV zTtTedWs+yHh697Osw6sC+DZp3yby6P9ZALc6L9oHN06T90ft1xXnz68NEvd HyY5MpH2/k/UA2Bk9jqGN/COegB6+644tsWpqAfi5AHJnEf5OPVAxN8ur9hI V1IXIqG19sQ2gwfUhcida+FNLDVTD0LGtMp2uDl9a15kQZBKW4eSSptofkRo Ukp//NaTUxehTnV/0qeezV8wYs5ncde5fdSDcaZ05rp0nM1vCPx5SW6SJJp3 WQiGRaWWmjJ2H0KxaO2feu0Y3R9ZKErc9awif7L7FAZ1knxD45SSehjimyrd p/rJnHAYmngsz3CGqYcjJS9ZIbYncyJgHDs1NuvJegTmuC91VHasiyHnCYz2 abAuhs90TY7pENs/EsWZUeLQTCX/L6H8NlY= "]]}, {RGBColor[0.6, 0.5470136627990908, 0.24], PointBox[CompressedData[" 1:eJw103ss1WEcx3FjNSPGnIxInKnc73Jccj4ux+W4NbUVU5GphrmsolyOE1E0 ZS7xB3MZR+oQHRxhnTPUWJmkY20xY9HmlqnR5Y8u831+27PfXr/3nt8fz76P xaX06CR1NTU1/Fv/33vPFv+2+peyhJZBBX0Aig6YC5OZNWG3eMtzncushzxB zU3r8RdkDt4UZelcOMNshN8bZpsXlQNkU8ReL7Ps0WA2RwDf8MY6R07mImnO LTx9tY9siWGf4Qzk9pKP4VPM0XoH5XOyFar8DfuM+7rJNohuWlDU9ErJdlip OzR6LqWd7AADWeEo2pvJTjDmV7vKTWvJzviubKs/LiwmuyBXVTjfjFj+nl2x UJmzuXOvjOyGnY0ak6GJWrI7Fg/rDaSnNJFPYPRtdmKUUSvZA9NZzrK0JQmZ B4eZOeWG6jHZE8sWFREm+k/IXljmXD7r9JLZG6cb9TqmtJ6SfSB7MJU6b8Z8 EgJ/mV0uh9kX3PZ1x/Fttp+P1IzslVwpM+AhTBuRcMliQKQ0CNa27aDuh2vz OSGv69qp+0F+p9zgvqCNuj+uNvdYSYxbqPvjXV5URCW3nnoAGhIcLMV5ldQD sA3eZIejiHog8gIfnp8wzdw7b3Eg3OOsRCXFFXT+AoT8UUjt5fXUBdAf+zgW w2ulHoS1g1pLnl0d1IPgPrDqIm/soh6MflFZbIqoh3ow4qW2o+PxMuohUCX+ 0nn0k+ZLHILQufyZ6sh+6qEQyO0HeFdoPsWh8N1faF0UxuZXiLXUxBIvDZp3 sRA/pFUBWTVs/sMwOVI6m6xJ90UcBonu3c64BHZ/wtFtd6RB0sp6OJ6FK7fG VaxH4EN/Zk/+LusRKND9PDutPUQ9Eu+NHHf9OGRxJCx9yydV+qxHofmrp5fr Ptaj0On9ar1pk/3/FEwKSm2+TQ4q/gJcTjX9 "]]}}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0.}, Method->{}, PlotRange->{{0, 64.}, {-0.712207169376192, 0.15804693606576672`}}, PlotRangeClipping->True, PlotRangePadding->{{1.28, 1.28}, {0.017405082108839173`, 0.017405082108839173`}}]], "Output", CellChangeTimes->{ 3.6007634945306745`*^9, 3.6007637002071667`*^9, {3.6007637385187535`*^9, 3.6007637800678473`*^9}, 3.600763954924234*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ SubscriptBox["g", "\[Mu]\[Nu]"], "=", " ", RowBox[{"e", ".", RowBox[{"Transpose", "[", "e", "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ SubscriptBox["M", "\[Mu]\[Nu]"], "=", " ", RowBox[{"e", " ", ".", " ", RowBox[{"Transpose", "[", RowBox[{"Table", "[", " ", RowBox[{ RowBox[{"Rt", "[", RowBox[{"R", "[", RowBox[{"e", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", "e", "]"}]}], "}"}]}], "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\<\!\(\*SubscriptBox[\(g\), \(\[Mu]\[Nu]\)]\) = \>\"", ",", RowBox[{"MatrixForm", "[", SubscriptBox["g", "\[Mu]\[Nu]"], "]"}], ",", "\"\<\\n\!\(\*SubscriptBox[\(M\), \(\[Mu]\[Nu]\)]\) = \>\"", ",", RowBox[{"MatrixForm", "[", SubscriptBox["M", "\[Mu]\[Nu]"], "]"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(g\\), \\(\[Mu]\[Nu]\\)]\\) = \"\>", "\[InvisibleSpace]", TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"19.3578902168881`", "6.95552094516942`", "1.3554213669348978`"}, {"6.95552094516942`", "10.151945186063003`", RowBox[{"-", "0.8705383581096324`"}]}, {"1.3554213669348978`", RowBox[{"-", "0.8705383581096324`"}], "0.3374834680041412`"} }, 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$]]], "\[InvisibleSpace]", "\<\"\\n\\!\\(\\*SubscriptBox[\\(M\\), \ \\(\[Mu]\[Nu]\\)]\\) = \"\>", "\[InvisibleSpace]", TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"11.785799316453526`", "4.234778329801587`", "0.825230644504391`"}, {"4.2347783298015855`", "5.828627329501755`", RowBox[{"-", "0.4730942818249515`"}]}, {"0.825230644504391`", RowBox[{"-", "0.47309428182495156`"}], "0.1957983155906643`"} }, 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$]]]}], SequenceForm["\!\(\*SubscriptBox[\(g\), \(\[Mu]\[Nu]\)]\) = ", MatrixForm[{{19.3578902168881, 6.95552094516942, 1.3554213669348978`}, { 6.95552094516942, 10.151945186063003`, -0.8705383581096324}, { 1.3554213669348978`, -0.8705383581096324, 0.3374834680041412}}], "\n\!\(\*SubscriptBox[\(M\), \(\[Mu]\[Nu]\)]\) = ", MatrixForm[{{11.785799316453526`, 4.234778329801587, 0.825230644504391}, { 4.2347783298015855`, 5.828627329501755, -0.4730942818249515}, { 0.825230644504391, -0.47309428182495156`, 0.1957983155906643}}]], Editable->False]], "Print", CellChangeTimes->{3.6007634946244793`*^9, 3.600763955111799*^9}] }, Open ]], Cell[TextData[{ "Este es un problema de diagonalizaci\[OAcute]n simult\[AAcute]nea de ", Cell[BoxData[ FormBox[ SubscriptBox["g", "\[Lambda]\[Nu]"], TraditionalForm]]], " sim\[EAcute]trica definida positiva y ", Cell[BoxData[ FormBox[ SubscriptBox["M", "\[Lambda]\[Nu]"], TraditionalForm]]], " sim\[EAcute]trica cuya soluci\[OAcute]n puede verse en \n\nPage 463 Matrix \ Computations (Third Edition). Gene H.Golub & Chales F. Van Loan \nThe Hohns \ Hopkins University Press\nBaltimore and London" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"SimultanousEigenvalues", "[", RowBox[{"ASym_", ",", "BSymDefPos_"}], "]"}], " ", ":=", " ", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ "G", ",", "Gi", ",", "H", ",", "QR", ",", "Q", ",", "D", ",", "X"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"G", "=", RowBox[{"Transpose", "[", RowBox[{"CholeskyDecomposition", "[", "BSymDefPos", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Gi", "=", RowBox[{"Inverse", "[", "G", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Git", "=", RowBox[{"Transpose", "[", "Gi", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"H", "=", RowBox[{"Gi", ".", "ASym", ".", "Git"}]}], ";", "\[IndentingNewLine]", RowBox[{"QR", "=", RowBox[{"SchurDecomposition", "[", "H", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Q", "=", RowBox[{"QR", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"D", "=", RowBox[{ RowBox[{"Transpose", "[", "Q", "]"}], ".", "H", ".", "Q"}]}], ";", "\[IndentingNewLine]", RowBox[{"X", "=", RowBox[{"Git", ".", "Q"}]}], ";", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"X", ",", "D"}], "}"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"SE", " ", "=", " ", RowBox[{"SimultanousEigenvalues", "[", RowBox[{ SubscriptBox["M", "\[Mu]\[Nu]"], ",", SubscriptBox["g", "\[Mu]\[Nu]"]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"NewBase", " ", "=", " ", RowBox[{"SE", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ SubscriptBox["\[Gamma]", "\[Mu]\[Nu]"], "=", RowBox[{"Chop", "[", RowBox[{"SE", "[", RowBox[{"[", "2", "]"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ SubscriptBox["\[Gamma]", "\[Mu]"], "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ SubscriptBox["\[Gamma]", "\[Mu]\[Nu]"], "[", RowBox[{"[", RowBox[{"i", ",", "i"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", SubscriptBox["\[Gamma]", "\[Mu]\[Nu]"], "]"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"MatrixForm", "[", "NewBase", "]"}], ",", "\"\<\\n\!\(\*SubscriptBox[\(\[Gamma]\), \(\[Mu]\[Nu]\)]\) = \>\"", ",", RowBox[{"MatrixForm", "[", SubscriptBox["\[Gamma]", "\[Mu]\[Nu]"], "]"}], ",", "\"\<\\n\!\(\*SubscriptBox[\(\[Gamma]\), \(\[Mu]\)]\) = \>\"", ",", RowBox[{"MatrixForm", "[", SubscriptBox["\[Gamma]", "\[Mu]"], "]"}]}], "]"}], ";"}], "\n"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NewBase = \"\>", "\[InvisibleSpace]", TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"-", "0.22728511196689066`"}], RowBox[{"-", "3.1702585156634466`"}], RowBox[{"-", "0.4054048485743056`"}]}, {"7.756343696797889`*^-14", "4.2360897845825285`", "0.3510899318995805`"}, {"4.60369719493342`*^-13", "23.539030637561257`", "3.988257311712885`"} }, 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$]]], "\[InvisibleSpace]", "\<\"\\n\\!\\(\\*SubscriptBox[\\(\[Gamma]\\), \ \\(\[Mu]\[Nu]\\)]\\) = \"\>", "\[InvisibleSpace]", TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0.6088369747118115`", "0", "0"}, {"0", "0.27928909474833336`", "0"}, {"0", "0", "0.5709464214200825`"} }, 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$]]], "\[InvisibleSpace]", "\<\"\\n\\!\\(\\*SubscriptBox[\\(\[Gamma]\\), \ \\(\[Mu]\\)]\\) = \"\>", "\[InvisibleSpace]", TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ {"0.6088369747118115`"}, {"0.27928909474833336`"}, {"0.5709464214200825`"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], Column], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}], SequenceForm["NewBase = ", MatrixForm[{{-0.22728511196689066`, -3.1702585156634466`, \ -0.4054048485743056}, {7.756343696797889*^-14, 4.2360897845825285`, 0.3510899318995805}, {4.60369719493342*^-13, 23.539030637561257`, 3.988257311712885}}], "\n\!\(\*SubscriptBox[\(\[Gamma]\), \(\[Mu]\[Nu]\)]\) = ", MatrixForm[{{0.6088369747118115, 0, 0}, {0, 0.27928909474833336`, 0}, { 0, 0, 0.5709464214200825}}], "\n\!\(\*SubscriptBox[\(\[Gamma]\), \(\[Mu]\)]\) = ", MatrixForm[{0.6088369747118115, 0.27928909474833336`, 0.5709464214200825}]], Editable->False]], "Print", CellChangeTimes->{3.6007634950775566`*^9, 3.600763955564928*^9}] }, Open ]], Cell["Las direcciones en la nueva base son ortogonales", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ OverscriptBox["e", "_"], "=", RowBox[{ RowBox[{"Transpose", "[", "NewBase", "]"}], " ", ".", " ", "e"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"MatrixForm", "[", RowBox[{"Chop", "[", RowBox[{ OverscriptBox["e", "_"], ".", RowBox[{"Transpose", "[", OverscriptBox["e", "_"], "]"}]}], "]"}], "]"}], "]"}]}], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0.999999999999999`", "0", "0"}, {"0", "0.999999999999759`", "0"}, {"0", "0", "0.9999999999999968`"} }, 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$]]]], "Print", CellChangeTimes->{3.6007634951400623`*^9, 3.600763955627387*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ OverscriptBox["e", "_"], "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ OverscriptBox["e", "_"], "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ OverscriptBox["e", "_"], "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", RowBox[{"Joined", "\[Rule]", "True"}]}], "]"}]], "Input", CellChangeTimes->{{3.6007637905210476`*^9, 3.6007638033803997`*^9}}], Cell[BoxData[ GraphicsBox[{{}, { {RGBColor[0.24720000000000014`, 0.24, 0.6], LineBox[CompressedData[" 1:eJxd0zsvQ2EcgPGGxGUghMEgLo1I0yDu9+pDq4rSYmBGB5e4DPazEYuvYBGD xTcgPgKDyWAzstlQRDx9k5M3v/Pk5P2fnJzW9f3lfEkoFOLrKuw/6zVe9VFY N/HfG9TJFTTLNfTJ9eTkBo7kRs7lFh7kMKXf++2v22iS2xmWIyzJUTbkDnbl LrbkblblHkbkXmrlPu41bz9n8gDb8iD78hAX8nDR+46QkEc5lse4k8d5kWO8 yRM8yXGuZTj87wA61Sd5/j9vMMmp5p9iU32KA/UEl+oJ3tWTRHR+kozOnyav Ps2eesrfM0ixpj4D6jOE1dOUqad51HyzXGn+WU7U59hRn2NFfZ6Y+jxR9QyN 6pmi/22BavUFKtUX7WCRMvUs5erZop77e/4TJ88usQ== "]]}, {RGBColor[0.6, 0.24, 0.4428931686004542], LineBox[CompressedData[" 1:eJw102lIFGEcx/HJDo99kS5GZoWriHkhrkempvvTXY/1aDU7JI00dbMSDCtJ MJxAURFvEl9UaqKUR5rHmtnieuSueRAdZFK2UkYqqe1GiYi04X8GhuEz3xfz PPM8j+3FzJNpRgzDwHD/f25fayKUVjoHdnQP0gswGa0Wrc+6yCbQevqODb/v JJtDZfZoU+bfSrYE1mNF55obyFZgIss+Ldrnkw+BXUh0GZ6oFm1bAFwP7d3Y 10K2Q9KE6bRXcC/ZHuZz2nJxoJLsgPR26++5ccNkR6yXdpnw0l6SnWGlnqo5 JVSTXcGGFeqdlsfIbhj9IDCrWBohu8MqsbIl6u8AWQhJaq/zVX4z2QPndW25 R5iHNH5P5AobflpFKMle+PHmYI2rSE32RkW7Ga+mdZJ8FJJsQXxT4xTZB5oW fpa8apx8DAzPRjkyoyL7AoVy+ZnKRrIfULJ/5bdxJ43HH44XbOz0Qdz4j6NW 6mesKXpFDoBJgn158eVJciBUet9dhXc1ZBFe59/qutk9SAaYnH7et2t122YN rsufVyqe0veDIKhLL0g1Hd02G4Q197MJc3kT1IOhVfgl9/2i+bHBKLIV9skF XBdDEN8b++XzCHUx2owq3CwauP0lAWaWfO4M5NH3Jdg8PW2qkHHrHwKNx/ST gY9D1ENww3trd+YWt96hSH+u7pA60XqzoVCFLe2cnuV6GJKKWX68mvYLa3BK /vieS9x+CkfS6tuxlf4e6uFgVZa6hazH1KVgHXIC5lxrqUvB3NMtK2Lk1CPA 9GfLxntKaX4Gz6cINzIe0PwiwepW/pSEN1OPBDQ7RhLlbdSjgHdhHtoyOk9s FJiCoYRaa+68RQNXqpkkPZmNhnZRrFtO5v7fCQii+fcPtJNZg19IbjfpuS4D E7cqanfpoS5D/de95fOBZCYG9VXDszGHewb/ASXxMXc= "]]}, {RGBColor[0.6, 0.5470136627990908, 0.24], LineBox[CompressedData[" 1:eJw10w0s1HEcx/FjRuWpEP/DDFOZXZ66zkX4OO7OwyHSMmlak+Zhxdposjpq NivVViY5Emoe8pAuocPh5CHpaMyqVUrIFFq2VqaH+f5/23+/vf7v3/7//ff7 /xyOn44+ocvhcPDv+j9vjGX/edGIrGGup5tuQK8lwzC6mvUmfCxycVRIWW+F mW32Fp0hNdkCa+kFa0J71gzUyVH1doVdZFtUDN4ZlZ5Ske3ByxB6c9+2kx2h Tfulrn3RSnZC12VB/W7hI/JOqAXni0qcG8nOYJjZoOfaarILnGWGq1/9i8g8 cJw5Zjvybfw37IoPxSf52YoisjsWZ7evq0T3yB4onm5/dWlXE9kTblXfU6/w lOQ9WMwSPTZNaSPzEVb5Uz4cqyLvRYqRUW6qgZosQEedqdZ8vIfsBV29B/rW N/vIQvx2mnlTGdZP3oesUJXV2JdnZG+cuXbUdfLgINkHCzE23MTkIfJ+1MVn pv0xHib7Qh1VMF/oztoPU6XCEpdRdr0/+rvHbpmssM8DOG4JHksNAxuWA/MP k9bz7Nn3B2Ay102ZkqShHgDJlMmTuVu91EWw1tzV5g3Q98pFSJbZM6sLndQD 0RdVplzhd1APRNy0xTvz8FbqQbjNt7o+sa2FehCa5yy9DN7XURej1T3HrHOi groYzZUOOQmmBdQlcOcy/M0zuRv7LZdAWx1g5f2jjPZfCm3+UqHdpxrqUihq lUHRq83Ug5G4MM4Tf1NSD0Z6rE6PJqKNegjOIkZXzDylHoJlX9WNl3Gd1ENR obD83OjaTT0UfvcluqWH2P8/DJ7H4vUFy2R5GF4Xep4TpbLnR4aO8ouHxwfJ chlSumRysWEv9XBU1GiMiwVkeTiESkvL6Ui2R4CZZcyLj7A9AoOBTV2G8WyP RGsjV3Ahiu2RcBqpGpH5sP0ANJlXx8uZ3u6/8kQucg== "]]}}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0.}, Method->{}, PlotRange->{{0, 64.}, {-0.23393909690405845`, 0.22197074403148798`}}, PlotRangeClipping->True, PlotRangePadding->{{1.28, 1.28}, {0.00911819681871093, 0.00911819681871093}}]], "Output", CellChangeTimes->{3.6007634951713057`*^9, 3.600763805364794*^9, 3.600763955658615*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ SubscriptBox["S", "\[Mu]"], "=", RowBox[{"Chop", "[", RowBox[{ OverscriptBox["e", "_"], ".", RowBox[{"\[Del]", "S"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ SubscriptBox["C", "\[Mu]"], "=", RowBox[{"Chop", "[", RowBox[{ OverscriptBox["e", "_"], ".", RowBox[{"\[Del]", "C"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"\"\<\!\(\*SubscriptBox[\(S\), \(\[Mu]\)]\) = \>\"", ",", RowBox[{"MatrixForm", "[", SubscriptBox["S", "\[Mu]"], " ", "]"}], ",", "\"\<\\n\!\(\*SubscriptBox[\(C\), \(\[Mu]\)]\) = \>\"", ",", RowBox[{"MatrixForm", "[", " ", SubscriptBox["C", "\[Mu]"], "]"}]}], "]"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(S\\), \\(\[Mu]\\)]\\) = \"\>", "\[InvisibleSpace]", TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{"-", "5.795603497930296`"}]}, {"0"}, {"0"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], Column], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]], "\[InvisibleSpace]", "\<\"\\n\\!\\(\\*SubscriptBox[\\(C\\), \ \\(\[Mu]\\)]\\) = \"\>", "\[InvisibleSpace]", TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{"-", "2.0824294935087453`"}]}, {"0.6087331635783053`"}, { RowBox[{"-", "3.592794413110383`"}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], Column], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]}], SequenceForm["\!\(\*SubscriptBox[\(S\), \(\[Mu]\)]\) = ", MatrixForm[{-5.795603497930296, 0, 0}], "\n\!\(\*SubscriptBox[\(C\), \(\[Mu]\)]\) = ", MatrixForm[{-2.0824294935087453`, 0.6087331635783053, -3.592794413110383}]], Editable->False]], "Print", CellChangeTimes->{3.6007634952025537`*^9, 3.600763955721116*^9}] }, Open ]], Cell["\<\ El problema se reduce ahora a la optmizaci\[OAcute]n de una funci\[OAcute]n \ cuadr\[AAcute]tica con restricciones cuadr\[AAcute]ticas de dimensi\[OAcute]n \ 3 y con derivadas cruzadas nulas\ \>", "Text"], Cell[BoxData[{ RowBox[{"Maximizar", " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{ OverscriptBox["S", "~"], RowBox[{"(", "x", ")"}]}], " ", "=", " ", RowBox[{ SubscriptBox["S", "0"], "+", RowBox[{ SubscriptBox["S", "\[Mu]"], " ", SubscriptBox["x", "\[Mu]"]}], "-", RowBox[{ FractionBox["1", "2"], SubscriptBox["x", "\[Mu]"], SubscriptBox["x", RowBox[{"\[Mu]", "\[IndentingNewLine]"}]]}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Sujeto", " ", "a"}], " "}], "\[IndentingNewLine]", RowBox[{" ", RowBox[{ RowBox[{ OverscriptBox["C", "~"], RowBox[{"(", "x", ")"}]}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["C", "0"], "+", RowBox[{ SubscriptBox["C", "\[Mu]"], " ", SubscriptBox["x", "\[Mu]"]}], "+", RowBox[{ FractionBox["1", "2"], SubscriptBox["\[Gamma]", "\[Mu]"], SubscriptBox["x", "\[Mu]"], SubscriptBox["x", "\[Mu]"]}]}], " ", "=", SubscriptBox["C", "aim"]}], "\[IndentingNewLine]"}]}]}]}], "Text", CellChangeTimes->{3.600763997357827*^9}], Cell["Que podemos reparametrizar como", "Text"], Cell[BoxData[{ RowBox[{"Maximizar", " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{ OverscriptBox["Q", "~"], RowBox[{"(", "x", ")"}]}], "=", RowBox[{ RowBox[{"\[Alpha]", OverscriptBox["S", "~"], RowBox[{"(", "x", ")"}]}], " ", "-", " ", RowBox[{ RowBox[{ OverscriptBox["C", "~"], RowBox[{"(", "x", ")"}]}], "\[IndentingNewLine]"}]}]}]}], "Text", CellChangeTimes->{3.600764010560989*^9}], Cell[TextData[{ "derivando respect a ", Cell[BoxData[ FormBox[ SubscriptBox["x", "\[Mu]"], TraditionalForm]]] }], "Text"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Alpha]", RowBox[{"(", RowBox[{ SubscriptBox["S", "\[Mu]"], "-", SubscriptBox["x", "\[Mu]"]}], ")"}]}], "-", " ", RowBox[{"(", RowBox[{ SubscriptBox["C", "\[Mu]"], " ", "+", RowBox[{ SubscriptBox["\[Gamma]", "\[Mu]"], SubscriptBox["x", "\[Mu]"]}]}], ")"}]}], "=", "0"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Alpha]S", "\[Mu]"], "-", SubscriptBox["\[Alpha]x", "\[Mu]"], "-", " ", SubscriptBox["C", "\[Mu]"], " ", "-", RowBox[{ SubscriptBox["\[Gamma]", "\[Mu]"], SubscriptBox["x", "\[Mu]"]}]}], "=", "0"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Alpha]S", "\[Mu]"], "-", " ", SubscriptBox["C", "\[Mu]"]}], " ", "=", RowBox[{ SubscriptBox["x", "\[Mu]"], RowBox[{"(", RowBox[{"\[Alpha]", "+", SubscriptBox["\[Gamma]", "\[Mu]"]}], ")"}]}]}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ SubscriptBox["x", "\[Mu]"], "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"\[Alpha]", " ", SubscriptBox["S", "\[Mu]"]}], "-", SubscriptBox["C", "\[Mu]"]}], ")"}], "/", RowBox[{"(", RowBox[{ SubscriptBox["\[Gamma]", "\[Mu]"], "+", "\[Alpha]"}], ")"}]}], "\[IndentingNewLine]"}]}]}], "Text", CellChangeTimes->{3.6007640051547365`*^9}], Cell["\<\ donde hay que elegir \[Alpha] para que se cumpla\ \>", "Text"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ SubscriptBox["C", "0"], "+", RowBox[{ SubscriptBox["C", "\[Mu]"], " ", SubscriptBox["x", "\[Mu]"]}], "+", RowBox[{ FractionBox["1", "2"], SubscriptBox["\[Gamma]", "\[Mu]"], SubscriptBox["x", "\[Mu]"], SubscriptBox["x", "\[Mu]"]}]}], " ", "=", SubscriptBox["C", "aim"]}], "\[IndentingNewLine]"}]], "Text", CellChangeTimes->{3.60076401807664*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"CreateX", " ", "[", "\[Alpha]_", "]"}], ":=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"\[Alpha]", " ", SubscriptBox["S", "\[Mu]"]}], "-", SubscriptBox["C", "\[Mu]"]}], ")"}], "/", RowBox[{"(", RowBox[{ SubscriptBox["\[Gamma]", "\[Mu]"], "+", "\[Alpha]"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ OverscriptBox["C", "~"], "[", "\[Alpha]_", "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", "x", "}"}], ",", RowBox[{ RowBox[{"x", "=", RowBox[{"CreateX", "[", "\[Alpha]", "]"}]}], ";", RowBox[{"C0", "+", RowBox[{ SubscriptBox["C", "\[Mu]"], " ", ".", "x"}], "+", RowBox[{ FractionBox["1", "2"], RowBox[{ RowBox[{"(", RowBox[{ SubscriptBox["\[Gamma]", "\[Mu]"], "x"}], ")"}], ".", " ", "x"}]}]}]}]}], " ", "]"}]}], ";", RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{ OverscriptBox["C", "~"], "[", "\[Alpha]", "]"}], "-", SubscriptBox["C", "aim"]}], ",", RowBox[{"{", RowBox[{"\[Alpha]", ",", "0", ",", "2"}], "}"}]}], "]"}]}]}], "Input"], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV0nk4lesaBnBDyk7JbGOtb5kzfZKxop7XlrS1JUMDyRASdUwnlaT2wjKE qBYhSYZKiMwZ8hrSbpAMbSlTQpRChpXxO+v88Vz39fvzvq9H4bifrScfDw9P KPf+n5aeY+143HNX6lOZ8W08DLzGTNZ/QGY3tMpHvyK51ruTkl4qYw8Cb6fi NnA9TFNmnlN1htsGanbvKALva+0uuS3jAe3WdWbRXCvH0XMphjfEWy+Gq3C9 bOl+003VFxz0R7YVrhL4a0RXs1tyACScQFEMrg97cOquy5yBJ9axNy6uELjZ TLayKf0seE7zhTUvE1hHaWfxLCMYNl7IOstZInA6n2ueSnYIhL9XSpLkWnAo LOuQ6iXIteLhZSwSuHqwo7XJ/m8wX/bvE18gsP8BD0unZCa0n/OzmuUQuG9L lFm8TAQ8vtsQeH6OwOw7v2OVoyxgJbWdps0S2HLTQ5On6ZEQOjMYXvCTwBU/ XhtMMmKgz+BTwrVJAp92di6JcrsCf7pnPx/9TmDFN5Pa8tmxIHlDJEttgsBX C8XUbFSvQujSusbAMQKb03Nyxr0SIHjk6t6QUQJfakkTabBPhIQtNTP+wwTW N9SKd0i+BjEB0UbKgwTOlLSOvCLDhqBox5HCfwkc0nqz+VA1G2IC+fiVugh8 iDXIp3Q0CT4LPdkV2U7gDXMBl2vSk0FAji9G+jWBz3exz08wUuEV3e3hSUxg u7i+iiqcCk0Sdks2dQTW3q06F+GWBhby/oc0qgk8XFoZQM++BZyks7cflBLY 5sYHn/2qGRBwUCiiP5fAWvuU8mSfZwDzstBOsyzu3vynv4x63QFrTv/AzQwC uzxZfVxvnwmap6ysRW8SWFI9UftI8l1YUD3aYBFFYKZg6eYYmRw4XNVT3+NB YA2La9El/jmgl3jiwqgLgTtYfuMfn+eAg0GB/pAj93/4tfK1z+aCZ4mr8v0D BG5ZydHu6rgHrUztvT7G3L4zyQaMuDzgZU0UTG4kcErfBbPy1UcQNJRWWHuf jk1pDjn9B4ugwduO9SmTjr86GgkIFhZBk29vz69UOiYmgt0LyWJ4NproRMXS sZ+3oQpPbzH03fcL8vWjYzGPoge520rAWuNex7geHR9xuFs0NV0GzEjzmZlS Gp78aL8+S7YcSnMNapQKaDjSWdDTzqwc+hYNNS1zaLjU3U+2gs21QeLuEDYN C/vuZIUYVcCSMYfwP0PDz8J6HNeGVkKpsLFeiS4N6xaIrqUJVkNnT8tPnjw5 LEQxnSyIeljJ2KPx4rIsVrH5MdYU2wyVXgzVLo40PvdmO6TMvoC8rrzMPb8k sHSZaD8vbxuwGf3dvOvFcMWQ2es2kw4YjnzGdyx1I270uS3nMNEFR0zXVXZu W4fnVHduKvxPN7CLt+dN2/FivoV8ETmVDyAkJvxrWHq+3pByrbB71AvCYdML i7tH6l9H6L9p0RgAAzaz26LqBjASvy3k/RiEtH9lvfnzvsA3zkD/O50hUHBc HV8xmYecrzUvha9+hn/snEuXr/EgYfr9yeHhYZB7lpPYXS2A4rzLT2QqjwL2 yjXxdlmP5lH+8MnLXyDopcI5px2b0Bnd5jKppjHgfRoQpflaFA08ao5tV/gK CjVFG/LEJdCST9ct9QvfQGrzxfL+HVJI03Td6omaCbAlElOZGr8jiywFr8fS P6CB4zDppyOLztFdFD85TUJCLjsqnI+GtCYOCnVbTYFJRI1zbQ0diW+aEmxe nILUxsYPN2IYaE39rSKl9GnQvZD2Wy1bHp3Sr+n5bd9P0CmcmUS+Cij/z+OZ T+d/gq2opJCRiSKa+yYrMpo8A/rFjc4f5hRRuZ0FP3v3LFQ0fTL1q1JCl/ey 7pSNzsJEgbC84EllpNBn+T6bNQeL11+KNwqroB1GEWwZ3XkwbfknaLBOBTkR bQu8nfPg5eHybpOTKpp7GT/kz+RAhF6rKYNvM0pZ43zGVvkXOEkU7zuQshmJ 95dn8TX/gpdbQxt0NqshgVCm/faABdCSjxwLqVNDPHnORUNiiyD3g+11xVwd BaYVxYs1LsKUs3JhUpc6qtOMFhT1XYLa66NSTU4aSNQyUCtMeBmWpd+1qbzX QL3k287g2mVosXaNZzhqouIMUB52XoH+GEMQfKuJjhUZitJXV2DGeKUtCmmh v1l3mb0PVyFj7dZ500otVDhGvyX8FwX9esdG1iuRKPo2y7V/PwXauubNI8ok 8rD9rvLIhgL1xdxTWJVEtNraYqvDFDzM4ogHaZDoSsLRlrjjFAys+fzq41YS nTRMnV4fTIFYr8inDEQixQiJvevuUaBrKTm0dIxEq9subux+QIGsnNFQmwuJ er5/7riXT8F81VJ+thuJEg+XHtvzmAJ5qevplp4kojRs/8uqpcCMd/Zn0mkS 9bYnZPB3UsC8FHtRMYREVZEc9853FFSMKsnNXCQR29hFPfs9BW8VpsObLpFo X6522R/9FDyw2l/mHkai6vOtL5jjFKytFrl2N4ZEyaRBgs0EBRsErML8Y0kU MJRurzBJQaDP2l0QTyK1v04P4FkKDoRLSPYmkmgNb1dOIocC/60hOg+vk2iw 3NjHdZEC48kwqfNsEtX6ZG/RWaHAsFz/uXkyiW4yhOYoitsvOPUP8RQS/Q8k CMJ3 "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{{0, 2}, {-7.750217096515962, 14.797021943316377`}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.6007634952494626`*^9, 3.600763955814881*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ OverscriptBox["\[Alpha]", "_"], "=", RowBox[{"\[Alpha]", "//.", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{ OverscriptBox["C", "~"], "[", "\[Alpha]", "]"}], "==", SubscriptBox["C", "aim"]}], ",", RowBox[{"{", RowBox[{"\[Alpha]", ",", "0", ",", "1"}], "}"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ OverscriptBox["x", "_"], "=", RowBox[{"CreateX", "[", "0", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"f", " ", "=", " ", RowBox[{"f0", "+", " ", RowBox[{ OverscriptBox["x", "_"], ".", OverscriptBox["e", "_"]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"recon1", " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"f", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"PlotStyle", " ", "->", " ", RowBox[{"RGBColor", "[", RowBox[{"0", ",", " ", "1", ",", " ", "0"}], "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"recon1", ",", "recon", ",", "ourhope"}], "]"}]}], "Input", CellChangeTimes->{{3.6007634699368563`*^9, 3.6007634705150533`*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {RGBColor[0, 1, 0], LineBox[CompressedData[" 1:eJw10n0s1HEcB3DPJ81qyJibOGoebh3Hebjk3vfj4u6cowebiBs9GEV5yGZZ R1tJij9qZKLNQ0papIs2TUbUrhVWVIhiHpan0aJQyvf73X777fV77/P5fvbb xz4u5cBxPR0dHWw8/96bZ14Ubdsg/xwxLCIf4GU+zTldRW2M1G+zW9bbqLfj U7rlE+l1aguoKle7WIbUVsgXn/C0YYaI2bgYJHyYfWmQ2A7dM32yESdqDuaV LhFt2dSOMH2xkhjeRet341xey6liuxFiJwwryi4I7b8Su4Alqu5Ldxon5uKW IWd2IGeKeA8SinhNvNg5Yjck8Ljs51GLxO64nWJ32KvsJzEf2z6sHOrYtUrs gZrkmSlb2f/ftnE8YZl1/56Znx6xALVxDg7GTfrEXjAauzsuyDcg9gZztYEd XU/tg9cu1hWLztS+8M7WHF1wpfVC3Mg9mMhf1iXeC3NtDUc5RO/3Q0FrrbvU jnofdgbifefAHzKvP6o0kRbucpqLUJ7ENGkdaT/gQWFaXZKEzK8GHrEk6yYn 6f1ipJdZlmwVkfnUYqgSbTSvrtF5GfgIOzW6XJoz6I6uz8gwo/UBUJWEIfMI uU8dgIbId47FyetkvkCkNi8xg70rm1YH4nz/zV+I+UFyCXb4uzrxUxZILsFC u1+lIP07yfejYE0sZ0omSL4f66GDHZPcMZIHYaLQU+G5NEryILSXxuuNVND9 CYbYXqXlPvtC8mBIzvLTVC/p/kqhdqiefMsiVkvRWvN0LiOe7qMMUbESVt00 2Ve1DL6NgnDrx3R/5ei4ImbetNBcjuaeHiuM0jwEbbNuZ3r1ST91CGQfvfJi TWh/BUp/zwo9lmm9ArlFjc7DY7Q+FOW8grDVSZqHwijHbHnNlNYrwdZertQo aH8luoODYgzKaR6GO1nHMvuXh0R/AVggGRY= "]]}, {}}, {{}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJw10gtIU2EYxnEVCgPzRplNi2XDtFTUNe+2Z/M2N2+ZUaKJClLJVEQrllLH GgtNTTJvRKmZYakQImLeho0wMZMUisqZttRZqSuzTfOS0fcdOBx+58/hvPC9 B1KzYtNMjIyMsHX/e/6/dPwrDwVRnHINn7xA6WR2y/oktSnuLLd22yxQW6Lz o61K1ke9C02rVUvXL1Hbgmf3mHM3htoem2dKPgzJqNmIS1bEXN72hdgBN0tX k5VT1BzIDQnGuTPTxI4oXA9eM/yYIXaCLjFLfmi7lvgwyqxbnEUn54hdIFTd 6GrkfSN2Q5Ri5z3BxjyxO6Dd2xQyryP2QELdrXRW3BKxJ4ZMRxzP6peJuXBm id/55hmIjyJEILLwy1sl5sFYbGey0fiH2AvNDbrBzDFqb2RXP+IpLah9IH3h WfCqeoXYFxdH49M0I3piPxR36ZuTOL+J/VGSIlWz134RB6C2KzAs+D51IEqe trYzB+m8x9BeURvqmkW/52M5sX2cxaXzA6rkdL3pJvk/A7zJKOxVxdP5BEjx LOVWThAzAhSrpx6wTtEuxOuGNvgP0u+FqBxzt/PPpfMHgVNR5ObRS+ZjgqAd zc9JtPpJejBsaiXKjbJF0oMRNFz+6bzDd9JDoJhK9Zt8T86TCYGrep9y9gI9 71BYr2j71VWzpIdu7cO5t3Nf6X6EwWxeVi8MIGbCwBq2l3pdpfskwoDLaqhT P9k3RgQzg9xQs4fuXzhO11X2cG+TfWXCcS3/eFUnj+6vGG3ZHeNlG59JFyOd 7f58YJrYSIIiQU2PUkO7BB38HEgXaY9ABKdKVL9GewRazOMsfddpj4Rswruo Zob2SPgcsdrBeUZ7FDozM7olctqjoFjwDnyZSns0npjnxjZl0h6NJA/bEwVt tMegb7+yejdbw/8LV+847g== "]]}, {}}, {{}, {RGBColor[0, 0, 1], LineBox[CompressedData[" 1:eJw10nsslXEcx/EjRCWT0RSSy0xCnHLp4JyP437cWrKwYloWJYcdMZfOnsyW jYn+oLlMWG6zEWZpSxe306Y0ly3LZJkhwsk1l6T5/p49z569nvc5v/2ePV/T W9KrsYd4PB72r//3g2NZ1F/SXNNgOSuiB4hY8bR8Hcysifw9wcqzLGYdONfU uWe8ZNZDsWZ89MwiswEyQyd6xp3myEZo75xSrOYxn0XYtq1LrZLZDGZDsQFV CT/JFojUmJy4oTpPtsT0uGOfsovZCpz/mCytaYFsjfq1C96u33+RbWDUGBMm TVki2+GpZGq3eWCZbA/FgNb6kuQ32QF50rD1uMRVMh9dO5Z/owo3yBfx6F2T 7LrVDvkSXM6b8MOFKjiwI8rzQ1MXWtXITtCtkw/aJh8mO6NF0ruRssq6C9qD Kupf9bL/XwaXxBdEjO7S+gLI1QLizbb/kF1hoyVfGlKw/bihtUlP1vJjnewO x+NJdrqGrAvRof5Rlv5wkyxC29KH8WNB22Rg4LF66dvZvQNzwHKa9EnqSVXa jwcMRkXy6Gn1A3MeyK3Q/6YwZO8jhs1waZxHDP2eE8OPJ9Oe6Kb1eJ7Ils6N RBhv0fqecCup0v58ju3XC4ap4t62jBXqXij5YlKtUa2k7o3UroS2yGv0vThv JJ7q2T8XqfsgN8s++XkkfW/OB+1F8tN3fNg8+ELljcC0O5HmhfOF7ljNC0E3 my8/9FdqZ+S4kDk/NJxIz+roY/Poj/l88/cWcWTOH7v6d60bDViXoHA3XOgw SvPOSZAt6lSLLWPzH4DNnMTYzHjWA6A0v2esI2Q9EAXlW+1HDVgPBMTz03nr M9SDUOVbcGbyK5kLwsSRteIHXawHQxgVujhSy3owpp3FpWVFrIdA0cZvSJOz HoLB25/kFvdZv4LK4VxJ3M0Z0T9LGC70 "]]}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0.4399948817005819}, Method->{}, PlotRange->{{0, 64.}, {0.4399948817005819, 2.5846483461076852`}}, PlotRangeClipping->True, PlotRangePadding->{{1.28, 1.28}, {0.042893069288142065`, 0.042893069288142065`}}]], "Output", CellChangeTimes->{3.600763495405743*^9, 3.6007639558929873`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ OverscriptBox["\[Alpha]", "_"]], "Input"], Cell[BoxData["0.4858771062576055`"], "Output", CellChangeTimes->{3.6007634954681826`*^9, 3.600763955908626*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"nEntropy", "[", "f0", "]"}], ",", RowBox[{"nEntropy", "[", "f", "]"}]}], "}"}]], "Input", CellChangeTimes->{{3.600763882005661*^9, 3.6007638870369425`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"13.387825579634745`", ",", RowBox[{"-", "31.20734739765913`"}]}], "}"}]], "Output", CellChangeTimes->{3.6007634954994864`*^9, 3.6007639559242687`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"ChiSqr", "[", "f0", "]"}], ",", RowBox[{"ChiSqr", "[", "f", "]"}]}], "}"}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"7.785071268521833`", ",", "7.785071268521717`"}], "}"}]], "Output", CellChangeTimes->{3.6007634955150585`*^9, 3.600763955955532*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"Apply", "[", RowBox[{"Plus", ",", "f0"}], "]"}], ",", RowBox[{"Apply", "[", RowBox[{"Plus", ",", "f"}], "]"}]}], "}"}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"48.585907542730936`", ",", "75.94862800818134`"}], "}"}]], "Output", CellChangeTimes->{3.6007634955306797`*^9, 3.6007639559867654`*^9}] }, Open ]] }, Open ]] }, Open ]] }, WindowToolbars->"EditBar", WindowSize->{2030, 986}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, Magnification->1.5, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)", StyleDefinitions->FrontEnd`FileName[{"Article"}, "Preprint.nb", CharacterEncoding -> "WindowsANSI"] ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1485, 35, 155, 3, 203, "Title"], Cell[CellGroupData[{ Cell[1665, 42, 379, 7, 57, "Author"], Cell[2047, 51, 849, 19, 66, "Text"], Cell[2899, 72, 168, 3, 41, "Text"], Cell[CellGroupData[{ Cell[3092, 79, 1492, 44, 128, "Input"], Cell[4587, 125, 1395, 29, 378, "Output"] }, Open ]], Cell[5997, 157, 148, 3, 41, "Text"], Cell[6148, 162, 78, 2, 43, "Input"], Cell[6229, 166, 641, 19, 48, "Input"], Cell[6873, 187, 264, 8, 43, "Input"], Cell[7140, 197, 415, 12, 71, "Input"], Cell[CellGroupData[{ Cell[7580, 213, 540, 13, 71, "Input"], Cell[8123, 228, 1387, 29, 385, "Output"] }, Open ]], Cell[9525, 260, 195, 4, 41, "Text"], Cell[9723, 266, 75, 2, 43, "Input"], Cell[9801, 270, 269, 8, 43, "Input"], Cell[10073, 280, 288, 8, 43, "Input"], Cell[CellGroupData[{ Cell[10386, 292, 690, 16, 98, "Input"], Cell[11079, 310, 1387, 29, 385, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12503, 344, 84, 2, 43, "Input"], Cell[12590, 348, 2407, 47, 385, "Output"] }, Open ]], Cell[15012, 398, 475, 9, 67, "Text"], Cell[15490, 409, 353, 11, 43, "Input"], Cell[CellGroupData[{ Cell[15868, 424, 567, 15, 71, "Input"], Cell[16438, 441, 1479, 30, 364, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17954, 476, 90, 2, 43, "Input"], Cell[18047, 480, 2505, 48, 364, "Output"] }, Open ]], Cell[20567, 531, 223, 4, 41, "Text"], Cell[20793, 537, 105, 4, 42, "Text"], Cell[20901, 543, 471, 18, 76, "Text"], Cell[21375, 563, 510, 15, 42, "Text"], Cell[21888, 580, 1069, 40, 58, "Text"], Cell[22960, 622, 392, 15, 76, "Text"], Cell[CellGroupData[{ Cell[23377, 641, 1745, 52, 179, "Input"], Cell[25125, 695, 6928, 121, 356, "Output"] }, Open ]], Cell[32068, 819, 94, 2, 41, "Text"], Cell[32165, 823, 1233, 40, 222, "NumberedEquation"], Cell[33401, 865, 553, 18, 75, "Input"], Cell[33957, 885, 204, 4, 41, "Text"], Cell[CellGroupData[{ Cell[34186, 893, 1104, 33, 207, "Input"], Cell[35293, 928, 252, 6, 42, "Output"] }, Open ]], Cell[35560, 937, 362, 10, 213, "Text"], Cell[CellGroupData[{ Cell[35947, 951, 1941, 66, 223, "Input"], Cell[37891, 1019, 2936, 56, 342, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[40864, 1080, 1037, 28, 98, "Input"], Cell[41904, 1110, 2571, 55, 140, "Print"] }, Open ]], Cell[44490, 1168, 523, 12, 139, "Text"], Cell[CellGroupData[{ Cell[45038, 1184, 2914, 78, 449, "Input"], Cell[47955, 1264, 3275, 78, 222, "Print"] }, Open ]], Cell[51245, 1345, 64, 0, 41, "Text"], Cell[CellGroupData[{ Cell[51334, 1349, 419, 13, 73, "Input"], Cell[51756, 1364, 734, 19, 70, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[52527, 1388, 506, 15, 46, "Input"], Cell[53036, 1405, 3054, 57, 342, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[56127, 1467, 756, 21, 102, "Input"], Cell[56886, 1490, 2034, 54, 140, "Print"] }, Open ]], Cell[58935, 1547, 214, 4, 41, "Text"], Cell[59152, 1553, 1122, 36, 209, "Text"], Cell[60277, 1591, 47, 0, 41, "Text"], Cell[60327, 1593, 434, 14, 92, "Text"], Cell[60764, 1609, 128, 5, 45, "Text"], Cell[60895, 1616, 1519, 50, 217, "Text"], Cell[62417, 1668, 72, 2, 41, "Text"], Cell[62492, 1672, 444, 14, 88, "Text"], Cell[CellGroupData[{ Cell[62961, 1690, 1246, 41, 98, "Input"], Cell[64210, 1733, 3158, 59, 354, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[67405, 1797, 1215, 35, 157, "Input"], Cell[68623, 1834, 3506, 65, 364, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[72166, 1904, 56, 1, 43, "Input"], Cell[72225, 1907, 112, 1, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[72374, 1913, 213, 5, 43, "Input"], Cell[72590, 1920, 195, 4, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[72822, 1929, 141, 4, 43, "Input"], Cell[72966, 1935, 173, 3, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[73176, 1943, 195, 6, 43, "Input"], Cell[73374, 1951, 176, 3, 42, "Output"] }, Open ]] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *) (* NotebookSignature @upckYQhFy@Y0Bw#WPBX5wzV *)