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: |
---|
10 | NotebookFileLineBreakTest |
---|
11 | NotebookFileLineBreakTest |
---|
12 | NotebookDataPosition[ 157, 7] |
---|
13 | NotebookDataLength[ 8731, 307] |
---|
14 | NotebookOptionsPosition[ 7215, 248] |
---|
15 | NotebookOutlinePosition[ 7559, 263] |
---|
16 | CellTagsIndexPosition[ 7516, 260] |
---|
17 | WindowFrame->Normal*) |
---|
18 | |
---|
19 | (* Beginning of Notebook Content *) |
---|
20 | Notebook[{ |
---|
21 | Cell[BoxData[ |
---|
22 | RowBox[{ |
---|
23 | RowBox[{"freq", "=", |
---|
24 | RowBox[{"{", |
---|
25 | RowBox[{"1280", ",", "1380", ",", "1410"}], "}"}]}], ";"}]], "Input", |
---|
26 | CellChangeTimes->{{3.510634896992296*^9, 3.5106349529288893`*^9}}], |
---|
27 | |
---|
28 | Cell[BoxData[""], "Input", |
---|
29 | CellChangeTimes->{{3.5106349560001965`*^9, 3.510635002484844*^9}}], |
---|
30 | |
---|
31 | Cell[BoxData[ |
---|
32 | RowBox[{"freq", "=."}]], "Input", |
---|
33 | CellChangeTimes->{{3.510635055560151*^9, 3.5106350577223673`*^9}}], |
---|
34 | |
---|
35 | Cell[CellGroupData[{ |
---|
36 | |
---|
37 | Cell[BoxData[ |
---|
38 | RowBox[{"freq", "=", |
---|
39 | RowBox[{"{", |
---|
40 | RowBox[{"f1", ",", "f2", ",", "f3"}], "}"}]}]], "Input", |
---|
41 | CellChangeTimes->{{3.510635061135709*^9, 3.5106350783984346`*^9}}], |
---|
42 | |
---|
43 | Cell[BoxData[ |
---|
44 | RowBox[{"{", |
---|
45 | RowBox[{"f1", ",", "f2", ",", "f3"}], "}"}]], "Output", |
---|
46 | CellChangeTimes->{3.510635079997595*^9, 3.510647760395031*^9, |
---|
47 | 3.5106569081720953`*^9}] |
---|
48 | }, Open ]], |
---|
49 | |
---|
50 | Cell[CellGroupData[{ |
---|
51 | |
---|
52 | Cell[BoxData[ |
---|
53 | RowBox[{"calRTEst", "=", |
---|
54 | RowBox[{"{", |
---|
55 | RowBox[{"c1E", ",", "c2E", ",", "c3E"}], "}"}]}]], "Input", |
---|
56 | CellChangeTimes->{{3.5106350832369184`*^9, 3.5106351113587303`*^9}}], |
---|
57 | |
---|
58 | Cell[BoxData[ |
---|
59 | RowBox[{"{", |
---|
60 | RowBox[{"c1E", ",", "c2E", ",", "c3E"}], "}"}]], "Output", |
---|
61 | CellChangeTimes->{3.5106351123658314`*^9, 3.510647760410631*^9, |
---|
62 | 3.510656908197098*^9}] |
---|
63 | }, Open ]], |
---|
64 | |
---|
65 | Cell[CellGroupData[{ |
---|
66 | |
---|
67 | Cell[BoxData[ |
---|
68 | RowBox[{"calRTWest", "=", |
---|
69 | RowBox[{"{", |
---|
70 | RowBox[{"c1W", ",", "c2W", ",", "c3W"}], "}"}]}]], "Input", |
---|
71 | CellChangeTimes->{{3.5106351149400883`*^9, 3.510635130544649*^9}}], |
---|
72 | |
---|
73 | Cell[BoxData[ |
---|
74 | RowBox[{"{", |
---|
75 | RowBox[{"c1W", ",", "c2W", ",", "c3W"}], "}"}]], "Output", |
---|
76 | CellChangeTimes->{3.5106351312737217`*^9, 3.5106477604418316`*^9, |
---|
77 | 3.5106569082221003`*^9}] |
---|
78 | }, Open ]], |
---|
79 | |
---|
80 | Cell[CellGroupData[{ |
---|
81 | |
---|
82 | Cell[BoxData[ |
---|
83 | RowBox[{"calBAOv0", "=", |
---|
84 | RowBox[{"{", |
---|
85 | RowBox[{"bao1v0", ",", "bao2v0", ",", "bao3v0"}], "}"}]}]], "Input", |
---|
86 | CellChangeTimes->{{3.5106351357401686`*^9, 3.5106351689704914`*^9}}], |
---|
87 | |
---|
88 | Cell[BoxData[ |
---|
89 | RowBox[{"{", |
---|
90 | RowBox[{"bao1v0", ",", "bao2v0", ",", "bao3v0"}], "}"}]], "Output", |
---|
91 | CellChangeTimes->{3.5106351700045943`*^9, 3.510647760457432*^9, |
---|
92 | 3.5106569082461033`*^9}] |
---|
93 | }, Open ]], |
---|
94 | |
---|
95 | Cell[CellGroupData[{ |
---|
96 | |
---|
97 | Cell[BoxData[ |
---|
98 | RowBox[{"calBAOv1", "=", |
---|
99 | RowBox[{"{", |
---|
100 | RowBox[{"bao1v1", ",", "bao2v1", ",", "bao3v1"}], "}"}]}]], "Input", |
---|
101 | CellChangeTimes->{{3.5106351809366875`*^9, 3.5106351868982835`*^9}}], |
---|
102 | |
---|
103 | Cell[BoxData[ |
---|
104 | RowBox[{"{", |
---|
105 | RowBox[{"bao1v1", ",", "bao2v1", ",", "bao3v1"}], "}"}]], "Output", |
---|
106 | CellChangeTimes->{3.5106351877323675`*^9, 3.510647760473032*^9, |
---|
107 | 3.510656908271106*^9}] |
---|
108 | }, Open ]], |
---|
109 | |
---|
110 | Cell[CellGroupData[{ |
---|
111 | |
---|
112 | Cell[BoxData[ |
---|
113 | RowBox[{"caliV0", "=", |
---|
114 | RowBox[{"calRTEst", "/", "calBAOv0"}]}]], "Input", |
---|
115 | CellChangeTimes->{{3.5106352159011836`*^9, 3.5106352792095137`*^9}}], |
---|
116 | |
---|
117 | Cell[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 | |
---|
127 | Cell[CellGroupData[{ |
---|
128 | |
---|
129 | Cell[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 | |
---|
135 | Cell[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 | |
---|
145 | Cell[CellGroupData[{ |
---|
146 | |
---|
147 | Cell[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 | |
---|
162 | Cell[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 | |
---|
178 | Cell[CellGroupData[{ |
---|
179 | |
---|
180 | Cell[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 | |
---|
195 | Cell[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 | |
---|
211 | Cell[CellGroupData[{ |
---|
212 | |
---|
213 | Cell[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 | |
---|
220 | Cell["\<\ |
---|
221 | c1E/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 | |
---|
231 | Cell[CellGroupData[{ |
---|
232 | |
---|
233 | Cell[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 | |
---|
240 | Cell["\<\ |
---|
241 | c1W/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 | }, |
---|
249 | WindowSize->{707, 761}, |
---|
250 | WindowMargins->{{19, Automatic}, {Automatic, 31}}, |
---|
251 | FrontEndVersion->"8.0 for Microsoft Windows (32-bit) (November 7, 2010)", |
---|
252 | StyleDefinitions->"Default.nb" |
---|
253 | ] |
---|
254 | (* End of Notebook Content *) |
---|
255 | |
---|
256 | (* Internal cache information *) |
---|
257 | (*CellTagsOutline |
---|
258 | CellTagsIndex->{} |
---|
259 | *) |
---|
260 | (*CellTagsIndex |
---|
261 | CellTagsIndex->{} |
---|
262 | *) |
---|
263 | (*NotebookFileOutline |
---|
264 | Notebook[{ |
---|
265 | Cell[557, 20, 206, 5, 31, "Input"], |
---|
266 | Cell[766, 27, 94, 1, 31, "Input"], |
---|
267 | Cell[863, 30, 116, 2, 31, "Input"], |
---|
268 | Cell[CellGroupData[{ |
---|
269 | Cell[1004, 36, 180, 4, 31, "Input"], |
---|
270 | Cell[1187, 42, 177, 4, 30, "Output"] |
---|
271 | }, Open ]], |
---|
272 | Cell[CellGroupData[{ |
---|
273 | Cell[1401, 51, 189, 4, 31, "Input"], |
---|
274 | Cell[1593, 57, 180, 4, 30, "Output"] |
---|
275 | }, Open ]], |
---|
276 | Cell[CellGroupData[{ |
---|
277 | Cell[1810, 66, 188, 4, 31, "Input"], |
---|
278 | Cell[2001, 72, 184, 4, 30, "Output"] |
---|
279 | }, Open ]], |
---|
280 | Cell[CellGroupData[{ |
---|
281 | Cell[2222, 81, 198, 4, 31, "Input"], |
---|
282 | Cell[2423, 87, 191, 4, 30, "Output"] |
---|
283 | }, Open ]], |
---|
284 | Cell[CellGroupData[{ |
---|
285 | Cell[2651, 96, 198, 4, 31, "Input"], |
---|
286 | Cell[2852, 102, 189, 4, 30, "Output"] |
---|
287 | }, Open ]], |
---|
288 | Cell[CellGroupData[{ |
---|
289 | Cell[3078, 111, 161, 3, 31, "Input"], |
---|
290 | Cell[3242, 116, 265, 7, 47, "Output"] |
---|
291 | }, Open ]], |
---|
292 | Cell[CellGroupData[{ |
---|
293 | Cell[3544, 128, 209, 4, 31, "Input"], |
---|
294 | Cell[3756, 134, 263, 7, 47, "Output"] |
---|
295 | }, Open ]], |
---|
296 | Cell[CellGroupData[{ |
---|
297 | Cell[4056, 146, 458, 13, 31, "Input"], |
---|
298 | Cell[4517, 161, 415, 13, 47, "Output"] |
---|
299 | }, Open ]], |
---|
300 | Cell[CellGroupData[{ |
---|
301 | Cell[4969, 179, 456, 13, 31, "Input"], |
---|
302 | Cell[5428, 194, 417, 13, 47, "Output"] |
---|
303 | }, Open ]], |
---|
304 | Cell[CellGroupData[{ |
---|
305 | Cell[5882, 212, 261, 5, 31, "Input"], |
---|
306 | Cell[6146, 219, 404, 8, 76, "Output"] |
---|
307 | }, Open ]], |
---|
308 | Cell[CellGroupData[{ |
---|
309 | Cell[6587, 232, 256, 5, 31, "Input"], |
---|
310 | Cell[6846, 239, 353, 6, 76, "Output"] |
---|
311 | }, Open ]] |
---|
312 | } |
---|
313 | ] |
---|
314 | *) |
---|
315 | |
---|
316 | (* End of internal cache information *) |
---|