(*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 0, 0] NotebookDataLength[ 83258, 1782] NotebookOptionsPosition[ 80191, 1721] NotebookOutlinePosition[ 80584, 1737] CellTagsIndexPosition[ 80541, 1734] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Undamped noisy oscillator (Problem 8.7)", "Section", CellChangeTimes->{{3.51722206395014*^9, 3.517222086707106*^9}, 3.517239902257428*^9, {3.762032770218287*^9, 3.762032771890544*^9}, { 3.793842408012704*^9, 3.793842412628067*^9}},ExpressionUUID->"0f17a414-3f5c-4602-9c7d-\ b49762641eab"], Cell["\<\ \[Bullet] revised, Oct. 25, 2022: \[Dash] use prediction estimator to conform to numbers given in text (and \ is simpler) \[Dash] compare explicit construction with LQGRegulator command for closed \ loop \[Dash] miscellaneous code cleanup\ \>", "Text", CellChangeTimes->{{3.716613652345023*^9, 3.716613661762876*^9}, { 3.716635650557143*^9, 3.716635673813285*^9}, {3.793904900379339*^9, 3.793904916165039*^9}, {3.843071915202323*^9, 3.843071961225065*^9}, { 3.844255613422617*^9, 3.844255654635604*^9}, {3.8457406313604803`*^9, 3.84574067346458*^9}, {3.875643242184225*^9, 3.87564326428683*^9}, { 3.875643298696796*^9, 3.875643322222205*^9}, {3.87575307507894*^9, 3.8757531548388357`*^9}}, FontSize->16, Background->GrayLevel[ 0.85],ExpressionUUID->"98bfe169-696b-4d7a-beae-415ad2a61460"], Cell["Define the (discretized) system, with noisy inputs", "Text", CellChangeTimes->{{3.517233046588011*^9, 3.517233070817132*^9}, 3.5172399022578583`*^9, {3.517337933977951*^9, 3.5173379479973783`*^9}, { 3.76203275466077*^9, 3.762032764866577*^9}, {3.762098697413294*^9, 3.76209872252246*^9}, {3.762133818830474*^9, 3.7621338200516977`*^9}, { 3.7621916003750763`*^9, 3.762191601916913*^9}, {3.762199057204894*^9, 3.7621990734109297`*^9}}, FontSize->18, Background->GrayLevel[ 0.85],ExpressionUUID->"428ad61b-4d61-4385-b9cc-a671d6e59502"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Ts", "=", "0.1"}], ";", " ", RowBox[{"Tend", " ", "=", " ", "15"}], ";", " ", RowBox[{"Nt", " ", "=", RowBox[{"Round", "[", RowBox[{"Tend", "/", "Ts"}], "]"}]}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"#", " ", "time", " ", "steps", " ", "to", " ", "simulate"}], " ", "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"G1ss", "=", RowBox[{"StateSpaceModel", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"x", "''"}], "[", "t", "]"}], "+", RowBox[{"x", "[", "t", "]"}]}], "\[Equal]", RowBox[{ RowBox[{"u", "[", "t", "]"}], "+", RowBox[{"w", "[", "t", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", "t", "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"x", "'"}], "[", "t", "]"}], ",", "0"}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"u", "[", "t", "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"w", "[", "t", "]"}], ",", "0"}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x", "[", "t", "]"}], "}"}], ",", "t", ",", RowBox[{"SystemsModelLabels", "->", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}], ",", RowBox[{"{", "\"\\"", "}"}], ",", RowBox[{"{", RowBox[{ "\"\\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(\[FilledSmallCircle]\)]\)\>\""}], "}"}]}], "}"}]}]}], "]"}]}], " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"undamped", " ", "oscillator"}], ";", " ", RowBox[{ "separate", " ", "inputs", " ", "for", " ", "feedback", " ", "and", " ", "process", " ", "noise"}]}], " ", "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"G1dss", " ", "=", " ", RowBox[{"ToDiscreteTimeModel", "[", RowBox[{"G1ss", ",", "Ts", ",", "z", ",", RowBox[{"Method", " ", "\[Rule]", " ", "\"\\""}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"zeroIn", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "Nt"}], "}"}]}], "]"}]}], " ", RowBox[{"(*", " ", RowBox[{"list", " ", "of", " ", "zero", " ", "inputs"}], " ", "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Nu]", "=", "0.3"}], ";", " ", RowBox[{"\[Xi]", "=", "0.3"}], " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"process", " ", "&"}], " ", "measurement", " ", "noise"}], " ", "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"SeedRandom", "[", "2", "]"}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"for", " ", "reproducible", " ", "plots"}], " ", "*)"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{"processNoise", "=", RowBox[{"RandomReal", "[", RowBox[{ RowBox[{"NormalDistribution", "[", RowBox[{"0", ",", "\[Nu]"}], "]"}], ",", "Nt"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"measurementNoise", "=", RowBox[{"RandomReal", "[", RowBox[{ RowBox[{"NormalDistribution", "[", RowBox[{"0", ",", "\[Xi]"}], "]"}], ",", "Nt"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{"G1ss", ",", "G1dss"}], "}"}]}], "Input", CellChangeTimes->{{3.76219198228271*^9, 3.76219202606069*^9}, { 3.762192059175441*^9, 3.762192073790366*^9}, {3.7621935228878736`*^9, 3.762193523283403*^9}, {3.762197171436796*^9, 3.762197172312241*^9}, { 3.762197692851165*^9, 3.762197693101293*^9}, {3.7621990873351593`*^9, 3.762199091178487*^9}, {3.76219919130022*^9, 3.7621991914236803`*^9}, { 3.7621994984190273`*^9, 3.7621995042319202`*^9}, {3.762199614393346*^9, 3.762199614485477*^9}, {3.762199837936109*^9, 3.762199885905244*^9}, { 3.762203351802896*^9, 3.7622033538982077`*^9}, {3.762204835642231*^9, 3.7622048498512363`*^9}, {3.7622133285475473`*^9, 3.762213331792344*^9}, { 3.762213476135552*^9, 3.762213478283379*^9}, {3.762213719023858*^9, 3.7622137786698437`*^9}, {3.762215413022902*^9, 3.762215416027739*^9}, { 3.762215653783928*^9, 3.762215677590703*^9}, {3.762215708261642*^9, 3.762215708658139*^9}, {3.762222835762754*^9, 3.7622228616625633`*^9}, 3.762222939407667*^9, {3.762222996006044*^9, 3.762223081503343*^9}, { 3.762223179648263*^9, 3.762223180091831*^9}, {3.7622234592628517`*^9, 3.7622234798986187`*^9}, {3.762223516437893*^9, 3.762223519537795*^9}, { 3.762224759870309*^9, 3.7622247805808783`*^9}, 3.76222764026654*^9, { 3.7622276744654217`*^9, 3.762227674661675*^9}, {3.762227977095234*^9, 3.762228000331336*^9}, 3.762228075684421*^9, {3.762228178600918*^9, 3.7622281798815947`*^9}, {3.762228224231744*^9, 3.7622282268359537`*^9}, 3.762228313043219*^9, {3.793842001397479*^9, 3.793842002298188*^9}, 3.875753176709236*^9}, CellLabel-> "In[3075]:=",ExpressionUUID->"87764a34-4e04-4711-b148-a77543ad8919"], Cell[BoxData[ RowBox[{"{", RowBox[{ TemplateBox[{"0", "1", "0", "0", RowBox[{"-", "1"}], "0", "1", "1", "1", "0", "0", "0", InterpretationBox["\[InvisibleSpace]", None], "\"u\"", "\"w\"", "\"y\"", "\"x\"", "\"\\!\\(\\*OverscriptBox[\\(x\\), \ \\(\[FilledSmallCircle]\\)]\\)\"", { StateSpaceModel, { True, True, { True, False}}, {$CellContext`stname1, $CellContext`stname2}, {{ "u", "w"}, {"y"}, { "x", "\!\(\*OverscriptBox[\(x\), \(\[FilledSmallCircle]\)]\)"}}, Identity, Automatic, 2, 1, 2, 1, {{1, 1}, {1, 2}, {2, 1}, {3, 1}, {3, 2}}, {True, True}, False, Automatic, {{"u", "w"}, {"y"}, { "x", "\!\(\*OverscriptBox[\(x\), \(\[FilledSmallCircle]\)]\)"}}, { ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", "t", "]"}], ",", "0"}], "}"}], ",", SubscriptBox["\[FormalX]", "1"]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"u", "[", "t", "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"w", "[", "t", "]"}], ",", "0"}], "}"}]}], "}"}], ",", "Automatic", ",", "t"}}, {}}, "StateSpaceModelFull", BaseStyle->{NumberMarks -> False}, DisplayFunction->(TagBox[ SubsuperscriptBox[ RowBox[{"(", GridBox[{{ "", "", "\[SpanFromLeft]", #14, #15}, {#17, #, #2, #3, #4}, {#18, \ #5, #6, #7, #8}, {#16, #9, #10, #11, #12}}, AllowedDimensions -> {5, 4}, GridBoxDividers -> {"Rows" -> {False, Directive[ GrayLevel[0], GrayLevel[0.6]], False, Directive[ GrayLevel[0], GrayLevel[0.6]], {False}, False}, "Columns" -> {False, Directive[ GrayLevel[0], GrayLevel[0.6]], False, Directive[ GrayLevel[0], GrayLevel[0.6]], {False}, False}}, DefaultBaseStyle -> "SystemsModelGrid11"], ")"}], #13, ButtonBox[ GraphicsBox[{ RectangleBox[{0, 0}, DefaultBaseStyle -> "SystemsModelSuperscript"], InsetBox[ StyleBox[ "\[ScriptCapitalS]", DefaultBaseStyle -> "SystemsModelSuperscript"], {0.5, 0.5}]}, ImageSize -> 20], ButtonFunction :> Control`Typesetting`SystemsModelTypesetToggle[ ButtonNotebook[], FE`CMObj], Evaluator -> Automatic, Appearance -> None], MultilineFunction -> None], EventHandlerTag[{"MouseEntered" :> {FE`CMObj = MathLink`CallFrontEnd[ FrontEnd`Value[ FEPrivate`Self[]]]}, PassEventsDown -> True, PassEventsUp -> True}]]& ), InterpretationFunction->(RowBox[{"StateSpaceModel", "[", RowBox[{"{", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{1, 2}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{5, 6}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{3, 4}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{7, 8}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{9, 10}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{11, 12}, ","], "}"}], "}"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", "t", "]"}], ",", "0"}], "}"}], ",", SubscriptBox["\[FormalX]", "1"]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"u", "[", "t", "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"w", "[", "t", "]"}], ",", "0"}], "}"}]}], "}"}], ",", "Automatic", ",", "t", ",", "SamplingPeriod", "->", RowBox[{#13}], ",", "SystemsModelLabels", "->", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{#14, ",", #15}], "}"}], ",", RowBox[{"{", #16, "}"}], ",", RowBox[{"{", RowBox[{#17, ",", #18}], "}"}]}], "}"}], "]"}]& ), Tooltip->StateSpaceModel], ",", TemplateBox[{ "0.9950041652780258`", "0.09983341664682817`", "0.004995834721974235`", "0.004995834721974235`", RowBox[{"-", "0.09983341664682817`"}], "0.9950041652780258`", "0.09983341664682815`", "0.09983341664682815`", "1", "0", "0", "0", "0.1`", "\"u\"", "\"w\"", "\"y\"", "\"x\"", "\"\\!\\(\\*OverscriptBox[\\(x\\), \\(\[FilledSmallCircle]\\)]\\)\"", { StateSpaceModel, { True, True, { True, False}}, {$CellContext`stname1, $CellContext`stname2}, {{ "u", "w"}, {"y"}, { "x", "\!\(\*OverscriptBox[\(x\), \(\[FilledSmallCircle]\)]\)"}}, Identity, Automatic, 2, 1, 2, 1, {{1, 1}, {1, 2}, {2, 1}, {3, 1}, {3, 2}}, {True, True}, False, Automatic, {{"u", "w"}, {"y"}, { "x", "\!\(\*OverscriptBox[\(x\), \(\[FilledSmallCircle]\)]\)"}}, { ",", "Automatic", ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"u", "[", "z", "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"w", "[", "z", "]"}], ",", "0"}], "}"}]}], "}"}], ",", "Automatic", ",", "z"}}, {}}, "StateSpaceModelFull", BaseStyle->{NumberMarks -> False}, DisplayFunction->(TagBox[ SubsuperscriptBox[ RowBox[{"(", GridBox[{{ "", "", "\[SpanFromLeft]", #14, #15}, {#17, #, #2, #3, #4}, {#18, \ #5, #6, #7, #8}, {#16, #9, #10, #11, #12}}, AllowedDimensions -> {5, 4}, GridBoxDividers -> {"Rows" -> {False, Directive[ GrayLevel[0], GrayLevel[0.6]], False, Directive[ GrayLevel[0], GrayLevel[0.6]], {False}, False}, "Columns" -> {False, Directive[ GrayLevel[0], GrayLevel[0.6]], False, Directive[ GrayLevel[0], GrayLevel[0.6]], {False}, False}}, DefaultBaseStyle -> "SystemsModelGrid11"], ")"}], #13, ButtonBox[ GraphicsBox[{ RectangleBox[{0, 0}, DefaultBaseStyle -> "SystemsModelSuperscript"], InsetBox[ StyleBox[ "\[ScriptCapitalS]", DefaultBaseStyle -> "SystemsModelSuperscript"], {0.5, 0.5}]}, ImageSize -> 20], ButtonFunction :> Control`Typesetting`SystemsModelTypesetToggle[ ButtonNotebook[], FE`CMObj], Evaluator -> Automatic, Appearance -> None], MultilineFunction -> None], EventHandlerTag[{"MouseEntered" :> {FE`CMObj = MathLink`CallFrontEnd[ FrontEnd`Value[ FEPrivate`Self[]]]}, PassEventsDown -> True, PassEventsUp -> True}]]& ), InterpretationFunction->(RowBox[{"StateSpaceModel", "[", RowBox[{"{", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{1, 2}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{5, 6}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{3, 4}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{7, 8}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{9, 10}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{11, 12}, ","], "}"}], "}"}], "}"}], ",", "Automatic", ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"u", "[", "z", "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"w", "[", "z", "]"}], ",", "0"}], "}"}]}], "}"}], ",", "Automatic", ",", "z", ",", "SamplingPeriod", "->", RowBox[{#13}], ",", "SystemsModelLabels", "->", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{#14, ",", #15}], "}"}], ",", RowBox[{"{", #16, "}"}], ",", RowBox[{"{", RowBox[{#17, ",", #18}], "}"}]}], "}"}], "]"}]& ), Tooltip->StateSpaceModel]}], "}"}]], "Output", CellChangeTimes->{3.823229881558577*^9, 3.875753178481382*^9, 3.875753422954742*^9}, CellLabel-> "Out[3084]=",ExpressionUUID->"075c041e-c778-4cdf-93de-822e9a9a39de"] }, Open ]], Cell["\<\ Open loop disturbance response, using Kalman filter (prediction observer)\ \>", "Text", CellChangeTimes->{{3.517233046588011*^9, 3.517233070817132*^9}, 3.5172399022578583`*^9, {3.517337933977951*^9, 3.5173379479973783`*^9}, { 3.76203275466077*^9, 3.762032764866577*^9}, {3.762098697413294*^9, 3.76209872252246*^9}, {3.762133818830474*^9, 3.7621338200516977`*^9}, { 3.7621916003750763`*^9, 3.762191601916913*^9}, {3.8757531851107388`*^9, 3.8757531860242033`*^9}}, FontSize->18, Background->GrayLevel[ 0.85],ExpressionUUID->"6ac269f8-a39a-4ca9-a823-bda5dd8fc0cc"], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"kalEst", "=", RowBox[{"KalmanEstimator", "[", RowBox[{ RowBox[{"{", RowBox[{"G1dss", ",", "1", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"(", "\[NoBreak]", GridBox[{ { SuperscriptBox["\[Nu]", "2"]} }], "\[NoBreak]", ")"}], ",", RowBox[{"(", "\[NoBreak]", GridBox[{ { SuperscriptBox["\[Xi]", "2"]} }], "\[NoBreak]", ")"}]}], "}"}], ",", RowBox[{"Method", "\[Rule]", "\"\\""}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ym", "=", RowBox[{ RowBox[{"Flatten", "[", RowBox[{"OutputResponse", "[", RowBox[{ RowBox[{"{", RowBox[{"G1dss", ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"zeroIn", ",", "processNoise"}], "}"}]}], "]"}], "]"}], "+", "measurementNoise"}]}], " ", RowBox[{"(*", " ", RowBox[{"noisy", " ", "position", " ", "measurement"}], " ", "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"xs", ",", "vs"}], "}"}], "=", RowBox[{"StateResponse", "[", RowBox[{ RowBox[{"{", RowBox[{"G1dss", ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"zeroIn", ",", "processNoise"}], "}"}]}], "]"}]}], " ", RowBox[{"(*", " ", RowBox[{"physical", " ", "states"}], " ", "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"xe", ",", "ve", ",", "ye"}], "}"}], "=", RowBox[{"OutputResponse", "[", RowBox[{"kalEst", ",", RowBox[{"{", RowBox[{"zeroIn", ",", "ym"}], "}"}]}], "]"}]}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"Kalman", " ", "estimate"}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"dr", "=", RowBox[{"DataRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "Tend"}], "}"}]}]}], ";", " ", RowBox[{"lRed", "=", RowBox[{"Lighter", "[", RowBox[{"Red", ",", "0.7"}], "]"}]}], ";", RowBox[{"plotmax", "=", "1.5"}], ";", RowBox[{"plot0a", "=", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"ym", ",", "xe", ",", "xs"}], "}"}], ",", "dr", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"InterpolationOrder", "\[Rule]", "0"}], ",", RowBox[{"PlotStyle", " ", "\[Rule]", RowBox[{"{", RowBox[{"lRed", ",", "Black", ",", "Blue"}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"-", "plotmax"}], ",", "plotmax"}], "}"}]}]}], "]"}]}], ";"}], " "}]}]], "Input", CellChangeTimes->{ 3.762192170481154*^9, 3.7621922320747843`*^9, {3.7621934894231863`*^9, 3.762193489785074*^9}, 3.762197405894991*^9, {3.76219893699751*^9, 3.762198954100754*^9}, 3.7621990430792847`*^9, 3.76219941050808*^9, { 3.762199644530333*^9, 3.762199645011698*^9}, {3.762222342387641*^9, 3.762222342799039*^9}, {3.762222405129356*^9, 3.762222406789*^9}, { 3.762225083602105*^9, 3.762225084135707*^9}, {3.762228061101931*^9, 3.762228063441267*^9}, {3.762228350259441*^9, 3.7622283601439047`*^9}, { 3.793842036599881*^9, 3.7938420724449043`*^9}, 3.8757532398780107`*^9}, CellLabel-> "In[3085]:=",ExpressionUUID->"eb36d036-fcff-46ab-86b1-8f0c73a6b48b"], Cell["\<\ LQG controller for disturbance response (explicit construction)\ \>", "Text", CellChangeTimes->{{3.513266788083189*^9, 3.5132667973147306`*^9}, { 3.513272093815999*^9, 3.513272098697585*^9}, {3.5172389588536386`*^9, 3.5172389629386263`*^9}, {3.517239692171343*^9, 3.5172397065137873`*^9}, 3.5172399022732773`*^9, {3.51723998911631*^9, 3.517239998600172*^9}, { 3.762184806607835*^9, 3.762184813214233*^9}, {3.762197448350705*^9, 3.762197451565967*^9}, {3.7621995236785192`*^9, 3.762199526597042*^9}, { 3.76221942397791*^9, 3.762219430125779*^9}, {3.7622248943547163`*^9, 3.7622249000126257`*^9}, 3.7935869580542803`*^9}, FontSize->16, Background->GrayLevel[ 0.85],ExpressionUUID->"e7b0b0ab-e691-482d-8404-e12dd17caddb"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"q", "=", RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "0"}, {"0", "1"} }], "\[NoBreak]", ")"}]}], ";", " ", RowBox[{"r", "=", RowBox[{"{", RowBox[{"{", "0.1", "}"}], "}"}]}], " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"weight", " ", "matrices", " ", "for", " ", "state"}], ",", " ", "input"}], " ", "*)"}], ";", RowBox[{"k", "=", RowBox[{"LQRegulatorGains", "[", RowBox[{ RowBox[{"{", RowBox[{"G1dss", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"q", ",", "r"}], "}"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.762212854907606*^9, 3.762212908304455*^9}, { 3.762213131502129*^9, 3.762213143440868*^9}, {3.762213174489119*^9, 3.7622132795375834`*^9}, 3.762215559075066*^9, 3.762225229633073*^9, { 3.762225309324801*^9, 3.762225309511846*^9}}, CellLabel-> "In[3090]:=",ExpressionUUID->"c860eed2-6467-43e3-9dea-e94056cf430b"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"1.7413559943513492`", ",", "3.2605789511795216`"}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.76221877461033*^9, 3.762218788512574*^9}, { 3.762219490100769*^9, 3.76221951192369*^9}, 3.762219743489787*^9, 3.762219846473543*^9, 3.762222166821191*^9, {3.762222205692368*^9, 3.7622222275481453`*^9}, {3.762222321917609*^9, 3.7622223492321367`*^9}, 3.762222412343783*^9, 3.762222568514304*^9, {3.762222708122336*^9, 3.762222728856097*^9}, {3.76222278850074*^9, 3.76222286713358*^9}, { 3.762222923839801*^9, 3.762222945116789*^9}, {3.762223001656773*^9, 3.762223086942349*^9}, {3.762223185914515*^9, 3.762223246140349*^9}, 3.7622233022502203`*^9, 3.762223364621974*^9, {3.762223464959053*^9, 3.7622234855712023`*^9}, 3.762223525597213*^9, 3.7622245642933598`*^9, { 3.762224639759614*^9, 3.762224656249084*^9}, {3.762224769371434*^9, 3.7622247865111837`*^9}, 3.762225185751719*^9, {3.76222523125709*^9, 3.7622252368652477`*^9}, 3.7622253140519857`*^9, {3.7622276163722878`*^9, 3.762227646055316*^9}, 3.762227679242899*^9, 3.762227905518566*^9, { 3.762227984516673*^9, 3.762228005151773*^9}, {3.762228068712532*^9, 3.7622280801783123`*^9}, 3.76222818435646*^9, 3.762228243432735*^9, 3.7622283176428213`*^9, 3.762228377232953*^9, 3.76222884773031*^9, 3.793586874910433*^9, 3.793841896441588*^9, 3.793842004574978*^9, 3.79384209318346*^9, 3.793842219752531*^9, 3.793842273499786*^9, 3.7938424163997507`*^9, 3.793908756174274*^9, 3.823229881815641*^9, 3.875753295390532*^9, 3.875753423339957*^9}, CellLabel-> "Out[3090]=",ExpressionUUID->"3b3b61db-adae-45a3-8946-eb775a5d318e"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"l", "=", RowBox[{"LQEstimatorGains", "[", RowBox[{ RowBox[{"{", RowBox[{"G1dss", ",", "All", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"(", "\[NoBreak]", GridBox[{ { SuperscriptBox["\[Nu]", "2"]} }], "\[NoBreak]", ")"}], ",", RowBox[{"(", "\[NoBreak]", GridBox[{ { SuperscriptBox["\[Xi]", "2"]} }], "\[NoBreak]", ")"}]}], "}"}], ",", RowBox[{"Method", "\[Rule]", "\"\\""}]}], "]"}]}]], "Input", CellChangeTimes->{{3.762213377243408*^9, 3.762213461462036*^9}, 3.7622155647857018`*^9, {3.762219622346703*^9, 3.762219623832375*^9}, 3.875753291438706*^9}, CellLabel-> "In[3091]:=",ExpressionUUID->"2861f37c-7d70-4363-b519-8f535afe8bb7"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", "0.09046726479597311`", "}"}], ",", RowBox[{"{", "0.03063492588555703`", "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.7622134053123837`*^9, 3.76221348454224*^9}, { 3.7622137281727123`*^9, 3.762213783363058*^9}, {3.762213820535761*^9, 3.762213833275033*^9}, 3.762213877070547*^9, {3.7622153975480213`*^9, 3.762215519343144*^9}, {3.7622155830998573`*^9, 3.7622156110515547`*^9}, { 3.762215657844756*^9, 3.7622156820404387`*^9}, 3.762215713086279*^9, { 3.762215772203617*^9, 3.762215831941369*^9}, 3.762218788653944*^9, { 3.7622194902514563`*^9, 3.762219512063903*^9}, 3.762219625010806*^9, 3.762219743635029*^9, 3.762219846619762*^9, 3.7622221669634447`*^9, { 3.762222205843301*^9, 3.762222227690982*^9}, {3.762222322063205*^9, 3.7622223493794928`*^9}, 3.7622224124819593`*^9, 3.7622225686608562`*^9, { 3.762222708272482*^9, 3.762222728996635*^9}, {3.762222788663588*^9, 3.762222867274719*^9}, {3.7622229239882803`*^9, 3.7622229452642*^9}, { 3.762223001801507*^9, 3.762223087083143*^9}, {3.7622231860593023`*^9, 3.76222324628251*^9}, 3.76222330239325*^9, 3.762223364762514*^9, { 3.762223465114674*^9, 3.762223485715006*^9}, 3.762223525741941*^9, 3.762224564436062*^9, {3.762224639902917*^9, 3.7622246563926067`*^9}, { 3.762224769512933*^9, 3.762224786653796*^9}, 3.762225185901537*^9, 3.762225237017243*^9, 3.7622253142015123`*^9, {3.7622276165137672`*^9, 3.762227646205687*^9}, 3.7622276793961143`*^9, 3.762227905660632*^9, { 3.762227984662586*^9, 3.762228005345522*^9}, {3.762228068856852*^9, 3.7622280803292217`*^9}, 3.7622281845098963`*^9, 3.762228243579829*^9, 3.762228317800507*^9, 3.762228377384074*^9, 3.762228847879817*^9, 3.793586874987342*^9, 3.79384189651998*^9, 3.793842004644559*^9, 3.793842093269507*^9, 3.793842219826549*^9, 3.793842273569686*^9, 3.793842416504443*^9, 3.7939087562491837`*^9, 3.8232298818911333`*^9, 3.875753298145361*^9, 3.875753423424055*^9}, CellLabel-> "Out[3091]=",ExpressionUUID->"7d2f1ec4-1b04-4d83-9f34-f3d58b53dc70"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c", ",", "d"}], "}"}], "=", RowBox[{"Normal", "[", RowBox[{"SystemsModelExtract", "[", RowBox[{"G1dss", ",", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Eigenvalues", "[", "a", "]"}], ",", RowBox[{"Eigenvalues", "[", RowBox[{"a", "-", RowBox[{"KroneckerProduct", "[", RowBox[{"b", ",", "k"}], "]"}]}], "]"}]}], "}"}]}], "Input", CellChangeTimes->{{3.762212919898875*^9, 3.762212948479576*^9}, { 3.762213036957861*^9, 3.762213052721231*^9}, {3.762215750263303*^9, 3.762215760737423*^9}, {3.762218283862629*^9, 3.762218292804352*^9}, { 3.762218342919735*^9, 3.7622184000908337`*^9}, {3.762218680946004*^9, 3.762218698268532*^9}, 3.762219444586433*^9, 3.762219518608609*^9}, CellLabel-> "In[3092]:=",ExpressionUUID->"fa5d7659-7981-4357-ae51-8483a7c99902"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"0.9950041652780258`", "\[VeryThinSpace]", "+", RowBox[{"0.09983341664682818`", " ", "\[ImaginaryI]"}]}], ",", RowBox[{"0.9950041652780258`", "\[VeryThinSpace]", "-", RowBox[{"0.09983341664682818`", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"0.8751062293223717`", ",", "0.7806878375507984`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.76221293488391*^9, 3.762212948978569*^9}, 3.762213053492734*^9, 3.76221317837819*^9, {3.762213305110923*^9, 3.762213336398535*^9}, 3.762213484671115*^9, {3.762213728303975*^9, 3.762213783494706*^9}, {3.762213820644436*^9, 3.7622138335036087`*^9}, 3.762213877186715*^9, {3.762215397660692*^9, 3.7622155194542427`*^9}, { 3.762215583204165*^9, 3.762215611129513*^9}, {3.762215657978653*^9, 3.762215682182475*^9}, 3.762215713204556*^9, {3.76221576134483*^9, 3.7622158320323753`*^9}, {3.762218329274063*^9, 3.7622184005202723`*^9}, { 3.762218684757825*^9, 3.762218699174477*^9}, 3.762218788770749*^9, 3.762219448946104*^9, {3.762219490375031*^9, 3.762219519469844*^9}, 3.762219743765257*^9, 3.7622198467458487`*^9, 3.7622221670669527`*^9, { 3.76222220596078*^9, 3.762222227816289*^9}, {3.762222322179907*^9, 3.7622223494925213`*^9}, 3.7622224126049557`*^9, 3.7622225687982483`*^9, { 3.762222708393509*^9, 3.762222729080247*^9}, {3.762222788791313*^9, 3.762222867398018*^9}, {3.762222924072154*^9, 3.762222945358169*^9}, { 3.762223001882638*^9, 3.7622230871925793`*^9}, {3.762223186171926*^9, 3.762223246395832*^9}, 3.762223302540559*^9, 3.7622233649162493`*^9, { 3.762223465208108*^9, 3.76222348586718*^9}, 3.762223525852454*^9, 3.762224564553309*^9, {3.7622246400306597`*^9, 3.762224656514076*^9}, { 3.7622247696299553`*^9, 3.762224786778933*^9}, 3.762225186738328*^9, 3.762225237103949*^9, 3.762225314284264*^9, {3.76222761663647*^9, 3.7622276463039007`*^9}, 3.762227679493078*^9, 3.762227905775669*^9, { 3.762227984820716*^9, 3.762228005505319*^9}, {3.762228068981227*^9, 3.762228080481154*^9}, 3.762228184674303*^9, 3.762228243699445*^9, 3.762228317882803*^9, 3.762228377477157*^9, 3.762228847960552*^9, 3.793586875012887*^9, 3.793841896543832*^9, 3.793842004660512*^9, 3.793842093275581*^9, 3.79384221985321*^9, 3.793842273599155*^9, 3.793842416527213*^9, 3.793908756277561*^9, 3.823229881907692*^9, 3.875753303541192*^9, 3.875753423438319*^9}, CellLabel-> "Out[3093]=",ExpressionUUID->"d9035ca0-8d6f-4ef4-8d9c-4f864f178299"] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"a1", "=", " ", RowBox[{"ArrayFlatten", "[", RowBox[{"(", "\[NoBreak]", GridBox[{ {"a", RowBox[{"-", RowBox[{"b", ".", "k"}]}]}, { RowBox[{"l", ".", "c"}], RowBox[{" ", RowBox[{"a", "-", RowBox[{"b", ".", "k"}], "-", RowBox[{"l", ".", "c"}]}]}]} }], "\[NoBreak]", ")"}], "]"}]}], " ", ";", RowBox[{"b1", "=", RowBox[{"ArrayFlatten", "[", RowBox[{"(", "\[NoBreak]", GridBox[{ {"b", "0"}, {"0", RowBox[{" ", "l"}]} }], "\[NoBreak]", ")"}], "]"}]}], " ", RowBox[{"(*", " ", RowBox[{"inputs", " ", "=", " ", RowBox[{"(", RowBox[{"\[Nu]", ",", "\[Xi]"}], ")"}]}], " ", "*)"}], ";", RowBox[{"c1", "=", RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "0", "0", "0"} }], "\[NoBreak]", ")"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"closedLoop2", "=", RowBox[{"StateSpaceModel", "[", RowBox[{ RowBox[{"{", RowBox[{"a1", ",", "b1", ",", "c1"}], "}"}], ",", RowBox[{"SamplingPeriod", "\[Rule]", "Ts"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"xs2", ",", "vs2", ",", "xe2", ",", "ve2"}], "}"}], "=", RowBox[{"StateResponse", "[", RowBox[{ RowBox[{"{", RowBox[{"closedLoop2", ",", RowBox[{"{", RowBox[{"0", ",", "1", ",", "0", ",", "0"}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"processNoise", ",", "measurementNoise"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ym2", "=", RowBox[{"xs2", "+", "measurementNoise"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"pclosed2", "=", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"ym2", ",", "xe2", ",", "xs2"}], "}"}], ",", "dr", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"InterpolationOrder", "\[Rule]", "0"}], ",", RowBox[{"PlotStyle", " ", "\[Rule]", RowBox[{"{", RowBox[{"lRed", ",", "Black", ",", "Blue"}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"-", "plotmax"}], ",", "plotmax"}], "}"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"pclosed3", "=", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xs", "-", "xe"}], ",", RowBox[{"xs2", "-", "xe2"}], ",", RowBox[{"xs", "-", "xe", "-", "xs2", "+", "xe2"}]}], " ", "}"}], ",", "dr", ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"InterpolationOrder", "\[Rule]", "0"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"-", "0.3"}], ",", "0.6"}], "}"}]}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.7622196017229548`*^9, 3.7622196300437183`*^9}, { 3.7622196708272743`*^9, 3.76221971977817*^9}, {3.762219757847953*^9, 3.762219760501882*^9}, {3.762222519365864*^9, 3.7622225464090157`*^9}, { 3.762224982047587*^9, 3.76222500916313*^9}, {3.762225050362632*^9, 3.7622250600240107`*^9}, {3.762227776567103*^9, 3.7622277910734367`*^9}, { 3.762228370132823*^9, 3.7622283726252127`*^9}, {3.762228825483676*^9, 3.762228834914062*^9}, {3.7938420858435097`*^9, 3.793842088091926*^9}, { 3.7938421762513638`*^9, 3.793842213249407*^9}, 3.875753332906352*^9}, CellLabel-> "In[3094]:=",ExpressionUUID->"69771da2-719e-4de9-8099-30433d7af0b1"], Cell["Graphs", "Text", CellChangeTimes->{{3.513266788083189*^9, 3.5132667973147306`*^9}, { 3.513272093815999*^9, 3.513272098697585*^9}, {3.5172389588536386`*^9, 3.5172389629386263`*^9}, {3.517239692171343*^9, 3.5172397065137873`*^9}, 3.5172399022732773`*^9, {3.51723998911631*^9, 3.517239998600172*^9}, { 3.762184806607835*^9, 3.762184813214233*^9}, {3.762197448350705*^9, 3.762197451565967*^9}, {3.76219972831061*^9, 3.762199729750489*^9}}, FontSize->16, Background->GrayLevel[ 0.85],ExpressionUUID->"9bcd96b6-baa5-4f4d-a7e0-c4737b8c9aba"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Grid", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"plot0a", ",", "pclosed2", ",", "pclosed3"}], "}"}], "}"}], ",", RowBox[{"Spacings", "->", RowBox[{"{", "4", "}"}]}]}], "]"}]], "Input", CellChangeTimes->{ 3.7621997408303537`*^9, {3.7622250966193933`*^9, 3.762225115551784*^9}, { 3.762227833420183*^9, 3.7622278487272663`*^9}, 3.875753352956853*^9}, CellLabel-> "In[3100]:=",ExpressionUUID->"13b557ae-8017-4685-b51a-0df23aa5523c"], Cell[BoxData[ TagBox[GridBox[{ { GraphicsBox[{{}, {{}, {}, {RGBColor[1., 0.7, 0.7], PointSize[0.007333333333333334], AbsoluteThickness[1.6], LineBox[CompressedData[" 1:eJxl2H04lPkaB/AhKe/CphcRNch7h6yzqZ+whY7N5bBKtF4uhxShFUmH1apl RamwaVtK0dZGw5JSPxdaZ4ZFm7znbRhjmDEIcRydmef57T97zz+uz3Ux8zy/ 577v7z30g055hsgyGAxTGQZD+vPPV+ikFoqLKcCeV5bP3nhVjf7quBUUqXKp F9FmA99xEh1r6inDe0+e32zl9QZYvVThj4tl5eT3u4HrHJvu5iS2I+uePZK/ GAT+zDzukF0kF9HvxwVOlLvMzZjpQy4rNZzhDB7wSn3Qcs/mDvJ548BF+4Jz nR7zkZ+B/Q6XlUlgr8XUW24uInI9YuAjBSNVKhuEyO7uu/qOqmngivzfGczS WXK9s8Bayxb/TJ8YQ067tpxOE78HTp75G67wGCH3Mw9so7bae4A/jdwbpVf8 AfhcZ+tOgZ2I3O8icPRXSf2VmiLkc+TWa1Hgf4HVjz6sTn8qIOexDPx7lb9c tGASBU70JhfeXAEuv87sc06cJufFcPirHXlFv4Q5ClDlS5UNP49D33d5zvqm fIicpwzw/pNZ7/P/N42eX/rl1KqLssDv94R/ZxnTQM57FXC1e/oFlTYBqvWQ npgcMP/C+EJsczt5HtAauhGpXm089GqjULfy5WrgGtcC01JPNnle8sAJBU9a 44uGkARn1I6uAVa+mtsQwxoiz3Mt8G4zebvJP4pw60OzlrAF6Gpvm8q6oF7S zwrAJSZeHBZqwG+/bmbWXVME1uRaP7bwacF0PSgB3zlj5v/ybi+mHqeVMrDO 3NvJGGM+pusF+tnGpeJU1jgekFeWfKQKcP3l6VkHVR6m60kVOIFXcC9J4hHq htWALXIXwtPTRZiuN3XgD0375nZljuBY6e2mQAd9chD/uIqP6XqEVj+rreDQ MIUXyuYiaj3WAXPMHgzNPRdiul6h5ef8Y0rs+TjhhrRiNIA5etkXEp0GMV3P 0GscRBX9PdNY+jQ3CqE9nBILipXnyTzWBD4s4yVz3HIeJx0Lid7/HPpyk+yP F1O5mO4HLeBX4nJmmbUQy1ANAb27SKksy38M0/0C/aj69kMj/gz+1mg1+yef T4CdNdzafDO5mO4n6IakbtXtBcOYKhfmemCbLjOrlGt8TPcbdOCKebDTKBen iYu2Ls5CO1gqf3TdNozpftQGFue5/JDpW4Woj6uDLhzzVvRC/Zju1w3AVnyz tHa5X3FmtXO85xVo+zVeG94N95J+hl4cOM3a8bYTr7uttT3p2EZgX+OkwGX9 btLv0NcTNNSTderQtRRpx2wCVo9aUziROYroeQAdKT9+wCeiC60PrTjXtQSt z9PjXJvrJHmyGbiSLS4OcZ5AeQelTxz6ftyItuXnA4ieJzrAg5rmSlej+pB0 2diZC21WnN89KPOE5BW01kz3E7ZiO6KOK2QL8FzzC2VDFg/R8wg6n/FExr9h EkmKxeR7a13glo5vCvl9YpKH0G42ripXC/mIjjM94IpR5tCBFg6i5xn0t/Wa DK+cKSSdViOt0DWveYbBHWMkb7cCx+O45pxJIZKe1tpD0Iad52r0TMWInofQ LzTrDW7YiBA1riqhGyoGPG2tB0meQ2u1Jd7+YksrktysZCToA3+Q692tFcFD 9DyFXl5xqC3cOYUkbyYpGegD8Stqlb0Csi9Acw469J407UZUOU5Ad5w3EkaN viHz2AD4SkCB4dy/uegm1bDQNnfiFLviMNlHoAcV1l1sSk/EfZJqyKyGVk50 yvIbHkD0PIeutWjv+f5dLZJO06Kt24APsTltYToDZN5DDz778tmJ2DeYGveX oGOabLc98+OSPIDWZwf02t7qxj9RDQ1dtmk2vMS7kexT24EPaHxgbu3j4yHX fTNj3tAnR0s9tIeaSZ5AhxbFyv38nzFMtXMNtPk/fv1yuUVA8gY6QxjNOp7P xsEW0onJBP57gb9so/EEySNoq8hjD0uejuN7VMND+1q0HrYcFZC8gpbl6zkr x09jnob0jqG/YPmE/SAjJnlmCPzI2k47q7wPU5dzBLqr9ETyBZ6I5B20AfvB 2lXXZ3EY1VDQVscN11fEDZM8hL7TqrbpX3odmGo3IyNgbZUgn9OaUyQvofvv 7s0JXTeGqXGYCa1t6ccrDhBgOk+hZYR+KXmccSxNw7I56GYTU4YwsInkrTHw +WzdLoEOD0dSBQmd5et7mNk0RvIYOm+teoKmAwuXsqQTBTo8wNAnNKQB0XkN vZDbMzWv2oap4zTZAfzdrrNVtrVPEZ3n0JkRaY9tXjRi6bRWzIbWNajem+HI Jd/3oO3dOfvtC9mIGheL0KUf5ZTL44fJPmACrFRsWeasMIikaW8bAF2mtJuX wuORfQG6tir2FXP+MZZOc/dG6GiLlSXHU6Nkn4C+OaooetTRgWwDSu4HW5gC rxbbVxm7t5F9A9rit8/dCm25KJ46EGh2ho953NERso9An3lQ+O7Te3xElcMy dHJHU75n5Xuyr5gBH6k43tz9egYtURcM/dlNhz0eiRNkn4F+NMr2zwlvR9S4 40B/SJhXDcqeIfsOtFqaPN/q0xlE/bdmpzlw1KWrvtdLJ8g+BO1adiK7EvER Fed50G7TNalvj86QfQm637Rz6euhefSRekFHRu19cWgPH/0furgcJw== "]]}, {GrayLevel[0], PointSize[0.007333333333333334], AbsoluteThickness[ 1.6], LineBox[CompressedData[" 1:eJxl2Hs0VXkbB3C5Nk0qjdyrUwk1x4lqmjTxk1Sa4XVCIZV7vDXIzBmjVNKI GHSRldcyKcUyJUxFifxcRiIpl3RcwnGce9ErlyGXOWfv31895x/rsxZn7/3b z/N8n2WFX5hLoLKSktLXc5SUFD8//7hcnDqeVluKPveHWL/0pKYgTLse2L1l cOR5XDay+fGUoYVbK3B+9Sq35vuZiP79DuADNuuda9uvow2d1vK/6AM+HWIu WrW7gHwfH9g2zCL4oko5cpgpb+hPEgIzXC6FbWNWk+tJgK3HrH4Kv/c3OrBy 6xqHmXfAmew3ZqNNz8n9fAB24My752HSgjbffFvT/vD/wDmOe9U9616T+/0I fHT00MkT3Z1o+zdLf074MALM/SXhauu1bvI8Y8CCP7UWGXq9RU51ijv+B3in zcksve/7yPNOAHs5ritQX8ND7p6ZzYO+n4Cr7Iaq5jfzyHlMAdezHuxnTfKQ r6zrzI2MGWCp2D88nNlPzkvJ9nP3FM1F5VX9qKRCU++2BDpkeuRG/6F+cp5z gHlOrolxOjxUFn83TCVOGbjIoEjl3WEeOW8V4ObQwwZmqW9RJVtxYqrA+7dP nlE72E3eB7RHSGNMpS4X1eq/X1ZSoQb80LipcXpdO3lf6sC60z5Jx12bkRwR C700gCsHa5K2nnpB3udc4IGJ4o1KOc/QyzvMpuBxaJstP2fqxD0m7/sL4G6X obDA2ST8mtO4ujp1HvDWR1VVOdalmK6HL4HLPk63+TrVYep1WswHdvCOLduc 3IzpeoG+qx/cFnKHi3vV58svqQk8omsovuHQi+l6WgDspM3MnvqmHw9QD7wQ mJmQefCg/wCm620RcFAbM3D5sBD/onjcs9Cbim29XzHFmK5H6Ma0wpcmXhI8 XjQaUsnWAtaeV7eooE+K6XqFdrW7tU7SJsMn0hQVsxj4sbWhstGQDNP1DP0V axqnq8iw4m3qv4e2m3nixrwuw3S9fwX8+4YvD6vde4ejDwWG7yyD5oly8+sj 3mO6H7SBMxOlAVP5MjyHaghol/o8ptYfUkz3C3Seesf2vi4xjjVVq89yXwIc 8XCmL9RZjOl+guZw7t5mjwkwVS6rdYDztc3YV2L4mO436KwBY/ecDB5O+HCL MfER+vCypTvr0now3Y+6wHuGrKQvajswdblq6OcWm1MYns2Y7lc94PWqu2LU 9tbjlFL7SJeL0DqhKxe+snuE6X6GbhiIzG03v461rmkbRx/SB65wXc7tPpZF +h16eLwwRvW3YpR6VtExBsBCdtcNy86niJ4H0E/sMz7uErcinaAHUdxJaBWL Y9PVbzpInhgCm4So/kc00oPSf1C8ceid2fXPj2rwET1PjIC9NvR7FC4QIMWy YXkV2rRzd5pDspDkFTRb+Vc/z0tCRB1X4FLg1Z7D8a9mhYieR9AlRakjpdvF SF4sa3/fsAxY5N9qlvFIQvIQOs74hSbXVoboOFsOzNGcWMC5IkP0PIM2937P 2v+PBCmm1cBL6LYw5fHCH6UkbxnArMWt3x4fkCDFac11ht4Xf/7CTJEE0fMQ +qYllz0aJUHUuCqBrvce8jbuFZM8h/bgm6G6V0Ikf1j5SFgBXG54xJj7jo/o eQrtvL6r93/yfUD+ZfKSgbbPij7+4DKP7AvQTKV63aOufYgqRxl0ZIhGtFZi F6Ln8Urgy0Yl51rPtqAMqmGhT1z3UddMbiL7CPS4875ep95q1C2vhpRS6CGH deMX1W4jep5DR0bseDK5PBoppuktxirgK0WsDm2rIjLvodNskn91DKjD1LiP hy5YJDt17VIzyQPoqKCxOMcYLs6iGhpaCVnXDHPekrwwBtas7Brb4tiDebu3 DYv2QjNyWRtvLekneQLN3fcs0kiFj6l2Loe+777x0x6GgOQNdM3VLTt6GCLs z1JMzNXATKfyPN8BIckj6LiaattYezHOoRoeesRjdNB2jYTkFTS3Se9p6GsJ Fi5WPDH01eACuzqZlOSZCXCA3vd+vwXLMHU7ntBvTp32CpiVkLyDzrIv85oR SXAw1VDQ7DMbzyxwk5I8hLZnB2yryRJjqt1MTYHbFmroHlwhJHkJfcSN85AT JMDUOEyBdhMIey3b+SRPoVXMw8/fjO3HijQsGoUudvbdn2LAI3lrBtzQw9KP knbhUKogoQu+zdvNz+CSPIb2M2otTSh+jQvvKSYK9OIdg38JdJtIXkMf8F0V lVvzBFPHuXYN8FRNWCaDn0PyHLoipeWSWkg+UkzreZeha8+dZAoSK0jeQ78x zB2pSX6JqHExAW2qk6GX2tRG9oG1wD/MHkkt9OhCirTf5AN9xUrNN2lPH9kX oO3s9L7jWvCRYpo71UG756jY2PvzyT4BrWH+U+uz/wrQJp+8XH/W18Cqk+nH zhoLyb4B3e1Sy1DfJkSR1IFAIyvd2DtLRGQfgTbIlgQ13BchqhymoCcizte2 B4nJvsIENo7gX/jDRYomqRuG1r7rleh4Tkb2GWizNP/4To4MUeOuAdrXpGry 9AUJ2Xeg2bZOmz6NSBD13xpLc+AZ52xB41wp2Yeg8/XnPDbUlSAqztOhn1b7 mE7ni8i+BN1dHNeivUuEZqkP9NPSigz1JBH6F2aODOo= "]]}, {RGBColor[0, 0, 1], PointSize[0.007333333333333334], AbsoluteThickness[1.6], LineBox[CompressedData[" 1:eJxl2Hs81OkeB3BRKimtlVzWpTCup1RyQjx1qEVrk02xW7ZpQ51NONK6bW1C 29ocNSu9HJ02l+ooUtLSyiNUiFza2cQyjDH3GWOHXEJn5vd79p++84/X+/Vi 5vc8z/f7/Txj1cHooHBNDQ0Nx3kaGuqf77+CsmcSc55Uo/ddHPr1ZtHzakS7 GbjUIziyjNuMvI5+a+q8+yXwvZ/WS596viS//xr4Utzv9hKX12hDj6fqLwaA j6/kGyj7+sn7DQHb1/xi98Kai3znalq4P/KBezeHna8r5pHPEwEnYu+Tn3Ty 0b7Vm+1956TApffD9/VWCsnzKIBr/JHI7zMx2lTY1/D7L6PAmw2YV9i3JeR5 lcA9PyQFfdMsRd4bzeLOKcaAzy8cSCkpkJH1vAFOHpsyxy5yFPBM/cSTwIrS nU2RP8jJeqeAlTyu/dsbcrQ3NL9TznwLzJQLw5P+JSf7MQOcXD80UtUjQ0xJ 73fX8uaAN/V8v2+cKyX7pbHlfQeYFOg2XZWgB7VLjUpE0CdNPR+mWP21n/OA M+qi9vsxhejXs6XRWhmawAVozEo3gE/2Wwt47J6XcsXsEKoLVO/YfODSdA9v x4RBch7QYT/bZOoF96EnxjLzB7ULgD/v1NT2X/OKnJc2cJj+Sb3qqx1IhRN6 XywE3tN7NPOKWyM5z0XAb99FyUqYpaj9ltOLwxPQcYcY7iGK65g+78XAHo2G 3H3Z9Zh9vNWmnqUD3H99p9O2Ze2YroclwDnx10Zn8tmYOk5nXeCQC51nXOJ7 MV0v0N8vc98/d3QAc7R1VR+5FDgxNiRXPMTFdD0tA+7/8h8sb6thzKMWrAds 3XDYv8pFgOl6Ww4c2lB3ydFOhOPVy02F3nPZNiVQLsZ0PUJPBR9ftz5WiifK x6PqAj8AFgzEtPlVyTBdr9A+ZtrfpD+R46QcdcXoAzP+jjmT2SOYrmfoZKVx cJGWAqtP01gGPWXmdmHOUIHpev8Q+OPMAoPp5Qp8Kiw8dvuv0FvOaK9qbx3B dD8YAPtrVSdEbB3B86iGgF6ROb9/2Wk5pvsFOueO808VCTKcZrug+ereFcCN E+xRoacU0/0EnVbbfee3RjGmysXGEFgXLTmTOCLEdL9BL/7Ze7StiY/PKYos p5TQKVEJ+l5HeJjux5XA3JfbWWtfDmLq4+qh3/bbvLO06sd0vxoBpw4t5rOP d+Osap+EoGzonfUHl++o7cJ0P0Obhj83jFU24w/+a2B9KswY2EjH5zbXpob0 O3TkDGPcWy8Ls1LVHWMCfIj56SIm+z6i5wF0o0PZPP3Qp8gw8n5y9zT00byg VkZdB8kTU+BvY6biORw2urxDfeLQNSFFYcfyehA9Tz4C3hrT6mdtxUHqy8a6 XGjbSXlZeeMgySvoYp8S16qmIURtV7gZsEnGf2oTsoYRPY+g/R3+nRThIECq YnHI3GAOnHhrJDqSJSR5CO1uxtp26LEI0XFmAaxtpzn59KEY0fMM+n8MC3wo WYLU04rXDm23xzvHYtFfeWsJzF8ytbvQR4rUu7VoJ/R3Zhl6bUZSRM9D6B3J dzu7r0kQNa4eQI/38hgXxGKS59BvxvsKCv8QIdViVSNhFbDB5duJvDAhoucp 9OOPe1tc0/hI9WaqkoGujjOpK4rhkfsC9Kmok9ud3bmIKkcJ9Nxhrv8rPgfR 83g18KPT53yZ83tRHtWw0PYdGU9ustjkPgL95Y3TmbvOtaM/VNWQVQ3t6Plj ycD9BkTPc2jDT7WK3iSVIfU0LbK0Ak5+zm9XzhaQeQ/ttKto4Ym9jzE17s9C xw2KS1372kgeQPc9K3uYK/gNX6UaGvomo3JiQVIPyQtrYNndlZURTRw86Lf1 T0Ew9MbPGFU3U7kkT6BX5xu8FkzzMNXONdAa5wPyHVwFJG+gH60vrLdxF+Gv 1qgnpg2wab1vmlJTQvII+hWrhbM6W4qLqYaHlloa3igUyUheQXsx4o8cmZBj vr56xdBZNcOxWypGSJ4xgMvtq96xPlRg6nFCoa+g2dbn5gqSd9AHzppkSqQj +DDVUNDTF/GjZao8p/MQeib92Odr/eSYajdbW2ChBZq8u0JG8hK6mtln1Fwi wdQ4zIKOrYjQYyhFJE+hF0d1uFfLBVidhuXj0M5FwRedLw2TvLUDzmWnN3Jm ufgYVZDQMV/EilmfDJA8hu7TEbdkG/fiO/fUEwX6XmUnW2cXm+Q1dJq8woip /QJT2+lgD/womhHpNllH8hx6bUaz5Q10Dauntc5F6Lw604vnDcrJ9z3obW7d G5kmjYgaF1PQw96js7aF7Yi+DzgAc4pLTPdHsJE67V0PQHedfth38GkPuS9A 165P0birMYDU0zzgGfRuEyfv/hNccp+ATuFr87Lyecj1wM3rX61xBB73Tfx6 lyoP6PsGtMWwJW7zEKIEakOg87oSrA2aReQ+Au1261JlkLkEUeUwA53Ntdi6 wUFK7itOwO5ej9fu/1OKpqkHhn4QvrtrTbiM3Geg/+lSETedJkPUuGuBDoxO Ty2KlJH7DvRwQ0+G34wUUf+tWfc34K648VyGm5Tch6BDOg4GhiMJouL8MrT/ lqXTgo/E5L4EPeZhsSm7VYjeUS9or3KTsKV+AvR/ahEuew== "]]}}, {{}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], Method->{ "OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ Identity[ Part[#, 1]], Identity[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ Identity[ Part[#, 1]], Identity[ Part[#, 2]]}& )}}, PlotRange->{{0, 14.999999999999998`}, {-1.5, 1.5}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, {0, 0}}, Ticks->{Automatic, Automatic}], GraphicsBox[{{}, {{}, {}, {RGBColor[1., 0.7, 0.7], PointSize[0.007333333333333334], AbsoluteThickness[1.6], LineBox[CompressedData[" 1:eJxl2Hs8lGkbB/ChVUSRLJFTBzHORb3qXd1sRJtdfexa2by2TaOD5LA6J+W4 JDkk3tYhpehgE+U01i1mtKNipnEYx8ZhB5EcopP0mue533/2mn98vp8PM89z P9f1u66xYk+AG0uWwWCYyDAY0p//f+0bUUPHgq9it8SZE6nccvRPH5tFhxfF dCDaPOCFpx0d1UWFePOh0OWWPwiBTyc7rr+VWkx+vw3YemB2x5s9Tciq3Xbu L8TA+uJc7cvufYh+vz5gydMyr5icTuQ8W1nfGy8Bbsm5abUyrZl83hCwvXzj rjsOg8hr5VdM59kR4JIg7sfGu6/I9YwB52lsvWNkN4JsrnfVtpSOA4vz92y8 qzNJrncSeOK266dP/hK0Zb3Or7Fjb4CNkh7ZtezvJfczDdwrtA3OqRhD3z6W XvE74Mij3sKsg8Pkft8D333Wgms7XyIPzwzB6C8fga39t+xrWj1AzmMGuG37 +tZkxSH0y3DH2Zwrs8DMt6wGT8koOS+G3T+tUZFTK9GUoJKqRctuD0GvaPh5 m6ZBBzlPGeDf5d/6sBljiB1TEDAvWhaYkTa+InJ+HjnvecDvih9srEgbQNU7 pCf2BXAUe9C93biRPA9oH0M/j+8+9SOu5ivdkio54Kbzo6Y5T7jkec0H/up3 36n/tPWgORxV3rUA+GRCWlK0Xh95nvLAErfF2mqeeajxjmnD/rfQHu2ViREL haSfFYC/7h5gpAmSUXPIU4OalIXAzFMnu0QPYzBdD4rA8cpspSPCWkw9Tksl YPuXvp0aFh2YrhfoMxKGw6Vn7fjFfKW5j1wErP1kd5iNAh/T9bQY2Jx/7DVT j4f7qRtWBn4em713UL8d0/WmAhzpEvLahZWHj0hvNxx6KvCwm/VeNqbrEXpn s9DPuUCI3xZO+VfvWAK8jWvmxSmsx3S9Qp/lGM50WpzDJ1OlFaMKPH5g8yiq 4pN6hv7aLHBvnUczlj5NzVfQmwIyXLa/6SV5vBT4eaTs1p/8+nCYNytoKxta X3dJ3hlJLekHNWBOeFn+iVgulqEaAtqSu2fNRNQFTPcLdI5gJEVfTYwjDeV4 2R5fAovN5AzUzNNJP0EzohT/VaW3A1HlYqAOXFi/0zVnVSOm+w3aq6A0Tauc g2PHcvXfT0JvYCYpHayuw3Q/agDbmJxqDvu7DVEfVwMtez6u7t8pDZju12XA 1zbF8vsC76OEcofjbonQo/F1zMTnrZjuZ2i/9rj020FteEmW2uowb03gyh+H R5ibxKTfoeM3mut9YN7DKeHSjtECXikf/yejVETyAPrqAnF0YXgsUt/34JTo A7Rgp+X6daIikhfLgZUiFPyvyzSj9O3SJw7d/WF6yarTbJIn2sAO3tyaa0uf YumysTYN2jIlP+ih7t8kb6APaV9hVSl2Y+q4WDrAKrbXzgQoYJJH0CK7pooT q7LRXLEYn7fSBbYPjFiTpcMh8xC6WMDWCi+uwfQ40wP27jbVdJ2QkDyDjta9 Jx/wZwGSplV/I/SXdutCnEvqSN7pA1+0Hf7+ZGgmkp6WvCv05OZ5fWah1YjO Q+jgXRIXJ+cyRMVVCXRZ37C7GUtI8hLat71/LChXjOdudi4SVgCrnJNMeqtF kjyFbqh89L6vtAnNvdlcyUDzjJjfsJz5ZF+ANlQ9+p1LSxWmynEYul4jQ2Ev Y5Dk8UpgWWbm+JWCRnSFaljoOpUeny3WHEznNTSvItbw8+J63DlXDQnl0LN8 nk/PExGi8xw61SJowfXrbCRN01z9VcAHKx19X2l0kLyH/lQW/H3CIg6m4j4G +rjjRMMDWRGZB9A2ZxeY+n8uwdlUQ0M7aA1dCLfvJfNiNfBlDVXOPK0m3LPN fmLAHXrptcVcQYYQ0fMEetrIMrPmYw2m2rkS+pGFR0LD0GMyb6D52yUhcdMv kI+5NDENgG/xQuSibrDJPIJ+92DQRDY3C9+gGh46Y4zLXJZ5A9PzCrpc3klO t1aAJarSO4aONfvL7b2wDtPzbA3wH7/5JJgPtiLqcjyhD4WlvWlwryXzDvrj kPNLO+MOvJ9qKOiSk8uPO7J4iJ6H0H8VlG7+4bc2RLWboSGw0NToHX+YT+Yl dLViXJ74v39gKg4ToMt0ZFxv3nqC6XkK/VNr9mObQD6WTsPCKeisjNLuhz8K yL5sBBy+OW3Q+pIAH6YKElp9wkllOqUN0/MYOqL3pquyMQ/dK5ImCrQoeqmJ 084Wso9DB3da5vXI1WDqOI2ZwNOiso8bv6hG9DyHttMo0vQsb8DStF6YDB05 XhIUVdRFvu9Bh/QHKrNOXcJUXLyHVr33TFHjrgDR+4Ax8NVJXpmtLEbSab9h N/ThAxe9RJl8si9AX12293l5ay+Wpvm3j6GjU5OcBluKyD4BLbC+JNRxEOIN u/Nv+pibANedvZ/+s30X2TegUz+9TknS4uDj1IFAq+sPbROacjC9j0DfPtNl 9br4FqbKYQZaZPFizAR1IXpfMQW+U33haaRhM/pAXTB04sYp9tEjcYjeZ6B1 o81Gp/LFmIq7emiV5a6RoXuayfd16Ba/Xy/f/6YZUf+tWWsGHPx0RmftxRJE 70PQ1p991lklJyJqnKdDBxtnrxvP7UT0vgS9JTTinO/AIPpMvaB93YZYKgc4 6H+RwBVz "]]}, {GrayLevel[0], PointSize[0.007333333333333334], AbsoluteThickness[ 1.6], LineBox[CompressedData[" 1:eJxt2Hs81OkeB3CFSsilVhwpFIZS1OFQzj560caoV+1Iaru4H9qoLSnJvkQk NSSyJEduI5dE436bR2ZlXYfjWuOSpRIqZNNY2mN+v+f80/fMP17v14uZ3/M8 3+/n+wwt17Msj6USEhKbl0hIiH9+/WJFz1+OqytHX3sy1DWB3eqJaTcAR7hY MkzZaehb75/VjQ51AB8/uHroB5skRP/+c+AJF0EZVyYF7Xjxz8W/eAncdX9F hQR+RN5vGPhl9OC+AFyBbL5UNf7Ofg3cJq82INB+Sj7vLXCB4ttrinJ8dFzb Qt/mywSw5vYuzkqXBvI8k8DDri0JURfbkFl6P7+7dArYPjh71C+4gzzvR2CG fUqby1APsjLR8I2YnAHWL8xMVBH0kvV8Ar5fNRO6Ma0X7a8XP/Fn4NSKb0IS Z1+Q9YqAt/nJb9AbESLHo0nt713+BG7OXjX4KKCP7Mc8cG7Y2mMppn3IZVx4 NTXxC7CTz8Kqa3/vI/slYfm1T8nI5lzY3o9KePKqOW+hi//RoSe80kf2cwmw uvye/9w0f4Eqw/POSl5fCrzP9lR9tJ+Q7LcksMOuyuvB2j2o5qB4x6SAzTSk i+Z1/nce0HNBv+kUnetEdWrv1pfwpIFr9/i1XPTuJOe1DPj+dGOPVXU7WsRF hWPLgXdOu7ysS2gn57kCuIvdv0mitx0Jcre0es1CD8Uy7kV8bCHnLQNsYXze /kPDr6jrQrNObexKYHOnK+iLRQ2pB1ng4Q70xG5/GaKO00gO2ExaelVOxiNS L9Catc0XfvtkjweXyS1+pDxwgcWzeX5ONqbraRXwBgVb/q6AQjxCLVgB2L2U o2fOLMV0vSkCT5xSnp1q52E/8XJDoJvPqrMa8qowXY/QAbsURxy9K/FswR8+ NQeVgFX8mMbMcB6m6xV6m1tpxJJRHg6IE1eMMvDI3gO75t0qMF3P0DwmO2z0 fj4Wn6baO+i/TBUkrZ4XYrreVwNfzzO085mpxkEnPc59VwmN7IPiIwf4mO6H NcCsTh07h34eXkI1BLR/cLF+46NqTPcLdLKelRR7oQyH6kk3PHD8BnjrgFug +gLGdD9BB1/iv3r1UwWmykVHBVjv8q3kjtJiTPcbdBFbV2CtUYojJjM0RR+h 5Vgm5qnVJZjux7XA7tlWE8+EJZj6uFrof/lZmNurp2O6X1WBdzTlf8/xysJR 5db+rGjoBpnn2WaSKZjuZ2jPMqOefqU8rJS8ZlPQSTXgKHQ4V9enmMxraJGa 3OPPTB6ODRF3zN+At/AOtX56XYnpPIBOluVdba3PwyqeRVd656DZrgOCseRM TOeFOnC32r8zpWIycIKd+MShndumWq2dHTCdJ+uAq5R6uO8tL2HxZcM4Hvpq 9Z7p9Sa/YDpvoL31VJkOgkJMbZeHBnCgcCFDPamK5BH0yBmTOXOLKrxYLAa3 dqwHVvS5EhrYWkryCvqS1hHjZu8CTI+zDcBZrLHVXrcKSJ5Be8YKpdYq8bA4 rUYE0O63wz0NKspI3mkCjyYaKEYNlWHxbq04AN1Sfj4sd1sRyUPoph+HQs+I cjAVVyXQ/j7uAxE/ppC8hGY84ca9ys/Gi4tdjAQt4Ak1qUTv0BKSp9Cq5/XW 5ZcV4sU3WywZ6DrbDJUP+9JI3kL7z73hDK7xwVQ5jkMHhzOGihghJI+1gZUs J3fq3OHiRKphoY9JT92MYHNIXkMLI2QLThtm477Faogqh/5eIdnd5vFjkufQ TUdvyXTGRWNxmmZobgQ+4v9tZHysFcl76M3LmUfl16VhKu7DoSXmft7xuTuT zAPoZNGJF/xYLn5ANTS0tVOM73JfLpkXm4BrXEsGW3AYHrLdPf3GAbqxsrjH 8UMKmSfQzpbyYXd8mYhq5yro5QMm0qzpEDJvoK0FJ+s6euKx21ZxYuoA/77T cX38p3uInkfQ8ZZqB1Se3kYcquGhVeNPCMzDb5D7ILS5F++pU2Qgeq0sXjE0 d1yr+wg3lswzXWDvdP6sTU4qph7nKLSRLqveM9eF3DehNbtlhrd9F4a9qIb6 Pw73rYk2yiLzEFp1u3JKbeA9TLWbnh7wSHG50jnZ2+Q+C72Pw+iJ6wrCVBxG QW86jVnD8+FknkLLXQyzeGedjMXTsOAP6GOnc/NDTmeTecsAPhx2W5Wz/y4+ QxUkNIsjtPX7NYPMY2hR+vuWmd1cnM8VJwr0IaeSprntD8m8hnZ+mJv60uwy prbTQB/Yz/DPcwnMSDLPoWtmRHaidQeQOK1XxkA/WcOXSVqIIfMeWjfax6l3 NBFRcSGC9jjRV2d5NwHR9wED4K6fpmTyTLKQeNqbOkP/or/+zg8Pcsn3Cehx xnHGMmsuEqf5/nroGo+xhZjSa4i+T0DLzTzM6lMNR6bOWZluWzcDOzTXj5se uUvuG9DJm+vffObnY39qQ6Aj+zr7JW9wyX0Eevfc471uHkWYKod56N7J+Wa+ fiG5r2wBltqr7XAhJx7PUQ8MbTlZvvRmERvR9xnoZ9c3Mg07YhAVd43QSZY3 NJuvJpH7DrSXwnSVqsgPUf+tMTYEHtBa9m7sSRqi70PQbXE3lL0xB1HjPAE6 /7CnfW5gBqLvS9AawoHpm8cL0V/UC7rp0phwPP0p+i8QFQdw "]]}, {RGBColor[0, 0, 1], PointSize[0.007333333333333334], AbsoluteThickness[1.6], LineBox[CompressedData[" 1:eJxl2H081NkeB/BBnhcrIkVCyEOoFlubPUWJRXdJW57aJEJtG3aLKGILt0jh 5uJmJKQsrvIwDGcaNMZ4bivKdstKxqS1aGWbWPP7nftP3/nH6/16MfM753y/ n+8ZBge/9wqWZjAYFlIMhuTnxy+vDHFMdhsLfexinyNbhAIWos0Hjikwm1EY 4aMvj55eaeP9ADgtTeqcmv0D8vuDwArLlR67rBxEG584LP7Fc+Ctp7dVVpQ/ I+/3G/B+pkewm/gFcplndwxfHAXmWTDMDgWNkM8TAmudr3S8e3gU+RtuMXOZ fw2stq648cCKMfI8k8BKXObQoVgh+rzo15ZHdX8AZ1pdUB1JHSfPOw1cvylP PO4mQk62elGpkzPAtiXH5gUNIrKeP4EfnT7i4zQgQh48yRO/A55zfxamyxaR 9c4BW+3erKwQIUJ7ffL73gS+Bz7HOKvR+myc7IcYOPmbRu6T90IUKHqaUJg7 DywuLYps6/3/fjG2fmzDN9cMnU++QrXNKstvCaGFrPXlydKjZD+lgG3mtmxJ OTyCGpN//l7mvDRwNv9d0hhzmOy3DPBGnbAybutzxPlasmNLgG1fnmV5sH4l 5wG9UG1gLVM5iNp0JlbVNssCh0yVybWEPyTnJQf8m/yYb5JiH1rECTU/eWDl vqP49QSfnKcCMK/gQG94EBf13LbsDp2F9mzsiO4dqyHnrQj8S0qYgZEJEz38 odOYm6kEPOmX9y0z8zKm60EZuNU261G7diWmjtPmE+CkHRtW5HqwMF0v0LZM qQDlcxz8P7lPFj9SBVjdOOCidlgLputJFbi55kPxpsg2PEItWA04rkVmZ8r4 fUzX26fApZMD98y2teMfJctNhHbfXqb70p+P6XqE9l61MsZHtQPPVr39jvO1 OnBjtruy5vUOTNcrtKNe0HCEnQCfypZUzFLg9FT7rAiWANP1DL32T4cO1lsB lpymzgS05ZV99SNGnZiudw3g5oNnqv9gdeL4/cERzo3Q5ZwPoUV2XZjuB03g gc7OGS9OF5aiGgKazXN94hfSjel+gY62T3DSmOrGP5nK8gv2LgN2nEyXOlHf g+l+ghZH8JK2berFVLkYawErbAzzwK69mO436AUd+dDrRr04dfLG6rlpaBND R9/NAz2Y7kdt4Iwd/Dey1j2Y+jgudE2MedVD325M9+ty4Lyazwt6KrpwOmt7 tFcGtFTGLSexXxem+xl6wWduONKyC6tf01wTv18H+Cu1PmdmXSem+x26Kdrm 8RRPgDMTJR2zAtg9/tlP4oMdmM4D6IIXNecihO1Y6/Dd2IG/oPtbNmhNP+KR vFgJ3I6aBjN28XCOm+TEobc3xr6/osEjeaILnFh8fkdcxX0suWysvwrNQndP ygjbSN5AH1FC3tZKbZjarmA94KlNT2aParSSPII277N8GZ7ExYvFYn5h4ypg t5xOY8U4DskraKHf2qjwUDamx5k+8JozhYLefXUkz6Az93okcd/9F0vSaqQH +iqLIfJtKiV5txq47jNR3kzCVSzZLYV/QM+q3quaNDlG8hBaE3cammSmISqu aqG/2jOd4ejOJPMc+nb0teLnKqVocbGLkWAAfOay/urX2rcQnafQrMjqk8fH bqHFN1ssGejDHxR2L7tym9wXoHfteTe5JK8cUeUogpZrSDOR8q9AdB4bAk/c bGAajVegXKphoflz3OSUoXJyH4HOelDIKDpVhoYWqyGdBb01P2hXmX4JovMc +suprJ0q7teRJE1vrDYC9rFPUuhxyyX3HeizrJB/rw9IQVTcJ0NXlKg1/T5x hMwD6OLXH25w72TgAqqhodNnecs89XLJvFgDXG7aFebpWYhfuG6berUHWjbN RpuJS8k8gU7Ur2p+tbQSU+3MhtYNOe6qpnKHzBvoSxfyh2w4tTjISpKYxsB1 A12awhMNZB5BR/n6WfUfasbFVMNDf+ec5qQdc4/MK+h94UGX4ota8OhSyYqh d6cVPrbOaiPzzARY8Etq9ZpgHqYexwd6it3vq1/dTuYddNjNnxu1Q/k4lGoo 6OOWribGHD6Zh9C6ehb1Lg/5mGo3U1Ng4Y3WqPznfDIvofefCvhsXrkDU3GY Dm1obTF9bZpP5il0f4VKxMxxPpZMw6q30PWW7Athru1k3q4FDk8VJDqU38fH qIKErnolCi1Z10bmMbTIyjhVQdiCK6sliQKt13y9/M03XDKvodXsvuWUDWFM bae5GbB2kHe67c0GMs+hY7VqQ3RzarEkrZWuQLv8sy+/2fYOmffQSoGc4EuW VZiKizno3ki98IH/3Cb3AXPgKF/1u6P+pVgy7e0OQMsHBl72di4i9wVonrt6 dfdMPpakuQcPOm6d+Onbqkxyn4D2rzDbUGB7EdsduFkSZGUBbDAUq1/JTiD3 Dej+nfKq0vGmOJraEOiBSGymmJuC6PsI9MLml/+Ka8xDVDmIoQWP4/S39hcj +r5iCTwc4PYAm1Sgv6gHhl4qo5dQlF2N6PsMdNVT2ZoO2xpExV0HtIM9w26O WUe+r0NfdPji0yW5LET9t2b9OuAEa13FlqcNiL4PQTsxvNjShmxEjfMc6BNX tVF5UxOi70vQg8pF+iEpGC1QL+hST1+hXA4H/Q2azQSj "]]}}, {{}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], Method->{ "OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ Identity[ Part[#, 1]], Identity[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ Identity[ Part[#, 1]], Identity[ Part[#, 2]]}& )}}, PlotRange->{{0, 14.999999999999998`}, {-1.5, 1.5}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, {0, 0}}, Ticks->{Automatic, Automatic}], GraphicsBox[{{}, {{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[ 0.007333333333333334], AbsoluteThickness[1.6], LineBox[CompressedData[" 1:eJxt2Hs41OkeAPBB95DEGnKZLF2Qy6noxPaOVid22SSRSIPY6aIbp2h5yArJ 6ThFO122FDG1qhGGibwzw5QRod3YymVQUi7p2EKJNb/f2z99m3/m+TyPmd97 +d7GguC9XqGqDAbDUoXBUL5//vJKH4vOlInQ564Nt760NfU2oi0H1kvW3/mB LUOrd8fOt/X+HTh582DTsQON5O8fA9dYDsZpGzSjZU++mfyEArhqB4P3QOcp +b4u4ISAEIdS11bkOl5e05nWDbyxJ9pOYKAgz3sJ/Mq5eIEgpwMFmDotcR3v A05S2GgbSzvIegaBeVqr0o3YHWhldmtlU8kb4Ffh9gZL5Aqy3iFgN9bpxZHC NvTtCqOIY4N/Aa85Yx98JLeN7Ocd8FzUuTREsx153FOueATYuz13w4X4FrLf UWBbux+Hez4+Qb5+5xsHgj4As32H1j1d85icxxjw9kXbLHQlzSio92n8pbPj wOc73Zqmhz8i58Vgf+62h/84p1feiIQVGsxrL6FdolhBXO16cp4qwBk1O9xn W9ahsuTre9WSVIGrT2kdiedVkvNWA7YeWMqb/qAKiT2VJzYFmMEa0nXxuE3u A5pZu93c6W0Rkun3GwsrpgJzMjIXJsv/S+5rGjDLd+v5OVFJeBIH5/hPB/5j 4Q3V1cm3MH2fM4CfR2+JznSX4PrfrB5wh6Gd7CsDFkxIMX3fM4E16vWOh4nF +FFkrbn01Cxgcc2e7w/slWI6HmYD7x120eT1VWLqOm3VgVnbz6YP11diOl6g 705j1hptFeP2aeqTj9QAHrH0NagbE2E6njSBUwJupTjsLMHPqA3PARbXB4U1 pwoxHW9awOmHLqYFCQrwv5XbTYBm6Du8YKgLMR2P0FyVhZr+glI8LHgbLvac C8wT/vC8tVyI6XiFTmnWdOxMLMGHM5URow3MuW/S3615B9PxDJ3W9YSPNauw 8jb1+6HF7xrsv3WQYjre5wEzmfX6B+PKcFxg6P5/lUGPPBm7bCm4hel80AHe J4u44JldjlWohIBmaral9ptgTOcLdFF2oLNHhBQnLpoqv+irC7xP9XZSQwTG dD5Bsxdf19FaVYmpcDH/Cjj+atr9jBIZpvMNmt9/mHu5tAofG8xhjQ5Bm0Xn 3dkRXUXyUQ+4e19elKmkElOPk0JzAp52MDLkJF+ZwA+ltd/lZtzDJ0QuUV7p 0Kc8lutb+8tJPkP3WC4N//BChude0DGLC9QHTnxt/kZd5VO+Q8c1ezP6XUT4 VIIyYwyAuWyuo2p7KakH0MLApOexhyT4qx+LfvrzPbS1pWPwSptP9WI+sNeo rmCOhRTzvlfeOHRWc+v+hz13ST0xBH7V97o+8i8ZVg4bdr9AM1J2ndlm+qne QG/uDw2cuVaIqeMKNQIW2IZnxujmkXoELQ5zzXLIvIQng8Xi+DJjYEWS4wnL jZdJvYIWR2b97DSch+l2ZgJ9vaXt8opzpJ59wYcM1a7a8pGyWj2rh2aHaj4b FV4l/ZYFLK6zeempUoiUpzVjPbRnQd9ExOpCRNdDaLZHzTRDmwJElSshNKfj N4cwfwHp59CDGXyjrsBSNLnZyZKwAJjzZuYmu2oJouspdPUmRn7EaTGa/LLJ kIHmZzctT3gtJPMCtO0jn16v3TcQFY690MyBXTZ6EQWIrsemwD41uUfnq0nQ WSphoQUTfjOPLReReQS64YbP5uneItQyGQ0nRNBtj+Ifj4TdRnQ9h153cYv7 8OlrSFlNc1hfA2+x7K/a1XOWzDvQOcXZ2G8ZH1HlPvkLnsh30GdeQXQ/gB4c cwrOrctHF6mEhhZL2QFhzXlknjIDZhjjrjLpL7jDzfn/LzZBs9f2fmQ5Hyf9 BDrLu2qKMOYmptK5HJp7TT5DdkNA+g10Q7NLr2ONAIdYKyumOXD573Z71h0W k34E7eJsEnPYWoKvUAkPrTD5WZGqUkn6FTRv6/7/vEquwt3ayh1Dt6xvLX9m 9KmfLQQe7OVFnnwvxdRy/KDTpg6O1o5Uk34HPWV27J2iedWYSyUUNP87q3eM okrSD6EVTYu3C67cxVS6LVoE7MlWJG5svE/6JbTTipyjGzTkmCqHJ6AjV07c aYioJv0UmsmIlfvwZVjZDQVvofs2NLanDEhIv10MXM0TtaeYVuE9VEB+wWMv rbRHMenH0H0ifxfD9GJ885ayokDnuO8c2qYjIv0aOv8nF+F8twpMHafFEuCY jzrxneNC0s+hD876wc+nSYiV1XrWSWifVMnym335pN9Dm8nOxcq6hZgqF6PQ DxuOho0NFZJ5wAJYEKrS7iYvwcpub8+B3ixD/2s+UELmBWhWjOKPfx4RYWU1 97gHzeG82YLTfiXzxBfcyO3/Ju8Mtufwc0OsLYHZ5fNKK4aSye9ZaFsz0fPo /JsoijoQ6D9r1NXWnitG9DwCPeLKDV/jXYaocBiDLg0oa+nSrkD0vGIFLGAP jI99LELvqQVDs1KPx5mNX0f0PAOt9T4meCJegKhyVwMdtc50tZWkgvxeh+Yb 7k48sKoUUf+tsVsKnGV85mRdQgGi5yFoTrFGofuaAkS1cx40qy8k0yeuENHz EnR8UIqf0JePJqgXNCNYdcn6inj8Nzn0BgQ= "]]}, {RGBColor[0.880722, 0.611041, 0.142051], PointSize[ 0.007333333333333334], AbsoluteThickness[1.6], LineBox[CompressedData[" 1:eJxl2Hs8lOkeAPDBlpQksS655dDFuGQTTtoeyp7YpYtcVi7lurqoLc5Gy4dK jUuntYVFu7VFqFWN25hBHjPMZiZCu3Eid5UwpaNIieN936d/+s0/8/l+Pmbe 5/K7jRVBh91D5VksFluOxaLeP325p03HZIgF6FM3Rlhc8U+pRIwlwJoc7f3v HcRo08G45Ws9/gLmfDvWlny0lfz9I2ApeyxeTacdrev4cu4TvcDifays++qd 5PsGgBP9gm35zl3IeaZa2n/2KfCuoRgrrk4ved5z4GHH8hXcvD7kZ7RxjfPM KDCn11JNX9RH1jMGnKW6IU3PoQ/Z5XbVtVW8Ah6OsNFZI+kl6x0HdjbMXB3F 60Zb1utFJo+9BnbMtgk6kd9N9jMBvBT1mwer9CC3u9SK3wLv7MnfeSnhMdnv FLC51XeTQx86kLfPr60vAt8Db/Ie39q5+RE5j2ng4FV7TDWE7ShwpDPhSs4M cHa/S5tixENyXiyHT9354IuLmtWtiFezWOvGc+gN0YaB4WrN5DzlgFOk+1wX sZtQFefmYYUz8sDVF1RPJGTVkfNWADZ6YZ6leL8e1e6gTuwz4NcG4xpObpXk PqDVG0NMNr4pQ2JtmT6vZh5wQnrGSo7kJ3Jf84FZ3v6/Lok+g+fwwxJfReAH K2/Jb+KUYOY+FwAPxOyOyXAV4uY/zO6HT0Lb2dT5rZgVYea+lYAVmzVTw2pr 8cOoRhPRhYXAFdJD3xw9LMJMPCwCDpp0UskarcP0da5VBlYNyUmbbK7DTLxA F8/XatTzr8U985XnHrkYWMr21mmaFmAmnlSA/f1Kkmz3V+BBesNLgNObA8Pa U3iYiTdV4KBjl88Gcovxv6ntnoRu1LJ9xlLmYSYeoS3kVqr4cvl4kvsmonbH UuC9vG1Puqp5mIlXaM92Ffv+xAp8PIOKGDVgu3sGsqcqdzATz9DRAx2FWKUe U7epLYMummix2WIrwky8LwOe0GzW/iG+CscHhB75VxV0R8f0VTa3BDP5oA4c II68tCO3GsvRCQHNUulOkRlgzOQL9PXcAEe3SBFOXDVPctlbA9hLvvJMSyTG TD5B262+qa66oQ7T4WLyOXDs9bP30ivEmMk36DzZ8fCr/HqcPJZnODUObRhT cGdfTD3JR03g7u8Loo2EdZh+nAjaz6+zj5UuIfmqBXxf1Ph1fvpdfE7gFO2e Bv2Tm7W2ha+E5DP0E7Z5xPtnYrz0krpxfIA28KmXJq+U5T7mO/SJdg+WzEmA L5ykMkYH+KBDuL18D5/UA2h+wJkncceE+PPvyn787ztoS7Z9kJ3lx3qxHNh9 SoO7xFSEs76hbhw6t73ryIOhP0k90QV+PfqyOeq1GFPDhtUv0EpJB7L3GH2s N9D7ZaEBSl/xMH1coXrALWsjMmI1Ckg9gu4Nc/7dNuMKngsW09R1+sDKHPtz 7F1XSb2CHoz6/dTGyQLMtDMD4LGbj7uvrr9I6hl04TFdhetrCxFVrQaboa1C VQaneNdJvzUErmyyfL5DrhRRp7VgO7Rn8ehs5KZSxNRDaGc36Xxdy2JElyse dGDfH7ZhvlzSz6GH0wv1BgL4aG6zcyVhBbDXKyVPqwYhYuopdKUnqygysxbN fdlcyEBn5rZZn3zJI/MC9MTfXiPuB28hOhxHoMdlByw1I4sRU4+NgJ2k+aeX KwhRDp2w0OdnfZSSrQVkHoHm3/L6VtFDgB7PRcM5AXTDw4RHb8MqEVPPoY0v 73adzLyBqGqaZ/gPYF22rP7AUA6Zd6BTy3Oxz7pCRJd7DnTGbJGtttY1xPQD 6OHpjUH5TUXoMp3Q0EKRg19YewGZp4yB3+rhgSrRL7jPxfF/zzyhR51GPhg6 ppJ+Ap3oUf8ZL/Y2ptO5Gtr3hmSB+BaX9BvomnanEXspFwdbUBXTBJj/l9Wh rcdrST+C3uJoEHvcQoiv0QkP3W1wqjdFro70K+gc/yP/GebU46dq1I6hu7Z3 VQ/qfexnK4HHRrKizr8TYXo5PtCp88amGt82kH4Hrbgo7k7ZsgYcTicUdOHX ZhOssjrSD6Efta0O4V77E9PptmoV8FaH3sRdrfdIv4S2Xp93eudiCabL4Tno 7+1m77RENpB+Cq3KipN4FYox1Q25b6D7drb2JL0Qkn67GhhnCXqSjOrxITog oSumn5upTWHSj6G7BL5Oumnl+HYJVVGg01z3j+9RF5B+DZ39oxNvuUsNpo/T dA1w+Af1hP4ZHunn0CELt/l4tfEwVa0Xnoe2SxFa3x4tIv0eWl18MU78lIfp cjEFLWw5HTY9XkrmAVPgvFC5HhdJBaa6vc1e6G1i9HP70QoyL0Drxfb+/c8T AkxVc7e70H57X+3GZ38j8wS0dWu47MuCbGyztzA/2IINnFy9jF8zziG/Z6G9 jAVPYopuo2j6QKCHpMoKX10sR8w8Ai3nEh6x2aMK0eEwDV3jV/V4QK0GMfOK GXCJw4uZ6Q9l6B29YGj1lNR445mbiJlnoGemYoNmE7iILndS6JCtRpvMhDXk 9zp0pu7BxKMb+Ij+b42VOfDP+tnnm04WI2YegvYtX1zqurkY0e08C1p/NDjD K74UMfMSdHRgkg/PuxDN0i9ohyD5NdtrEvD/AapSB2Q= "]]}, {RGBColor[0.560181, 0.691569, 0.194885], PointSize[ 0.007333333333333334], AbsoluteThickness[1.6], LineBox[CompressedData[" 1:eJx12H9o1VUYx/EzTVPzR/7InKZNc02niSKJfxxLUUQJSUYhkog/GPmHP4p+ YEV4k1DEGBdFGhJ6wh8MRXcaY93hg1zWxCmiiBRWSqWiiBorM5uN2fec5+yf 2Pv+M16w3fv9nu/zPJ/nbvyaTVXVvYwxU0qMCT///6rKd36051Tzqz37Has+ AzbmlfWfjpn+xiW0/v6P6Jk/zcn+4ldwXvT9roOdLOqSs9e+uAnu/rzb6BUT 7ORFXXfBTvR62sHGzD5w9bsfvv0Drdd7Hz3/5bHv7Wj/C+ys3s/fYGOWnA5X /A/Yp+fVAXZ22fKvLv6++l+0nkcnevWdn3Nf7+0Ce6vnZeb2bGebTg4adeQ2 2Vs9zxL0ie3HNvXe1gtcsHrevdHFpeHEngA7q8+DnLenSu+NazrZh5zqsS86 w4dD3noSnBN9nv3Au+TC0ann1z0k7xN93v3BdfL9++fKW3YPADeI1sNTYJH4 OKcPBBdF64XcKr/0HZh95CBwm2g9DQa3yo14w0PAbalfn0Z/EG53K1vrkXxO HvoHG4pLh6K1Xslt8vGeUDHDwK2i9Uxuk/A0S++xtd6Hg4uyZWX1uwtPsLUf RoC9lMSGIBdE+4X9eUWfM/uXPYPWfiJ7ieVSPhLcnQfsHe0Hyzrus7UfnwV7 iR/XQnapX0eha5oXbK7Ks7WfybUydN+IiVtWloJzqd/Jn9ndW0PHjAbnrc4D 9si3Gz+5/IhsjM6LMeja18ITJzur8+Q5cKMNy8aML8k+5RU7Hlf1WLTOI3ZW LJU7Z44DF1Iekr3VOHserfOM7GyYVjcusHXelaHDafV7HZ3mITuOqyayMTov yc5mN5uNhPFgb3Wekgs2e7OsZMjd+wK5aGM53iG3WJ3HE9B7Y8OSJe0j5IK9 klVDTTO5yeo8B+e+sWGaHix7oWeb+rTvkOtsHPfbyc5qHpDzdn9saHJ2ojEv JoLz8tvieX/eepPsUp6QvcR2FrbmDXvttDAxy8Eu5RHZmEOx4dmaV+ybw8Id szXPXkTHy1nO1rwjO7suNhTZGM1DspfYbhUVaM1LdhyHNWSX8pTsJaShf0Au pLydBPayMRYkW/OYXd8QJgq5MeU1Ox5n5WS05jk41yBhWg/YBTbHU96TvcRx 0UFuEN0HKsF1EtJ+1iqyS98nyMaEab7kNFv3CbKTWavqDq+dNgWt+wbZy+Z4 IGzdR8hOYjl0snVfmQo25lG8YLJL+wzZ2zjuzrJ132HH/9bMeAmt+xDZ2Rjn tWRjdF9iP44vsrP/AXKxQuU= "]]}}, {{}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], Method->{ "OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ Identity[ Part[#, 1]], Identity[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ Identity[ Part[#, 1]], Identity[ Part[#, 2]]}& )}}, PlotRange->{{0, 14.999999999999998`}, {-0.3, 0.6}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, {0, 0}}, Ticks->{Automatic, Automatic}]} }, AutoDelete->False, GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, GridBoxSpacings->{"Columns" -> {{4}}}], "Grid"]], "Output", CellChangeTimes->{ 3.762197212818427*^9, {3.7621972625466633`*^9, 3.762197278810548*^9}, 3.762197504955035*^9, 3.762197642392365*^9, {3.762197673280938*^9, 3.762197698088641*^9}, {3.76219775171242*^9, 3.762197763405024*^9}, 3.7621988649554453`*^9, {3.762198916518241*^9, 3.7621989586372337`*^9}, 3.762199014540296*^9, 3.762199098210876*^9, 3.762199195864786*^9, { 3.762199249587709*^9, 3.762199268891148*^9}, {3.76219936662719*^9, 3.7621993862570677`*^9}, {3.762199510334671*^9, 3.762199531081738*^9}, 3.762199569559883*^9, 3.762199618541822*^9, {3.762199703792391*^9, 3.7621997463205633`*^9}, {3.762199841900545*^9, 3.7621998902185907`*^9}, { 3.76220333538874*^9, 3.7622033770009127`*^9}, {3.762203409007957*^9, 3.762203481928541*^9}, {3.762204792468906*^9, 3.762204854855832*^9}, { 3.762213304799882*^9, 3.76221333611662*^9}, 3.762213484238714*^9, { 3.762213727788582*^9, 3.762213783069066*^9}, {3.7622138202410316`*^9, 3.762213832858261*^9}, 3.762213876762684*^9, {3.762215397244705*^9, 3.7622155190608807`*^9}, {3.762215582719809*^9, 3.762215610745611*^9}, { 3.762215657549638*^9, 3.762215681750853*^9}, 3.762215712717291*^9, { 3.762215771857085*^9, 3.762215831624331*^9}, 3.762218787075348*^9, { 3.762219488693756*^9, 3.7622195097558603`*^9}, 3.762219741991323*^9, 3.762219844656371*^9, 3.762222165375514*^9, {3.762222198691024*^9, 3.762222226160482*^9}, {3.762222320567809*^9, 3.7622223478128843`*^9}, 3.762222411040936*^9, 3.762222567382531*^9, {3.7622227068749447`*^9, 3.762222727234448*^9}, {3.762222786238966*^9, 3.7622228660510397`*^9}, { 3.7622229225041847`*^9, 3.762222943804055*^9}, {3.762223000553392*^9, 3.762223085792268*^9}, {3.762223184780827*^9, 3.7622232448946657`*^9}, 3.762223301126424*^9, 3.762223363445265*^9, {3.7622234635902443`*^9, 3.762223484416829*^9}, 3.762223524452116*^9, 3.762224562735259*^9, { 3.7622246385605793`*^9, 3.762224655126856*^9}, {3.762224768256753*^9, 3.7622247853046837`*^9}, 3.7622252373533297`*^9, 3.762225314461473*^9, { 3.762227616906666*^9, 3.762227646511689*^9}, 3.762227679720969*^9, 3.762227796179632*^9, {3.7622278413286753`*^9, 3.7622278491960907`*^9}, 3.762227906024839*^9, {3.762227985114732*^9, 3.7622280058835382`*^9}, { 3.762228069216201*^9, 3.762228080787154*^9}, 3.76222818497937*^9, 3.762228244020208*^9, 3.762228318181788*^9, 3.762228377855871*^9, 3.762228848186655*^9, 3.793586875249881*^9, 3.7938418967835712`*^9, 3.7938420048078527`*^9, 3.7938420934865932`*^9, {3.793842187309124*^9, 3.793842220047513*^9}, 3.793842273792713*^9, 3.793842416703652*^9, 3.793908756447013*^9, 3.823229882079051*^9, 3.875753353477861*^9, 3.8757534235089693`*^9}, CellLabel-> "Out[3100]=",ExpressionUUID->"d2acd1ad-043f-41ee-897a-aec6ba6bafc1"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"StandardDeviation", "/@", RowBox[{"{", RowBox[{ RowBox[{"Take", "[", RowBox[{"xs2", ",", RowBox[{"-", "50"}]}], "]"}], ",", RowBox[{"Take", "[", RowBox[{"xe2", ",", RowBox[{"-", "50"}]}], "]"}]}], "}"}], " ", RowBox[{"(*", " ", RowBox[{"sd", " ", "of", " ", "last", " ", "50", " ", "elements"}], " ", "*)"}]}]], "Input", CellChangeTimes->{{3.7622237347202578`*^9, 3.762223743650178*^9}, 3.762223775375497*^9, {3.762223992192687*^9, 3.762224046754156*^9}, 3.762224627849821*^9, {3.7622276883692636`*^9, 3.762227688508553*^9}, { 3.762228099812374*^9, 3.7622281654473143`*^9}, {3.762228237534542*^9, 3.762228280905727*^9}, {3.793842257950529*^9, 3.7938422688784037`*^9}}, CellLabel-> "In[3101]:=",ExpressionUUID->"ef7cc903-78d3-48d2-b31b-cb31a2386605"], Cell[BoxData[ RowBox[{"{", RowBox[{"0.10717310039386918`", ",", "0.04157494717510941`"}], "}"}]], "Output", CellChangeTimes->{{3.762223738646949*^9, 3.762223744235549*^9}, 3.762223775953066*^9, {3.762223996761142*^9, 3.762224048096767*^9}, 3.7622245657693663`*^9, {3.762224628609942*^9, 3.762224657546302*^9}, { 3.762224770648658*^9, 3.762224787834373*^9}, 3.762225237496768*^9, 3.76222531458786*^9, {3.7622276169749517`*^9, 3.76222764670891*^9}, { 3.76222767986157*^9, 3.762227689134707*^9}, 3.762227906177816*^9, { 3.762227985214151*^9, 3.762228006269342*^9}, {3.762228069376869*^9, 3.762228114592658*^9}, {3.7622281661494627`*^9, 3.762228185112812*^9}, { 3.762228244155208*^9, 3.762228281406764*^9}, 3.762228318377733*^9, 3.76222837798374*^9, 3.76222884831153*^9, 3.7935868752632313`*^9, 3.793841896797511*^9, 3.793842004841667*^9, 3.7938420935130653`*^9, 3.7938422200740013`*^9, {3.793842261491123*^9, 3.793842273806065*^9}, 3.7938424167186947`*^9, 3.793908756461687*^9, 3.823229882093197*^9, 3.87575336011583*^9, 3.875753423548922*^9}, CellLabel-> "Out[3101]=",ExpressionUUID->"163e09ad-cf9b-4974-9241-312cb66fff5f"] }, Open ]], Cell["\<\ Compare the above closed-loop system with the one-step routine LQGRegulator\ \>", "Text", CellChangeTimes->{{3.513266788083189*^9, 3.5132667973147306`*^9}, { 3.513272093815999*^9, 3.513272098697585*^9}, {3.5172389588536386`*^9, 3.5172389629386263`*^9}, {3.517239692171343*^9, 3.5172397065137873`*^9}, 3.5172399022732773`*^9, {3.51723998911631*^9, 3.517239998600172*^9}, { 3.762184806607835*^9, 3.762184813214233*^9}, {3.762197448350705*^9, 3.762197451565967*^9}, {3.7621995236785192`*^9, 3.762199526597042*^9}, { 3.875744646747241*^9, 3.875744668172064*^9}}, FontSize->16, Background->GrayLevel[ 0.85],ExpressionUUID->"997155b3-7240-42e1-8ec5-f64619ae0c91"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"sspec", "=", RowBox[{"<|", RowBox[{ RowBox[{"\"\\"", "\[Rule]", "G1dss"}], ",", RowBox[{"\"\\"", "\[Rule]", " ", "1"}], ",", RowBox[{"\"\\"", "\[Rule]", "1"}], ",", RowBox[{"SystemsModelLabels", "->", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}], ",", RowBox[{"{", "\"\\"", "}"}], ",", RowBox[{"{", RowBox[{ "\"\\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(\[FilledSmallCircle]\)]\)\>\""}], "}"}]}], "}"}]}]}], "|>"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"closedLoop1", "=", RowBox[{"LQGRegulator", "[", RowBox[{"sspec", ",", RowBox[{"{", RowBox[{ RowBox[{"(", "\[NoBreak]", GridBox[{ { SuperscriptBox["\[Nu]", "2"]} }], "\[NoBreak]", ")"}], ",", RowBox[{"(", "\[NoBreak]", GridBox[{ { SuperscriptBox["\[Xi]", "2"]} }], "\[NoBreak]", ")"}]}], "}"}], ",", RowBox[{"{", RowBox[{"q", ",", "r"}], "}"}], ",", "\"\\"", ",", RowBox[{"Method", "\[Rule]", "\"\\""}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{"closedLoop1", ",", "closedLoop2"}], "}"}]}], "Input", CellLabel-> "In[3102]:=",ExpressionUUID->"3bcffffd-1d9f-43b9-b5dd-473337bae631"], Cell[BoxData[ RowBox[{"{", RowBox[{ TemplateBox[{"0.9950041652780258`", "0.09983341664682817`", RowBox[{"-", "0.00869952673989844`"}], RowBox[{"-", "0.016289313538040986`"}], "0.004995834721974235`", "0.004995834721974235`", "0.`", RowBox[{"-", "0.09983341664682817`"}], "0.9950041652780258`", RowBox[{"-", "0.17384551851452998`"}], RowBox[{"-", "0.3255147369429831`"}], "0.09983341664682815`", "0.09983341664682815`", "0.`", "0.09046726479597311`", "0.`", "0.8958373737421542`", "0.08354410310878718`", "0.004995834721974235`", "0.`", "0.09046726479597311`", "0.03063492588555703`", "0.`", RowBox[{"-", "0.30431386104691516`"}], "0.6694894283350428`", "0.09983341664682815`", "0.`", "0.03063492588555703`", "1.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.1`", "\"w\"", "\"y\"", { StateSpaceModel, { True, True, { False, False}}, {$CellContext`stname1, $CellContext`stname2, \ $CellContext`stname3, $CellContext`stname4}, {{ Control`CommonDump`$DUMMY$, "w", Control`CommonDump`$DUMMY$}, {"y"}, { Control`CommonDump`$DUMMY$, Control`CommonDump`$DUMMY$, Control`CommonDump`$DUMMY$, Control`CommonDump`$DUMMY$}}, Identity, Automatic, 3, 1, 4, 1, {{1, 2}, {2, 1}}, {True, True}, False, Automatic, {{None, "w"}, {"y"}}, {",", "Automatic", ",", RowBox[{"{", RowBox[{"Automatic", ",", "Automatic", ",", "Automatic"}], "}"}], ",", "Automatic", ",", "z"}}, {}}, "StateSpaceModelFull", BaseStyle->{NumberMarks -> False}, DisplayFunction->(TagBox[ SubsuperscriptBox[ RowBox[{"(", GridBox[{{ "", "", "\[SpanFromLeft]", "\[SpanFromLeft]", "\[SpanFromLeft]", "", #37, ""}, {"", #, #2, #3, #4, #5, #6, #7}, { "", #8, #9, #10, #11, #12, #13, #14}, { "", #15, #16, #17, #18, #19, #20, #21}, { "", #22, #23, #24, #25, #26, #27, #28}, {#38, #29, #30, #31, #32, \ #33, #34, #35}}, AllowedDimensions -> {8, 6}, GridBoxDividers -> {"Rows" -> {False, Directive[ GrayLevel[0], GrayLevel[0.6]], False, False, False, Directive[ GrayLevel[0], GrayLevel[0.6]], {False}, False}, "Columns" -> {False, Directive[ GrayLevel[0], GrayLevel[0.6]], False, False, False, Directive[ GrayLevel[0], GrayLevel[0.6]], {False}, False}}, DefaultBaseStyle -> "SystemsModelGrid11"], ")"}], #36, ButtonBox[ GraphicsBox[{ RectangleBox[{0, 0}, DefaultBaseStyle -> "SystemsModelSuperscript"], InsetBox[ StyleBox[ "\[ScriptCapitalS]", DefaultBaseStyle -> "SystemsModelSuperscript"], {0.5, 0.5}]}, ImageSize -> 20], ButtonFunction :> Control`Typesetting`SystemsModelTypesetToggle[ ButtonNotebook[], FE`CMObj], Evaluator -> Automatic, Appearance -> None], MultilineFunction -> None], EventHandlerTag[{"MouseEntered" :> {FE`CMObj = MathLink`CallFrontEnd[ FrontEnd`Value[ FEPrivate`Self[]]]}, PassEventsDown -> True, PassEventsUp -> True}]]& ), InterpretationFunction->(RowBox[{"StateSpaceModel", "[", RowBox[{"{", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{1, 4}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{8, 11}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{15, 18}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{22, 25}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{5, 7}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{12, 14}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{19, 21}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{26, 28}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{29, 32}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{33, 35}, ","], "}"}], "}"}], "}"}], ",", "Automatic", ",", RowBox[{"{", RowBox[{"Automatic", ",", "Automatic", ",", "Automatic"}], "}"}], ",", "Automatic", ",", "z", ",", "SamplingPeriod", "->", RowBox[{#36}], ",", "SystemsModelLabels", "->", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"None", ",", #37}], "}"}], ",", RowBox[{"{", #38, "}"}]}], "}"}], "]"}]& ), Tooltip->StateSpaceModel], ",", TemplateBox[{"0.9950041652780258`", "0.09983341664682817`", RowBox[{"-", "0.00869952673989844`"}], RowBox[{"-", "0.016289313538040986`"}], "0.004995834721974235`", "0", RowBox[{"-", "0.09983341664682817`"}], "0.9950041652780258`", RowBox[{"-", "0.17384551851452998`"}], RowBox[{"-", "0.3255147369429831`"}], "0.09983341664682815`", "0", "0.09046726479597311`", "0.`", "0.8958373737421543`", "0.08354410310878718`", "0", "0.09046726479597311`", "0.03063492588555703`", "0.`", RowBox[{"-", "0.30431386104691516`"}], "0.6694894283350428`", "0", "0.03063492588555703`", "1", "0", "0", "0", "0", "0", "0.1`", { StateSpaceModel, { False, False, { False, False}}, {$CellContext`stname1, $CellContext`stname2, \ $CellContext`stname3, $CellContext`stname4}, {{ Control`CommonDump`$DUMMY$, Control`CommonDump`$DUMMY$}, { Control`CommonDump`$DUMMY$}, { Control`CommonDump`$DUMMY$, Control`CommonDump`$DUMMY$, Control`CommonDump`$DUMMY$, Control`CommonDump`$DUMMY$}}, Identity, Automatic, 2, 1, 4, 1, {}, {False, False}, False, Automatic, None, Automatic}, {}}, "StateSpaceModelFull", BaseStyle->{NumberMarks -> False}, DisplayFunction->(TagBox[ SubsuperscriptBox[ RowBox[{"(", GridBox[{{#, #2, #3, #4, #5, #6}, {#7, #8, #9, #10, #11, #12}, {#13, \ #14, #15, #16, #17, #18}, {#19, #20, #21, #22, #23, #24}, {#25, #26, #27, \ #28, #29, #30}}, AllowedDimensions -> {6, 5}, GridBoxDividers -> {"Rows" -> {False, False, False, False, Directive[ GrayLevel[0], GrayLevel[0.6]], {False}, False}, "Columns" -> {False, False, False, False, Directive[ GrayLevel[0], GrayLevel[0.6]], {False}, False}}, DefaultBaseStyle -> "SystemsModelGrid00"], ")"}], #31, ButtonBox[ GraphicsBox[{ RectangleBox[{0, 0}, DefaultBaseStyle -> "SystemsModelSuperscript"], InsetBox[ StyleBox[ "\[ScriptCapitalS]", DefaultBaseStyle -> "SystemsModelSuperscript"], {0.5, 0.5}]}, ImageSize -> 20], ButtonFunction :> Control`Typesetting`SystemsModelTypesetToggle[ ButtonNotebook[], FE`CMObj], Evaluator -> Automatic, Appearance -> None], MultilineFunction -> None], EventHandlerTag[{"MouseEntered" :> {FE`CMObj = MathLink`CallFrontEnd[ FrontEnd`Value[ FEPrivate`Self[]]]}, PassEventsDown -> True, PassEventsUp -> True}]]& ), InterpretationFunction->(RowBox[{"StateSpaceModel", "[", RowBox[{"{", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{1, 4}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{7, 10}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{13, 16}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{19, 22}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{5, 6}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{11, 12}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{17, 18}, ","], "}"}], ",", RowBox[{"{", TemplateSlotSequence[{23, 24}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{25, 28}, ","], "}"}], "}"}], ",", RowBox[{"{", RowBox[{"{", TemplateSlotSequence[{29, 30}, ","], "}"}], "}"}], "}"}], ",", "SamplingPeriod", "->", RowBox[{#31}], ",", "SystemsModelLabels", "->", "None", "]"}]& ), Tooltip->StateSpaceModel]}], "}"}]], "Output", CellChangeTimes->{{3.875753382114778*^9, 3.875753424221181*^9}}, CellLabel-> "Out[3104]=",ExpressionUUID->"b9db6498-5ddf-4fd7-948a-c0560cd077cb"] }, Open ]], Cell["Export data", "Text", CellChangeTimes->{{3.513266788083189*^9, 3.5132667973147306`*^9}, { 3.513272093815999*^9, 3.513272098697585*^9}, {3.5172389588536386`*^9, 3.5172389629386263`*^9}, {3.517239692171343*^9, 3.5172397065137873`*^9}, 3.5172399022732773`*^9, {3.51723998911631*^9, 3.517239998600172*^9}, { 3.762184806607835*^9, 3.762184813214233*^9}, {3.762197448350705*^9, 3.762197451565967*^9}, {3.762225124743793*^9, 3.762225160180443*^9}, 3.762227896642788*^9, {3.76222863689554*^9, 3.762228638861969*^9}}, FontSize->16, Background->GrayLevel[ 0.85],ExpressionUUID->"b2802ce5-6413-4518-9174-a0c45efebe69"], Cell[BoxData[ RowBox[{"(*", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"SetDirectory", "[", RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"Export", "[", RowBox[{"\"\\"", ",", " ", RowBox[{ RowBox[{"{", RowBox[{ "ym", ",", "xe", ",", "xs", ",", "ym2", ",", "xe2", ",", "xs2"}], "}"}], "\[Transpose]"}]}], "]"}], ";"}], "\[IndentingNewLine]", "*)"}]], "Input", CellChangeTimes->{{3.762228654906502*^9, 3.762228662400475*^9}, { 3.762228758646339*^9, 3.762228799971019*^9}, {3.762228839185369*^9, 3.762228880638526*^9}, {3.762231986597589*^9, 3.762231990761389*^9}, { 3.793586853678639*^9, 3.7935868612305737`*^9}}, CellLabel-> "In[3105]:=",ExpressionUUID->"304bddec-c63d-4a1f-ac3f-d4788b0e1f5e"] }, Open ]] }, WindowSize->{1393, 652}, WindowMargins->{{Automatic, 39}, {Automatic, 0}}, FrontEndVersion->"13.1 for Mac OS X x86 (64-bit) (June 16, 2022)", StyleDefinitions->"Default.nb", ExpressionUUID->"f72d2bcd-59f3-435a-8473-741ebf035891" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[422, 15, 306, 5, 67, "Section",ExpressionUUID->"0f17a414-3f5c-4602-9c7d-b49762641eab"], Cell[731, 22, 824, 17, 124, "Text",ExpressionUUID->"98bfe169-696b-4d7a-beae-415ad2a61460"], Cell[1558, 41, 561, 9, 54, "Text",ExpressionUUID->"428ad61b-4d61-4385-b9cc-a671d6e59502"], Cell[CellGroupData[{ Cell[2144, 54, 5316, 131, 272, "Input",ExpressionUUID->"87764a34-4e04-4711-b148-a77543ad8919"], Cell[7463, 187, 8783, 213, 113, "Output",ExpressionUUID->"075c041e-c778-4cdf-93de-822e9a9a39de"] }, Open ]], Cell[16261, 403, 594, 11, 54, "Text",ExpressionUUID->"6ac269f8-a39a-4ca9-a823-bda5dd8fc0cc"], Cell[16858, 416, 3701, 97, 201, "Input",ExpressionUUID->"eb36d036-fcff-46ab-86b1-8f0c73a6b48b"], Cell[20562, 515, 759, 13, 52, "Text",ExpressionUUID->"e7b0b0ab-e691-482d-8404-e12dd17caddb"], Cell[CellGroupData[{ Cell[21346, 532, 956, 26, 45, "Input",ExpressionUUID->"c860eed2-6467-43e3-9dea-e94056cf430b"], Cell[22305, 560, 1710, 27, 48, "Output",ExpressionUUID->"3b3b61db-adae-45a3-8946-eb775a5d318e"] }, Open ]], Cell[CellGroupData[{ Cell[24052, 592, 806, 22, 33, "Input",ExpressionUUID->"2861f37c-7d70-4363-b519-8f535afe8bb7"], Cell[24861, 616, 2113, 32, 48, "Output",ExpressionUUID->"7d2f1ec4-1b04-4d83-9f34-f3d58b53dc70"] }, Open ]], Cell[CellGroupData[{ Cell[27011, 653, 944, 22, 52, "Input",ExpressionUUID->"fa5d7659-7981-4357-ae51-8483a7c99902"], Cell[27958, 677, 2613, 43, 48, "Output",ExpressionUUID->"d9035ca0-8d6f-4ef4-8d9c-4f864f178299"] }, Open ]], Cell[30586, 723, 3597, 97, 175, "Input",ExpressionUUID->"69771da2-719e-4de9-8099-30433d7af0b1"], Cell[34186, 822, 563, 9, 52, "Text",ExpressionUUID->"9bcd96b6-baa5-4f4d-a7e0-c4737b8c9aba"], Cell[CellGroupData[{ Cell[34774, 835, 498, 13, 30, "Input",ExpressionUUID->"13b557ae-8017-4685-b51a-0df23aa5523c"], Cell[35275, 850, 30287, 540, 145, "Output",ExpressionUUID->"d2acd1ad-043f-41ee-897a-aec6ba6bafc1"] }, Open ]], Cell[CellGroupData[{ Cell[65599, 1395, 846, 19, 30, "Input",ExpressionUUID->"ef7cc903-78d3-48d2-b31b-cb31a2386605"], Cell[66448, 1416, 1180, 19, 48, "Output",ExpressionUUID->"163e09ad-cf9b-4974-9241-312cb66fff5f"] }, Open ]], Cell[67643, 1438, 693, 12, 52, "Text",ExpressionUUID->"997155b3-7240-42e1-8ec5-f64619ae0c91"], Cell[CellGroupData[{ Cell[68361, 1454, 1492, 40, 85, "Input",ExpressionUUID->"3bcffffd-1d9f-43b9-b5dd-473337bae631"], Cell[69856, 1496, 8825, 188, 140, "Output",ExpressionUUID->"b9db6498-5ddf-4fd7-948a-c0560cd077cb"] }, Open ]], Cell[78696, 1687, 640, 10, 52, "Text",ExpressionUUID->"b2802ce5-6413-4518-9174-a0c45efebe69"], Cell[79339, 1699, 836, 19, 94, "Input",ExpressionUUID->"304bddec-c63d-4a1f-ac3f-d4788b0e1f5e"] }, Open ]] } ] *)