(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 10.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 8963, 270] NotebookOptionsPosition[ 8062, 233] NotebookOutlinePosition[ 8406, 248] CellTagsIndexPosition[ 8363, 245] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"f", "[", RowBox[{"t_", ",", "a_", ",", "\[Omega]_", ",", "\[Gamma]_"}], "]"}], ":=", RowBox[{"a", " ", RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "\[Gamma]"}], " ", "t"}], "]"}], RowBox[{"Cos", "[", RowBox[{"\[Omega]", " ", "t"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.6709450914214063`*^9, 3.670945152907493*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"NotebookDirectory", "[", "]"}]], "Input"], Cell[BoxData["\<\"C:\\\\Users\\\\vanin\\\\Desktop\\\\\"\>"], "Output", CellChangeTimes->{3.6709451694175167`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"t", "=", RowBox[{"Range", "[", RowBox[{"0", ",", "5", ",", "0.25"}], "]"}]}]], "Input", CellChangeTimes->{{3.6709453008157053`*^9, 3.6709453372347574`*^9}, { 3.6709454153538675`*^9, 3.670945457252927*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "0.`", ",", "0.25`", ",", "0.5`", ",", "0.75`", ",", "1.`", ",", "1.25`", ",", "1.5`", ",", "1.75`", ",", "2.`", ",", "2.25`", ",", "2.5`", ",", "2.75`", ",", "3.`", ",", "3.25`", ",", "3.5`", ",", "3.75`", ",", "4.`", ",", "4.25`", ",", "4.5`", ",", "4.75`", ",", "5.`"}], "}"}]], "Output", CellChangeTimes->{ 3.670945338174759*^9, {3.6709454211238756`*^9, 3.6709454577629275`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"y0", "=", RowBox[{"Table", "[", " ", RowBox[{ RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", "10.", ",", "1.5", ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "t", "]"}]}], " ", "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.6709453431747656`*^9, 3.670945410563861*^9}, { 3.6709454638929367`*^9, 3.6709454682729425`*^9}, {3.670945506541997*^9, 3.6709455084919996`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "10.`", ",", "9.075333068689924`", ",", "6.9600398165242305`", ",", "4.000212049177964`", ",", "0.6400566691609326`", ",", RowBox[{"-", "2.643373914326029`"}], ",", RowBox[{"-", "5.406740476572865`"}], ",", RowBox[{"-", "7.299139080908658`"}], ",", RowBox[{"-", "8.10537302283234`"}], ",", RowBox[{"-", "7.768636098352096`"}], ",", RowBox[{"-", "6.390522700526187`"}], ",", RowBox[{"-", "4.209469153068849`"}], ",", RowBox[{"-", "1.561613690614905`"}], ",", "1.1697380465359173`", ",", "3.6086053673024843`", ",", "5.4371698755234865`", ",", "6.436213907495263`", ",", "6.510161010763925`", ",", "5.694059849506122`", ",", "4.142441217788538`", ",", "2.1024494800617655`"}], "}"}]], "Output", CellChangeTimes->{{3.6709453904638324`*^9, 3.670945468612943*^9}, 3.6709455087019997`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"s", "=", " ", RowBox[{"Round", "[", " ", RowBox[{ RowBox[{"Sqrt", "[", " ", RowBox[{ RowBox[{"0.003", RowBox[{"y0", "^", "2"}]}], "+", RowBox[{"0.2", "^", "2"}]}], "]"}], ",", "0.01"}], "]"}]}]], "Input", CellChangeTimes->{{3.670945197366556*^9, 3.6709452102265744`*^9}, { 3.6709452474846296`*^9, 3.6709452970257*^9}, {3.670945489122972*^9, 3.6709455684710846`*^9}, {3.67094561441015*^9, 3.6709456433701906`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "0.58`", ",", "0.54`", ",", "0.43`", ",", "0.3`", ",", "0.2`", ",", "0.25`", ",", "0.36`", ",", "0.45`", ",", "0.49`", ",", "0.47000000000000003`", ",", "0.4`", ",", "0.31`", ",", "0.22`", ",", "0.21`", ",", "0.28`", ",", "0.36`", ",", "0.41000000000000003`", ",", "0.41000000000000003`", ",", "0.37`", ",", "0.3`", ",", "0.23`"}], "}"}]], "Output", CellChangeTimes->{{3.670945481972962*^9, 3.6709454984519854`*^9}, { 3.6709455647910795`*^9, 3.670945569491086*^9}, {3.6709456189401565`*^9, 3.670945644050192*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", " ", RowBox[{"Round", "[", " ", RowBox[{ RowBox[{"RandomVariate", "[", " ", RowBox[{"NormalDistribution", "[", " ", RowBox[{ RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", "10.", ",", "1.5", ",", "0.1"}], "]"}], ",", RowBox[{"s", "[", RowBox[{"[", "i", "]"}], "]"}]}], " ", "]"}], " ", "]"}], ",", "0.01"}], "]"}], ",", RowBox[{"s", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "t", "]"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.6709455748810935`*^9, 3.6709456019711323`*^9}, { 3.67094565650021*^9, 3.670945702859275*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0.`", ",", "11.06`", ",", "0.58`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.25`", ",", "9.85`", ",", "0.54`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.5`", ",", "7.5200000000000005`", ",", "0.43`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.75`", ",", "4.3`", ",", "0.3`"}], "}"}], ",", RowBox[{"{", RowBox[{"1.`", ",", "0.32`", ",", "0.2`"}], "}"}], ",", RowBox[{"{", RowBox[{"1.25`", ",", RowBox[{"-", "2.86`"}], ",", "0.25`"}], "}"}], ",", RowBox[{"{", RowBox[{"1.5`", ",", RowBox[{"-", "5.41`"}], ",", "0.36`"}], "}"}], ",", RowBox[{"{", RowBox[{"1.75`", ",", RowBox[{"-", "7.09`"}], ",", "0.45`"}], "}"}], ",", RowBox[{"{", RowBox[{"2.`", ",", RowBox[{"-", "8.34`"}], ",", "0.49`"}], "}"}], ",", RowBox[{"{", RowBox[{"2.25`", ",", RowBox[{"-", "8.52`"}], ",", "0.47000000000000003`"}], "}"}], ",", RowBox[{"{", RowBox[{"2.5`", ",", RowBox[{"-", "6.09`"}], ",", "0.4`"}], "}"}], ",", RowBox[{"{", RowBox[{"2.75`", ",", RowBox[{"-", "4.07`"}], ",", "0.31`"}], "}"}], ",", RowBox[{"{", RowBox[{"3.`", ",", RowBox[{"-", "1.1400000000000001`"}], ",", "0.22`"}], "}"}], ",", RowBox[{"{", RowBox[{"3.25`", ",", "1.1300000000000001`", ",", "0.21`"}], "}"}], ",", RowBox[{"{", RowBox[{"3.5`", ",", "3.91`", ",", "0.28`"}], "}"}], ",", RowBox[{"{", RowBox[{"3.75`", ",", "5.54`", ",", "0.36`"}], "}"}], ",", RowBox[{"{", RowBox[{"4.`", ",", "7.05`", ",", "0.41000000000000003`"}], "}"}], ",", RowBox[{"{", RowBox[{"4.25`", ",", "6.36`", ",", "0.41000000000000003`"}], "}"}], ",", RowBox[{"{", RowBox[{"4.5`", ",", "5.15`", ",", "0.37`"}], "}"}], ",", RowBox[{"{", RowBox[{"4.75`", ",", "4.48`", ",", "0.3`"}], "}"}], ",", RowBox[{"{", RowBox[{"5.`", ",", "2.`", ",", "0.23`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.670945602721133*^9, {3.6709456502902007`*^9, 3.670945703569276*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Export", "[", RowBox[{ RowBox[{"FileNameJoin", "[", RowBox[{"{", RowBox[{ RowBox[{"NotebookDirectory", "[", "]"}], ",", "\"\\""}], "}"}], "]"}], ",", "d"}], "]"}]], "Input", CellChangeTimes->{{3.6709456953292646`*^9, 3.6709457186792974`*^9}, { 3.670945799667414*^9, 3.6709458395574703`*^9}}], Cell[BoxData["\<\"C:\\\\Users\\\\vanin\\\\Desktop\\\\OHamortecido.dat\"\>"], \ "Output", CellChangeTimes->{3.6709458402274714`*^9}] }, Open ]] }, WindowSize->{767, 738}, WindowMargins->{{Automatic, 52}, {4, Automatic}}, FrontEndVersion->"10.0 for Microsoft Windows (64-bit) (December 4, 2014)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 383, 10, 31, "Input"], Cell[CellGroupData[{ Cell[966, 34, 65, 1, 31, "Input"], Cell[1034, 37, 114, 1, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[1185, 43, 239, 5, 31, "Input"], Cell[1427, 50, 440, 8, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[1904, 63, 524, 14, 31, "Input"], Cell[2431, 79, 878, 17, 72, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3346, 101, 481, 11, 31, "Input"], Cell[3830, 114, 577, 10, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4444, 129, 981, 27, 112, "Input"], Cell[5428, 158, 2079, 54, 132, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7544, 217, 367, 9, 31, "Input"], Cell[7914, 228, 132, 2, 88, "Output"] }, Open ]] } ] *) (* End of internal cache information *)