source: BAORadio/AmasArchive/Amas 21cm/fonctionCalib.nb@ 689

Last change on this file since 689 was 672, checked in by campagne, 14 years ago
File size: 9.0 KB
Line 
1(* Content-type: application/vnd.wolfram.mathematica *)
2
3(*** Wolfram Notebook File ***)
4(* http://www.wolfram.com/nb *)
5
6(* CreatedBy='Mathematica 8.0' *)
7
8(*CacheID: 234*)
9(* Internal cache information:
10NotebookFileLineBreakTest
11NotebookFileLineBreakTest
12NotebookDataPosition[ 157, 7]
13NotebookDataLength[ 8731, 307]
14NotebookOptionsPosition[ 7215, 248]
15NotebookOutlinePosition[ 7559, 263]
16CellTagsIndexPosition[ 7516, 260]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[
22 RowBox[{
23 RowBox[{"freq", "=",
24 RowBox[{"{",
25 RowBox[{"1280", ",", "1380", ",", "1410"}], "}"}]}], ";"}]], "Input",
26 CellChangeTimes->{{3.510634896992296*^9, 3.5106349529288893`*^9}}],
27
28Cell[BoxData[""], "Input",
29 CellChangeTimes->{{3.5106349560001965`*^9, 3.510635002484844*^9}}],
30
31Cell[BoxData[
32 RowBox[{"freq", "=."}]], "Input",
33 CellChangeTimes->{{3.510635055560151*^9, 3.5106350577223673`*^9}}],
34
35Cell[CellGroupData[{
36
37Cell[BoxData[
38 RowBox[{"freq", "=",
39 RowBox[{"{",
40 RowBox[{"f1", ",", "f2", ",", "f3"}], "}"}]}]], "Input",
41 CellChangeTimes->{{3.510635061135709*^9, 3.5106350783984346`*^9}}],
42
43Cell[BoxData[
44 RowBox[{"{",
45 RowBox[{"f1", ",", "f2", ",", "f3"}], "}"}]], "Output",
46 CellChangeTimes->{3.510635079997595*^9, 3.510647760395031*^9,
47 3.5106569081720953`*^9}]
48}, Open ]],
49
50Cell[CellGroupData[{
51
52Cell[BoxData[
53 RowBox[{"calRTEst", "=",
54 RowBox[{"{",
55 RowBox[{"c1E", ",", "c2E", ",", "c3E"}], "}"}]}]], "Input",
56 CellChangeTimes->{{3.5106350832369184`*^9, 3.5106351113587303`*^9}}],
57
58Cell[BoxData[
59 RowBox[{"{",
60 RowBox[{"c1E", ",", "c2E", ",", "c3E"}], "}"}]], "Output",
61 CellChangeTimes->{3.5106351123658314`*^9, 3.510647760410631*^9,
62 3.510656908197098*^9}]
63}, Open ]],
64
65Cell[CellGroupData[{
66
67Cell[BoxData[
68 RowBox[{"calRTWest", "=",
69 RowBox[{"{",
70 RowBox[{"c1W", ",", "c2W", ",", "c3W"}], "}"}]}]], "Input",
71 CellChangeTimes->{{3.5106351149400883`*^9, 3.510635130544649*^9}}],
72
73Cell[BoxData[
74 RowBox[{"{",
75 RowBox[{"c1W", ",", "c2W", ",", "c3W"}], "}"}]], "Output",
76 CellChangeTimes->{3.5106351312737217`*^9, 3.5106477604418316`*^9,
77 3.5106569082221003`*^9}]
78}, Open ]],
79
80Cell[CellGroupData[{
81
82Cell[BoxData[
83 RowBox[{"calBAOv0", "=",
84 RowBox[{"{",
85 RowBox[{"bao1v0", ",", "bao2v0", ",", "bao3v0"}], "}"}]}]], "Input",
86 CellChangeTimes->{{3.5106351357401686`*^9, 3.5106351689704914`*^9}}],
87
88Cell[BoxData[
89 RowBox[{"{",
90 RowBox[{"bao1v0", ",", "bao2v0", ",", "bao3v0"}], "}"}]], "Output",
91 CellChangeTimes->{3.5106351700045943`*^9, 3.510647760457432*^9,
92 3.5106569082461033`*^9}]
93}, Open ]],
94
95Cell[CellGroupData[{
96
97Cell[BoxData[
98 RowBox[{"calBAOv1", "=",
99 RowBox[{"{",
100 RowBox[{"bao1v1", ",", "bao2v1", ",", "bao3v1"}], "}"}]}]], "Input",
101 CellChangeTimes->{{3.5106351809366875`*^9, 3.5106351868982835`*^9}}],
102
103Cell[BoxData[
104 RowBox[{"{",
105 RowBox[{"bao1v1", ",", "bao2v1", ",", "bao3v1"}], "}"}]], "Output",
106 CellChangeTimes->{3.5106351877323675`*^9, 3.510647760473032*^9,
107 3.510656908271106*^9}]
108}, Open ]],
109
110Cell[CellGroupData[{
111
112Cell[BoxData[
113 RowBox[{"caliV0", "=",
114 RowBox[{"calRTEst", "/", "calBAOv0"}]}]], "Input",
115 CellChangeTimes->{{3.5106352159011836`*^9, 3.5106352792095137`*^9}}],
116
117Cell[BoxData[
118 RowBox[{"{",
119 RowBox[{
120 FractionBox["c1E", "bao1v0"], ",",
121 FractionBox["c2E", "bao2v0"], ",",
122 FractionBox["c3E", "bao3v0"]}], "}"}]], "Output",
123 CellChangeTimes->{3.5106352804596386`*^9, 3.5106477605042324`*^9,
124 3.5106569082961082`*^9}]
125}, Open ]],
126
127Cell[CellGroupData[{
128
129Cell[BoxData[
130 RowBox[{"caliV1", "=",
131 RowBox[{"calRTWest", "/", "calBAOv1"}]}]], "Input",
132 CellChangeTimes->{{3.5106352883114243`*^9, 3.510635290746668*^9}, {
133 3.5106568922185*^9, 3.5106568933526134`*^9}}],
134
135Cell[BoxData[
136 RowBox[{"{",
137 RowBox[{
138 FractionBox["c1W", "bao1v1"], ",",
139 FractionBox["c2W", "bao2v1"], ",",
140 FractionBox["c3W", "bao3v1"]}], "}"}]], "Output",
141 CellChangeTimes->{3.510635291263719*^9, 3.5106477605198326`*^9,
142 3.5106569083211107`*^9}]
143}, Open ]],
144
145Cell[CellGroupData[{
146
147Cell[BoxData[
148 RowBox[{"dataV0", "=",
149 RowBox[{"Table", "[",
150 RowBox[{
151 RowBox[{"{",
152 RowBox[{
153 RowBox[{"freq", "[",
154 RowBox[{"[", "n", "]"}], "]"}], ",",
155 RowBox[{"caliV0", "[",
156 RowBox[{"[", "n", "]"}], "]"}]}], "}"}], ",",
157 RowBox[{"{",
158 RowBox[{"n", ",", "1", ",",
159 RowBox[{"Length", "[", "freq", "]"}]}], "}"}]}], "]"}]}]], "Input",
160 CellChangeTimes->{{3.5106353155721498`*^9, 3.510635383830975*^9}}],
161
162Cell[BoxData[
163 RowBox[{"{",
164 RowBox[{
165 RowBox[{"{",
166 RowBox[{"f1", ",",
167 FractionBox["c1E", "bao1v0"]}], "}"}], ",",
168 RowBox[{"{",
169 RowBox[{"f2", ",",
170 FractionBox["c2E", "bao2v0"]}], "}"}], ",",
171 RowBox[{"{",
172 RowBox[{"f3", ",",
173 FractionBox["c3E", "bao3v0"]}], "}"}]}], "}"}]], "Output",
174 CellChangeTimes->{3.510635386760268*^9, 3.510647760551033*^9,
175 3.510656908347113*^9}]
176}, Open ]],
177
178Cell[CellGroupData[{
179
180Cell[BoxData[
181 RowBox[{"dataV1", "=",
182 RowBox[{"Table", "[",
183 RowBox[{
184 RowBox[{"{",
185 RowBox[{
186 RowBox[{"freq", "[",
187 RowBox[{"[", "n", "]"}], "]"}], ",",
188 RowBox[{"caliV1", "[",
189 RowBox[{"[", "n", "]"}], "]"}]}], "}"}], ",",
190 RowBox[{"{",
191 RowBox[{"n", ",", "1", ",",
192 RowBox[{"Length", "[", "freq", "]"}]}], "}"}]}], "]"}]}]], "Input",
193 CellChangeTimes->{{3.510635395911183*^9, 3.510635400627654*^9}}],
194
195Cell[BoxData[
196 RowBox[{"{",
197 RowBox[{
198 RowBox[{"{",
199 RowBox[{"f1", ",",
200 FractionBox["c1W", "bao1v1"]}], "}"}], ",",
201 RowBox[{"{",
202 RowBox[{"f2", ",",
203 FractionBox["c2W", "bao2v1"]}], "}"}], ",",
204 RowBox[{"{",
205 RowBox[{"f3", ",",
206 FractionBox["c3W", "bao3v1"]}], "}"}]}], "}"}]], "Output",
207 CellChangeTimes->{3.510635401669759*^9, 3.510647760566633*^9,
208 3.5106569083731155`*^9}]
209}, Open ]],
210
211Cell[CellGroupData[{
212
213Cell[BoxData[
214 RowBox[{
215 RowBox[{"InterpolatingPolynomial", "[",
216 RowBox[{"dataV0", ",", "f"}], "]"}], " ", "//", "InputForm"}]], "Input",
217 CellChangeTimes->{{3.5106355722628164`*^9, 3.510635578758466*^9}, {
218 3.5106479258039513`*^9, 3.510647928377984*^9}}],
219
220Cell["\<\
221c1E/bao1v0 + (f - f1)*((-(c1E/bao1v0) + c2E/bao2v0)/(-f1 + f2) +
222 ((f - f2)*(-((-(c1E/bao1v0) + c2E/bao2v0)/(-f1 + f2)) +
223 (-(c2E/bao2v0) + c3E/bao3v0)/(-f2 + f3)))/(-f1 + f3))\
224\>", "Output",
225 CellChangeTimes->{
226 3.5106355795115414`*^9, {3.510635611380728*^9, 3.510635617842374*^9},
227 3.5106477605978336`*^9, 3.510647864994372*^9, 3.5106479288927913`*^9,
228 3.510656908397118*^9}]
229}, Open ]],
230
231Cell[CellGroupData[{
232
233Cell[BoxData[
234 RowBox[{
235 RowBox[{"InterpolatingPolynomial", "[",
236 RowBox[{"dataV1", ",", "f"}], "]"}], "//", "InputForm"}]], "Input",
237 CellChangeTimes->{{3.5106477465576534`*^9, 3.5106477467916565`*^9}, {
238 3.510647934368461*^9, 3.510647937145297*^9}}],
239
240Cell["\<\
241c1W/bao1v1 + (f - f1)*((-(c1W/bao1v1) + c2W/bao2v1)/(-f1 + f2) +
242 ((f - f2)*(-((-(c1W/bao1v1) + c2W/bao2v1)/(-f1 + f2)) +
243 (-(c2W/bao2v1) + c3W/bao3v1)/(-f2 + f3)))/(-f1 + f3))\
244\>", "Output",
245 CellChangeTimes->{{3.5106477473064632`*^9, 3.510647760629034*^9},
246 3.510647901592441*^9, 3.5106479376757035`*^9, 3.5106569084201202`*^9}]
247}, Open ]]
248},
249WindowSize->{707, 761},
250WindowMargins->{{19, Automatic}, {Automatic, 31}},
251FrontEndVersion->"8.0 for Microsoft Windows (32-bit) (November 7, 2010)",
252StyleDefinitions->"Default.nb"
253]
254(* End of Notebook Content *)
255
256(* Internal cache information *)
257(*CellTagsOutline
258CellTagsIndex->{}
259*)
260(*CellTagsIndex
261CellTagsIndex->{}
262*)
263(*NotebookFileOutline
264Notebook[{
265Cell[557, 20, 206, 5, 31, "Input"],
266Cell[766, 27, 94, 1, 31, "Input"],
267Cell[863, 30, 116, 2, 31, "Input"],
268Cell[CellGroupData[{
269Cell[1004, 36, 180, 4, 31, "Input"],
270Cell[1187, 42, 177, 4, 30, "Output"]
271}, Open ]],
272Cell[CellGroupData[{
273Cell[1401, 51, 189, 4, 31, "Input"],
274Cell[1593, 57, 180, 4, 30, "Output"]
275}, Open ]],
276Cell[CellGroupData[{
277Cell[1810, 66, 188, 4, 31, "Input"],
278Cell[2001, 72, 184, 4, 30, "Output"]
279}, Open ]],
280Cell[CellGroupData[{
281Cell[2222, 81, 198, 4, 31, "Input"],
282Cell[2423, 87, 191, 4, 30, "Output"]
283}, Open ]],
284Cell[CellGroupData[{
285Cell[2651, 96, 198, 4, 31, "Input"],
286Cell[2852, 102, 189, 4, 30, "Output"]
287}, Open ]],
288Cell[CellGroupData[{
289Cell[3078, 111, 161, 3, 31, "Input"],
290Cell[3242, 116, 265, 7, 47, "Output"]
291}, Open ]],
292Cell[CellGroupData[{
293Cell[3544, 128, 209, 4, 31, "Input"],
294Cell[3756, 134, 263, 7, 47, "Output"]
295}, Open ]],
296Cell[CellGroupData[{
297Cell[4056, 146, 458, 13, 31, "Input"],
298Cell[4517, 161, 415, 13, 47, "Output"]
299}, Open ]],
300Cell[CellGroupData[{
301Cell[4969, 179, 456, 13, 31, "Input"],
302Cell[5428, 194, 417, 13, 47, "Output"]
303}, Open ]],
304Cell[CellGroupData[{
305Cell[5882, 212, 261, 5, 31, "Input"],
306Cell[6146, 219, 404, 8, 76, "Output"]
307}, Open ]],
308Cell[CellGroupData[{
309Cell[6587, 232, 256, 5, 31, "Input"],
310Cell[6846, 239, 353, 6, 76, "Output"]
311}, Open ]]
312}
313]
314*)
315
316(* End of internal cache information *)
Note: See TracBrowser for help on using the repository browser.